|
%I #19 Feb 04 2024 03:29:52
%S 1,4,5,15,13,31,26,57,42,91,66,139,95,209,129,283,171,365,216,463,272,
%T 573,333,697,401,825,468,993,545,1139,629,1315,725,1509,815,1689,920,
%U 1921,1030,2139,1147,2367,1261,2619,1391,2861,1521,3135,1659,3409,1802,3703,1952
%N Number of different coefficient values in expansion of Product_{k=1..n} (1-x^k)^3.
%H Seiichi Manyama, <a href="/A369790/b369790.txt">Table of n, a(n) for n = 0..1000</a>
%o (PARI) a(n) = #Set(Vec(prod(k=1, n, (1-x^k)^3)));
%o (Python)
%o from collections import Counter
%o def A369790(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 len(set(c.values()))+int(max(c)+1>len(c)) # _Chai Wah Wu_, Feb 01 2024
%Y Cf. A010816, A039822.
%K nonn
%O 0,2
%A _Seiichi Manyama_, Feb 01 2024
|
