login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A158947 Inverse Euler transform of A156217. 2
1, 4, 21, 598, 1555, 497652, 299593, 320361028, 1178277701, 357046721884, 67546215517, 19351522090518670, 61054982558011, 1502551437369035044, 33657152197069739919, 45463945109198918808616, 128583032925805678351 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

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

FORMULA

a(n) = (1/n)*Sum_{d|n} sigma(d)^d*moebius(n/d).

MAPLE

A158947 := proc(n) add(numtheory[sigma](d)^d*numtheory[mobius](n/d), d=numtheory[divisors](n))/n ; end: seq( A158947(n), n=1..40) ; [From R. J. Mathar, Apr 02 2009]

MATHEMATICA

f[n_] := Block[{d = Divisors@n}, Plus @@ (DivisorSigma[1, d]^d*MoebiusMu[n/d])/n]; Array[f, 17] [From Robert G. Wilson v, May 04 2009]

CROSSREFS

Sequence in context: A048164 A118909 A225157 * A000868 A000875 A094046

Adjacent sequences:  A158944 A158945 A158946 * A158948 A158949 A158950

KEYWORD

easy,nonn

AUTHOR

Vladeta Jovovic (vladeta(AT)eunet.yu), Mar 31 2009

EXTENSIONS

Extended by R. J. Mathar, Apr 02 2009

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified May 21 07:04 EDT 2013. Contains 225477 sequences.