|
|
A307100
|
|
a(n) = Sum_{d|n} (tau(d)*pod(d)) where tau(k) = the number of divisors of k (A000005) and pod(k) = the product of the divisors of k (A007955).
|
|
1
|
|
|
1, 5, 7, 29, 11, 155, 15, 285, 88, 415, 23, 10547, 27, 803, 917, 5405, 35, 35228, 39, 48439, 1785, 1963, 47, 2665011, 386, 2735, 3004, 132539, 59, 6481465, 63, 202013, 4385, 4663, 4925, 90744884, 75, 5819, 6117, 20528695, 83, 24896285, 87, 513091, 547748, 8515
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
n divides a(n) for n = 1, 21, 333592, ...
|
|
LINKS
|
|
|
FORMULA
|
a(p) = 2*(p+1) for p = primes (A000040).
|
|
EXAMPLE
|
a(6) = tau(1)*pod(1) + tau(2)*pod(2) + tau(3)*pod(3) + tau(6)*pod(6) = (1*1) + (2*2) + (2*3) + (4*36) = 155.
|
|
PROG
|
(Magma) [&+[# [c: c in Divisors(d)] * &*[c: c in Divisors(d)]: d in Divisors(n)]: n in [1..100]]
(PARI) a(n) = sumdiv(n, d, my(dd=divisors(d)); #dd*vecprod(dd)); \\ Michel Marcus, Apr 25 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|