OFFSET
1,1
PROG
(PARI) sopf(f) = sum(j=1, #f~, f[j, 1]); \\ A008472
lista(nn) = {/* nn should be > 10^7 */ my(nmax = 43, v = vector(nmax, k, List())); for (n=2, nn, my(f=factor(n), q); if (denominator(q=sigma(f)/(numdiv(f)*sopf(f))) == 1, if (q <= nmax, listput(v[q], n)); ); ); for (i=1, nmax, if (#v[i] == 0, break); print1(vecmax(Vec(v[i])), ", "); ); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Oct 06 2019
STATUS
approved