OFFSET
1,2
COMMENTS
2*10^9 < a(31) <= 2938669883. a(32) <= 7511549827. a(33) <= 11754740251. a(34) <= 30050593523. - Donovan Johnson, Dec 22 2010
LINKS
Jason Earls, Smarandache iterations of the first kind on functions involving divisors and prime factors, in Smarandache Notions Journal (2004), Vol. 14.1, page 259.
PROG
(PARI) {m=25; z=11000000; v=listcreate(m); for(i=1, m, listinsert(v, -1, i)); for(n=1, z, c=1; b=1; k=n; while(b&&c<=m, d=divisors(k); i=matsize(d)[2]-1; p=if(i>0, d[i], 1); q=if(k==1, 1, vecmax(component(factor(k), 1))); a=abs(p-q); if(a==0, b=0, k=a; c++)); if(a==0, if(v[c]<0, v[c]=n; print1([c, n])))); print(); for(i=1, m, print1(v[i], ", "))}
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Jason Earls, Sep 23 2002
EXTENSIONS
Four more terms from Klaus Brockhaus, Oct 01 2002
a(24)-a(30) from Donovan Johnson, Dec 22 2010
STATUS
approved