OFFSET
1,1
REFERENCES
S. Kahan, "Divisor Advisory", Journal of Recreational Mathematics 30(1) 41-4 1999-2000 Baywood NY.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Donovan Johnson)
EXAMPLE
504 = 2^3*3^2*7 is in the sequence because d(504) = A000005(504) = (3+1)*(2+1)*(1+1) = 24 = 2^3 + 3^2 + 7. Similarly for 32340 = 2^2*3*5*7^2*11, where d(32340) = 2^2 + 3 + 5 + 7^2 + 11 = 72.
MATHEMATICA
ndppdQ[n_]:=DivisorSigma[0, n]==Total[#[[1]]^#[[2]]&/@FactorInteger[n]]; Select[Range[2, 120000], ndppdQ] (* Harvey P. Dale, Nov 22 2013 *)
PROG
(PARI) isok(n) = my(f = factor(n)); numdiv(n) == sum(i=1, #f~, f[i, 1]^f[i, 2]); \\ Michel Marcus, Oct 26 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Oct 13 2003
EXTENSIONS
a(1) = 2 prepended by Michel Marcus, Oct 26 2013
STATUS
approved