default(realprecision,38);
factor(x^2+I)
factor(x^2-1.)
factor(x^2+I+1.)
factor(x^2+Mod(1,5))
factor(x^2+Mod(1,3)+I)
factor(x^2+Mod(1,5)*I)
factor(x^2+(1+O(5))*I)
factor(x^2+(1+O(3))+I)
factor(x^2+(1+O(5)))
factor(x^2+quadgen(-3))
factor(x^2+quadgen(-3)*Mod(1,3))
factor(x^2+quadgen(-3)*Mod(1,5))
factor(x^2+quadgen(-3)*(1+O(5)))
factor(x^2+quadgen(-3)*(1+O(3)))
factor(x^2+Mod(y,y^2+1))
factor(x^2+Mod(y*Mod(1,3),y^2+1))
factor(x^2+Mod(y*Mod(1,5),y^2+1))
factor(x^2+Mod(y*(1+O(3)),y^2+1))
factor(5/2 + I/3)
factor(5^4 + 5^4*I)
factor(5 + I)
factor(15+3*I)
factor(4+4*I)
factor((3+4*I)/25)
factor(15,3)
q=31271192761826143388782348951/31274945109847936339856761591;
factor(q,2)
factor(q,500100)
factor(257377849,100)
factor((x-2.)*(x^2+1))
factor((x+1)^2/(x^2-4))
factor((x^2+1)/x)
factor(x^2+a^3*x-Mod(6,a^2-2))

p1 = (289*x^4 + 1)^8; p2 = (19*x^2 + 1)^16;
factor(p1)
factor(p2)

factor(x^2-1/y^2)
factor(y/x)
factor((x^3+y^3)^2*(x^4+y^4))
factor(x^3+y^3)
factor((x^7-y^7)*Mod(1,7))
factor((x^7-y)*Mod(1,7))
factor(x^5+y*x^4+y^2*x^3+(y^2-1)*x^2+(y^3-y)*x+(y^4-y^2))
factor((x^2+y^2) * Mod(1,13))
factor(x^3+I*y*x^2+(y^2+I)*x+(I*y^3-y))
factor((x^3+y^3+z^3-3*x*y*z)*Mod(1,7))
centerlift(factor((x^3+y^3+1-3*x*y)*(1+O(7^3))))
factor(x*y*(x+y))
E(n)=ellinit([2,1-n^2,(1-n^2)*(1+n),0,0]);
factor(elldivpol(ellinit(E(t)),4))
factor(0.*x^2+x+1)
factor(0.*x^2+x+I)

factor(x^2+1, I)
factor(x^2-1, 1.)
factor(x^2+I+1, 1.)
factor(x^2+1,Mod(1,5))
factor(x^2+Mod(1,3),I)
factor(x^2+1,O(5))
factor(x^2+O(5^2))
factor(x^2+(1+O(3)),I)
factor(x^2+x+1,quadgen(-3))
factor(x^2+quadgen(-3),Mod(1,3))
factor(x^2+quadgen(-3),Mod(1,5))
factor(x^2+quadgen(-3),O(5))
factor(x^2+quadgen(-3),O(3))
factor(x^2+y,Mod(1,y^2+1))
factor(x^2+Mod(y,y^2+1),Mod(1,3))
factor(x^2+Mod(y,y^2+1),Mod(1,5))
factor(x^2+Mod(y,y^2+1),O(3))
factor(5/2,I)
factor(x^2+a^3*x-6,Mod(1,a^2-2))
factor((x*y+1)*(x*y+2),Mod(1,2))
factor((x*y+1)*x*(y^2+y+1)*(z^3+z+1),Mod(1,2))
factor(1/(x^5+3*y*x^4+3*y^2*x^3+y^3*x^2))
factor((x-1)*(x-y)^2)

factor(x^3+y^3,quadgen(-3))
factor(x^7-y^7,Mod(1,7))
factor(x^7-y,Mod(1,7))
factor(x^2+y^2,Mod(1,13))
factor(x^3+y^3+z^3-3*x*y*z,Mod(1,7))
sum(i=1,10,setrand(i);polisirreducible((Mod(1,2)*((t^2+1)*x^4+(t^4+1)*x^2+1))))
lift(factor(Mod(1,2)*x^15+(Mod(1,2)*a^3+Mod(1,2)*a^2+Mod(1,2)*a+Mod(1,2))*x^13))
