OFFSET
1,2
PROG
(PARI) { /* replace M below by a <cr> (<ctrl+v><ctrl+m>) */ a=vector(1000000); c=0; s=Set(); q=0; for(p=1, 2000000, if(isprime(p), q+=p; v=divisors(q); for(i=1, length(v), if(setsearch(s, v[i])==0, c++; if(v[i]<1000000, a[v[i]]=c); print1(p, "M"); s=setunion(s, [v[i]]); break)))); m=0; i=1; while(a[i], if(a[i]>m, m=a[i]; print1(i, ", ")); i++) } (Klasen)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 03 2005
EXTENSIONS
More terms from Lambert Klasen (lambert.klasen(AT)gmx.net), Nov 08 2005
STATUS
approved