OFFSET
0,3
MATHEMATICA
Table[Max[CoefficientList[Product[(1 + Sum[x^Prime[j], {j, 1, i}]), {i, 1, n}], x]], {n, 0, 22}]
PROG
(PARI) a(n) = vecmax(Vec(prod(k=1, n, 1 + sum(i=1, k, x^prime(i))))); \\ Michel Marcus, Feb 01 2024
(Python)
from collections import Counter
from sympy import prime, primerange
def A369775(n):
if n == 0: return 1
c, p = {0:1}, list(primerange(prime(n)+1))
for k in range(1, n+1):
d = Counter(c)
for j in c:
a = c[j]
for i in p[:k]:
d[j+i] += a
c = d
return max(c.values()) # Chai Wah Wu, Feb 01 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 31 2024
STATUS
approved