- Timestamp:
- Apr 22, 1998, 11:13:27 AM (26 years ago)
- Branches:
- (u'spielwiese', 'a719bcf0b8dbc648b128303a49777a094b57592c')
- Children:
- 03483470cb89c43677d7dc4875a95a0fb672f9a9
- Parents:
- c0971eb773a055d9b3183e7672eb903762e398c9
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Tst/Long/gcdp_l.tst
rc0971e rcc8350 1 // $Id: gcdp_l.tst,v 1. 1.1.1 1998-04-17 15:07:35 obachmanExp $1 // $Id: gcdp_l.tst,v 1.2 1998-04-22 09:13:27 schmidt Exp $ 2 2 3 3 // 4 4 // gcdp_l.tst - long tests for gcd calculations mod p. 5 5 // 6 // To Do: gcd calculations over transcendental extensions of 7 // finite fields. 6 // All univariate non-trivial examples come from gcdUnivP.fex 7 // or gcdUnivPAlpha.fex. 8 // 9 // The example in variables `a', `c', `e', ... came from `untitled2(6/p/1)' 10 // in stdMultivPGcd.in. 11 // 12 // To Do: 13 // 14 // o multivariate gcd calculations with parameters 8 15 // 9 16 10 17 LIB "tst.lib"; 11 18 tst_init(); 12 tst_ignore("CVS ID : $Id: gcdp_l.tst,v 1. 1.1.1 1998-04-17 15:07:35 obachmanExp $");19 tst_ignore("CVS ID : $Id: gcdp_l.tst,v 1.2 1998-04-22 09:13:27 schmidt Exp $"); 13 20 14 21 // … … 16 23 // 17 24 25 tst_ignore( "ring r1=32003,x,dp;" ); 18 26 ring r1=32003,x,dp; 19 27 20 28 poly f; 21 29 poly g; 30 31 f=106*x^13-583*x^12-1060*x^11+8056*x^10+1696*x^9-327*x^8+12508*x^7+8277*x^6+5609*x^5-12879*x^4+13144*x^3-2544*x^2; 32 g=11118*x^14-5404*x^13-908*x^12-13908*x^11+3188*x^10-8818*x^9+10439*x^8-14811*x^7+15530*x^6-4891*x^5+6322*x^4+15829*x^3-13686*x^2; 33 34 gcd(f, g); 35 36 f=-14613*x^17+2235*x^16-298*x^15+4540*x^14+7214*x^13+5494*x^12-3122*x^11-4720*x^10+8300*x^9-6582*x^8-9908*x^7-15983*x^6-5802*x^5-8634*x^4+7899*x^3+10556*x^2+6931*x+11063; 37 g=-10168*x^20+6674*x^19+3004*x^18+13113*x^17+9449*x^16+9097*x^15-6420*x^14+50*x^13+434*x^12-15226*x^11+3727*x^10-14065*x^9-9751*x^8-15792*x^7+6004*x^6-5059*x^5+2479*x^4-12504*x^3-11328*x^2-11338*x-8280; 38 39 gcd(f, g); 40 41 f=3812*x^22-6874*x^21+12586*x^20+3003*x^19-9568*x^18+11117*x^17 +7524*x^16+11138*x^15-9743*x^14+1892*x^13+12485*x^12-569*x^11-8265*x^10-5991*x^9+13701*x^8+2644*x^7-3936*x^6-15875*x^5+1289*x^4+3956*x^3-10099*x^2-6616*x+5401; 42 g=10652*x^20-4144*x^19-11810*x^18+8237*x^17-8675*x^16+6545*x^15-3601*x^14+14559*x^13+8090*x^12-8378*x^11+14255*x^10+8767*x^9-13932*x^8+11602*x^7-10751*x^6-4899*x^5+8637*x^4+14084*x^3-11583*x^2+5882*x+885; 43 44 gcd(f, g); 22 45 23 46 f=3812*x^22-6874*x^21+12586*x^20+3003*x^19-9568*x^18+11117*x^17+7524*x^16+11138*x^15-9743*x^14+1892*x^13+12485*x^12-569*x^11-8265*x^10-5991*x^9+13701*x^8+2644*x^7-3936*x^6-15875*x^5+1289*x^4+3956*x^3-10099*x^2-6616*x+5401; … … 46 69 47 70 // 48 // - ring r2=(32003,a),x,dp .71 // - ring r2=(32003,a),x,dp; minpoly=a^4+8734*a^3+a^2+11817*a+1. 49 72 // 50 73 74 tst_ignore( "ring r2=(32003,a),x,dp;" ); 51 75 ring r2=(32003,a),x,dp; 52 76 minpoly=a^4+8734*a^3+a^2+11817*a+1; … … 55 79 poly g; 56 80 57 // first, some of the aboveexamples81 // first, some of the r1 examples 58 82 59 83 f=3812*x^22-6874*x^21+12586*x^20+3003*x^19-9568*x^18+11117*x^17+7524*x^16+11138*x^15-9743*x^14+1892*x^13+12485*x^12-569*x^11-8265*x^10-5991*x^9+13701*x^8+2644*x^7-3936*x^6-15875*x^5+1289*x^4+3956*x^3-10099*x^2-6616*x+5401; … … 100 124 gcd(f,g); 101 125 126 // 127 // - ring r3=(32003,a,c,e),(f,h,k),dp. 128 // 129 130 tst_ignore( "ring r3=(32003,a,c,e),(f,h,k),dp;" ); 131 ring r3=(32003,a,c,e),(f,h,k),dp; 132 133 poly fPoly; 134 poly gPoly; 135 136 fPoly=(10501*c^11*e^3*f^3*k^4+500*c^10*e^4*f^2*k^5+12001*c^10*e^3*f^4*h*k^3-500*c^9*e^5*f*k^6-14501*c^9*e^4*f^3*h*k^4-3500*c^9*e^3*f^5*h^2*k^2-10501*c^8*e^6*k^7+4000*c^8*e^5*f^2*h*k^5+7000*c^8*e^4*f^4*h^2*k^3-7001*c^8*e^3*f^6*h^3*k-12501*c^7*e^6*f*h*k^6-13501*c^7*e^5*f^3*h^2*k^4-12001*c^7*e^3*f^7*h^4+11001*c^6*e^7*h*k^7+9001*c^6*e^6*f^2*h^2*k^5+13002*c^6*e^5*f^4*h^3*k^3+3001*c^6*e^4*f^6*h^4*k+1000*c^5*e^7*f*h^2*k^6+14002*c^5*e^6*f^3*h^3*k^4-7001*c^5*e^5*f^5*h^4*k^2+4000*c^5*e^4*f^7*h^5-4003*c^4*e^7*f^2*h^3*k^5+7997*c^4*e^6*f^4*h^4*k^3-12002*c^4*e^5*f^6*h^5*k+2*c^3*e^8*f*h^3*k^6-7995*c^3*e^7*f^3*h^4*k^4-7997*c^3*e^6*f^5*h^5*k^2-16001*c^3*e^5*f^7*h^6+16001*c^2*e^9*h^3*k^7-3*c^2*e^8*f^2*h^4*k^5+15998*c^2*e^7*f^4*h^5*k^3-c^2*e^6*f^6*h^6*k-16001*c*e^9*f*h^4*k^6+c*e^8*f^3*h^5*k^4-16001*c*e^7*f^5*h^6*k^2); 137 138 gPoly=(10501*c^9*f^8*k^4+500*c^8*e*f^7*k^5-9001*c^8*f^9*h*k^3-500*c^7*e^2*f^6*k^6+15002*c^7*e*f^8*h*k^4-5500*c^7*f^10*h^2*k^2-10501*c^6*e^3*f^5*k^7+9000*c^6*e^2*f^7*h*k^5-7503*c^6*e*f^9*h^2*k^3-4001*c^6*f^11*h^3*k+5001*c^5*e^3*f^6*h*k^6-496*c^5*e^2*f^8*h^2*k^4+3*c^5*e*f^10*h^3*k^2+8001*c^5*f^12*h^4+12001*c^4*e^4*f^5*h*k^7+5497*c^4*e^3*f^7*h^2*k^5+7995*c^4*e^2*f^9*h^3*k^3-8002*c^4*e*f^11*h^4*k-16000*c^3*e^4*f^6*h^2*k^6-15996*c^3*e^3*f^8*h^3*k^4-15999*c^3*e^2*f^10*h^4*k^2-8001*c^2*e^5*f^5*h^2*k^7-4003*c^2*e^4*f^7*h^3*k^5+15999*c^2*e^3*f^9*h^4*k^3-16001*c*e^5*f^6*h^3*k^6+8002*c*e^4*f^8*h^4*k^4-8001*e^5*f^7*h^4*k^5); 139 140 gcd(fPoly, gPoly); 141 142 fPoly=(8001*a*c^10*e^2*f^4*k^4+16001*a*c^9*e^3*f^3*k^5+15999*a*c^9*e^2*f^5*h*k^3+8001*a*c^8*e^4*f^2*k^6+7*a*c^8*e^3*f^4*h*k^4+8009*a*c^8*e^2*f^6*h^2*k^2+15995*a*c^7*e^4*f^3*h*k^5+15973*a*c^7*e^3*f^5*h^2*k^3-10*a*c^7*e^2*f^7*h^3*k+2*a*c^6*e^5*f^2*h*k^6+8037*a*c^6*e^4*f^4*h^2*k^4+38*a*c^6*e^3*f^6*h^3*k^2+4*a*c^6*e^2*f^8*h^4-20*a*c^5*e^5*f^3*h^2*k^5-54*a*c^5*e^4*f^5*h^3*k^3-16*a*c^5*e^3*f^7*h^4*k+4*a*c^4*e^6*f^2*h^2*k^6+34*a*c^4*e^5*f^4*h^3*k^4+24*a*c^4*e^4*f^6*h^4*k^2-8*a*c^3*e^6*f^3*h^3*k^5-16*a*c^3*e^5*f^5*h^4*k^3+4*a*c^2*e^6*f^4*h^4*k^4); 143 144 gPoly=(8001*c^11*f^5*k^4+8000*c^10*e*f^4*k^5+15999*c^10*f^6*h*k^3-8000*c^9*e^2*f^3*k^6-15992*c^9*e*f^5*h*k^4+8009*c^9*f^7*h^2*k^2-8001*c^8*e^3*f^2*k^7+15988*c^8*e^2*f^4*h*k^5+7964*c^8*e*f^6*h^2*k^3-10*c^8*f^8*h^3*k-15993*c^7*e^3*f^3*h*k^6-7936*c^7*e^2*f^5*h^2*k^4+48*c^7*e*f^7*h^3*k^2+4*c^7*f^9*h^4-2*c^6*e^4*f^2*h*k^7-8057*c^6*e^3*f^4*h^2*k^5-92*c^6*e^2*f^6*h^3*k^3-20*c^6*e*f^8*h^4*k+24*c^5*e^4*f^3*h^2*k^6+88*c^5*e^3*f^5*h^3*k^4+40*c^5*e^2*f^7*h^4*k^2-4*c^4*e^5*f^2*h^2*k^7-42*c^4*e^4*f^4*h^3*k^5-40*c^4*e^3*f^6*h^4*k^3+8*c^3*e^5*f^3*h^3*k^6+20*c^3*e^4*f^5*h^4*k^4-4*c^2*e^5*f^4*h^4*k^5); 145 146 gcd(fPoly, gPoly); 147 148 fPoly=(-8001*a^3*c^12*e*f^3*k^6-16000*a^3*c^11*e^2*f^2*k^7+3*a^3*c^11*e*f^4*h*k^5-2*a^3*c^10*e^3*f*k^8-8022*a^3*c^10*e^2*f^3*h*k^6+15988*a^3*c^10*e*f^5*h^2*k^4-15963*a^3*c^9*e^3*f^2*h*k^7-15889*a^3*c^9*e^2*f^4*h^2*k^5+29*a^3*c^9*e*f^6*h^3*k^3-16*a^3*c^8*e^4*f*h*k^8-250*a^3*c^8*e^3*f^3*h^2*k^6-8291*a^3*c^8*e^2*f^5*h^3*k^4-8033*a^3*c^8*e*f^7*h^4*k^2+174*a^3*c^7*e^4*f^2*h^2*k^7-15251*a^3*c^7*e^3*f^4*h^3*k^5-15593*a^3*c^7*e^2*f^6*h^4*k^3+18*a^3*c^7*e*f^8*h^5*k-32*a^3*c^6*e^5*f*h^2*k^8-664*a^3*c^6*e^4*f^3*h^3*k^6-1211*a^3*c^6*e^3*f^5*h^4*k^4-315*a^3*c^6*e^2*f^7*h^5*k^2-4*a^3*c^6*e*f^9*h^6+184*a^3*c^5*e^5*f^2*h^3*k^7+1278*a^3*c^5*e^4*f^4*h^4*k^5+1054*a^3*c^5*e^3*f^6*h^5*k^3+120*a^3*c^5*e^2*f^8*h^6*k-480*a^3*c^4*e^5*f^3*h^4*k^6-1260*a^3*c^4*e^4*f^5*h^5*k^4-444*a^3*c^4*e^3*f^7*h^6*k^2-16*a^3*c^4*e^2*f^9*h^7+32*a^3*c^3*e^6*f^2*h^4*k^7+568*a^3*c^3*e^5*f^4*h^5*k^5+568*a^3*c^3*e^4*f^6*h^6*k^3+64*a^3*c^3*e^3*f^8*h^7*k-64*a^3*c^2*e^6*f^3*h^5*k^6-272*a^3*c^2*e^5*f^5*h^6*k^4-80*a^3*c^2*e^4*f^7*h^7*k^2+32*a^3*c*e^6*f^4*h^6*k^5+32*a^3*c*e^5*f^6*h^7*k^3); 149 150 gPoly=(12001*c^14*k^8-15999*c^13*f*h*k^7-2*c^12*e*h*k^8+7980*c^12*f^2*h^2*k^6+32*c^11*e*f*h^2*k^7-15909*c^11*f^3*h^3*k^5-12*c^10*e^2*h^2*k^8-204*c^10*e*f^2*h^3*k^6+11761*c^10*f^4*h^4*k^4+144*c^9*e^2*f*h^3*k^7+664*c^9*e*f^3*h^4*k^5+370*c^9*f^5*h^5*k^3-32*c^8*e^3*h^3*k^8-648*c^8*e^2*f^2*h^4*k^6-1186*c^8*e*f^4*h^5*k^4-332*c^8*f^6*h^6*k^2+256*c^7*e^3*f*h^4*k^7+1392*c^7*e^2*f^3*h^5*k^5+1176*c^7*e*f^5*h^6*k^3+160*c^7*f^7*h^7*k-32*c^6*e^4*h^4*k^8-704*c^6*e^3*f^2*h^5*k^6-1548*c^6*e^2*f^4*h^6*k^4-608*c^6*e*f^6*h^7*k^2-32*c^6*f^8*h^8+128*c^5*e^4*f*h^5*k^7+896*c^5*e^3*f^3*h^6*k^5+864*c^5*e^2*f^5*h^7*k^3+128*c^5*e*f^7*h^8*k-192*c^4*e^4*f^2*h^6*k^6-544*c^4*e^3*f^4*h^7*k^4-192*c^4*e^2*f^6*h^8*k^2+128*c^3*e^4*f^3*h^7*k^5+128*c^3*e^3*f^5*h^8*k^3-32*c^2*e^4*f^4*h^8*k^4); 151 152 gcd(fPoly, gPoly); 153 154 fPoly=(-15939*a^2*c^18*e*f^5*k^10+15314*a^2*c^17*e^2*f^4*k^11-5813*a^2*c^17*e*f^6*h*k^9+12626*a^2*c^16*e^3*f^3*k^12-6189*a^2*c^16*e^2*f^5*h*k^10-13501*a^2*c^16*e*f^7*h^2*k^8+2938*a^2*c^15*e^4*f^2*k^13+6254*a^2*c^15*e^3*f^4*h*k^11-809*a^2*c^15*e^2*f^6*h^2*k^9+14497*a^2*c^15*e*f^8*h^3*k^7+6688*a^2*c^14*e^5*f*k^14+10933*a^2*c^14*e^4*f^3*h*k^12+2667*a^2*c^14*e^3*f^5*h^2*k^10-15807*a^2*c^14*e^2*f^7*h^3*k^8+12710*a^2*c^14*e*f^9*h^4*k^6+10376*a^2*c^13*e^6*k^15-11809*a^2*c^13*e^5*f^2*h*k^13+15983*a^2*c^13*e^4*f^4*h^2*k^11-7531*a^2*c^13*e^3*f^6*h^3*k^9-11595*a^2*c^13*e^2*f^8*h^4*k^7+7759*a^2*c^13*e*f^10*h^5*k^5-377*a^2*c^12*e^6*f*h*k^14-8237*a^2*c^12*e^5*f^3*h^2*k^12+12315*a^2*c^12*e^4*f^5*h^3*k^10+14837*a^2*c^12*e^3*f^7*h^4*k^8-2945*a^2*c^12*e^2*f^9*h^5*k^6-1669*a^2*c^12*e*f^11*h^6*k^4+7001*a^2*c^11*e^7*h*k^15-5095*a^2*c^11*e^6*f^2*h^2*k^13+15751*a^2*c^11*e^5*f^4*h^3*k^11-6433*a^2*c^11*e^4*f^6*h^4*k^9-1178*a^2*c^11*e^3*f^8*h^5*k^7+13741*a^2*c^11*e^2*f^10*h^6*k^5-6075*a^2*c^11*e*f^12*h^7*k^3-7011*a^2*c^10*e^7*f*h^2*k^14+9719*a^2*c^10*e^6*f^3*h^3*k^12+2122*a^2*c^10*e^5*f^5*h^4*k^10-7622*a^2*c^10*e^4*f^7*h^5*k^8-6837*a^2*c^10*e^3*f^9*h^6*k^6+496*a^2*c^10*e^2*f^11*h^7*k^4-7959*a^2*c^10*e*f^13*h^8*k^2-16000*a^2*c^9*e^8*h^2*k^15-4923*a^2*c^9*e^7*f^2*h^3*k^13+9471*a^2*c^9*e^6*f^4*h^4*k^11+9695*a^2*c^9*e^5*f^6*h^5*k^9-5510*a^2*c^9*e^4*f^8*h^6*k^7+10543*a^2*c^9*e^3*f^10*h^7*k^5+7714*a^2*c^9*e^2*f^12*h^8*k^3+15988*a^2*c^9*e*f^14*h^9*k-8022*a^2*c^8*e^8*f*h^3*k^14+2836*a^2*c^8*e^7*f^3*h^4*k^12+341*a^2*c^8*e^6*f^5*h^5*k^10+1267*a^2*c^8*e^5*f^7*h^6*k^8-5523*a^2*c^8*e^4*f^9*h^7*k^6-15137*a^2*c^8*e^3*f^11*h^8*k^4+94*a^2*c^8*e^2*f^13*h^9*k^2+2*a^2*c^8*e*f^15*h^10-15999*a^2*c^7*e^9*h^3*k^15+8068*a^2*c^7*e^8*f^2*h^4*k^13+1994*a^2*c^7*e^7*f^4*h^5*k^11-9422*a^2*c^7*e^6*f^6*h^6*k^9+7341*a^2*c^7*e^5*f^8*h^7*k^7+14503*a^2*c^7*e^4*f^10*h^8*k^5+15717*a^2*c^7*e^3*f^12*h^9*k^3-14*a^2*c^7*e^2*f^14*h^10*k-14*a^2*c^6*e^9*f*h^4*k^14-8072*a^2*c^6*e^8*f^3*h^5*k^12+487*a^2*c^6*e^7*f^5*h^6*k^10+9827*a^2*c^6*e^6*f^7*h^7*k^8-6364*a^2*c^6*e^5*f^9*h^8*k^6+489*a^2*c^6*e^4*f^11*h^9*k^4+42*a^2*c^6*e^3*f^13*h^10*k^2+a^2*c^5*e^10*h^4*k^15+33*a^2*c^5*e^9*f^2*h^5*k^13+7970*a^2*c^5*e^8*f^4*h^6*k^11+15252*a^2*c^5*e^7*f^6*h^7*k^9+6849*a^2*c^5*e^6*f^8*h^8*k^7+15479*a^2*c^5*e^5*f^10*h^9*k^5-70*a^2*c^5*e^4*f^12*h^10*k^3-5*a^2*c^4*e^10*f*h^5*k^14-36*a^2*c^4*e^9*f^3*h^6*k^12-15865*a^2*c^4*e^8*f^5*h^7*k^10+503*a^2*c^4*e^7*f^7*h^8*k^8+356*a^2*c^4*e^6*f^9*h^9*k^6+70*a^2*c^4*e^5*f^11*h^10*k^4+9*a^2*c^3*e^10*f^2*h^6*k^13-15989*a^2*c^3*e^9*f^4*h^7*k^11-117*a^2*c^3*e^8*f^6*h^8*k^9+15850*a^2*c^3*e^7*f^8*h^9*k^7-42*a^2*c^3*e^6*f^10*h^10*k^5-7*a^2*c^2*e^10*f^3*h^7*k^12+6*a^2*c^2*e^9*f^5*h^8*k^10+37*a^2*c^2*e^8*f^7*h^9*k^8+14*a^2*c^2*e^7*f^9*h^10*k^6+2*a^2*c*e^10*f^4*h^8*k^11-4*a^2*c*e^9*f^6*h^9*k^9-2*a^2*c*e^8*f^8*h^10*k^7); 155 156 gPoly=(-5313*c^17*e*f^3*k^13+5313*c^17*f^5*h*k^11+15939*c^16*e^2*f^2*k^14-5688*c^16*e*f^4*h*k^12-10251*c^16*f^6*h^2*k^10-15939*c^15*e^3*f*k^15-3813*c^15*e^2*f^3*h*k^13+1875*c^15*e*f^5*h^2*k^11-14126*c^15*f^7*h^3*k^9+5313*c^14*e^4*k^16-7563*c^14*e^3*f^2*h*k^14+7877*c^14*e^2*f^4*h^2*k^12-10376*c^14*e*f^6*h^3*k^10+4749*c^14*f^8*h^4*k^8-9251*c^13*e^4*f*h*k^15+5623*c^13*e^3*f^3*h^2*k^13-12131*c^13*e^2*f^5*h^3*k^11-13560*c^13*e*f^7*h^4*k^9-2684*c^13*f^9*h^5*k^7-11001*c^12*e^5*h*k^16+6878*c^12*e^4*f^2*h^2*k^14-8112*c^12*e^3*f^4*h^3*k^12+193*c^12*e^2*f^6*h^4*k^10-6954*c^12*e*f^8*h^5*k^8-13007*c^12*f^10*h^6*k^6-10002*c^11*e^5*f*h^2*k^15+12734*c^11*e^4*f^3*h^3*k^13-2722*c^11*e^3*f^5*h^4*k^11-12047*c^11*e^2*f^7*h^5*k^9+14032*c^11*e*f^9*h^6*k^7-1995*c^11*f^11*h^7*k^5-2000*c^10*e^6*h^2*k^16+8012*c^10*e^5*f^2*h^3*k^14+7371*c^10*e^4*f^4*h^4*k^12+713*c^10*e^3*f^6*h^5*k^10+7940*c^10*e^2*f^8*h^6*k^8-6032*c^10*e*f^10*h^7*k^6+15999*c^10*f^12*h^8*k^4+7997*c^9*e^6*f*h^3*k^15+15952*c^9*e^5*f^3*h^4*k^13+3916*c^9*e^4*f^5*h^5*k^11-7958*c^9*e^3*f^7*h^6*k^9+12079*c^9*e^2*f^9*h^7*k^7-15985*c^9*e*f^11*h^8*k^5-16001*c^9*f^13*h^9*k^3-16001*c^8*e^7*h^3*k^16-11978*c^8*e^6*f^2*h^4*k^14+9096*c^8*e^5*f^4*h^5*k^12-972*c^8*e^4*f^6*h^6*k^10+3903*c^8*e^3*f^8*h^7*k^8-46*c^8*e^2*f^10*h^8*k^6+15998*c^8*e*f^12*h^9*k^4+15996*c^7*e^7*f*h^4*k^15-8056*c^7*e^6*f^3*h^5*k^13+1925*c^7*e^5*f^5*h^6*k^11-9945*c^7*e^4*f^7*h^7*k^9-15932*c^7*e^3*f^9*h^8*k^7-15991*c^7*e^2*f^11*h^9*k^5-16001*c^6*e^8*h^4*k^16-15985*c^6*e^7*f^2*h^5*k^14-1942*c^6*e^6*f^4*h^6*k^12-13999*c^6*e^5*f^6*h^7*k^10-60*c^6*e^4*f^8*h^8*k^8+15984*c^6*e^3*f^10*h^9*k^6-2*c^5*e^8*f*h^5*k^15-21*c^5*e^7*f^3*h^6*k^13-22*c^5*e^6*f^5*h^7*k^11-15974*c^5*e^5*f^7*h^8*k^9-15984*c^5*e^4*f^9*h^9*k^7+3*c^4*e^8*f^2*h^6*k^14-15990*c^4*e^7*f^4*h^7*k^12-4*c^4*e^6*f^6*h^8*k^10+15991*c^4*e^5*f^8*h^9*k^8-2*c^3*e^8*f^3*h^7*k^13+16000*c^3*e^7*f^5*h^8*k^11-15998*c^3*e^6*f^7*h^9*k^9-16001*c^2*e^8*f^4*h^8*k^12+16001*c^2*e^7*f^6*h^9*k^10); 157 158 gcd(fPoly, gPoly); 159 160 // 161 // - ring r4=32003,(a,c,e,f,h,k),dp. 162 // 163 // The examples from r4 are those from r3 with parameters 164 // changed to ring variables. 165 // 166 167 tst_ignore( "ring r4=32003,(a,c,e,f,h,k),dp;" ); 168 ring r4=32003,(a,c,e,f,h,k),dp; 169 170 poly fPoly; 171 poly gPoly; 172 173 fPoly=(10501*c^11*e^3*f^3*k^4+500*c^10*e^4*f^2*k^5+12001*c^10*e^3*f^4*h*k^3-500*c^9*e^5*f*k^6-14501*c^9*e^4*f^3*h*k^4-3500*c^9*e^3*f^5*h^2*k^2-10501*c^8*e^6*k^7+4000*c^8*e^5*f^2*h*k^5+7000*c^8*e^4*f^4*h^2*k^3-7001*c^8*e^3*f^6*h^3*k-12501*c^7*e^6*f*h*k^6-13501*c^7*e^5*f^3*h^2*k^4-12001*c^7*e^3*f^7*h^4+11001*c^6*e^7*h*k^7+9001*c^6*e^6*f^2*h^2*k^5+13002*c^6*e^5*f^4*h^3*k^3+3001*c^6*e^4*f^6*h^4*k+1000*c^5*e^7*f*h^2*k^6+14002*c^5*e^6*f^3*h^3*k^4-7001*c^5*e^5*f^5*h^4*k^2+4000*c^5*e^4*f^7*h^5-4003*c^4*e^7*f^2*h^3*k^5+7997*c^4*e^6*f^4*h^4*k^3-12002*c^4*e^5*f^6*h^5*k+2*c^3*e^8*f*h^3*k^6-7995*c^3*e^7*f^3*h^4*k^4-7997*c^3*e^6*f^5*h^5*k^2-16001*c^3*e^5*f^7*h^6+16001*c^2*e^9*h^3*k^7-3*c^2*e^8*f^2*h^4*k^5+15998*c^2*e^7*f^4*h^5*k^3-c^2*e^6*f^6*h^6*k-16001*c*e^9*f*h^4*k^6+c*e^8*f^3*h^5*k^4-16001*c*e^7*f^5*h^6*k^2); 174 175 gPoly=(10501*c^9*f^8*k^4+500*c^8*e*f^7*k^5-9001*c^8*f^9*h*k^3-500*c^7*e^2*f^6*k^6+15002*c^7*e*f^8*h*k^4-5500*c^7*f^10*h^2*k^2-10501*c^6*e^3*f^5*k^7+9000*c^6*e^2*f^7*h*k^5-7503*c^6*e*f^9*h^2*k^3-4001*c^6*f^11*h^3*k+5001*c^5*e^3*f^6*h*k^6-496*c^5*e^2*f^8*h^2*k^4+3*c^5*e*f^10*h^3*k^2+8001*c^5*f^12*h^4+12001*c^4*e^4*f^5*h*k^7+5497*c^4*e^3*f^7*h^2*k^5+7995*c^4*e^2*f^9*h^3*k^3-8002*c^4*e*f^11*h^4*k-16000*c^3*e^4*f^6*h^2*k^6-15996*c^3*e^3*f^8*h^3*k^4-15999*c^3*e^2*f^10*h^4*k^2-8001*c^2*e^5*f^5*h^2*k^7-4003*c^2*e^4*f^7*h^3*k^5+15999*c^2*e^3*f^9*h^4*k^3-16001*c*e^5*f^6*h^3*k^6+8002*c*e^4*f^8*h^4*k^4-8001*e^5*f^7*h^4*k^5); 176 177 gcd(fPoly, gPoly); 178 179 fPoly=(8001*a*c^10*e^2*f^4*k^4+16001*a*c^9*e^3*f^3*k^5+15999*a*c^9*e^2*f^5*h*k^3+8001*a*c^8*e^4*f^2*k^6+7*a*c^8*e^3*f^4*h*k^4+8009*a*c^8*e^2*f^6*h^2*k^2+15995*a*c^7*e^4*f^3*h*k^5+15973*a*c^7*e^3*f^5*h^2*k^3-10*a*c^7*e^2*f^7*h^3*k+2*a*c^6*e^5*f^2*h*k^6+8037*a*c^6*e^4*f^4*h^2*k^4+38*a*c^6*e^3*f^6*h^3*k^2+4*a*c^6*e^2*f^8*h^4-20*a*c^5*e^5*f^3*h^2*k^5-54*a*c^5*e^4*f^5*h^3*k^3-16*a*c^5*e^3*f^7*h^4*k+4*a*c^4*e^6*f^2*h^2*k^6+34*a*c^4*e^5*f^4*h^3*k^4+24*a*c^4*e^4*f^6*h^4*k^2-8*a*c^3*e^6*f^3*h^3*k^5-16*a*c^3*e^5*f^5*h^4*k^3+4*a*c^2*e^6*f^4*h^4*k^4); 180 181 gPoly=(8001*c^11*f^5*k^4+8000*c^10*e*f^4*k^5+15999*c^10*f^6*h*k^3-8000*c^9*e^2*f^3*k^6-15992*c^9*e*f^5*h*k^4+8009*c^9*f^7*h^2*k^2-8001*c^8*e^3*f^2*k^7+15988*c^8*e^2*f^4*h*k^5+7964*c^8*e*f^6*h^2*k^3-10*c^8*f^8*h^3*k-15993*c^7*e^3*f^3*h*k^6-7936*c^7*e^2*f^5*h^2*k^4+48*c^7*e*f^7*h^3*k^2+4*c^7*f^9*h^4-2*c^6*e^4*f^2*h*k^7-8057*c^6*e^3*f^4*h^2*k^5-92*c^6*e^2*f^6*h^3*k^3-20*c^6*e*f^8*h^4*k+24*c^5*e^4*f^3*h^2*k^6+88*c^5*e^3*f^5*h^3*k^4+40*c^5*e^2*f^7*h^4*k^2-4*c^4*e^5*f^2*h^2*k^7-42*c^4*e^4*f^4*h^3*k^5-40*c^4*e^3*f^6*h^4*k^3+8*c^3*e^5*f^3*h^3*k^6+20*c^3*e^4*f^5*h^4*k^4-4*c^2*e^5*f^4*h^4*k^5); 182 183 gcd(fPoly, gPoly); 184 185 fPoly=(-8001*a^3*c^12*e*f^3*k^6-16000*a^3*c^11*e^2*f^2*k^7+3*a^3*c^11*e*f^4*h*k^5-2*a^3*c^10*e^3*f*k^8-8022*a^3*c^10*e^2*f^3*h*k^6+15988*a^3*c^10*e*f^5*h^2*k^4-15963*a^3*c^9*e^3*f^2*h*k^7-15889*a^3*c^9*e^2*f^4*h^2*k^5+29*a^3*c^9*e*f^6*h^3*k^3-16*a^3*c^8*e^4*f*h*k^8-250*a^3*c^8*e^3*f^3*h^2*k^6-8291*a^3*c^8*e^2*f^5*h^3*k^4-8033*a^3*c^8*e*f^7*h^4*k^2+174*a^3*c^7*e^4*f^2*h^2*k^7-15251*a^3*c^7*e^3*f^4*h^3*k^5-15593*a^3*c^7*e^2*f^6*h^4*k^3+18*a^3*c^7*e*f^8*h^5*k-32*a^3*c^6*e^5*f*h^2*k^8-664*a^3*c^6*e^4*f^3*h^3*k^6-1211*a^3*c^6*e^3*f^5*h^4*k^4-315*a^3*c^6*e^2*f^7*h^5*k^2-4*a^3*c^6*e*f^9*h^6+184*a^3*c^5*e^5*f^2*h^3*k^7+1278*a^3*c^5*e^4*f^4*h^4*k^5+1054*a^3*c^5*e^3*f^6*h^5*k^3+120*a^3*c^5*e^2*f^8*h^6*k-480*a^3*c^4*e^5*f^3*h^4*k^6-1260*a^3*c^4*e^4*f^5*h^5*k^4-444*a^3*c^4*e^3*f^7*h^6*k^2-16*a^3*c^4*e^2*f^9*h^7+32*a^3*c^3*e^6*f^2*h^4*k^7+568*a^3*c^3*e^5*f^4*h^5*k^5+568*a^3*c^3*e^4*f^6*h^6*k^3+64*a^3*c^3*e^3*f^8*h^7*k-64*a^3*c^2*e^6*f^3*h^5*k^6-272*a^3*c^2*e^5*f^5*h^6*k^4-80*a^3*c^2*e^4*f^7*h^7*k^2+32*a^3*c*e^6*f^4*h^6*k^5+32*a^3*c*e^5*f^6*h^7*k^3); 186 187 gPoly=(12001*c^14*k^8-15999*c^13*f*h*k^7-2*c^12*e*h*k^8+7980*c^12*f^2*h^2*k^6+32*c^11*e*f*h^2*k^7-15909*c^11*f^3*h^3*k^5-12*c^10*e^2*h^2*k^8-204*c^10*e*f^2*h^3*k^6+11761*c^10*f^4*h^4*k^4+144*c^9*e^2*f*h^3*k^7+664*c^9*e*f^3*h^4*k^5+370*c^9*f^5*h^5*k^3-32*c^8*e^3*h^3*k^8-648*c^8*e^2*f^2*h^4*k^6-1186*c^8*e*f^4*h^5*k^4-332*c^8*f^6*h^6*k^2+256*c^7*e^3*f*h^4*k^7+1392*c^7*e^2*f^3*h^5*k^5+1176*c^7*e*f^5*h^6*k^3+160*c^7*f^7*h^7*k-32*c^6*e^4*h^4*k^8-704*c^6*e^3*f^2*h^5*k^6-1548*c^6*e^2*f^4*h^6*k^4-608*c^6*e*f^6*h^7*k^2-32*c^6*f^8*h^8+128*c^5*e^4*f*h^5*k^7+896*c^5*e^3*f^3*h^6*k^5+864*c^5*e^2*f^5*h^7*k^3+128*c^5*e*f^7*h^8*k-192*c^4*e^4*f^2*h^6*k^6-544*c^4*e^3*f^4*h^7*k^4-192*c^4*e^2*f^6*h^8*k^2+128*c^3*e^4*f^3*h^7*k^5+128*c^3*e^3*f^5*h^8*k^3-32*c^2*e^4*f^4*h^8*k^4); 188 189 gcd(fPoly, gPoly); 190 191 fPoly=(-15939*a^2*c^18*e*f^5*k^10+15314*a^2*c^17*e^2*f^4*k^11-5813*a^2*c^17*e*f^6*h*k^9+12626*a^2*c^16*e^3*f^3*k^12-6189*a^2*c^16*e^2*f^5*h*k^10-13501*a^2*c^16*e*f^7*h^2*k^8+2938*a^2*c^15*e^4*f^2*k^13+6254*a^2*c^15*e^3*f^4*h*k^11-809*a^2*c^15*e^2*f^6*h^2*k^9+14497*a^2*c^15*e*f^8*h^3*k^7+6688*a^2*c^14*e^5*f*k^14+10933*a^2*c^14*e^4*f^3*h*k^12+2667*a^2*c^14*e^3*f^5*h^2*k^10-15807*a^2*c^14*e^2*f^7*h^3*k^8+12710*a^2*c^14*e*f^9*h^4*k^6+10376*a^2*c^13*e^6*k^15-11809*a^2*c^13*e^5*f^2*h*k^13+15983*a^2*c^13*e^4*f^4*h^2*k^11-7531*a^2*c^13*e^3*f^6*h^3*k^9-11595*a^2*c^13*e^2*f^8*h^4*k^7+7759*a^2*c^13*e*f^10*h^5*k^5-377*a^2*c^12*e^6*f*h*k^14-8237*a^2*c^12*e^5*f^3*h^2*k^12+12315*a^2*c^12*e^4*f^5*h^3*k^10+14837*a^2*c^12*e^3*f^7*h^4*k^8-2945*a^2*c^12*e^2*f^9*h^5*k^6-1669*a^2*c^12*e*f^11*h^6*k^4+7001*a^2*c^11*e^7*h*k^15-5095*a^2*c^11*e^6*f^2*h^2*k^13+15751*a^2*c^11*e^5*f^4*h^3*k^11-6433*a^2*c^11*e^4*f^6*h^4*k^9-1178*a^2*c^11*e^3*f^8*h^5*k^7+13741*a^2*c^11*e^2*f^10*h^6*k^5-6075*a^2*c^11*e*f^12*h^7*k^3-7011*a^2*c^10*e^7*f*h^2*k^14+9719*a^2*c^10*e^6*f^3*h^3*k^12+2122*a^2*c^10*e^5*f^5*h^4*k^10-7622*a^2*c^10*e^4*f^7*h^5*k^8-6837*a^2*c^10*e^3*f^9*h^6*k^6+496*a^2*c^10*e^2*f^11*h^7*k^4-7959*a^2*c^10*e*f^13*h^8*k^2-16000*a^2*c^9*e^8*h^2*k^15-4923*a^2*c^9*e^7*f^2*h^3*k^13+9471*a^2*c^9*e^6*f^4*h^4*k^11+9695*a^2*c^9*e^5*f^6*h^5*k^9-5510*a^2*c^9*e^4*f^8*h^6*k^7+10543*a^2*c^9*e^3*f^10*h^7*k^5+7714*a^2*c^9*e^2*f^12*h^8*k^3+15988*a^2*c^9*e*f^14*h^9*k-8022*a^2*c^8*e^8*f*h^3*k^14+2836*a^2*c^8*e^7*f^3*h^4*k^12+341*a^2*c^8*e^6*f^5*h^5*k^10+1267*a^2*c^8*e^5*f^7*h^6*k^8-5523*a^2*c^8*e^4*f^9*h^7*k^6-15137*a^2*c^8*e^3*f^11*h^8*k^4+94*a^2*c^8*e^2*f^13*h^9*k^2+2*a^2*c^8*e*f^15*h^10-15999*a^2*c^7*e^9*h^3*k^15+8068*a^2*c^7*e^8*f^2*h^4*k^13+1994*a^2*c^7*e^7*f^4*h^5*k^11-9422*a^2*c^7*e^6*f^6*h^6*k^9+7341*a^2*c^7*e^5*f^8*h^7*k^7+14503*a^2*c^7*e^4*f^10*h^8*k^5+15717*a^2*c^7*e^3*f^12*h^9*k^3-14*a^2*c^7*e^2*f^14*h^10*k-14*a^2*c^6*e^9*f*h^4*k^14-8072*a^2*c^6*e^8*f^3*h^5*k^12+487*a^2*c^6*e^7*f^5*h^6*k^10+9827*a^2*c^6*e^6*f^7*h^7*k^8-6364*a^2*c^6*e^5*f^9*h^8*k^6+489*a^2*c^6*e^4*f^11*h^9*k^4+42*a^2*c^6*e^3*f^13*h^10*k^2+a^2*c^5*e^10*h^4*k^15+33*a^2*c^5*e^9*f^2*h^5*k^13+7970*a^2*c^5*e^8*f^4*h^6*k^11+15252*a^2*c^5*e^7*f^6*h^7*k^9+6849*a^2*c^5*e^6*f^8*h^8*k^7+15479*a^2*c^5*e^5*f^10*h^9*k^5-70*a^2*c^5*e^4*f^12*h^10*k^3-5*a^2*c^4*e^10*f*h^5*k^14-36*a^2*c^4*e^9*f^3*h^6*k^12-15865*a^2*c^4*e^8*f^5*h^7*k^10+503*a^2*c^4*e^7*f^7*h^8*k^8+356*a^2*c^4*e^6*f^9*h^9*k^6+70*a^2*c^4*e^5*f^11*h^10*k^4+9*a^2*c^3*e^10*f^2*h^6*k^13-15989*a^2*c^3*e^9*f^4*h^7*k^11-117*a^2*c^3*e^8*f^6*h^8*k^9+15850*a^2*c^3*e^7*f^8*h^9*k^7-42*a^2*c^3*e^6*f^10*h^10*k^5-7*a^2*c^2*e^10*f^3*h^7*k^12+6*a^2*c^2*e^9*f^5*h^8*k^10+37*a^2*c^2*e^8*f^7*h^9*k^8+14*a^2*c^2*e^7*f^9*h^10*k^6+2*a^2*c*e^10*f^4*h^8*k^11-4*a^2*c*e^9*f^6*h^9*k^9-2*a^2*c*e^8*f^8*h^10*k^7); 192 193 gPoly=(-5313*c^17*e*f^3*k^13+5313*c^17*f^5*h*k^11+15939*c^16*e^2*f^2*k^14-5688*c^16*e*f^4*h*k^12-10251*c^16*f^6*h^2*k^10-15939*c^15*e^3*f*k^15-3813*c^15*e^2*f^3*h*k^13+1875*c^15*e*f^5*h^2*k^11-14126*c^15*f^7*h^3*k^9+5313*c^14*e^4*k^16-7563*c^14*e^3*f^2*h*k^14+7877*c^14*e^2*f^4*h^2*k^12-10376*c^14*e*f^6*h^3*k^10+4749*c^14*f^8*h^4*k^8-9251*c^13*e^4*f*h*k^15+5623*c^13*e^3*f^3*h^2*k^13-12131*c^13*e^2*f^5*h^3*k^11-13560*c^13*e*f^7*h^4*k^9-2684*c^13*f^9*h^5*k^7-11001*c^12*e^5*h*k^16+6878*c^12*e^4*f^2*h^2*k^14-8112*c^12*e^3*f^4*h^3*k^12+193*c^12*e^2*f^6*h^4*k^10-6954*c^12*e*f^8*h^5*k^8-13007*c^12*f^10*h^6*k^6-10002*c^11*e^5*f*h^2*k^15+12734*c^11*e^4*f^3*h^3*k^13-2722*c^11*e^3*f^5*h^4*k^11-12047*c^11*e^2*f^7*h^5*k^9+14032*c^11*e*f^9*h^6*k^7-1995*c^11*f^11*h^7*k^5-2000*c^10*e^6*h^2*k^16+8012*c^10*e^5*f^2*h^3*k^14+7371*c^10*e^4*f^4*h^4*k^12+713*c^10*e^3*f^6*h^5*k^10+7940*c^10*e^2*f^8*h^6*k^8-6032*c^10*e*f^10*h^7*k^6+15999*c^10*f^12*h^8*k^4+7997*c^9*e^6*f*h^3*k^15+15952*c^9*e^5*f^3*h^4*k^13+3916*c^9*e^4*f^5*h^5*k^11-7958*c^9*e^3*f^7*h^6*k^9+12079*c^9*e^2*f^9*h^7*k^7-15985*c^9*e*f^11*h^8*k^5-16001*c^9*f^13*h^9*k^3-16001*c^8*e^7*h^3*k^16-11978*c^8*e^6*f^2*h^4*k^14+9096*c^8*e^5*f^4*h^5*k^12-972*c^8*e^4*f^6*h^6*k^10+3903*c^8*e^3*f^8*h^7*k^8-46*c^8*e^2*f^10*h^8*k^6+15998*c^8*e*f^12*h^9*k^4+15996*c^7*e^7*f*h^4*k^15-8056*c^7*e^6*f^3*h^5*k^13+1925*c^7*e^5*f^5*h^6*k^11-9945*c^7*e^4*f^7*h^7*k^9-15932*c^7*e^3*f^9*h^8*k^7-15991*c^7*e^2*f^11*h^9*k^5-16001*c^6*e^8*h^4*k^16-15985*c^6*e^7*f^2*h^5*k^14-1942*c^6*e^6*f^4*h^6*k^12-13999*c^6*e^5*f^6*h^7*k^10-60*c^6*e^4*f^8*h^8*k^8+15984*c^6*e^3*f^10*h^9*k^6-2*c^5*e^8*f*h^5*k^15-21*c^5*e^7*f^3*h^6*k^13-22*c^5*e^6*f^5*h^7*k^11-15974*c^5*e^5*f^7*h^8*k^9-15984*c^5*e^4*f^9*h^9*k^7+3*c^4*e^8*f^2*h^6*k^14-15990*c^4*e^7*f^4*h^7*k^12-4*c^4*e^6*f^6*h^8*k^10+15991*c^4*e^5*f^8*h^9*k^8-2*c^3*e^8*f^3*h^7*k^13+16000*c^3*e^7*f^5*h^8*k^11-15998*c^3*e^6*f^7*h^9*k^9-16001*c^2*e^8*f^4*h^8*k^12+16001*c^2*e^7*f^6*h^9*k^10); 194 195 gcd(fPoly, gPoly); 102 196 $
Note: See TracChangeset
for help on using the changeset viewer.