A342395 -id:A342395 - cp's OEIS frontend

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

1, 5, 30, 264, 3135, 46709, 823564, 16777528, 387420759, 10000003265, 285311670666, 8916100500148, 302875106592331, 11112006826381965, 437893890380965260, 18446744073726350224, 827240261886336764313, 39346408075296928032645

Offset: 1

Seiichi Manyama, Mar 10 2021

Cf. A000217, A228640, A342394, A342395.

a[n_] := Sum[k^GCD[k, n], {k, 1, n}]; Array[a, 18] (* Amiram Eldar, Mar 10 2021 *)

PARI
```
a(n) = sum(k=1, n, k^gcd(k, n));
```

If p is prime, a(p) = A000217(p-1) + p^p = (p-1)*p/2 + p^p.

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

Original entry on oeis.org

1, 2, 6, 31, 355, 3150, 67172, 904085, 22998481, 427799450, 14914341926, 287337926355, 13421957361111, 339940911160914, 15434209582905140, 493467700905592777, 28101527071305611529, 836396358233559195382

Offset: 1

Seiichi Manyama, Mar 10 2021

Cf. A031971, A332621, A342394, A342395.

a[n_] := Sum[k^(n/GCD[k, n] - 1), {k, 1, n}]; Array[a, 18] (* Amiram Eldar, Mar 10 2021 *)

PARI
```
a(n) = sum(k=1, n, k^(n/gcd(k, n)-1));
```

If p is prime, a(p) = A031971(p-1) + 1.

cp's OEIS Frontend

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

Views

Author

Keywords

Crossrefs

Programs

Mathematica

PARI

Formula