OFFSET
1,6
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = Sum_{k=5..n} C(n,k) * A079260(k).
EXAMPLE
a(13) = [1287,1716,1716,1287,715,286,78,13,1] * [1,0,0,0,0,0,0,0,1] = 1287+1 = 1288.
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)=1 then 1 else 0 fi end: a:= bintrans(f): seq (a(n), n=1..40);
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Sep 06 2008
STATUS
approved