cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A301624 G.f. A(x) satisfies: A(x) = Product_{k>=1} (1 - x^k*A(x)^k)^k.

Original entry on oeis.org

1, -1, -1, 4, 1, -17, -6, 118, -8, -876, 625, 5966, -7486, -41937, 75969, 306312, -768637, -2164992, 7487063, 14461466, -70259884, -89410774, 646971980, 459817892, -5861484630, -1128608133, 52082250637, -15894742662, -453574650852, 366848121166, 3866670213663, -5215687717614
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 24 2018

Keywords

Examples

			G.f. A(x) = 1 - x - x^2 + 4*x^3 + x^4 - 17*x^5 - 6*x^6 + 118*x^7 - 8*x^8 - 876*x^9 + 625*x^10 + ...
G.f. A(x) satisfies: A(x) = (1 - x*A(x)) * (1 - x^2*A(x)^2)^2 * (1 - x^3*A(x)^3)^3 * (1 - x^4*A(x)^4)^4 * ...
log(A(x)) = -x - 3*x^2/2 + 8*x^3/3 + 13*x^4/4 - 51*x^5/5 - 120*x^6/6 + 538*x^7/7 + 781*x^8/8 - 5419*x^9/9 - 3053*x^10/10 + ... + A281267(n)*x^n/n + ...
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    Order := 33:
    Gser := solve(series(x*exp(add(sigma[2](n)*x^n/n, n = 1..32)), x) = y, x):
    seq(coeff(Gser, y^k), k = 1..32); # Peter Bala, Feb 09 2020

Formula

From Peter Bala, Feb 09 2020: (Start)
A(x) = 1/x * series reversion of ( exp( Sum_{n >= 1} sigma_2(n)*x^n/n ) ), where sigma_2(n) = A001157(n).
Equivalently, the o.g.f. A(x) satisfies [x^n](1/A(x))^n = sigma_2(n) for n >= 1. Cf. A066398. (End)
A(x) equals (1/x) * series reversion of (x * the o.g.f. for the sequence of planar partitions A000219). - Peter Bala, Feb 11 2020

A301625 G.f. A(x) satisfies: A(x) = Product_{k>=1} ((1 + x^k*A(x)^k)/(1 - x^k*A(x)^k))^k.

Original entry on oeis.org

1, 2, 10, 60, 398, 2820, 20892, 159868, 1253758, 10024070, 81400672, 669532924, 5566386324, 46701736772, 394910202608, 3362210548344, 28797181196766, 247955463799812, 2145088563952510, 18636002388075260, 162523319555310664, 1422259430668179592, 12485554521209720492, 109922263517662775292
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 24 2018

Keywords

Examples

			G.f. A(x) = 1 + 2*x + 10*x^2 + 60*x^3 + 398*x^4 + 2820*x^5 + 20892*x^6 + 159868*x^7 + 1253758*x^8 + ...
G.f. A(x) satisfies: A(x) = ((1 + x*A(x)) * (1 + x^2*A(x)^2)^2 * (1 + x^3*A(x)^3)^3 * ...)/((1 - x*A(x)) * (1 - x^2*A(x)^2)^2 * (1 - x^3*A(x)^3)^3 * ...).
log(A(x)) = 2*x + 16*x^2/2 + 128*x^3/3 + 1056*x^4/4 + 8952*x^5/5 + 77200*x^6/6 + 673948*x^7/7 + 5937792*x^8/8 + ... + A270924(n)*x^n/n + ...
		

Crossrefs

A291418 Expansion of the series reversion of Sum_{k>=1} x^(k*(k+1)/2).

Original entry on oeis.org

1, 0, -1, 0, 3, -1, -12, 9, 55, -67, -267, 468, 1323, -3180, -6513, 21267, 30969, -140581, -135995, 919698, 494361, -5954217, -829116, 38113425, -9433359, -240844482, 154219912, 1499076989, -1585801575, -9161079266, 13958031252, 54710928759, -113373461193, -317030478360, 875491422246
Offset: 1

Views

Author

Ilya Gutkovskiy, Aug 23 2017

Keywords

Comments

Reversion of g.f. (with constant term omitted) for A010054.

Crossrefs

Programs

  • Mathematica
    nmax = 35; Rest[CoefficientList[InverseSeries[Series[Sum[x^(k (k + 1)/2), {k, 1, nmax}], {x, 0, nmax}], x], x]]
    nmax = 35; Rest[CoefficientList[InverseSeries[Series[(-2 x^(1/8) + EllipticTheta[2, 0, Sqrt[x]])/(2 x^(1/8)), {x, 0, nmax}], x], x]]

Formula

G.f. A(x) satisfies: Sum_{k>=1} A(x)^(k*(k+1)/2) = x.
Showing 1-3 of 3 results.