OFFSET
1,5
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
K. Cattell, C. R. Miers, F. Ruskey, J. Sawada and M. Serra, The Number of Irreducible Polynomials over GF(2) with Given Trace and Subtrace, J. Comb. Math. and Comb. Comp., 47 (2003) 31-64.
FORMULA
a(n) = (1/n) * Sum_{ L(n, k) : n+k = 3 mod 4}, where L(n, k) = Sum_{ mu(d)*binomial(n/d, k/d) : d|gcd(n, k)}.
MATHEMATICA
L[n_, k_] := Sum[ MoebiusMu[d]*Binomial[n/d, k/d], {d, Divisors[GCD[n, k]]}]/n; a[n_] := Sum[ If[ Mod[n+k, 4] == 3, L[n, k], 0], {k, 0, n}]; Table[a[n], {n, 1, 32}] (* Jean-François Alcover, Jun 28 2012, from formula *)
PROG
(PARI)
L(n, k) = sumdiv(gcd(n, k), d, moebius(d) * binomial(n/d, k/d) );
a(n) = sum(k=0, n, if( (n+k)%4==3, L(n, k), 0 ) ) / n;
vector(33, n, a(n))
/* Joerg Arndt, Jun 28 2012 */
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved