OFFSET
1,2
COMMENTS
The corresponding number of divisors are: 1, 2, 6, 12, 40.
LINKS
Michel Marcus, solveBA PARI script
PROG
(PARI) isok(k) = sigma(k, -1) == sigma(numdiv(k), -1);
(PARI) lista(nn, lim) = for (n=1, nn, my(ad = sigma(n, -1), v = solveBA(numerator(ad), denominator(ad), lim)); if (#v, my(vn = apply(numdiv, v), vs = select(x->x==n, vn, 1)); if (#vs, print1(v[vs[1]], ", ")); ); ); \\ using solveBA script
lista(10000, 10^30)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Michel Marcus, Jan 15 2026
STATUS
approved
