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,3

Links

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

Mathematica

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

Prog

(PARI) a(n) = sum(k=1, n, (gcd(n, k)==1)*k^n) \\ Felix Fröhlich, May 30 2019
(Python)
from math import gcd
def A308481(n): return sum(k**n for k in range(1, n+1) if gcd(n, k)==1) # Chai Wah Wu, Apr 11 2026

Crossrefs

First superdiagonal of A308477.

Cf. A031971.

Sequence in context: A294645 A338663 A308668 * A041146 A320991 A015497

Adjacent sequences: A308478 A308479 A308480 * A308482 A308483 A308484

Keyword

nonn

Author

Ilya Gutkovskiy, May 30 2019

Status

approved