A288391 -id:A288391 - cp's OEIS frontend

A288392 Expansion of Product_{k>=1} (1 - x^k)^(sigma_3(k)).

1, -1, -9, -19, -9, 163, 573, 1127, 109, -7198, -27159, -58611, -50378, 157532, 892986, 2431694, 4040909, 1605559, -16109148, -68261139, -167737209, -263590908, -109589779, 934422499, 3976197701, 9922490735, 16765911071, 13022553978, -33008232762

Offset: 0

Seiichi Manyama, Jun 08 2017

sign

Seiichi Manyama, Table of n, a(n) for n = 0..1000

Cf. A027848, A288391.

Product_{k>=1} (1 - x^k)^sigma_m(k): A288098 (m=0), A288385 (m=1), A288389 (m=2), this sequence (m=3).

m:=30; R:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[(1-q^k)^DivisorSigma(3,k): k in [1..m]]) )); // G. C. Greubel, Oct 30 2018

with(numtheory):
b:= proc(n) option remember; `if`(n=0, 1, add(add(
      d*sigma[3](d), d=divisors(j))*b(n-j), j=1..n)/n)
    end:
a:= proc(n) option remember; `if`(n=0, 1,
      -add(b(n-i)*a(i), i=0..n-1))
    end:
seq(a(n), n=0..30);  # Alois P. Heinz, Jun 08 2017

nmax = 30; CoefficientList[Series[Product[(1-x^k)^DivisorSigma[3, k], {k, 1, nmax}], {x, 0, nmax}], x] (* G. C. Greubel, Oct 30 2018 *)

m=30; x='x+O('x^m); Vec(prod(k=1, m, (1-x^k)^sigma(k,3))) \\ G. C. Greubel, Oct 30 2018

Convolution inverse of A288391.
a(0) = 1, a(n) = -(1/n)*Sum_{k=1..n} A027848(k)*a(n-k) for n > 0.
G.f.: exp(-Sum_{k>=1} sigma_4(k)*x^k/(k*(1 - x^k))). - Ilya Gutkovskiy, Oct 29 2018

A321876 Square array A(n,k), n >= 0, k >= 0, read by antidiagonals, where column k is the expansion of Product_{j>=1} 1/(1 - x^j)^sigma_k(j).

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 1, 4, 5, 1, 1, 6, 8, 11, 1, 1, 10, 16, 21, 17, 1, 1, 18, 38, 52, 39, 34, 1, 1, 34, 100, 156, 128, 92, 52, 1, 1, 66, 278, 526, 534, 373, 170, 94, 1, 1, 130, 796, 1896, 2546, 2014, 913, 360, 145, 1, 1, 258, 2318, 7102, 13074, 12953, 6796, 2399, 667, 244

Offset: 0

Ilya Gutkovskiy, Nov 20 2018

			Square array begins:
   1,   1,    1,    1,     1,      1,  ...
   1,   1,    1,    1,     1,      1,  ...
   3,   4,    6,   10,    18,     34,  ...
   5,   8,   16,   38,   100,    278,  ...
  11,  21,   52,  156,   526,   1896,  ...
  17,  39,  128,  534,  2546,  13074,  ...

Seiichi Manyama, Antidiagonals n = 0..139, flattened

Columns k=0..9 give A006171, A061256, A275585, A288391, A301542, A301543, A301544, A301545, A301546, A301547.
Main diagonal gives A319647.
Cf. A321877.

Table[Function[k, SeriesCoefficient[Product[1/(1 - x^j)^DivisorSigma[k, j], {j, 1, n}], {x, 0, n}]][i - n], {i, 0, 10}, {n, 0, i}] // Flatten
Table[Function[k, SeriesCoefficient[Exp[Sum[DivisorSigma[k + 1, j] x^j/(j (1 - x^j)), {j, 1, n}]], {x, 0, n}]][i - n], {i, 0, 10}, {n, 0, i}] // Flatten

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

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

cp's OEIS Frontend

A288392 Expansion of Product_{k>=1} (1 - x^k)^(sigma_3(k)).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula