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.

A320590 Expansion of Product_{k>=1} 1/(1 - x^k/(1 + x)^k).

Original entry on oeis.org

1, 1, 1, 0, 1, -2, 5, -12, 28, -63, 137, -290, 604, -1253, 2617, -5537, 11870, -25666, 55617, -120103, 257582, -548119, 1158437, -2437114, 5117165, -10748530, 22621055, -47728657, 100932549, -213750621, 452855190, -958925784, 2028187595, -4283531490, 9033779224
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 16 2018

Keywords

Comments

The zero-based binomial transform of this sequence is A000070, and if we remove first terms it becomes A000041.

Links

Crossrefs

Cf. A000203, A103446, A218482, A320568.
Row n=1 of A175804 (except first term). Row n=0 is A281425.
The version for strict partitions is A320591, row n=1 of A378622, first column A293467.
A000009 counts strict integer partitions, differences A087897, A378972.
A000041 counts integer partitions, differences A002865.
Cf. A047966, A008284, A377056, A378621.

Programs

  • Magma
    m:=50; R:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[1/(1 - x^k/(1 + x)^k): k in [1..(m+2)]]) )); // G. C. Greubel, Oct 29 2018
  • Maple
    seq(coeff(series(mul(1/(1-x^k/(1+x)^k),k=1..n),x,n+1), x, n), n = 0 .. 35); # Muniru A Asiru, Oct 16 2018
  • Mathematica
    nmax = 34; CoefficientList[Series[Product[1/(1 - x^k/(1 + x)^k), {k, 1, nmax}], {x, 0, nmax}], x]
nmax = 34; CoefficientList[Series[Exp[Sum[DivisorSigma[1, k] x^k/(k (1 + x)^k), {k, 1, nmax}]], {x, 0, nmax}], x]
  • PARI
    m=50; x='x+O('x^m); Vec(prod(k=1, m+2, 1/(1 - x^k/(1 + x)^k))) \\ G. C. Greubel, Oct 29 2018

Formula

G.f.: exp(Sum_{k>=1} x^k/(k*((1 + x)^k - x^k))).
G.f.: exp(Sum_{k>=1} sigma(k)*x^k/(k*(1 + x)^k)).

A320591 Expansion of Product_{k>=1} (1 + x^k/(1 + x)^k).

Original entry on oeis.org

1, 1, 0, 1, -2, 4, -7, 11, -16, 23, -36, 65, -129, 256, -473, 772, -1028, 835, 776, -5755, 17562, -41750, 86678, -165145, 299949, -541837, 1020029, -2068203, 4509512, -10252952, 23465297, -52762788, 115160832, -243018459, 496094524, -982431070, 1894710043, -3574095362
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 16 2018

Keywords

Comments

After the first term, this is the second term of the n-th differences of A000009, or column n=1 of A378622. - Gus Wiseman, Feb 03 2025

Links

Crossrefs

Cf. A000593, A320589, A307548.
The version for non-strict partitions is A320590, row n=1 of A175804.
Column n=1 (except first term) of A378622. See also A293467, A377285, A378970, A378971, A380412 (column n=0).
A000009 counts strict integer partitions, differences A087897, A378972.
A266232 gives zero-based binomial transform of strict partitions, differences A129519.
Cf. A000041, A007442, A002865, A008284, A095195, A103446, A116608, A218482, A281425, A320568.

Programs

  • Magma
    m:=50; R:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[(1 + x^k/(1 + x)^k): k in [1..(m+2)]]) )); // G. C. Greubel, Oct 29 2018
  • Maple
    seq(coeff(series(mul((1+x^k/(1+x)^k),k=1..n),x,n+1), x, n), n = 0 .. 37); # Muniru A Asiru, Oct 16 2018
  • Mathematica
    nmax = 37; CoefficientList[Series[Product[(1 + x^k/(1 + x)^k), {k, 1, nmax}], {x, 0, nmax}], x]
