A308481 - OEIS

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

(list; graph; refs; listen; history; text; internal format)

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

CROSSREFS

First superdiagonal of A308477.

Cf. A031971.