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.
%I A128577 #10 Mar 19 2016 06:31:35
%S A128577 1,2,9,64,624,7736,116416,2060808,41952600,965497440,24786054816,
%T A128577 702201877920,21761251764672,732269872931712,26589359234860560,
%U A128577 1036241806935453696,43142510740036313088,1911022260200150482944,89737455913330610995200,4452805047268938247981056,232806644343118618035904512,12791828071344703747110764544,736928909474399720669590216704,44416721474748725213260027514880
%N A128577 Self-convolution of A128318.
%C A128577 A128318 equals row 0 of table A128570.
%H A128577 Paul D. Hanna, <a href="/A128577/b128577.txt">Table of n, a(n) for n = 0..200</a>
%F A128577 G.f.: A(x) = [1 + x*R(x,1)^2]^2, where R(x,1) = 1 + 2*x*R(x,2)^2, R(x,2) = 1 + 3*x*R(x,3)^2, ..., R(x,n) = 1 + (n+1)*x*R(x,n+1)^2, ... and R(x,n) is the g.f. of row n of table A128570.
%F A128577 a(n) ~ 2*A128318(n). - _Vaclav Kotesovec_, Mar 19 2016
%o A128577 (PARI) {a(n)=local(A=1+(n+1)*x);for(j=0,n,A=1+(n+1-j)*x*A^2 +x*O(x^n)); polcoeff(A^2,n)}
%o A128577 for(n=0, 25, print1(a(n), ", "))
%Y A128577 Cf. A128570 (triangle), rows: A128318, A128571, A128572, A128573, A128574, A128575, A128576; A128578 (main diagonal), A128579 (antidiagonal sums).
%K A128577 nonn
%O A128577 0,2
%A A128577 _Paul D. Hanna_, Mar 11 2007