OFFSET
0,1
PROG
(PARI) p=7; print1(p, ", "); for(n=1, 50, v=divisors(p); r=sqrt(p); t=0; for(k=1, matsize(v)[2], if(v[k]>=r, t=k; break)); if(v[t]^2==p, u=t, u=t-1); if(v[t]-v[u]<2, u=u-1; t=t+1); print1(v[t]-v[u]", "); p=p*(v[t]-v[u]))
CROSSREFS
KEYWORD
nonn,hard
AUTHOR
Ralf Stephan, Apr 04 2003
STATUS
approved