%I #10 Mar 30 2012 18:49:34
%S 0,2,6,17,30,177,417,1857,4302,19457,47731,203777,509769,2134017,
%T 5462701,22347777,58104062,234029057,616919457,2450784257,6533317815,
%U 25664946177,69085604341,268766806017,729558799305,2814562533377
%N 2*A203579 - A204449. Difference between the exponential convolution of A000032 (Lucas) with itself and the corresponding exponential half-convolution.
%C See A204449 for the exponential half-convolution of A000032 (Lucas). The present sequence gives the numbers to be added to A204449 to obtain the corresponding (full) exponential convolution A203579.
%F a(n) = sum(binomial(n,k)*L(k)*L(n-k),k=floor(n/2)+1..n), n>=0,
%F with the Lucas numbers L(n)=A000032(n). For n=0 this is 0.
%F E.g.f.: exp(x)*(cosh(sqrt(5)*x)+1) - (BesselI(0,2*phi*x) + BesselI(0,2*(phi-1)*x) + 2*BesselI(0,2*I*x))/2. Compare this with the e.g.f. of A204449, where phi and BesselI are explained.
%F Bisection: a(2*k) = (2^(2*k)-binomial(2*k,k))*L(2*k)/2 + 1 - ((-1)^k)*binomial(2*k,k), a(2*k+1) = 2^(2*k)*L(2*k+1)+ 1 = A204449(2*k+1), k>=0.
%e With A000032 = {2, 1, 3, 4, 7, 11, ...}
%e a(4) = 4*4*1 + 1*7*2 = 30.
%e a(5) = 10*4*3 + 5*7*1 + 1*11*2 = 177.
%Y Cf. A000032, A203579, A204449.
%K nonn
%O 0,2
%A _Wolfdieter Lang_, Jan 16 2012