A332652 -id:A332652 - cp's OEIS frontend

A332653 a(n) = (1/n) * Sum_{k=1..n} n^(k/gcd(n, k)).

1, 2, 5, 19, 157, 1306, 19609, 266372, 5321721, 101001214, 2593742461, 61920391842, 1941507093541, 56984643437138, 2076518238897649, 72340172854919941, 3041324492229179281, 121440691499123469858, 5784852794328402307381, 262799364106291328009626

Offset: 1

Ilya Gutkovskiy, Feb 18 2020

nonn

Cf. A023037, A056665, A057661, A226561, A228640, A332620, A332621, A332652, A332655.

[(1/n)*&+[n^(k div Gcd(n,k)):k in [1..n]]:n in [1..21]]; // Marius A. Burtea, Feb 18 2020

Table[(1/n) Sum[n^(k/GCD[n, k]), {k, 1, n}], {n, 1, 20}]
Table[Sum[Sum[If[GCD[k, d] == 1, n^(k - 1), 0], {k, 1, d}], {d, Divisors[n]}], {n, 1, 20}]

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

A332655 a(n) = Sum_{k=1..n} (k/gcd(n, k))^n.

Original entry on oeis.org

1, 2, 10, 84, 1301, 15693, 376762, 6168552, 176787631, 3770427352, 142364319626, 3152758480715, 154718778284149, 4340093860950619, 210971170836848270, 7281694486114555088, 435659030617933827137, 14181121059071691716406, 1052864393300587929716722, 41673907052879908244100770

Offset: 1

Ilya Gutkovskiy, Feb 18 2020

nonn

Cf. A031971, A057661, A226561, A332517, A332621, A332652, A332653, A332654.

[&+[(k div Gcd(n,k))^n:k in [1..n]]:n in [1..20]]; // Marius A. Burtea, Feb 18 2020

Table[Sum[(k/GCD[n, k])^n, {k, 1, n}], {n, 1, 20}]
Table[Sum[Sum[If[GCD[k, d] == 1, k^n, 0], {k, 1, d}], {d, Divisors[n]}], {n, 1, 20}]

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

a(n) = Sum_{d|n} Sum_{k=1..d, gcd(k, d) = 1} k^n.

cp's OEIS Frontend

A332653 a(n) = (1/n) * Sum_{k=1..n} n^(k/gcd(n, k)).

Views

Author

Keywords

Crossrefs

Programs

Magma

Mathematica

Formula