login
A356238
a(n) = Sum_{k=1..n} (k * floor(n/k))^n.
3
1, 8, 62, 849, 8541, 206345, 2581403, 76623522, 1617299079, 49463993875, 952905453423, 59000021366675, 1198427462876421, 54128102218676115, 2321105129608323165, 117387839988330848902, 3205342976298888473968, 268263812478494295219717
OFFSET
1,2
LINKS
FORMULA
a(n) = Sum_{k=1..n} k^n * Sum_{d|k} (1 - (1 - 1/d)^n).
MATHEMATICA
a[n_] := Sum[(k * Floor[n/k])^n, {k, 1, n}]; Array[a, 18] (* Amiram Eldar, Jul 30 2022 *)
PROG
(PARI) a(n) = sum(k=1, n, (k*(n\k))^n);
(PARI) a(n) = sum(k=1, n, k^n*sumdiv(k, d, 1-(1-1/d)^n));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jul 30 2022
STATUS
approved