A387017 -id:A387017 - cp's OEIS frontend

A329156 Expansion of Product_{k>=1} 1 / (1 - Sum_{j>=1} j * x^(k*j)).

1, 1, 4, 10, 29, 72, 200, 510, 1364, 3546, 9348, 24400, 64090, 167562, 439200, 1149360, 3010349, 7879832, 20633304, 54014950, 141422328, 370239300, 969323000, 2537696160, 6643839400, 17393731933, 45537549048, 119218684970, 312119004990, 817137724392, 2139295489200, 5600747143950

Offset: 0

Ilya Gutkovskiy, Nov 06 2019

nonn

Euler transform of A032198.

Alois P. Heinz, Table of n, a(n) for n = 0..2392
Christian Kassel and Christophe Reutenauer, Pairs of intertwined integer sequences, arXiv:2507.15780 [math.NT], 2025. See p. 13.

Cf. A006951, A032198, A088305, A162891, A258210, A329157, A329163, A386706, A387017.

Convolution inverse of A329157.

b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i>1, b(n, i-1), 0)-
      add(b(n-i*j, min(n-i*j, i-1))*j, j=`if`(i=1, n, 1..n/i)))
    end:
a:= proc(n) option remember; `if`(n=0, 1,
      -add(a(j)*b(n-j$2), j=0..n-1))
    end:
seq(a(n), n=0..31);  # Alois P. Heinz, Jul 25 2025

nmax = 31; CoefficientList[Series[Product[1/(1 - Sum[j x^(k j), {j, 1, nmax}]), {k, 1, nmax}], {x, 0, nmax}], x]
nmax = 31; CoefficientList[Series[Product[1/(1 - x^k/(1 - x^k)^2), {k, 1, nmax}], {x, 0, nmax}], x]

G.f.: Product_{k>=1} 1 / (1 - x^k / (1 - x^k)^2).
G.f.: exp(Sum_{k>=1} ( Sum_{d|k} 1 / (d * (1 - x^(k/d))^(2*d)) ) * x^k).
G.f.: Product_{k>=1} 1 / (1 - x^k)^A032198(k).
G.f.: A(x) = Product_{k>=1} B(x^k), where B(x) = g.f. of A088305.
a(n) ~ phi^(2*n-1), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Nov 07 2019
a(2^k) = A002878(2^k-1) for all nonnegative integers k. Follows from Cor. 4.5 on page 11 of Kassel-Reutenauer paper. - Michael De Vlieger, Jul 28 2025

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

Original entry on oeis.org

0, 1, 5, 18, 71, 260, 990, 3672, 13775, 51343, 191860, 715770, 2672298, 9972092, 37220040, 138903480, 518408351, 1934712530, 7220497115, 26947209762, 100568547820, 375326739216, 1400739172470, 5227629044040, 19509779871450, 72811487038701, 271736178975820, 1014133216234068

Offset: 0

Christian Kassel, Jul 30 2025

nonn

a(n) is the value at q = 2 + sqrt(3) of C_n(q)/(q^{n-1}(q - 1)^2), where C_n(q) is the number of codimension n ideals of the algebra of two-variable Laurent polynomials over a finite field of order q. The number C_n(q) is a palindromic polynomial of degree 2n with integer coefficients in the variable q and it is divisible by (q-1)^2.

Vaclav Kotesovec, Table of n, a(n) for n = 0..1000
Christian Kassel and Christophe Reutenauer, Counting the ideals of given codimension of the algebra of Laurent polynomials in two variables, arXiv:1505.07229 [math.AG], 2015-2016; Michigan Math. J. 67 (2018), 715-741.
Christian Kassel and Christophe Reutenauer, Complete determination of the zeta function of the Hilbert scheme of n points on a two-dimensional torus, arXiv:1610.07793 [math.NT], 2016; The Ramanujan Journal 46 (2018), 633-655.
Christian Kassel and Christophe Reutenauer, Pairs of intertwined integer sequences, arXiv:2507.15780 [math.NT], 2025.

Cf. A001834, A329156, A387017.

nmax = 30; CoefficientList[Series[(Product[(1 - x^k)^2/(1 - 4*x^k + x^(2*k)), {k, 1, nmax}] - 1)/2, {x, 0, nmax}], x] (* Vaclav Kotesovec, Jul 30 2025 *)

G.f.: ((Product_{k>=1} (1 - x^k)^2/(1 - 4*x^k + x^(2k))) - 1)/2.
a(2^k) = A001834(2^k-1) for all nonnegative integers k. Follows from Cor. 4.5 of Kassel-Reutenauer paper "Pairs of intertwined integer sequences".
a(n) ~ (1 + sqrt(3))^(2*n-1) / 2^n. - Vaclav Kotesovec, Jul 30 2025

a(0)=0 added, offset changed to 0, a(7) corrected and more terms added by Vaclav Kotesovec, Jul 30 2025

cp's OEIS Frontend

A329156 Expansion of Product_{k>=1} 1 / (1 - Sum_{j>=1} j * x^(k*j)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula