|
%I #12 Dec 10 2015 16:50:47
%S 3,5,12,34,104,292,916,2791,8660,26538,81584,248554
%N Number of (-1,0,1) polynomials of degree-n irreducible over the integers.
%C A (-1,0,1) polynomial is defined as a monic polynomial whose remaining coefficients are either -1, 0, or 1. For each n, there are 3^n polynomials to consider.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/IrreduciblePolynomial.html">Irreducible Polynomial</a>
%e a(2) = 5 because 1+x+x^2, 1+x^2, 1-x+x^2, -1+x+x^2, -1-x+x^2 are irreducible over the integers.
%p F:= proc(n) local T, count, t,x,p;
%p if n::odd then
%p T:= combinat:-cartprod([[-1,0,1]$(n-1),[1]])
%p else
%p T:= combinat:-cartprod([[-1,0,1]$(n-1),[-1,1]])
%p fi;
%p count:= 0;
%p while not T[finished] do
%p t:= T[nextvalue]();
%p p:= x^n + add(t[i]*x^(n-i),i=1..n);
%p if irreduc(p) then count:= count+1 fi;
%p od;
%p if n::odd then 2*count else count fi;
%p end proc:
%p 3, seq(F(n),n=2..11); # _Robert Israel_, Dec 10 2015
%t Irreducible[p_, n_] := Module[{f}, f=FactorList[p, Modulus->n]; Length[f]==1 || Simplify[p-f[[2, 1]]]===0]; Table[xx=x^Range[0, n-1]; cnt=0; Do[p=x^n+xx.(IntegerDigits[i, 3, n]-1); If[Irreducible[p, 0], cnt++ ], {i, 0, 3^n-1}]; cnt, {n, 10}]
%Y Cf. A087481 (irreducible polynomials of the form x^n +- x^(n-1) +- x^(n-2) +- ... +- 1), A087482 (irreducible binary polynomials).
%K nonn,more
%O 1,1
%A _T. D. Noe_, Sep 11 2003
%E a(11) and a(12) from _Robert Israel_, Dec 10 2015
|
