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.

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

Original entry on oeis.org

1, 0, 4, 4, 10, 16, 26, 44, 63, 100, 144, 212, 297, 420, 584, 796, 1081, 1452, 1940, 2556, 3355, 4372, 5668, 7288, 9327, 11892, 15076, 19012, 23884, 29904, 37276, 46284, 57276, 70680, 86918, 106528, 130220, 158784, 193054, 234076, 283178, 341824, 411616, 494512, 592933
Offset: 4

Views

Author

Ilya Gutkovskiy, Feb 07 2021

Keywords

Links

Crossrefs

Cf. A000700, A001482, A022599, A112160, A327382, A338463, A341222, A341241, A341244, A341245, A341246, A341247, A341251.
Column k=4 of A341279.

Programs

  • Maple
    g:= proc(n) option remember; `if`(n=0, 1, add(add([0, d, -d, d]
      [1+irem(d, 4)], d=numtheory[divisors](j))*g(n-j), j=1..n)/n)
    end:
b:= proc(n, k) option remember; `if`(k<2, `if`(n=0, 1-k, g(n)),
      (q-> add(b(j, q)*b(n-j, k-q), j=0..n))(iquo(k, 2)))
    end:
a:= n-> b(n, 4):
seq(a(n), n=4..48);  # Alois P. Heinz, Feb 07 2021
  • Mathematica
    nmax = 48; CoefficientList[Series[(-1 + Product[1/(1 + (-x)^k), {k, 1, nmax}])^4, {x, 0, nmax}], x] // Drop[#, 4] &

Formula

G.f.: (-1 + Product_{k>=1} (1 + x^(2*k - 1)))^4.
a(n) ~ A112160(n). - Vaclav Kotesovec, Feb 20 2021

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

Original entry on oeis.org

1, 1, 1, 2, 4, 7, 12, 21, 38, 68, 120, 212, 377, 670, 1188, 2107, 3740, 6638, 11778, 20898, 37084, 65808, 116775, 207212, 367696, 652478, 1157815, 2054524, 3645730, 6469316, 11479734, 20370656, 36147506, 64143372, 113821732, 201975429, 358403220, 635982680, 1128544452, 2002589998
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 21 2019

Keywords

Comments

Invert transform of A000700.

Links

Crossrefs

Cf. A000700, A055887, A302017, A304969, A307059.
Row sums of A341279.

Programs

  • Magma
    m:=80;
R:=PowerSeriesRing(Integers(), m);
Coefficients(R!( 1/(2 - (&*[1 + x^(2*j-1): j in [1..m+2]])) )); // G. C. Greubel, Jan 24 2024
  • Maple
    g:= proc(n) option remember; `if`(n=0, 1, add(add([0, d, -d, d]
      [1+irem(d, 4)], d=numtheory[divisors](j))*g(n-j), j=1..n)/n)
    end:
a:= proc(n) option remember; `if`(n=0, 1,
      add(a(n-i)*g(i), i=1..n))
    end:
seq(a(n), n=0..39);  # Alois P. Heinz, Feb 09 2021
  • Mathematica
    nmax = 39; CoefficientList[Series[1/(2 - Product[(1 + x^(2 k - 1)), {k, 1, nmax}]), {x, 0, nmax}], x]
  • SageMath
    m=80;
def f(x): return 1/(2 - product(1+x^(2*j-1) for j in range(1,m+3)))
def A307058_list(prec):
    P. = PowerSeriesRing(QQ, prec)
    return P( f(x) ).list()
A307058_list(m) # G. C. Greubel, Jan 24 2024

Formula

a(0) = 1; a(n) = Sum_{k=1..n} A000700(k)*a(n-k).
From G. C. Greubel, Jan 24 2024: (Start)
G.f.: (1+x)/(2*(1+x) - x*QPochhammer(-1/x; x^2)).
G.f.: 1/( 2 - x^(1/24)*etx(x^2)^2/(eta(x^4)*eta(x)) ), where eta(x) is the Dedekind eta function. (End)

A341265 Coefficient of x^(2*n) in (-1 + Product_{k>=1} 1 / (1 + x^k))^n.

Original entry on oeis.org

1, 0, 2, 3, 10, 25, 71, 203, 562, 1650, 4667, 13673, 39427, 115440, 336639, 987628, 2898658, 8529257, 25134200, 74173606, 219207815, 648546314, 1921045953, 5695642513, 16902924883, 50203798050, 149229323544, 443895849894, 1321292939459, 3935377071154, 11728037768186
Offset: 0

Views

Author

Ilya Gutkovskiy, Feb 07 2021

Keywords

Links

Crossrefs

Cf. A000700, A081362, A255526, A324595, A338463, A340987, A341241, A341243, A341244, A341245, A341246, A341247, A341251, A341253, A341263, A341279.

Programs

  • Maple
    g:= proc(n) option remember; `if`(n=0, 1, add(add([0, d, -d, d]
      [1+irem(d, 4)], d=numtheory[divisors](j))*g(n-j), j=1..n)/n)
    end:
b:= proc(n, k) option remember; `if`(k=0, 1, `if`(k=1, g(n+1),
      (q-> add(b(j, q)*b(n-j, k-q), j=0..n))(iquo(k, 2))))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..30);  # Alois P. Heinz, Feb 07 2021
  • Mathematica
    Table[SeriesCoefficient[(-1 + 1/QPochhammer[-x, x])^n, {x, 0, 2 n}], {n, 0, 30}]
A[n_, k_] := A[n, k] = If[n == 0, 1, -k Sum[A[n - j, k] Sum[Mod[d, 2] d, {d, Divisors[j]}], {j, 1, n}]/n]; T[n_, k_] := Sum[(-1)^i Binomial[k, i] A[n, k - i], {i, 0, k}]; Table[T[2 n, n], {n, 0, 30}]

Formula

a(n) = A341279(2n,n).
a(n) ~ c * d^n / sqrt(n), where d = 3.03044218957412050685579849718626198523346... and c = 0.2319377657497495246637662111041144... - Vaclav Kotesovec, Feb 20 2021
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.