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!)
A299021 G.f.: A(x) = Sum_{n>=1} a(n)*x^n = x*Product_{n>=1} (1 + a(n)*x^n)/(1 - a(n)*x^n). 0
1, 2, 6, 22, 86, 358, 1558, 6966, 31894, 148918, 705062, 3380054, 16381158, 80056550, 394266950, 1955139942, 9749771926, 48873487942, 246160229782, 1244801094742, 6318514387638, 32184084454166, 164425969781062, 842429440124854, 4327629345403078, 22283328480744070 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
LINKS
FORMULA
G.f.: A(x) = Sum_{n>=1} a(n)*x^n = x*exp(Sum_{k>=1} Sum_{n>=1} (1 + (-1)^(k+1))*a(n)^k*x^(n*k)/k).
EXAMPLE
G.f.: A(x) = x + 2*x^2 + 6*x^3 + 22*x^4 + 86*x^5 + ... = x * ((1 + x) * (1 + 2*x^2) * (1 + 6*x^3) * (1 + 22*x^4) * (1 + 86*x^5) * ...) / ((1 - x) * (1 - 2*x^2) * (1 - 6*x^3) * (1 - 22*x^4) * (1 - 86*x^5) * ...).
MATHEMATICA
a[n_] := a[n] = SeriesCoefficient[x Product[(1 + a[k] x^k)/(1 - a[k] x^k), {k, 1, n - 1}], {x, 0, n}]; a[1] = 1; Table[a[n], {n, 26}]
a[n_] := a[n] = SeriesCoefficient[x Exp[Sum[Sum[(1 + (-1)^(k + 1)) a[j]^k x^(j k)/k, {j, 1, n - 1}], {k, 1, n - 1}]], {x, 0, n}]; a[1] = 1; Table[a[n], {n, 26}]
CROSSREFS
Sequence in context: A150257 A150258 A073075 * A153475 A150259 A165531
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jun 18 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 29 06:15 EDT 2024. Contains 371265 sequences. (Running on oeis4.)