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)
q=31271192761826143388782348951/31274945109847936339856761591;
factor(q,2)
factor(q,500100)
factor((x-2.)*(x^2+1))
factor((x+1)^2/(x^2-4))

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))