nmax = 37; CoefficientList[Series[Exp[Sum[Sum[(-1)^(k/d + 1) d, {d, Divisors[k]}] x^k/(k (1 + x)^k), {k, 1, nmax}]], {x, 0, nmax}], x]
Prepend[Table[Differences[PartitionsQ/@Range[0,k+1],k][[2]],{k,0,30}],1] (* Gus Wiseman, Jan 29 2025 *)
  • PARI
    m=50; x='x+O('x^m); Vec(prod(k=1, m+2, (1 + x^k/(1 + x)^k))) \\ G. C. Greubel, Oct 29 2018

Formula

G.f.: exp(Sum_{k>=1} (-1)^(k+1)*x^k/(k*((1 + x)^k - x^k))).
G.f.: exp(Sum_{k>=1} A000593(k)*x^k/(k*(1 + x)^k)).
From Peter Bala, Dec 22 2020: (Start)
O.g.f.: Sum_{n >= 0} x^(n*(n+1)/2)/Product_{k = 1..n} ((1 + x)^k - x^k). Cf. A307548.
Conjectural o.g.f.: (1/2) * Sum_{n >= 0} x^(n*(n-1)/2)*(1 + x)^n/( Product_{k = 1..n} ( (1 + x)^k - x^k ) ). (End)
a(n+1) = Sum_{k=0..n} (-1)^(n-k) binomial(n,k) A000009(k+1). - Gus Wiseman, Feb 03 2025

A320589 Expansion of (1/(1 + x)) * Sum_{k>=1} k*x^k/(x^k + (1 + x)^k).

Original entry on oeis.org

1, -1, 4, -13, 36, -88, 197, -421, 895, -1946, 4346, -9832, 22140, -49043, 106389, -226213, 473366, -980413, 2022418, -4179198, 8687753, -18201140, 38398455, -81343408, 172383461, -364158198, 764854519, -1595107695, 3302884966, -6796646603, 13921482698, -28437025029, 58034908034
Offset: 1

Views

Author

Ilya Gutkovskiy, Oct 16 2018

Keywords

Comments

Inverse binomial transform of A000593.

Links

Crossrefs

Cf. A000593, A320568, A320586, A320591.

Programs

  • Magma
    m:=50; R:=PowerSeriesRing(Rationals(), m); Coefficients(R!( (1/(1 + x))*(&+[k*x^k/(x^k + (1 + x)^k): k in [1..(m+2)]]) )); // G. C. Greubel, Oct 29 2018
  • Maple
    seq(coeff(series((1/(1+x))*add(k*x^k/(x^k+(1+x)^k),k=1..n),x,n+1), x, n), n = 1 .. 35); # Muniru A Asiru, Oct 16 2018
  • Mathematica
    nmax = 33; Rest[CoefficientList[Series[1/(1 + x) Sum[k x^k/(x^k + (1 + x)^k), {k, 1, nmax}], {x, 0, nmax}], x]]
nmax = 33; Rest[CoefficientList[Series[(EllipticTheta[3, 0, x/(1 + x)]^4 + EllipticTheta[2, 0, x/(1 + x)]^4 - 1)/(24 (1 + x)), {x, 0, nmax}], x]]
Table[Sum[(-1)^(n - k) Binomial[n, k] Sum[(-1)^(k/d + 1) d, {d, Divisors[k]}], {k, n}], {n, 33}]
  • PARI
    m=50; x='x+O('x^m); Vec((1/(1 + x))*sum(k=1, m+2, k*x^k/(x^k + (1 + x)^k))) \\ G. C. Greubel, Oct 29 2018

Formula

G.f.: (theta_3(x/(1 + x))^4 + theta_2(x/(1 + x))^4 - 1)/(24*(1 + x)), where theta_() is the Jacobi theta function.
L.g.f.: Sum_{k>=1} A000593(k)*x^k/(k*(1 + x)^k) = Sum_{n>=1} a(n)*x^n/n.
a(n) = Sum_{k=1..n} (-1)^(n-k)*binomial(n,k)*A000593(k).
Conjecture: a(n) ~ -(-1)^n * c * 2^n * n, where c = Pi^2/48 = 0.205616758356... - Vaclav Kotesovec, Jun 26 2019
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

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