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.

A116364 Row squared sums of Catalan triangle A033184.

This page as a plain text file.
%I A116364 #18 Jun 21 2025 16:17:39
%S A116364 1,2,9,60,490,4534,45689,489920,5508000,64276492,773029466,9531003552,
%T A116364 119990158054,1537695160070,20009930706137,263883333450760,
%U A116364 3521003563829212,47470845904561648,645960472314074400
%N A116364 Row squared sums of Catalan triangle A033184.
%C A116364 Number of 321-avoiding permutations in which the length of the longest increasing subsequence is n. Example: a(2)=9 because we have 12, 132, 312, 213, 231, 3142, 3412, 2143 and 2413. Column sums of triangle in A126217 (n >= 1). - _Emeric Deutsch_, Sep 07 2007
%H A116364 G. C. Greubel, <a href="/A116364/b116364.txt">Table of n, a(n) for n = 0..830</a>
%F A116364 a(n) = Sum_{k=0..n} (C(2*n-k+1,n-k)*(k+1)/(2*n-k+1))^2.
%e A116364 The dot product of Catalan row 4 with itself equals
%e A116364   a(4) = [14,14,9,4,1]*[14,14,9,4,1] = 490
%e A116364 which is equivalent to obtaining the final term in these repeated partial sums of Catalan row 4:
%e A116364   14,   14,    9,    4,    1
%e A116364      28,   37,   41,   42
%e A116364         65,  106,  148
%e A116364           171,  319
%e A116364              490
%p A116364 a:=proc(k) options operator, arrow: sum((2*k-n+1)^2*binomial(n+1,k+1)^2/(n+1)^2,n=k..2*k) end proc: 1,seq(a(k),k=1..17); # _Emeric Deutsch_, Sep 07 2007
%t A116364 Table[Sum[(Binomial[2*n-j+1, n-j]*(j+1)/(2*n-j+1))^2, {j, 0, n}], {n, 0, 30}] (* _G. C. Greubel_, May 12 2019 *)
%o A116364 (PARI) a(n)=sum(k=0,n,((k+1)*binomial(2*n-k+1,n-k)/(2*n-k+1))^2)
%o A116364 (Magma) [(&+[(Binomial(2*n-j+1, n-j)*(j+1)/(2*n-j+1))^2: j in [0..n]]): n in [0..30]]; // _G. C. Greubel_, May 12 2019
%o A116364 (Sage) [sum(( binomial(2*n-j+1, n-j)*(j+1)/(2*n-j+1) )^2 for j in (0..n)) for n in (0..30)] # _G. C. Greubel_, May 12 2019
%o A116364 (GAP) List([0..30], n-> Sum([0..n], j-> (Binomial(2*n-j+1, n-j)* (j+1)/(2*n-j+1))^2 )); # _G. C. Greubel_, May 12 2019
%Y A116364 Cf. A033184, A116363.
%Y A116364 Cf. A126217.
%K A116364 nonn
%O A116364 0,2
%A A116364 _Paul D. Hanna_, Feb 04 2006