%I #13 Feb 07 2024 11:51:14
%S 1,3,12,62,332,1974,12345,80006,531524,3602358,24836850,173607568,
%T 1226700784,8748861828,62922343566,455805857978,3321800235936,
%U 24338840717799,179217603427200,1325490660318216,9841000101286172,73319407735938570,548051770664957631,4108826483323392880
%N Maximal coefficient of (1 + x)^3 * (1 + x^2)^3 * (1 + x^3)^3 * ... * (1 + x^n)^3.
%t Table[Max[CoefficientList[Product[(1 + x^k)^3, {k, 1, n}], x]], {n, 0, 23}]
%o (PARI) a(n) = vecmax(Vec(prod(k=1, n, (1+x^k)^3))); \\ _Michel Marcus_, Jan 30 2024
%o (Python)
%o from collections import Counter
%o def A369709(n):
%o c = {0:1}
%o for k in range(1,n+1):
%o d = Counter(c)
%o for j in c:
%o a = c[j]
%o d[j+k] += 3*a
%o d[j+2*k] += 3*a
%o d[j+3*k] += a
%o c = d
%o return max(c.values()) # _Chai Wah Wu_, Feb 07 2024
%Y Cf. A022568, A025591, A047653.
%K nonn
%O 0,2
%A _Ilya Gutkovskiy_, Jan 29 2024