OFFSET
0,11
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..500
MAPLE
b:= proc(n) option remember; `if`(n=0, 1,
expand((1-x^ithprime(n))*b(n-1)))
end:
a:= n-> max(coeffs(b(n))):
seq(a(n), n=0..60);
MATHEMATICA
a[n_] := Times@@(1-x^Prime[Range[n]])//Expand//CoefficientList[#, x]&//Max;
Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Jun 02 2022 *)
PROG
(Python)
from sympy.abc import x
from sympy import prime, prod
def A350514(n): return 1 if n == 0 else max(prod(1-x**prime(i) for i in range(1, n+1)).as_poly().coeffs()) # Chai Wah Wu, Jan 04 2022
(PARI) a(n) = vecmax(Vec(prod(j=1, n, 1-'x^prime(j)))); \\ Michel Marcus, Jan 04 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jan 02 2022
STATUS
approved