login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A345263 a(n) = Sum_{d|n} d^rad(d). 0
1, 5, 28, 21, 3126, 46688, 823544, 85, 757, 10000003130, 285311670612, 3032688, 302875106592254, 11112006826381564, 437893890380862528, 341, 827240261886336764178, 34059641, 1978419655660313589123980, 10250000003146, 5842587018385982521381947992, 341427877364219557681958394200 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

If p is prime, a(p) = Sum_{d|p} d^rad(d) = 1^1 + p^p = p^p + 1.

LINKS

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

FORMULA

a(prime(n)) = A125137(n). - Michel Marcus, Jun 12 2021

EXAMPLE

a(4) = Sum_{d|4} d^rad(d) = 1^1 + 2^2 + 4^2 = 21.

a(6) = Sum_{d|6} d^rad(d) = 1^1 + 2^2 + 3^3 + 6^6 = 46688.

MATHEMATICA

Table[Sum[(1 - Ceiling[n/i] + Floor[n/i]) i^Product[k^((PrimePi[k] - PrimePi[k - 1]) (1 - Ceiling[i/k] + Floor[i/k])), {k, i}], {i, n}], {n, 30}]

PROG

(PARI) rad(n) = factorback(factorint(n)[, 1]);

a(n) = sumdiv(n, d, d^rad(d)); \\ Michel Marcus, Jun 12 2021

CROSSREFS

Cf. A007947 (rad), A125137.

Sequence in context: A321236 A351774 A344579 * A218346 A106679 A341061

Adjacent sequences:  A345260 A345261 A345262 * A345264 A345265 A345266

KEYWORD

nonn

AUTHOR

Wesley Ivan Hurt, Jun 12 2021

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 20 10:33 EDT 2022. Contains 353871 sequences. (Running on oeis4.)