A137636
a(n) = Sum_{k=0..n} C(2k+1,k)*C(2k+1,n-k) ; equals row 1 of square array A137634; also equals the convolution of A137635 and A073157.
Original entry on oeis.org
1, 4, 19, 94, 474, 2431, 12609, 65972, 347524, 1840680, 9792986, 52296799, 280163091, 1504969409, 8103433329, 43722788132, 236340999038, 1279602656590, 6938126362948, 37668424608552, 204751452911832, 1114151447523038
Offset: 0
-
{a(n)=sum(k=0,n,binomial(2*k+1,k)*binomial(2*k+1,n-k))} /* Using the g.f.: */ {a(n)=local(R=1/sqrt(1-4*x*(1+x +x*O(x^n))^2), G=(1-sqrt(1-4*x*(1+x)^2+x^2*O(x^n)))/(2*x*(1+x+x*O(x^n)))); polcoeff(R*G,n,x)}
A137637
a(n) = Sum_{k=0..n} C(2k+2,k)*C(2k+2,n-k) ; equals row 2 of square array A137634 ; also equals the convolution of A137635 and the self-convolution of A073157.
Original entry on oeis.org
1, 6, 32, 170, 899, 4764, 25318, 134964, 721562, 3868024, 20785035, 111931154, 603938905, 3264309644, 17671408012, 95800342628, 520022296700, 2826089180652, 15374990077568, 83727902852188, 456370687687082
Offset: 0
-
{a(n)=sum(k=0,n,binomial(2*k+2,k)*binomial(2*k+2,n-k))} /* Using the g.f.: */ {a(n)=local(R=1/sqrt(1-4*x*(1+x +x*O(x^n))^2), G=(1-sqrt(1-4*x*(1+x)^2+x^2*O(x^n)))/(2*x*(1+x+x*O(x^n)))); polcoeff(R*G^2,n,x)}
A073154
Triangle of numbers relating two sequences (A073157 and A073155).
Original entry on oeis.org
1, 2, 4, 5, 9, 14, 18, 28, 38, 56, 70, 106, 131, 167, 237, 293, 433, 523, 613, 753, 1046, 1283, 1869, 2219, 2543, 2893, 3479, 4762, 5808, 8374, 9839, 11099, 12359, 13824, 16390, 22198
Offset: 0
a(4,0)=a(3,3)+a(2,2)=56+14=70.
a(5,2)=A073157(0)*A073157(5)+A073157(1)*A073157(4)+A073157(2)*A073157(3)= 1*293+2*70+5*18=523.
Rows:
{1};
{2,4};
{5,9,14};
{18,28,38,56};
{70,106,131,167,237};
{293,433,523,613,753,1046};
{1283,1869,2219,2543,2893,3479,4762};
...
A198953
G.f. satisfies A(x) = (1 + x*A(x)) * (1 + x*A(x)^3).
Original entry on oeis.org
1, 2, 9, 56, 400, 3095, 25240, 213633, 1859006, 16527544, 149472480, 1370794835, 12718060947, 119158146283, 1125816405458, 10714275588727, 102615375322564, 988302823695146, 9565859385140272, 93000625498797314, 907782305262566776, 8892941663606408172
Offset: 0
G.f.: A(x) = 1 + 2*x + 9*x^2 + 56*x^3 + 400*x^4 + 3095*x^5 + 25240*x^6 +...
Related expansions.
A(x)^2 = 1 + 4*x + 22*x^2 + 148*x^3 + 1105*x^4 + 8798*x^5 + 73196*x^6 +...
A(x)^3 = 1 + 6*x + 39*x^2 + 284*x^3 + 2223*x^4 + 18267*x^5 + 155445*x^6 +...
A(x)^4 = 1 + 8*x + 60*x^2 + 472*x^3 + 3878*x^4 + 32948*x^5 + 287300*x^6 +...
where A(x) = 1 + x*(A(x) + A(x)^3) + x^2*A(x)^4.
The logarithm of the g.f. equals the series:
log(A(x)) = (1 + A(x)^2)*x + (1 + 2^2*A(x)^2 + A(x)^4)*x^2/2 +
(1 + 3^2*A(x)^2 + 3^2*A(x)^4 + A(x)^6)*x^3/3 +
(1 + 4^2*A(x)^2 + 6^2*A(x)^4 + 4^2*A(x)^6 + A(x)^8)*x^4/4 +
(1 + 5^2*A(x)^2 + 10^2*A(x)^4 + 10^2*A(x)^6 + 5^2*A(x)^8 + A(x)^10)*x^5/5 +...
more explicitly,
log(A(x)) = 2*x + 14*x^2/2 + 122*x^3/3 + 1118*x^4/4 + 10557*x^5/5 + 101642*x^6/6 + 991916*x^7/7 +...
-
nmax=20; aa=ConstantArray[0,nmax]; aa[[1]]=2; Do[AGF=1+Sum[aa[[n]]*x^n,{n,1,j-1}]+koef*x^j; sol=Solve[Coefficient[(1 + x*AGF)*(1 + x*AGF^3) - AGF,x,j]==0,koef][[1]];aa[[j]]=koef/.sol[[1]],{j,2,nmax}]; Flatten[{1,aa}] (* Vaclav Kotesovec, Sep 19 2013 *)
-
a(n):=sum((sum((binomial(2*n+2*k+2,j-k)*binomial(n+2*k,k))/(k+n+1),k,0,j))*(-1)^(n-j)*binomial(2*n-j,n-j),j,0,n); /* Vladimir Kruchinin, Mar 13 2016 */
-
{a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*(A+x*O(x^n))^(2*j))*x^m/m))); polcoeff(A, n)}
-
{a(n)=polcoeff((1/x)*serreverse( 2*x^2*(1+x) / (1 - sqrt(1 - 4*x*(1+x)^2 +x^3*O(x^n)))),n)}
-
{a(n)=local(A=1+x);for(i=1,n,A=(1 + x*A)*(1 + x*(A+x*O(x^n))^3));polcoeff(A,n)}
A137635
a(n) = Sum_{k=0..n} C(2k,k)*C(2k,n-k); equals row 0 of square array A137634.
Original entry on oeis.org
1, 2, 10, 46, 226, 1136, 5810, 30080, 157162, 826992, 4376408, 23267332, 124179570, 664919780, 3570265000, 19216805476, 103652442922, 560127574340, 3031887311256, 16435458039076, 89213101943000, 484839755040768, 2637805800869740, 14365506336197816
Offset: 0
-
CoefficientList[Series[1/Sqrt[1 - 4*x*(1 + x)^2],{x,0,50}],x] (* Stefano Spezia, Sep 01 2018 *)
Table[Sum[Binomial[2k,k]Binomial[2k,n-k],{k,0,n}],{n,0,30}] (* Harvey P. Dale, Dec 31 2018 *)
a[n_]:=Binomial[2n, n]HypergeometricPFQ[{(1-2*n)/3, 2(1-n)/3, -2n/3}, {1/2-n, 1/2-n}, -3^3/2^4]; Array[a,24,0] (* Stefano Spezia, Jul 11 2024 *)
-
a(n)=sum(k=0,n,binomial(2*k,k)*binomial(2*k,n-k));
-
a(n)=polcoeff(1/sqrt(1-4*x*(1+x +x*O(x^n))^2),n,x); /* Using the g.f.: */
A073155
Leftmost column sequence of triangle A073153.
Original entry on oeis.org
1, 1, 4, 14, 56, 237, 1046, 4762, 22198, 105430, 508384, 2482297, 12248416, 60980875, 305955356, 1545397464, 7852100294, 40105277640, 205798130604, 1060467961508, 5485199090812, 28469067353686, 148220323891460
Offset: 0
a(3)=a(0)*[a(2)+a(1)]+[a(1)+a(0)]*[a(1)+a(0)]+[a(2)+a(1)]*a(0) =1*[4+1] + [1+1]*[1+1] + [4+1]*1 = 5 + 2*2 + 5 = 14.
A364336
G.f. satisfies A(x) = (1 + x) * (1 + x*A(x)^3).
Original entry on oeis.org
1, 2, 7, 39, 242, 1634, 11631, 85957, 653245, 5072862, 40077807, 321106623, 2602911282, 21308131235, 175909559897, 1462846379247, 12242600576066, 103035285071630, 871490142773640, 7404121610615520, 63157400073057627, 540689217572662413, 4644083121177225292
Offset: 0
-
A364336 := proc(n)
add( binomial(3*k+1,k) * binomial(3*k+1,n-k)/(3*k+1),k=0..n) ;
end proc:
seq(A364336(n),n=0..80); # R. J. Mathar, Jul 25 2023
-
nmax = 80; A[_] = 1;
Do[A[x_] = (1 + x)*(1 + x*A[x]^3) + O[x]^(nmax+1) // Normal, {nmax+1}];
CoefficientList[A[x], x] (* Jean-François Alcover, Mar 03 2024 *)
-
a(n) = sum(k=0, n, binomial(3*k+1, k)*binomial(3*k+1, n-k)/(3*k+1));
A364337
G.f. satisfies A(x) = (1 + x) * (1 + x*A(x)^4).
Original entry on oeis.org
1, 2, 9, 68, 580, 5406, 53270, 545844, 5757332, 62094217, 681653493, 7591431752, 85558696024, 974024788280, 11184192097016, 129378232148016, 1506363564912368, 17639001584452320, 207593804132718948, 2454236122156830254, 29132714097692056954, 347086786035103983446
Offset: 0
-
terms = 22; A[] = 0; Do[A[x] = (1+x)(1+x*A[x]^4) + O[x]^terms // Normal, terms]; CoefficientList[A[x], x] (* Stefano Spezia, Mar 24 2025 *)
-
a(n) = sum(k=0, n, binomial(4*k+1, k)*binomial(4*k+1, n-k)/(4*k+1));
A143330
G.f. satisfies: A(x) = (1 + x*A(x)^2)/(1 - x^2).
Original entry on oeis.org
1, 1, 3, 8, 25, 83, 289, 1041, 3847, 14504, 55569, 215727, 846761, 3354858, 13398965, 53888063, 218053915, 887107888, 3626373205, 14887942624, 61358959587, 253771944529, 1052917272543, 4381374717994, 18280470530047, 76459765772375
Offset: 0
G.f. = 1 + x + 3*x^2 + 8*x^3 + 25*x^4 + 83*x^5 + 289*x^6 + 1041*x^7 + ...
- Vincenzo Librandi, Table of n, a(n) for n = 0..1000
- Paul Barry, Notes on Riordan arrays and lattice paths, arXiv:2504.09719 [math.CO], 2025. See pp. 25, 29.
- James East and Nicholas Ham, Lattice paths and submonoids of Z^2, arXiv:1811.05735 [math.CO], 2018.
- Joseph P. S. Kung and Anna de Mier, Catalan lattice paths with rook, bishop and spider steps, Journal of Combinatorial Theory, Series A 120 (2013) 379-389. - From _N. J. A. Sloane_, Dec 27 2012
-
a := proc(n) option remember; if n <= 3 then return [1, 1, 3, 8][n + 1] fi;
((5 - n)*a(n - 4) + (2*n - 4)*a(n - 2) + (4*n - 2)*a(n - 1))/(n + 1) end:
seq(a(n), n = 0..25); # Peter Luschny, Jan 25 2023
-
CoefficientList[Series[(1 - x^2 - Sqrt[1 - 4 x - 2 x^2 + x^4])/(2 x), {x, 0, 30}], x] (* Vaclav Kotesovec, Sep 17 2013 *)
-
{a(n)=polcoeff((1-x^2-sqrt((1-x^2)^2-4*x+x^2*O(x^n)))/(2*x),n)}
A364338
G.f. satisfies A(x) = (1 + x) * (1 + x*A(x)^5).
Original entry on oeis.org
1, 2, 11, 105, 1140, 13555, 170637, 2235472, 30161255, 416248640, 5848462880, 83378361111, 1203100853951, 17537182300140, 257858115407535, 3819894878557990, 56958234329850060, 854192593184162160, 12875579347191388830, 194963091634569681550, 2964229359714424159370, 45234864131654311730160
Offset: 0
-
terms = 22; A[] = 0; Do[A[x] = (1+x)(1+x*A[x]^5) + O[x]^terms // Normal, terms]; CoefficientList[A[x], x] (* Stefano Spezia, Mar 24 2025 *)
-
a(n) = sum(k=0, n, binomial(5*k+1, k)*binomial(5*k+1, n-k)/(5*k+1));
Showing 1-10 of 29 results.
Comments