naif(n)=if(n<4,dd=[1,n],\
if(n%2==0,dd=[2,n/2],\
if(n%3==0,dd=[3,n/3],\
sq=sqrt(n);e=2;d=5;dd=[1,n];\
while(d<sq,if(n%d==0,\
dd=[d,n/d];d=n,d=d+e;e=6-e)))));\
return(dd)

q=nextprime(10^40+987879798796876876876087987);\
nn=vector(18);for(i=1,18,nn[i]=nextprime(10^i)*q);\
for(i=1,18,print([i,nn[i]]));print();\
for(i=1,18,t=getabstime();pp=naif(nn[i]);\
print(i,": ",pp[1]," * ",pp[2],". time:",strtime(getabstime()-t)))




 
 
pollard(n)=x=Mod(random,n);y=x;d=1;\
while(d==1,\
for(j=1,10,x=x^2+2;y=y^2+2;y=y^2+2;d=d*(y-x));\
d=gcd(lift(d),n));\
return([d,n/d])

q=nextprime(10^40+987879798796876876876087987);\
nn=vector(18);for(i=1,18,nn[i]=nextprime(10^i)*q);\
for(i=1,18,print([i,nn[i]]));print();\
for(i=1,18,t=getabstime();pp=pollard(nn[i]);\
print(i,": ",pp[1]," * ",pp[2],". time:",strtime(getabstime()-t)))