OFFSET
1,1
PROG
(PARI) sopf(f) = sum(j=1, #f~, f[j, 1]); \\ A008472
isok(k, n) = my(fk=factor(k)); n*numdiv(fk)*sopf(fk) == sigma(fk);
a(n) = {my(k=1); while (!isok(k, n), k++); k; }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Oct 06 2019
STATUS
approved