A341265 -id:A341265 - cp's OEIS frontend

A341279 Triangle read by rows: T(n,k) = coefficient of x^n in expansion of (-1 + Product_{j>=1} 1 / (1 + (-x)^j))^k, n >= 0, 0 <= k <= n.

1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 2, 0, 1, 0, 1, 2, 3, 0, 1, 0, 1, 3, 3, 4, 0, 1, 0, 1, 4, 6, 4, 5, 0, 1, 0, 2, 5, 9, 10, 5, 6, 0, 1, 0, 2, 8, 13, 16, 15, 6, 7, 0, 1, 0, 2, 9, 21, 26, 25, 21, 7, 8, 0, 1, 0, 2, 12, 27, 44, 45, 36, 28, 8, 9, 0, 1, 0, 3, 15, 40, 63, 80, 71, 49, 36, 9, 10, 0, 1

Offset: 0

Ilya Gutkovskiy, Feb 08 2021

			Triangle T(n,k) begins:
  1;
  0,  1;
  0,  0,  1;
  0,  1,  0,  1;
  0,  1,  2,  0,  1;
  0,  1,  2,  3,  0,  1;
  0,  1,  3,  3,  4,  0,  1;
  0,  1,  4,  6,  4,  5,  0,  1;
  0,  2,  5,  9, 10,  5,  6,  0,  1;
  0,  2,  8, 13, 16, 15,  6,  7,  0,  1;
  0,  2,  9, 21, 26, 25, 21,  7,  8,  0,  1;
  0,  2, 12, 27, 44, 45, 36, 28,  8,  9,  0,  1;
  ...

Alois P. Heinz, Rows n = 0..200, flattened

Columns k=0-10 give A000007, A000700, A338463, A341241, A341243, A341244, A341245, A341246, A341247, A341251, A341253.
Main diagonal and lower diagonals give A000012, A000004, A001477, A000217, A000290.
Row sums give A307058.
T(2n,n) gives A341265.
Cf. A000009, A047265, A060642, A286352, A308680.

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:
T:= proc(n, k) option remember;
      `if`(k=0, `if`(n=0, 1, 0), `if`(k=1, `if`(n=0, 0, g(n)),
      (q-> add(T(j, q)*T(n-j, k-q), j=0..n))(iquo(k, 2))))
    end:
seq(seq(T(n, k), k=0..n), n=0..12);  # Alois P. Heinz, Feb 09 2021

T[n_, k_] := SeriesCoefficient[(-1 + 2/QPochhammer[-1, -x])^k, {x, 0, n}]; Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten

G.f. of column k: (-1 + Product_{j>=1} (1 + x^(2*j-1)))^k.

Sum_{k=0..n} (-1)^(n-k) * T(n,k) = A000009(n).

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

Original entry on oeis.org

1, -1, 1, -1, -3, 19, -65, 181, -419, 755, -749, -1530, 12255, -47477, 141065, -343526, 660941, -770917, -911369, 9721976, -40135713, 124134772, -313463842, 631382751, -824406065, -492101356, 8192253811, -35948431288, 115087580857, -299576625051, 627027769120, -894734468883

Offset: 0

Ilya Gutkovskiy, Feb 07 2021

sign

Cf. A001482, A001483, A001484, A001485, A001486, A001487, A001488, A008705, A010815, A047654, A047655, A324595, A340987, A341265.

g:= proc(n) option remember; `if`(n=0, 1, add(add(
     -d, 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..31);  # Alois P. Heinz, Feb 07 2021

Table[SeriesCoefficient[(-1 + QPochhammer[x, x])^n, {x, 0, 2 n}], {n, 0, 31}]
A[n_, k_] := A[n, k] = If[n == 0, 1, -k Sum[A[n - j, k] DivisorSigma[1, 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, 31}]

cp's OEIS Frontend

A341279 Triangle read by rows: T(n,k) = coefficient of x^n in expansion of (-1 + Product_{j>=1} 1 / (1 + (-x)^j))^k, n >= 0, 0 <= k <= n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula