Changeset 0777ea in git for Tst/Short/factorizep_s.tst
 Timestamp:
 Jun 17, 2011, 12:08:12 PM (12 years ago)
 Branches:
 (u'spielwiese', 'ec94ef7a30b928574c0c3daf41f6804dff5f6b69')
 Children:
 96992c61ea7928aadcfccb69aecf86ad5f3bd13c
 Parents:
 d40294dfa77c784d669bfb161359a9b40c51af81
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

Tst/Short/factorizep_s.tst
rd40294 r0777ea 2 2 tst_init(); 3 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 16 17 // polys from L. Bernardin's thesis 18 ring r= 17,(x,y),dp; 19 list l; 20 poly f; 21 for (int n=2; n <= 700; n++) 22 { 23 f= x^n*y^n+x^((n div 2)+1)*y^(n div 2)*(y+1)+x^2*y+(n+1)*x*y+(n^2+3)*x2; 24 f; 25 l= factorize (f); testfactors (l, f); 26 } 27 tst_status(); 28 29 kill r; 30 4 31 ring r= 31991,(z,w),dp; 32 33 list l; 34 5 35 poly f= 10582*w^3*z^4+21325*z^4+29620*w^4*z^3+23697*w*z^3+12439*w^6*z^2+3572*w^3*z^2+5463*z^2+16590*w*z+24885*w^3+31963; 6 36 7 factorize (f); 37 l= factorize (f); testfactors (l, f); 38 39 kill r; 40 41 ring r= 3, (s,t),dp; 42 list l; 43 44 poly f= s^85+(t)*s^84+(t^2)*s^83+(1*t^18)*s^67+(1*t^19)*s^66+(t^20)*s^65 45 +(t^21)*s^64+(t^22)*s^63+(1*t^24)*s^61+(1*t^25)*s^60+(t^29)*s^56 46 +(t^38)*s^47+(1*t^39)*s^46+(1*t^40)*s^45+(1*t^42)*s^43+(1*t^43)*s^42 47 +(1*t^45)*s^40+(1*t^46)*s^39+(t^47)*s^38+(t^56)*s^29+(1*t^60)*s^25 48 +(1*t^61)*s^24+(t^63)*s^22+(t^64)*s^21+(t^65)*s^20+(1*t^66)*s^19 49 +(1*t^67)*s^18+(t^83)*s^2+(t^84)*s+(t^85); 50 51 l= factorize (f); testfactors (l, f); 52 53 kill r; 54 55 ring r= 32003, (f,g,v,y,u,x), dp; 56 list l; 57 poly h= f^3*u^34*f^2*g*v*y*u^2 + (8*f^3*v^2*y+6*f^2*g*v^2*x+4*f*g^2*v^2*y^2)*u 58 8*f^3*v^3*x4*f*g^2*v^3*x*y + g^3*v^3*x^2; 59 60 l= factorize (h); testfactors (l, h); 61 62 kill r; 63 ring r= 32003, (a,b,c,f), dp; 64 list l; 65 poly h= a^3*b^3*f^4+186*a^2*b^4*c*f^3+11532*a*b^5*c^2*f^2+14307*b^6*c^3*f; 66 67 l= factorize (h); testfactors (l, h); 68 69 kill r; 70 ring r= 32003,(x,y,z),dp; 71 list l; 72 poly f= x^25+y^23+z^18+x^4*y^10+x^8*y^5+2*x^5*y^8+x^5*y^5*z^3+x^6*y^6; 73 74 l= factorize (f); testfactors (l, f); 75 76 kill r; 77 ring r= 32003, (x,y,z,w),dp; 78 list l; 79 poly f= 80 573002*x^10+240517*x^9*y+396389*x^8*y^2+535082*x^7*y^3+288429*x^6*y^4 81 +763919*x^5*y^5+3766*x^4*y^6649376*x^3*y^759851*x^2*y^8542799*x*y^9+20979*y^10 82 461020*x^9*z650459*x^8*y*z944957*x^7*y^2*z+1677170*x^6*y^3*z+1666154*x^5*y^4*z 83 1232710*x^4*y^5*z+862729*x^3*y^6*z618762*x^2*y^7*z705768*x*y^8*z187036*y^9*z741032*x^8*z^2 84 +1446654*x^7*y*z^2643320*x^6*y^2*z^21635676*x^5*y^3*z^2189420*x^4*y^4*z^2+682529*x^3*y^5*z^2 85 450284*x^2*y^6*z^2327970*x*y^7*z^2+66106*y^8*z^2+386688*x^7*z^3+198264*x^6*y*z^3757313*x^5*y^2*z^3+2008508*x^4*y^3*z^3675589*x^3*y^4*z^32338517*x^2*y^5*z^3+413997*x*y^6*z^3+150157*y^7*z^3 86 +115299*x^6*z^4+1838075*x^5*y*z^4821489*x^4*y^2*z^4772496*x^3*y^3*z^4+94982*x^2*y^4*z^4 87 +552427*x*y^5*z^4+894534*y^6*z^4173749*x^5*z^51322895*x^4*y*z^5+545850*x^3*y^2*z^5 88 +1506535*x^2*y^3*z^5+31641*x*y^4*z^5+858761*y^5*z^5+188464*x^4*z^6+676365*x^3*y*z^6 89 1175321*x^2*y^2*z^694611*x*y^3*z^6+365391*y^4*z^6300517*x^3*z^774262*x^2*y*z^7 90 +756469*x*y^2*z^7+702099*y^3*z^7+651058*x^2*z^8+261997*x*y*z^8+70610*y^2*z^8493574*x*z^9 91 +129784*y*z^9+265065*z^10532141*x^9*w285474*x^8*y*w1300132*x^7*y^2*w+908974*x^6*y^3*w 92 1966687*x^5*y^4*w127968*x^4*y^5*w+58220*x^3*y^6*w+41679*x^2*y^7*w+479035*x*y^8*w+466832*y^9*w 93 +1566035*x^8*z*w+1817709*x^7*y*z*w688806*x^6*y^2*z*w748313*x^5*y^3*z*w1930577*x^4*y^4*z*w 94 957911*x^3*y^5*z*w+81915*x^2*y^6*z*w+133232*x*y^7*z*w+980198*y^8*z*w110745*x^7*z^2*w 95 1370738*x^6*y*z^2*w1182195*x^5*y^2*z^2*w+3215028*x^4*y^3*z^2*w2006230*x^3*y^4*z^2*w 96 +2364130*x^2*y^5*z^2*w+1457270*x*y^6*z^2*w+1484168*y^7*z^2*w+1648545*x^6*z^3*w+772215*x^5*y*z^3*w 97 +554722*x^4*y^2*z^3*w+667926*x^3*y^3*z^3*w+375356*x^2*y^4*z^3*w+3601160*x*y^5*z^3*w 98 +625131*y^6*z^3*w2510503*x^5*z^4*w2729919*x^4*y*z^4*w+1243946*x^3*y^2*z^4*w 99 +2503349*x^2*y^3*z^4*w+1909123*x*y^4*z^4*w+588366*y^5*z^4*w+94641*x^4*z^5*w+1274333*x^3*y*z^5*w 100 +1615834*x^2*y^2*z^5*w1286036*x*y^3*z^5*w491438*y^4*z^5*w864129*x^3*z^6*w+96287*x^2*y*z^6*w 101 1641596*x*y^2*z^6*w479951*y^3*z^6*w+865513*x^2*z^7*w+37697*x*y*z^7*w425879*y^2*z^7*w 102 246899*x*z^8*w1043521*y*z^8*w615277*z^9*w353183*x^8*w^2+1903124*x^7*y*w^2+307847*x^6*y^2*w^2 103 +153254*x^5*y^3*w^2+4494*x^4*y^4*w^2+775116*x^3*y^5*w^2+84747*x^2*y^6*w^2+1970329*x*y^7*w^2 104 +248796*y^8*w^2200676*x^7*z*w^21283772*x^6*y*z*w^2+312236*x^5*y^2*z*w^2+258350*x^4*y^3*z*w^2 105 +899318*x^3*y^4*z*w^2+879951*x^2*y^5*z*w^2197099*x*y^6*z*w^2+1027373*y^7*z*w^2 106 939408*x^6*z^2*w^2+1231368*x^5*y*z^2*w^2+811932*x^4*y^2*z^2*w^2+136753*x^3*y^3*z^2*w^2 107 1060314*x^2*y^4*z^2*w^2373787*x*y^5*z^2*w^2419717*y^6*z^2*w^2602528*x^5*z^3*w^2 108 279785*x^4*y*z^3*w^2+140803*x^3*y^2*z^3*w^2+1245653*x^2*y^3*z^3*w^2+611906*x*y^4*z^3*w^2 109 1166551*y^5*z^3*w^2878201*x^4*z^4*w^2+597628*x^3*y*z^4*w^2220458*x^2*y^2*z^4*w^2 110 3038116*x*y^3*z^4*w^22785259*y^4*z^4*w^2425903*x^3*z^5*w^2316206*x^2*y*z^5*w^2 111 3764903*x*y^2*z^5*w^22194636*y^3*z^5*w^2185521*x^2*z^6*w^2+419827*x*y*z^6*w^2 112 70946*y^2*z^6*w^2+561958*x*z^7*w^2513594*y*z^7*w^2487215*z^8*w^2902070*x^7*w^3 113 +260239*x^6*y*w^31618965*x^5*y^2*w^3+1123916*x^4*y^3*w^31107669*x^3*y^4*w^31433250*x^2*y^5*w^3 114 +665446*x*y^6*w^31463775*y^7*w^3+778712*x^6*z*w^380174*x^5*y*z*w^393922*x^4*y^2*z*w^3 115 +1249290*x^3*y^3*z*w^3+1272307*x^2*y^4*z*w^3+513673*x*y^5*z*w^3322676*y^6*z*w^3 116 634533*x^5*z^2*w^31875741*x^4*y*z^2*w^32698209*x^3*y^2*z^2*w^3+937176*x^2*y^3*z^2*w^3 117 2883756*x*y^4*z^2*w^31430757*y^5*z^2*w^3+485806*x^4*z^3*w^3+731831*x^3*y*z^3*w^3 118 +2196178*x^2*y^2*z^3*w^31491104*x*y^3*z^3*w^31456373*y^4*z^3*w^3206031*x^3*z^4*w^3+314576*x^2*y*z^4*w^3225867*x*y^2*z^4*w^32750466*y^3*z^4*w^31087247*x^2*z^5*w^3+2172840*x*y*z^5*w^3 119 541045*y^2*z^5*w^3+624945*x*z^6*w^3+1085162*y*z^6*w^3+565094*z^7*w^3457315*x^6*w^4 120 +342040*x^5*y*w^4681000*x^4*y^2*w^4+297164*x^3*y^3*w^4+689*x^2*y^4*w^4+529745*x*y^5*w^4 121 1310613*y^6*w^4786298*x^5*z*w^41023635*x^4*y*z*w^42098764*x^3*y^2*z*w^4+1733995*x^2*y^3*z*w^4 122 +1561806*x*y^4*z*w^4+133720*y^5*z*w^4240204*x^4*z^2*w^4+1124674*x^3*y*z^2*w^4 123 +2411690*x^2*y^2*z^2*w^4600532*x*y^3*z^2*w^471275*y^4*z^2*w^4+128215*x^3*z^3*w^4 124 +567626*x^2*y*z^3*w^4+973911*x*y^2*z^3*w^4562780*y^3*z^3*w^4894759*x^2*z^4*w^4 125 +1265582*x*y*z^4*w^4+509433*y^2*z^4*w^4671007*x*z^5*w^4+345729*y*z^5*w^4+80717*z^6*w^4 126 +1144765*x^5*w^5+137320*x^4*y*w^5+390809*x^3*y^2*w^5+207350*x^2*y^3*w^5+65943*x*y^4*w^5 127 1283191*y^5*w^5109804*x^4*z*w^5+951743*x^3*y*z*w^5154107*x^2*y^2*z*w^5394628*x*y^3*z*w^5 128 3201737*y^4*z*w^5+2307230*x^3*z^2*w^51082386*x^2*y*z^2*w^5240669*x*y^2*z^2*w^5 129 1274548*y^3*z^2*w^51768397*x^2*z^3*w^52788826*x*y*z^3*w^5+1689141*y^2*z^3*w^5 130 1051474*x*z^4*w^5+126201*y*z^4*w^5789876*z^5*w^5288558*x^4*w^6697110*x^3*y*w^6 131 +659535*x^2*y^2*w^61659535*x*y^3*w^6+371290*y^4*w^61605798*x^3*z*w^6+63599*x^2*y*z*w^6 132 +273476*x*y^2*z*w^6774681*y^3*z*w^6113162*x^2*z^2*w^6+861980*x*y*z^2*w^6 133 +1824646*y^2*z^2*w^6275200*x*z^3*w^6+1083821*y*z^3*w^6+803980*z^4*w^6 134 +912560*x^3*w^7373953*x^2*y*w^7+635559*x*y^2*w^7+2023871*y^3*w^7252773*x^2*z*w^7 135 1013176*x*y*z*w^7356255*y^2*z*w^7+412544*x*z^2*w^7199680*y*z^2*w^7+190561*z^3*w^7 136 +33346*x^2*w^81189055*x*y*w^8+246983*y^2*w^8220965*x*z*w^8273745*y*z*w^8493392*z^2*w^8 137 +411602*x*w^9+853400*y*w^9+1110886*z*w^9580548*w^10; 138 139 l= factorize (f); testfactors (l, f); 140 141 kill r; 142 ring r= 31991,(x,y,z,w),dp; 143 list l; 144 poly f= (15*y^2*z^16+29*w^4*x^12*y^12*z^3+21*x^3*z^2+3*w^15*y^20)*(1*z^31w^12*z^20+y^18y^14+x^2*y^2+x^21+w^2); 145 146 l= factorize (f); testfactors (l, f); 147 148 f= (w^4*z^3x*y^2*z^2w^4*x^5*y^6w^2*x^3*y)*(1*x^5z^3+y*z+x^2*y^3)*(w^4*z^6+y^2*z^3w^2*x^2*y^2*z^2+x^5*zx^4*y^2w^3*x^3*y); 149 150 l= factorize (f); testfactors (l, f); 151 152 kill r; 153 ring r= 2,(x,y),dp; 154 list l; 155 poly f= x*(x+y)*(x^2+x+y^2)*(x^2+x+y+1)*(x^2+x+y^2+1)*(x^3+x+y^3)*(x^3+x+y+1)*(x^3+x+y^3+1)*(x^3+x+y^3+y+1)*(x^3+x+y^3+y^2+1); 156 157 l= factorize (f); testfactors (l, f); 158 159 f= x*(x+y)*(x^2+x+y^2)*(x^2+x+y+1)*(x^2+x+y^2+1); 160 161 l= factorize (f); testfactors (l, f); 162 163 f= (x+y^2)*(x^2+y)*(x+y^4)*(x^4+y); 164 165 l= factorize (f); testfactors (l, f); 166 167 f= (x+y^2)*(x^2+y); 168 169 l= factorize (f); testfactors (l, f); 170 171 kill r; 172 173 ring r= 3,(x,y),dp; 174 list l; 175 poly f= (x+y^2)*(x^2+y)*(x+y^4)*(x^4+y); 176 177 l= factorize (f); testfactors (l, f); 178 179 f= x^3  y^3; 180 181 l= factorize (f); testfactors (l, f); 182 183 kill r; 184 185 ring r= 5,(x,y,z),dp; 186 187 list l; 188 189 poly f= (x*z+y+30)*(y*z+x+20)*(z+x*y+10); 190 191 l= factorize (f); testfactors (l, f); 192 193 f= (x+y^2)*(x^2+y)*(x+y^4)*(x^4+y); 194 195 l= factorize (f); testfactors (l, f); 196 197 f= x^10  y^4; 198 199 l= factorize (f); testfactors (l, f); 200 201 kill r; 202 ring r= 101, (a,b,c,d), dp; 203 list l; 204 poly f= (ab)*(a^47*b* c4)*(a^5b+c11*d^3+1); 205 206 l= factorize (f); testfactors (l, f); 207 208 kill r; 209 ring r= 31991, (x,y,z,w,u,v,s,t), dp; 210 list l; 211 poly f= x^6*y^3*z^2*(3*z^3+2*w*z8*x*y^2+14*w^2*y^2y^2+18*x^3*y)*(w^2*z^312*w^2*x*y*z^3+3*x*y^2+29*xw^2); 212 213 l= factorize (f); testfactors (l, f); 214 215 f= (22*y^518*x^4*y^526*x^3*y^438*x^2*y^4+29*x^2*y^341*x^4*y^2+37*x^4)*(33*x^5*y^6+11*y^2+35*x^3*y22*x^4); 216 217 l= factorize (f); testfactors (l, f); 218 219 f= (3*z^3+2*w*z9*y^3y^2+45*x^3)*(w^2*z^3+47*x*yw^2); 220 221 l= factorize (f); testfactors (l, f); 222 223 f= (z+y+x3)*(z+y+x3)*(z+y+x3); 224 225 l= factorize (f); testfactors (l, f); 226 227 f= (z+y+x3)*(z+y+x3)*(z+y+x3)*(z+y+x2)*(z+y+x2); 228 229 l= factorize (f); testfactors (l, f); 230 231 f= (z^2+x^3*y^4+u^2)*((y^2+x)*z^2 + 3*u^2*x^3*y^4*z+19*y^2)*(u^2*y^4*z^2+x^2*z+5); 232 233 l= factorize (f); testfactors (l, f); 234 235 f= (z^2x^3*y+3)*(z^2+x^3*y^4)*(y^4*z^2+x^2*z+5); 236 237 l= factorize (f); testfactors (l, f); 238 239 f= (y*z^3+x*y*z+y^2+x^3)*(x*(z^4+1)+z+x^3*y^2); 240 241 l= factorize (f); testfactors (l, f); 242 243 f= (x^3*(z+y)+z11)*(x^2*(z^2+y^2)+y+90); 244 245 l= factorize (f); testfactors (l, f); 246 247 f= x^10  y^4; 248 249 l= factorize (f); testfactors (l, f); 250 251 f= (x*z^2+2*x*z*w+x*w^2+z^2+w+2*z*w)*(x*z^2+2*x*z*w+x*w^2+z); 252 253 l= factorize (f); testfactors (l, f); 254 255 f= (x^2x*z+z*y+2)*(x^2+3*xy^215); 256 257 l= factorize (f); testfactors (l, f); 258 259 f= (z^3+x*y*z+y^2+x^3)*(x*(z^4+1)+z+x^3*y^2); 260 261 l= factorize (f); testfactors (l, f); 262 263 f= (x^2*y+x*y+1)*(x^2*y+2*x*y+3)*(2*x^2*y+3*x*y+5)*(3*x^2*y+4*x*y+7); 264 265 l= factorize (f); testfactors (l, f); 266 267 f= w^2*x^7*y^6*12*z^5+31962*z^2*y^5*x^7+87*z^5*y^3*x^7+31967*z^2*y^7*x^8+31759*z^2*y^5*x^8+31988*z^2*y^7*x^7+54*z^2*y^6*x^10+522*z^2*y^4*x^10+9*z^5*y^5*x^7+z^2*y^5*x^6*w^2+3*z^8*y^3*x^6*w^2+2*z^6*y^3*x^6*w^3+31989*z^3*y^3*x^6*w^3+31977*z^2*y^5*x^6*w^4+31988*z^5*y^3*x^6*w^2+414*z^2*y^5*x^7*w^2+42*z^2*y^7*x^7*w^2+14*z^5*y^5*x^6*w^4+31990*z^5*y^5*x^6*w^2+31973*z^2*y^4*x^9*w^2+58*z^3*y^3*x^7*w+31955*z^8*y^4*x^7*w^2+6*z^3*y^5*x^7*w+31967*z^6*y^4*x^7*w^3+96*z^5*y^6*x^8*w^2+31983*z^5*y^5*x^7*w^2+31823*z^5*y^6*x^7*w^4+31775*z^5*y^5*x^10*w^2+18*z^5*y^4*x^9*w^2; 268 269 l= factorize (f); testfactors (l, f); 270 271 f= (z+x*y+10)*(x*z+y+30)*(y*z+x+20); 272 273 l= factorize (f); testfactors (l, f); 274 275 f= (x^3*y+x^3*z+z11)*(x^2*z^2+x^2*y^2+y+90); 276 277 l= factorize (f); testfactors (l, f); 278 279 f= (xw)*(yw)*(xy)*(zw)*(xz)*(yz); 280 281 l= factorize (f); testfactors (l, f); 282 283 f= (10*y+x^3+9)*(6*x*y^2+x+y+x^3+15); 284 285 l= factorize (f); testfactors (l, f); 286 287 f= (3*z^3+2*w*z10*y^2+45*x^3)*(w^2*z^3+47*x*yw^2); 288 289 l= factorize (f); testfactors (l, f); 290 291 f= (x^210748*y*x+y^2)*(y^2+x^2)*(x^2+10748*y*x+y^2); 292 293 l= factorize (f); testfactors (l, f); 294 295 f= (y*zx*z+x^2+2)*(z^3y^2+x^2+3*x15); 296 297 l= factorize (f); testfactors (l, f); 298 299 f= (x+y+z)*(zy+x)*(xzy); 300 301 l= factorize (f); testfactors (l, f); 302 303 f= (y*z+x*y^2+9)*(y*z+x*z+x^2*y); 304 305 l= factorize (f); testfactors (l, f); 306 307 f= (3*z^2+2*y^2+x^2)*(4*z^23*y^2+2*x^2); 308 309 l= factorize (f); testfactors (l, f); 310 311 f= x^4 + y^3 + z^5; 312 313 l= factorize (f); testfactors (l, f); 314 315 f= 9944*(x^25265*x+14343)*(x^2+5943*x3555)*(x8995)*(x15035); 316 317 l= factorize (f); testfactors (l, f); 318 319 f= x^8 + x +1 ; 320 321 l= factorize (f); testfactors (l, f); 322 323 f= x^141; 324 325 l= factorize (f); testfactors (l, f); 8 326 9 327 tst_status(1);$
Note: See TracChangeset
for help on using the changeset viewer.