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.

A125275 Eigensequence of triangle A039599: a(n) = Sum_{k=0..n-1} A039599(n-1,k)*a(k) for n > 0 with a(0) = 1.

This page as a plain text file.
%I A125275 #23 Jun 05 2021 01:37:48
%S A125275 1,1,2,7,31,162,968,6481,47893,386098,3364562,31460324,313743665,
%T A125275 3320211313,37124987124,436985496790,5397178181290,69748452377058,
%U A125275 940762812167126,13213888481979449,192891251215160017
%N A125275 Eigensequence of triangle A039599: a(n) = Sum_{k=0..n-1} A039599(n-1,k)*a(k) for n > 0 with a(0) = 1.
%C A125275 Starting with offset 1, these are the row sums of triangle A147294. - _Gary W. Adamson_, Nov 05 2008
%H A125275 Vaclav Kotesovec, <a href="/A125275/b125275.txt">Table of n, a(n) for n = 0..500</a>
%H A125275 Guo-Niu Han, <a href="/A196265/a196265.pdf">Enumeration of Standard Puzzles</a>, 2011. [Cached copy]
%H A125275 Guo-Niu Han, <a href="https://arxiv.org/abs/2006.14070">Enumeration of Standard Puzzles</a>, arXiv:2006.14070 [math.CO], 2020.
%F A125275 a(n) = Sum_{k=0..n-1} a(k) * C(2*n-1, n-k-1) * (2*k + 1)/(2*n - 1) for n > 0 with a(0) = 1.
%e A125275 a(3) = 2*(1) + 3*(1) + 1*(2) = 7;
%e A125275 a(4) = 5*(1) + 9*(1) + 5*(2) + 1*(7) = 31;
%e A125275 a(5) = 14*(1) + 28*(1) + 20*(2) + 7*(7) + 1*(31) = 162.
%e A125275 Triangle A039599(n,k) = C(2*n+1, n-k)*(2*k+1)/(2*n+1) (with rows n >= 0 and columns k = 0..n) begins:
%e A125275    1;
%e A125275    1,  1;
%e A125275    2,  3,  1;
%e A125275    5,  9,  5,  1;
%e A125275   14, 28, 20,  7, 1;
%e A125275   42, 90, 75, 35, 9, 1;
%e A125275   ...
%e A125275 where the g.f. of column k is G000108(x)^(2*k+1)
%e A125275 and G000108(x) = (1 - sqrt(1 - 4*x))/(2*x) is the Catalan g.f. function.
%t A125275 A125275=ConstantArray[0,20]; A125275[[1]]=1; Do[A125275[[n]]=Binomial[2*n-1,n-1]/(2*n-1)+Sum[A125275[[k]]*Binomial[2*n-1,n-k-1]*(2*k+1)/(2*n-1),{k,1,n-1}];,{n,2,20}]; Flatten[{1,A125275}] (* _Vaclav Kotesovec_, Dec 09 2013 *)
%o A125275 (PARI) a(n)=if(n==0,1,sum(k=0,n-1, a(k)*binomial(2*n-1, n-k-1)*(2*k+1)/(2*n-1)))
%Y A125275 Cf. A000108, A039599, A125276 (variant), A147294.
%K A125275 nonn
%O A125275 0,3
%A A125275 _Paul D. Hanna_, Nov 26 2006