login
A050214
a(n) = Product_{d|n, d^2<=n} (d+n/d); a(1)=1.
0
1, 3, 4, 20, 6, 35, 8, 54, 60, 77, 12, 728, 14, 135, 128, 1360, 18, 1881, 20, 2268, 220, 299, 24, 38500, 260, 405, 336, 5104, 30, 75361, 32, 7128, 476, 665, 432, 1731600, 38, 819, 640, 164164, 42, 218569, 44, 16200, 11592, 1175, 48, 5422144, 700, 20655, 1040
OFFSET
1,2
COMMENTS
a(n) = n+1 iff n is prime.
EXAMPLE
a(12) = (1+12)*(2+6)*(3+4) = 13*8*7 = 728.
MATHEMATICA
pdn[n_]:=Module[{d=Select[Divisors[n], #^2<=n&]}, Times@@ (Total/@ ({#, n/#}&/@d ))]; Join[{1}, Array[pdn, 60, 2]] (* Harvey P. Dale, May 03 2013 *)
CROSSREFS
Sequence in context: A041989 A196198 A041561 * A256605 A237884 A256532
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, May 08 2003
EXTENSIONS
Corrected and extended by Harvey P. Dale, May 03 2013
STATUS
approved