OFFSET
1,1
COMMENTS
LINKS
Joerg Arndt, Mar 31 2008, Table of n, a(n) for n = 1..201
Joerg Arndt, Matters Computational (The Fxtbook), section 42.9 "Gaussian normal bases", pp.914-920
EXAMPLE
12 is in the list because 3*12+1=37 is prime and the index of 2 mod 37 (==36/ord(2 mod 37)==1, 2 is a generator mod 37) is coprime to 12.
PROG
(PARI)
gauss_test(n, t)=
{ /* test whether a type-t Gaussian normal basis exists for GF(2^n) */
local( p, r2, g, d );
p = t*n + 1;
if ( !isprime(p), return( 0 ) );
if ( p<=2, return( 0 ) );
r2 = znorder( Mod(2, p) );
d = (p-1)/r2;
g = gcd(d, n);
return ( if ( 1==g, 1, 0) );
}
/* generate this sequence: */
t=3; ct=1; for(n=1, 10^7, if(gauss_test(n, t), print1(n, ", "); ct+=1; if(ct>200, break())))
CROSSREFS
KEYWORD
nonn
AUTHOR
Joerg Arndt, Mar 31 2008
STATUS
approved