OFFSET
1,2
COMMENTS
EXAMPLE
PROG
(PARI) upto(n) = my(u = List(), t, res=List()); for(i=1, n, c=sigma(i); if(c<=n, listput(u, c))); listsort(u, 1); u=Vec(u); for(i=1, #u, t=1; d=divisors(u[i]); for(j=2, (#d + 1)\2, if(vecsearch(u, d[j]) > 0 && vecsearch(u, u[i]/d[j]) > 0, t=0; next(1))); if(t==1, listput(res, u[i]))); res
CROSSREFS
KEYWORD
nonn
AUTHOR
David A. Corneth, Jan 24 2019
STATUS
approved