A356238 a(n) = Sum_{k=1..n} (k * floor(n/k))^n.

1, 8, 62, 849, 8541, 206345, 2581403, 76623522, 1617299079, 49463993875, 952905453423, 59000021366675, 1198427462876421, 54128102218676115, 2321105129608323165, 117387839988330848902, 3205342976298888473968, 268263812478494295219717

Offset: 1

Author

Seiichi Manyama, Jul 30 2022

Keywords

• nonn

Links

• Seiichi Manyama, Table of n, a(n) for n = 1..386

Crossrefs

Cf. A007778, A350109, A350125, A356239, A356240.

Programs

Mathematica

a[n_] := Sum[(k * Floor[n/k])^n, {k, 1, n}]; Array[a, 18] (* Amiram Eldar, Jul 30 2022 *)

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));

Formula

a(n) = Sum_{k=1..n} k^n * Sum_{d|k} (1 - (1 - 1/d)^n).