OFFSET
1,1
COMMENTS
Except for the first three, the members are all powers of 3. Proved by Luke Pebody, pers. comm.
LINKS
Sean A. Irvine, Table of n, a(n) for n = 1..21
PROG
(PARI) p=19; 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
AUTHOR
Ralf Stephan, Apr 04 2003
EXTENSIONS
Offset changed and more terms from Sean A. Irvine, Nov 24 2025
STATUS
approved
