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!)
A321261 a(n) = [x^n] Product_{k>=1} (1 + x^k)^(sigma_n(k)-k^n). 2
1, 0, 1, 1, 17, 2, 859, 131, 105508, 40907, 72916903, 6834168, 228239366293, 27616985835, 2050004858009336, 352807044193881, 87173272463714343166, 6798224808203572198, 18318379579349549499397403, 1187836799227050499295342, 11258903016282277676462826232428 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,5
LINKS
FORMULA
a(n) = [x^n] exp(Sum_{k>=1} ( Sum_{d|k} (-1)^(k/d+1)*d*(sigma_n(d) - d^n) ) * x^k/k).
MATHEMATICA
Table[SeriesCoefficient[Product[(1 + x^k)^(DivisorSigma[n, k] - k^n), {k, 1, n}], {x, 0, n}], {n, 0, 20}]
Table[SeriesCoefficient[Exp[Sum[Sum[(-1)^(k/d + 1) d (DivisorSigma[n, d] - d^n), {d, Divisors[k]}] x^k/k, {k, 1, n}]], {x, 0, n}], {n, 0, 20}]
CROSSREFS
Sequence in context: A303733 A040288 A036171 * A367604 A040280 A134907
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 01 2018
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 March 28 17:42 EDT 2024. Contains 371254 sequences. (Running on oeis4.)