spielwiese
|
Last change
on this file since 09f10e was
d2d057e,
checked in by Martin Lee <martinlee84@…>, 12 years ago
|
|
chg: input factorized poly to avoid diff problems on 64bit machines
|
-
Property mode set to
100644
|
|
File size:
2.3 KB
|
| Line | |
|---|
| 1 | LIB "tst.lib"; |
|---|
| 2 | tst_init(); |
|---|
| 3 | |
|---|
| 4 | proc testfactors (list l, poly f) |
|---|
| 5 | { |
|---|
| 6 | poly g= 1; |
|---|
| 7 | for (int i= 1; i <= size (l[1]); i++) |
|---|
| 8 | { |
|---|
| 9 | g= g*(l[1][i]^l[2][i]); |
|---|
| 10 | } |
|---|
| 11 | g == f; |
|---|
| 12 | l; |
|---|
| 13 | } |
|---|
| 14 | |
|---|
| 15 | ring r=(0,a),(x,y),lp; |
|---|
| 16 | minpoly= 1000*a8+a5+233*a+1; |
|---|
| 17 | |
|---|
| 18 | list l; |
|---|
| 19 | poly g= (-106671/379000000a5+24854343/379000a4+106671/379000a3-24854343/379000000a-106671/379000000)*(x4y7+(17872400039344000/999a7-76705579568000/999a6+329208496000/999a5+17870987127344/999a4-76699515568/999a3+329208496/999a2-1412912/999a+462696578797024/111))*(x6y2+(-228/961a2)*x5y8+(295976258105814000/961a7-1270284369551000/961a6+5451864247000/961a5+9546866436994/31a4-1270183946551/961a3+5451433247/961a2-23398559/961a+2224595746411451/31)*x4y8+(797000/961a7+797/961a4+185701/961)*x4y3-642/961*x3y3+(-100/961a3)*x3y2+(-753/961a2)*x3y+(-171954000/961a7+738000/961a6-171954/961a4+738/961a3-40065282/961)*x-831/961*y2)*(x+y)^2*(x2-xy+y2)^2*y^5*x^10; |
|---|
| 20 | |
|---|
| 21 | l= factorize (g); testfactors (l,g); |
|---|
| 22 | |
|---|
| 23 | poly f= (77/108000a6-17941/108a5-77/108a4+17941/108000a2+77/108000a)*(x11y6+(-27/192500a7-27/192500000000a6+6291/96250000a5+27/192500000a4-6291/192500000000a2+1465802973/192500000000a+6291/192500000))*(x19y8+(-11420770143875a7+49016180875a6-210369875a5-91358938151/8a4+392098447/8a3-1682959/8a2+7223/8a-10644157774107/4)*x18y18+(1140a7+57/50a4+13281/50)*x18y9+(-13/25a7+431/200a)*x17y19+(-50407607945a7+216341665a6-928505a5-10080724589/200a4+43268333/200a3-185701/200a2+797/200a-2348994530237/200)*x17y4+(-369/100a5)*x14y+(-831/200a7)*x13y3+(-41/25a3)*x10y13+(161/40a)*x10y3+(961/200a7)*x9y3+(-854966221210530a7+3669382923645a6-15748424565a5-85489863140553/100a4+733818567729/200a3-3149435913/200a2+13517961/200a-7968285181693743/40)*x8y15+(6384614902795584474915a7-27401780698653596340a6+117604209006925725a5+1276822032740226976983/200a4-1369980721822283067/50a3+4703796469769229/40a2-50469919203459/100a+148761527235353727360813/100)*x8-377/200*x7y8+(174267690a7-747930a6+3210a5+17426769/100a4-74793/100a3+321/100a2+4060437177/100)*x6y4+(1/2000a7+233/2000a3+1/2000a2)*x6y3+(753/200000a6+175449/200000a2+753/200000a)*x6y2+(253/200a2)*x4y9+(358412827449678400a7-1538252478322560a6+6601941967040a5+1791922464673992/5a4-38453271774064/25a3+165035501176/25a2-708306872/25a+2087754719897416864/25)*xy15+(443/100a7)*y9)*(x+y)^2*(x2-xy+y2)^2*x^5*y^5; |
|---|
| 24 | |
|---|
| 25 | l= factorize (f); testfactors (l,f); |
|---|
| 26 | |
|---|
| 27 | tst_status(1); $ |
|---|
Note: See
TracBrowser
for help on using the repository browser.
