# Omega Calculator v1.2 (based on Omega Library 1.2, August, 2000): # symbolic m; # # I := {[p1,p2,c2,a1,a2,b1,b2] : # 1 <= p1 <= 4 && 1 <= p2 <= 2 && 17p1<=16+2m # && 1<=2c2+p2 <= 6 # && 0 <= 10c2+5p2-a2 <= 4 # && 0 <= 17p1-2a1 <= 16 && a1 <= 30 && a2-m+a1 <= 1 # && 17 <= 2b1-2a1+17p1 <= 18 # && b2 = 5 + a2-5p2-5c2}; # # I; {[p1,p2,c2,a1,a2,b1,a2-5p2-5c2+5]: 1, -2c2+1 <= p2 <= 2, -2c2+6 && a1 <= 30 && 1 <= p1 && 17p1+2b1 <= 18+2a1 && 17+2a1 <= 17p1+2b1 && a2 <= 5p2+10c2 && 5p2+10c2 <= 4+a2 && 2a1 <= 17p1 && 17p1 <= 16+2a1 && a1+a2 <= 1+m} # # codegen I; for(t1 = 1; t1 <= min(intDiv(2*m+16,17),4); t1++) { for(t2 = 1; t2 <= min(intDiv(2*m-17*t1+26,10),2); t2++) { for(t3 = 0; t3 <= min(intDiv(2*m-10*t2-17*t1+26,20),2); t3++) { for(t4 = intDiv(17*t1-16+1,2); t4 <= min(intDiv(17*t1,2),m-5*t2-10*t3+5,30); t4++) { for(t5 = 5*t2+10*t3-4; t5 <= min(m-t4+1,5*t2+10*t3); t5++) { for(t6 = intDiv(-17*t1+2*t4+17+1,2); t6 <= intDiv(-17*t1+2*t4+18,2); t6++) { s1(t1,t2,t3,t4,t5,t6,t5+-5*t2+-5*t3+5); } } } } } } #