login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A332652 a(n) = Sum_{k=1..n} n^(k/gcd(n, k)). 2
1, 4, 15, 76, 785, 7836, 137263, 2130976, 47895489, 1010012140, 28531167071, 743044702104, 25239592216033, 797785008119932, 31147773583464735, 1157442765678719056, 51702516367896047777, 2185932446984222457444, 109912203092239643840239, 5255987282125826560192520 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Table of n, a(n) for n=1..20.

FORMULA

a(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.

a(n) = n * A332653(n).

MATHEMATICA

Table[Sum[n^(k/GCD[n, k]), {k, 1, n}], {n, 1, 20}]

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

PROG

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

CROSSREFS

Cf. A031972, A056665, A057661, A226561, A228640, A332620, A332621, A332653, A332655.

Sequence in context: A002750 A178887 A002467 * A243327 A179511 A111726

Adjacent sequences:  A332649 A332650 A332651 * A332653 A332654 A332655

KEYWORD

nonn

AUTHOR

Ilya Gutkovskiy, Feb 18 2020

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 31 08:46 EDT 2020. Contains 338101 sequences. (Running on oeis4.)