login
G.f. A(x) satisfies: A(x) = A(x^3) / (1 - x - x^2 - x^3).
1

%I #4 Feb 15 2026 22:59:28

%S 1,1,2,5,8,15,30,53,98,186,337,621,1152,2110,3883,7160,13153,24196,

%T 44539,81888,150623,277103,509614,937340,1724155,3171109,5832604,

%U 10728054,19731767,36292425,66752583,122776775,225821783,415351762,763950320,1405123865,2584427099

%N G.f. A(x) satisfies: A(x) = A(x^3) / (1 - x - x^2 - x^3).

%F G.f.: Product_{k>=0} 1 / (1 - x^(3^k) - x^(2*3^k) - x^(3^(k+1))).

%F a(0) = 1; a(n) = Sum_{k=0..floor(n/3)} A000073(n-3*k+2) * a(k).

%t nmax = 36; A[_] = 1; Do[A[x_] = A[x^3]/(1 - x - x^2 - x^3) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]

%t nmax = 36; CoefficientList[Series[Product[1/(1 - x^(3^k) - x^(2 3^k) - x^(3^(k+1))), {k, 0, Floor[Log[3, nmax]] + 1}], {x, 0, nmax}], x]

%Y Cf. A000073, A054390, A062051, A309702, A345007, A393315.

%K nonn

%O 0,3

%A _Ilya Gutkovskiy_, Feb 10 2026