OFFSET
0,2
REFERENCES
Roman E. Maeder, Programming in Mathematica, Addison-Wesley, 1990, page 105.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10
Eric Weisstein's World of Mathematics, Swinnerton-Dyer Polynomial
EXAMPLE
SwinnertonDyerP(3, x) = x^8 - 40*x^6 + 352*x^4 - 960*x^2 + 576, hence a(3) = 960.
MAPLE
p:= proc(n) option remember; expand(`if`(n=0, x, mul(
subs(x=x+i*sqrt(ithprime(n)), p(n-1)), i=[1, -1])))
end:
a:= n-> max(map(abs, [coeffs(p(n))])):
seq(a(n), n=0..8); # Alois P. Heinz, Jun 12 2022
MATHEMATICA
SwinnertonDyerP[0, x_] := x; SwinnertonDyerP[n_, x_] := Module[{sd, srp = Sqrt[Prime[n]] }, sd[y_] = SwinnertonDyerP[n - 1, y]; Expand[sd[x + srp] sd[x - srp]]]; a[n_] := CoefficientList[SwinnertonDyerP[n, x], x^2] // Abs // Max; Table[a[n], {n, 1, 10}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Jean-François Alcover, Feb 27 2015
EXTENSIONS
a(0)=1 prepended by Alois P. Heinz, Jun 12 2022
STATUS
approved