cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A287041 Column 1 of A286783.

This page as a plain text file.
%I A287041 #8 May 20 2017 21:55:19
%S A287041 5,77,1044,14784,227877,3862305,71983440,1469813400,32718512925,
%T A287041 789901955325,20578796752500,575836554270600,17232413940017325,
%U A287041 549370878062313825,18591830334684129600,665771181527890746000,25154357611638416671125,1000094581801108086418125,41741166856778766269392500
%N A287041 Column 1 of A286783.
%H A287041 Gheorghe Coserea, <a href="/A287041/b287041.txt">Table of n, a(n) for n = 2..211</a>
%o A287041 (PARI)
%o A287041 A286781_ser(N, t='t) = {
%o A287041   my(x='x+O('x^N), y0=1+O('x^N), y1=0, n=1);
%o A287041   while(n++,
%o A287041     y1 = (1 + x*y0 + 2*x^2*y0')*(1 - x*y0*(1-t))/(1-x*y0)^2;
%o A287041     if (y1 == y0, break()); y0 = y1;);
%o A287041   y0;
%o A287041 };
%o A287041 Kol(K, N=20) = {
%o A287041   my(s = A286781_ser(N+K+1, 't+O('t^(K+1))),
%o A287041      p = (1 + x*s + 2*x^2*s')/(1-x*s)^2);
%o A287041   vector(N, n, polcoeff(polcoeff(p, K+n), K));
%o A287041 };
%o A287041 Kol(1)
%Y A287041 Cf. A286783.
%K A287041 nonn
%O A287041 2,1
%A A287041 _Gheorghe Coserea_, May 19 2017