OFFSET
0,10
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..100
MAPLE
f:= proc(n) local i; max(coeffs(expand(mul(1+x^(i^4), i=1..n)))) end proc:
map(f, [$1..50]); # Robert Israel, Dec 26 2022
PROG
(PARI) a(n) = vecmax(Vec(prod(k=1, n, 1+x^(k^4)))); \\ Michel Marcus, Dec 26 2022
(Python)
from collections import Counter
def A359320(n):
c = {0:1, 1:1}
for i in range(2, n+1):
j, d = i**4, Counter(c)
for k in c:
d[k+j] += c[k]
c = d
return max(c.values()) # Chai Wah Wu, Jan 31 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 25 2022
EXTENSIONS
a(38)-a(50) from Seiichi Manyama, Dec 26 2022
STATUS
approved