|
%I #9 May 20 2017 21:54:16
%S 3,45,637,9567,156123,2781389,54043365,1141864959,26137086451,
%T 645573913005,17138071687053,487130581207775,14771502665168715,
%U 476244289169954253,16274079501768450421,587711952558035789055,22370368544174736836835,895269564785985107345453,37585737112233443277746589
%N Column 1 of A286782.
%H Gheorghe Coserea, <a href="/A287031/b287031.txt">Table of n, a(n) for n = 2..211</a>
%F a(n) = (2*n-1)*A286786(n-1).
%o (PARI)
%o A286781_ser(N, t='t) = {
%o my(x='x+O('x^N), y0=1+O('x^N), y1=0, n=1);
%o while(n++,
%o y1 = (1 + x*y0 + 2*x^2*y0')*(1 - x*y0*(1-t))/(1-x*y0)^2;
%o if (y1 == y0, break()); y0 = y1; );
%o y0;
%o };
%o A286782_ser(N, t='t) = my(s=A286781_ser(N, t)); 1 + x*s + 2*x^2 * deriv(s,'x);
%o Kol(K, N=20) = {
%o my(v = A286782_ser(N+K+1, 't+O('t^(K+1))));
%o vector(N, n, polcoeff(polcoeff(v, K+n), K));
%o };
%o Kol(1)
%Y Cf. A286782.
%K nonn
%O 2,1
%A _Gheorghe Coserea_, May 18 2017
|
