|
%I #12 May 17 2017 22:53:12
%S 1,23,416,7344,134613,2620379,54636792,1223392968,29409134545,
%T 757686550455,20870680635528,612964613117960,19140704352872949,
%U 633710701752022635,22185391759982205904,819180275431111135536,31826528430233802890049,1298154677953792936043447,55473817874169725876166480
%N Column 2 of A286781.
%H Gheorghe Coserea, <a href="/A286787/b286787.txt">Table of n, a(n) for n = 2..211</a>
%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 Kol(K,N=20) = {
%o my(s = A286781_ser(N+K, 't+O('t^(K+1))));
%o vector(N, n, polcoeff(polcoeff(s, K+n-1), K));
%o };
%o Kol(2)
%Y Cf. A286781.
%K nonn
%O 2,2
%A _Gheorghe Coserea_, May 16 2017
|
