%I #14 Oct 24 2022 07:49:11
%S 1,7,9,5,30941,63,5229043,32,27,11,23,45,53,29,4651,64,103,19,
%T 12865927,25,43,161,1381,288,701,371,81,145,1973,31,311,128,207,721,
%U 211,37,1481,90061489,79,41,6740847065723,261,119627,115,181,47,183959,576,1667,101
%N Minimal order of degree-n irreducible polynomials over GF(13).
%C a(n) < 13^n.
%H Max Alekseyev, <a href="/A218360/b218360.txt">Table of n, a(n) for n = 1..416</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/IrreduciblePolynomial.html">Irreducible Polynomial</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PolynomialOrder.html">Polynomial Order</a>
%F a(n) = min(M(n)) with M(n) = {d : d|(13^n-1)} \ U(n-1) and U(n) = M(n) union U(n-1) for n>0, U(0) = {}.
%F a(n) = A218337(n,1) = A213224(n,6).
%p with(numtheory):
%p M:= proc(n) M(n):= divisors(13^n-1) minus U(n-1) end:
%p U:= proc(n) U(n):= `if`(n=0, {}, M(n) union U(n-1)) end:
%p a:= n-> min(M(n)[]):
%p seq(a(n), n=1..33);
%t M[n_] := M[n] = Divisors[13^n - 1]~Complement~U[n - 1];
%t U[n_] := U[n] = If[n == 0, {}, M[n]~Union~U[n - 1]];
%t a[n_] := Min[M[n]];
%t Table[a[n], {n, 1, 50}] (* _Jean-François Alcover_, Oct 24 2022, after _Alois P. Heinz_ *)
%Y Cf. A213224, A218337.
%K nonn
%O 1,2
%A _Alois P. Heinz_, Oct 27 2012
|