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 A102087 #7 Dec 15 2014 05:17:54 %S A102087 0,2,4,20,156,1664,22684,378572,7504640,172785512,4540756252, %T A102087 134330010172,4423176368332,160596613105384,6378859853838480, %U A102087 275308217428662672,12836003750434047344,643227594173121801096 %N A102087 Column 1 of triangular matrix A102086, which shifts upward to exclude the main diagonal under matrix square. %t A102087 T[n_, n_] := n+1; T[n_, k_] /; k>n = 0; T[n_, k_] /; k == n-1 := n^2; T[n_, k_] := T[n, k] = Coefficient[1-Sum[T[i, k]*x^i*Product[1-(j+k)*x, {j, 1, i-k+1}], {i, k, n-1}], x, n]; a[n_] := T[n, 1]; Table[a[n], {n, 0, 17} ] (* _Jean-François Alcover_, Dec 15 2014 *) %o A102087 (PARI) {a(n)=local(A=matrix(2,2),B);A[1,1]=1; for(m=2,n+1,B=matrix(m,m);for(i=1,m, for(j=1,i, if(j==i,B[i,j]=j,if(j==1,B[i,j]=(A^2)[i-1,1], B[i,j]=(A^2)[i-1,j]));));A=B); return(A[n+1,2])} %Y A102087 Cf. A102086, A102088. %K A102087 nonn %O A102087 0,2 %A A102087 _Paul D. Hanna_, Dec 29 2004