A104312 Prime coefficient of x^n in (x^3+x^2+x+1)^n for n in A104311.

3, 31, 101, 3823, 2266366724843687556556015073508073201681

Offset: 1

Author

T. D. Noe, Mar 01 2005

Keywords

• nonn

Comments

a(6), which corresponds to n=649, is too large to include.

Crossrefs

Cf. A005725 (quadrinomial coefficients), A104314 (nontrivial prime pentanomial coefficients).

Programs

Mathematica

f=1; Do[f=Expand[f*(x^3+x^2+x+1)]; s=Coefficient[f, x, n]; If[PrimeQ[s], Print[{n, s}]], {n, 1000}]