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
33
34
35
36
37
38
39
40
41
42
|
# Omega Calculator v1.2 (based on Omega Library 1.2, August, 2000):
# # From Fabien Coelho
#
# R := { [i, j, k, l, m, n, o, p] :
# 0 <= i <= 3 &&
# 0 <= j <= 3 &&
# 0 <= k <= 3 &&
# 0 <= l <= 3 &&
# 0 <= m <= 3 &&
# 0 <= n <= 3 &&
# 1 <= o <= 1000 && o <= p <= 2 o +1 &&
# 0 <= -1000 i + o + p <= 999 &&
# 0 <= -1000 j - o + 2p <= 999 &&
# -1 <= -1000 k + 3 o <= 998 &&
# 2 <= -1000 l + 2 o + p <= 1001 &&
# 0 <= -1000 m + o <= 999 &&
# 3 <= -1000 n + 2 p<= 1002
# } ;
#
#
# codegen R;
for(t1 = 0; t1 <= 3; t1++) {
for(t2 = max(2*t1-3,0); t2 <= min(intDiv(3*t1+2,2),3); t2++) {
for(t3 = t1; t3 <= min(-t2+2*t1+1,intDiv(3*t1+2,2),3); t3++) {
for(t4 = max(intDiv(2*t3-2+2,3),intDiv(-t2+5*t1-4+2,3),intDiv(3*t2+5*t3-6+5,6),intDiv(4*t2-3+2,3)); t4 <= min(intDiv(t2+5,2),intDiv(4*t3+3,3),intDiv(-t2+5*t1+4,3),3); t4++) {
for(t5 = max(intDiv(t2-2+2,3),intDiv(t3-2+2,3)); t5 <= intDiv(t3,3); t5++) {
for(t6 = max(intDiv(6*t4-4*t3-6+2,3),intDiv(2*t3-3+2,3),2*t4-4,intDiv(4*t2-3+2,3),intDiv(t3+3*t2-3+2,3),0); t6 <= min(intDiv(6*t4-4*t3+5,3),t2+1,t4); t6++) {
if (t4 >= t1-1) {
for(t7 = max(intDiv(2000*t1-1000*t2-999+2,3),1000*t1-500*t6-501,250*t4+1,400*t4-200*t2-199,intDiv(1000*t3-1+2,3),intDiv(1000*t2-2+2,3)); t7 <= min(1000*t1-500*t6+997,intDiv(1000*t4+1001,3),400*t4-200*t2+400,500*t1+499,500*t6+501,intDiv(1000*t3+998,3),1000); t7++) {
for(t8 = max(1000*t1-t7,500*t6+2,1000*t4-2*t7+2,t7,intDiv(t7+1000*t2+1,2)); t8 <= min(1000*t1-t7+999,500*t6+501,1000*t4-2*t7+1001,2*t7+1,intDiv(t7+1000*t2+999,2)); t8++) {
s1(t1,t2,t3,t4,t5,t6,t7,t8);
}
}
}
}
}
}
}
}
}
#
|
