1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
# 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);
}
}
}
}
}
}
#
|