login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

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