OFFSET
1,1
PROG
(PARI) p=21*4*5*13; 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(round(log(v[t]-v[u])/log(2))", "); p=p*(v[t]-v[u]))
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Ralf Stephan, Apr 04 2003
EXTENSIONS
Revised and more terms from Sean A. Irvine, Nov 25 2025
STATUS
approved
