|
%I #5 Oct 17 2012 20:31:35
%S 0,0,0,0,1,6,21,56,126,252,462,792,1288,2016,3108,4928,8569,17154,
%T 38931,94164,229824,550088,1278662,2884752,6335005,13590930,28575315,
%U 59014620,119878606,239662236,471605976,913302656,1740247806,3262146492,6015853242
%N Binomial transform of A079260.
%H Alois P. Heinz, <a href="/A143980/b143980.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) = Sum_{k=5..n} C(n,k) * A079260(k).
%e a(13) = [1287,1716,1716,1287,715,286,78,13,1] * [1,0,0,0,0,0,0,0,1] = 1287+1 = 1288.
%p 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);
%Y Cf. A007318, A079260.
%K nonn
%O 1,6
%A _Alois P. Heinz_, Sep 06 2008
|
