A344509 a(n) = (1/n) * Sum_{k=1..n} k * lcm(k,n).

1, 3, 8, 16, 35, 45, 98, 116, 183, 215, 396, 316, 663, 609, 750, 912, 1513, 1107, 2128, 1560, 2135, 2453, 3818, 2428, 4325, 4095, 4914, 4452, 7743, 4425, 9486, 7280, 8613, 9299, 10640, 8088, 16243, 13053, 14378, 12100, 22181, 12663, 25628, 17996, 19305, 23345, 33558, 19200, 34251, 26075

Offset

1,2

Links

Table of n, a(n) for n=1..50.

Formula

a(n) = A344508(n)/n.

Sum_{k=1..n} a(k) ~ zeta(3) * n^4 / (2*Pi^2). - Vaclav Kotesovec, May 29 2021

Mathematica

a[n_] := Sum[k * LCM[k, n], {k, 1, n}]/n; Array[a, 50] (* Amiram Eldar, May 21 2021 *)

Prog

(PARI) a(n) = sum(k=1, n, k*lcm(k, n))/n;

Crossrefs

Cf. A057661, A344508.

Sequence in context: A171497 A024623 A337118 * A201250 A196373 A027291

Adjacent sequences: A344506 A344507 A344508 * A344510 A344511 A344512

Keyword

nonn

Author

Seiichi Manyama, May 21 2021

Status

approved