OFFSET
1,4
PROG
(PARI) factz(n, minn) = {my(v=[]); fordiv(n, d, if ((d>=minn) && (d<=sqrtint(n)), w = factz(n/d, d); for (i=1, #w, w[i] = concat([d], w[i]); ); v = concat(v, w); ); ); concat(v, [[n]]); }
factorz(n) = factz(n, 2);
lista(nn) = {for (n=1, nn, my(vf = factorz(n)); my(vs = apply(x->vecsum(x), vf)); my(d = #vs - #Set(vs)); if (d>0, print1(d, ", ")); ); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Aug 15 2020
STATUS
approved