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

1, 1, 9, 82, 1300, 15626, 376761, 6161988, 176787117, 3769318700, 142364319625, 3152513804548, 154718778284148, 4340009120036086, 210971169748692000, 7281661100510001416, 435659030617933827136, 14181101408561469791694, 1052864393300587929716721, 41673894815421072916530408

Offset: 1

Author

Ilya Gutkovskiy, May 30 2019

Keywords

• nonn

Links

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

Crossrefs

First superdiagonal of A308477.

Cf. A031971.

Programs

Mathematica

a[n_] := Sum[If[GCD[n, k] == 1, k^n, 0], {k, 1, n}]; Table[a[n], {n, 1, 20}]

PARI

a(n) = sum(k=1, n, (gcd(n,k)==1)*k^n) \\ Felix Fröhlich, May 30 2019