1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | LIB "primdecint.lib"; |
---|
5 | |
---|
6 | ring R1 = integer,(a,b,c,d,e,f,g),dp; |
---|
7 | ideal I = 2*3,a2+2de+2cf+2bg+a, 2ab+e2+2df+2cg+b, |
---|
8 | b2+2ac+2ef+2dg+c, 2bc+2ad+f2+2eg+d, |
---|
9 | c2+2bd+2ae+2fg+e, 2cd+2be+2af+g2+f, d2+2ce+2bf+2ag+g; |
---|
10 | primdecZ(I,1); |
---|
11 | primdecZ(I,3); |
---|
12 | |
---|
13 | |
---|
14 | ring R3 = integer,(w,z,y,x),dp; |
---|
15 | ideal I = 2*3*181*32003, xzw+(-y^2+y)*z^2, (-x^2+x)*w^2+yzw, |
---|
16 | ((y^4-2*y^3+y^2)*x-y^4+y^3)*z^3,y2z2w+(-y*4+2*y^3-y^2)*z3; |
---|
17 | primdecZ(I,1); |
---|
18 | primdecZ(I,3); |
---|
19 | |
---|
20 | ring R5 = integer,(x,y,z),dp; |
---|
21 | ideal I = 2*13*181*32003, x2-y2-z2, xy-z2, y3+xz2-yz2+2z3+xy-z2, |
---|
22 | -y2z2+2z4+x2-y2+z2, y3z9+3y2z10+3yz11+z12-y2z2+2z4; |
---|
23 | primdecZ(I,1); |
---|
24 | primdecZ(I,3); |
---|
25 | |
---|
26 | ring R7 = integer,(x,y,z),dp; |
---|
27 | ideal I = 2^4*3*5*7*11*13*17, x2-y2-(z+2)^2, xy-(z+2)^2, |
---|
28 | y3+x*(z+2)^2-y*(z+2)^2+2*(z+2)^3+xy-(z+2)^2, |
---|
29 | -y^2*(z+2)^2+2*(z+2)^4+x2-y2+(z+2)^2, |
---|
30 | y3z9+3y2z10+3yz11+z12-y2z2+2z4; |
---|
31 | primdecZ(I,1); |
---|
32 | primdecZ(I,2); |
---|
33 | primdecZ(I,3); |
---|
34 | primdecZ(I,4); |
---|
35 | |
---|
36 | ring R8 = integer,(x,y,z),dp; |
---|
37 | ideal I = x2-y2-(z+2)^2, xy-(z+2)^2, |
---|
38 | y3+x*(z+2)^2-y*(z+2)^2+2*(z+2)^3+xy-(z+2)^2, |
---|
39 | -y^2*(z+2)^2+2*(z+2)^4+x2-y2+(z+2)^2, |
---|
40 | y3z9+3y2z10+3yz11+z12-y2z2+2z4; |
---|
41 | primdecZ(I,1); |
---|
42 | primdecZ(I,3); |
---|
43 | |
---|
44 | ring R9 = integer,(w,z,y,x),dp; |
---|
45 | ideal I = 2*3*5*630, ((y^2-y)*x-y^3+y^2)*z^2, (x-y)*zw, |
---|
46 | (x-y^2)*zw+(-y^2+y)*z^2, (-x^2+x)*w^2+(-yx+y)*zw; |
---|
47 | primdecZ(I,1); |
---|
48 | primdecZ(I,3); |
---|
49 | |
---|
50 | ring R10 = integer,(w,z,y,x),dp; |
---|
51 | ideal I = 2*3*5*1260, -yxzw+(-y^2+y)*z^2, (-x^2+x)*w^2-yxzw, |
---|
52 | ((-y^2+y)*x-y^3+2*y^2-y)*z^3, |
---|
53 | (y^2-y)*z^2*w+(-y^2+y)*z^2*w+(-y^2+y)*z^3; |
---|
54 | primdecZ(I,1); |
---|
55 | primdecZ(I,3); |
---|
56 | |
---|
57 | ring R11 = integer,(w,z,y,x),dp; |
---|
58 | ideal I = (4*y^2*x^2+(4*y^3+4*y^2-y)*x-y^2-y)*z^2, |
---|
59 | (x+y+1)*zw+(-4*y^2*x-4*y^3-4*y^2)*z^2, |
---|
60 | (-x-2*y^2 - 2*y - 1)*zw + (8*y^3*x + 8*y^4 + 8*y^3 + 2*y^2+y)*z^2, |
---|
61 | ((y^3 + y^2)*x - y^2 - y)*z^2, |
---|
62 | (y +1)*zw + (-y^3 -y^2)*z^2, |
---|
63 | (x + 1)*zw +(- y^2 -y)*z^2, |
---|
64 | (x^2 +x)*w^2 + (-yx - y)*zw; |
---|
65 | primdecZ(I,1); |
---|
66 | primdecZ(I,2); |
---|
67 | primdecZ(I,3); |
---|
68 | primdecZ(I,4); |
---|
69 | |
---|
70 | ring R12 = integer,(w,z,y,x),dp; |
---|
71 | ideal I = 2*3*181*72, |
---|
72 | ((y^3 + y^2)*x - y^2 - y)*z^2, |
---|
73 | (y + 1)*zw + (-y^3 -y^2)*z^2, |
---|
74 | (x + 1)*zw + (-y^2 -y)*z^2, (x^2 + x)*w^2 + (-yx - y)*zw; |
---|
75 | primdecZ(I,1); |
---|
76 | primdecZ(I,3); |
---|
77 | |
---|
78 | ring R14 = integer,(x(1),x(2),x(3),x(4)),dp; |
---|
79 | ideal I = 181*49^2, x(4)^4, x(1)*x(4)^3, x(1)*x(2)*x(4)^2, |
---|
80 | x(2)^2*x(4)^2, x(2)^2*x(3)*x(4), x(1)*x(2)*x(3)*x(4), |
---|
81 | x(1)*x(3)^2*x(4), x(3)^3*x(4); |
---|
82 | primdecZ(I,1); |
---|
83 | primdecZ(I,3); |
---|
84 | |
---|
85 | ring R15 = integer,(x,y,z),dp; |
---|
86 | ideal I = 32003*181*64, |
---|
87 | ((z^2-z)*y^2 + (z^2 -z)*y)*x, |
---|
88 | (z*y^3 + z*y^2)*x, |
---|
89 | (y^4 - y^2)*x, |
---|
90 | (z^2 - z)*y*x^2, |
---|
91 | (y^3 - y^2)*x^2, |
---|
92 | (z^3 - z^2)*x^4 + (2*z^3 -2*z^2)*x^3 + (z^3 -z^2)*x^2, |
---|
93 | z*y^2*x^2; |
---|
94 | primdecZ(I,1); |
---|
95 | primdecZ(I,3); |
---|
96 | |
---|
97 | ring R16 = integer,(x(1),x(2),x(3),x(4),x(5)),dp; |
---|
98 | ideal I = x(5)^5, x(1)*x(5)^4, x(1)*x(2)*x(5)^3, x(2)^2*x(5)^3, |
---|
99 | x(2)^2*x(3)*x(5)^2, x(1)*x(2)*x(3)*x(5)^2, x(1)*x(3)^2*x(5)^2, |
---|
100 | x(3)^3*x(5)^2, x(3)^3*x(4)*x(5), x(1)*x(3)^2*x(4)*x(5), |
---|
101 | x(1)*x(2)*x(3)*x(4)*x(5), x(2)^2*x(3)*x(4)*x(5), |
---|
102 | x(2)^2*x(4)^2*x(5), x(1)*x(2)*x(4)^2*x(5), x(1)*x(4)^3*x(5), |
---|
103 | x(4)^4*x(5); |
---|
104 | I = intersectZ(I,ideal(64*181,x(1)^2)); |
---|
105 | I; |
---|
106 | primdecZ(I,1); |
---|
107 | primdecZ(I,2); |
---|
108 | primdecZ(I,3); |
---|
109 | primdecZ(I,4); |
---|
110 | |
---|
111 | ring R17 = integer,(x,y,z),dp; |
---|
112 | ideal I = 374, (z+2)^8-140z6+2622*(z+2)^4-1820*(z+2)^2+169, |
---|
113 | 17y*(z+2)^4-374*y*(z+2)^2+221y+2z7-281z5+5240z3-3081z, |
---|
114 | 204y2+136yz3-3128yz+z6-149z4+2739z2+117, |
---|
115 | 17xz4-374xz2+221x+2z7-281z5+5240z3-3081z, |
---|
116 | 136xy-136xz-136yz+2z6-281z4+5376z2-3081, |
---|
117 | 204x2+136xz3-3128xz+z6-149z4+2739z2+117; |
---|
118 | primdecZ(I,1); |
---|
119 | primdecZ(I,3); |
---|
120 | |
---|
121 | ring R24 = integer,(a,b,c,d,e,f,g,h,j),dp; |
---|
122 | ideal I = 2*3*5*7*11*13*17, a+2b+c-d+g, f2gh-a, efg-c, fg2j-b, |
---|
123 | a+b+c+f+g-1, 3ad+3bd+2cd+df+dg-a-2b-c-g; |
---|
124 | primdecZ(I,1); |
---|
125 | primdecZ(I,3); |
---|
126 | |
---|
127 | ring R25 = integer,(s,r,b,u,v,w,x,y,z),dp; |
---|
128 | ideal I = 2*3*5*7*11*13*17*19*23, su+bv, ru+bw, rv+sw, sx+by, rx+bz, |
---|
129 | ry+sz, vx+uy, wx+uz, wy+vz; |
---|
130 | primdecZ(I,1); |
---|
131 | primdecZ(I,3); |
---|
132 | |
---|
133 | tst_status(1);$ |
---|