OFFSET
1,4
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
N. J. A. Sloane, Transforms
FORMULA
a(n) = Sum_{k=3..n} C(n,k)*A079261(k).
EXAMPLE
a(11) = [165,330,462,462,330,165,55,11,1] * [1,0,0,0,1,0,0,0,1] = 165+330+1 = 496.
MAPLE
bintrans:= proc(p) proc(n) add(p(k) *binomial(n, k), k=0..n) end end:
f:= proc(n) if isprime(n) and modp(n, 4)=3 then 1 else 0 fi end:
a:= bintrans(f):
seq(a(n), n=1..40);
MATHEMATICA
a[n_] := Sum[Binomial[n, k] Boole[PrimeQ[k] && Mod[k, 4] == 3], {k, 3, n}];
Array[a, 40] (* Jean-François Alcover, May 27 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Sep 06 2008
STATUS
approved