|
|
A283535
|
|
a(n) = Sum_{d|n} d^(3*d + 1).
|
|
7
|
|
|
1, 129, 59050, 67108993, 152587890626, 609359740069674, 3909821048582988050, 37778931862957228818561, 523347633027360537213570571, 10000000000000000000152587890754, 255476698618765889551019445759400442, 8505622499821102144576132293474637113130
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
L.g.f.: -log(Product_{k>=1} (1 - x^k)^(k^(3*k))) = Sum_{k>=1} a(k)*x^k/k. - Seiichi Manyama, Jun 18 2019
|
|
EXAMPLE
|
a(6) = 1^(3+1) + 2^(6+1) + 3^(9+1) + 6^(18+1) = 609359740069674.
|
|
MATHEMATICA
|
f[n_] := Block[{d = Divisors[n]}, Total[d^(3 d + 1)]]; Array[f, 12] (* Robert G. Wilson v, Mar 10 2017 *)
|
|
PROG
|
(PARI) a(n) = sumdiv(n, d, d^(3*d+1)); \\ Michel Marcus, Mar 11 2017
(PARI) N=20; x='x+O('x^N); Vec(x*deriv(-log(prod(k=1, N, (1-x^k)^k^(3*k))))) \\ Seiichi Manyama, Jun 18 2019
|
|
CROSSREFS
|
Cf. Sum_{d|n} d^(k*d+1): A283498 (k=1), A283533 (k=2), this sequence (k=3).
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|