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.

Previous Showing 41-50 of 146 results. Next

A122366 Triangle read by rows: T(n,k) = binomial(2*n+1,k), 0 <= k <= n.

Original entry on oeis.org

1, 1, 3, 1, 5, 10, 1, 7, 21, 35, 1, 9, 36, 84, 126, 1, 11, 55, 165, 330, 462, 1, 13, 78, 286, 715, 1287, 1716, 1, 15, 105, 455, 1365, 3003, 5005, 6435, 1, 17, 136, 680, 2380, 6188, 12376, 19448, 24310, 1, 19, 171, 969, 3876, 11628, 27132, 50388, 75582, 92378, 1, 21
Offset: 0

Views

Author

Reinhard Zumkeller, Aug 30 2006

Keywords

Comments

Sum of n-th row = A000302(n) = 4^n.
Central terms give A052203.
Reversal of A111418. - Philippe Deléham, Mar 22 2007
Coefficient triangle for the expansion of one half of odd powers of 2*x in terms of Chebyshev's T-polynomials: ((2*x)^(2*n+1))/2 = Sum_{k=0..n} a(n,k)*T(2*(n-k)+1,x) with Chebyshev's T-polynomials. See A053120. - Wolfdieter Lang, Mar 07 2007
The signed triangle T(n,k)*(-1)^(n-k) appears in the formula ((2*sin(phi))^(2*n+1))/2 = Sum_{k=0..n} ((-1)^(n-k))*a(n,k)*sin((2*(n-k)+1)*phi). - Wolfdieter Lang, Mar 07 2007
The signed triangle T(n,k)*(-1)^(n-k) appears therefore in the formula (4-x^2)^n = Sum_{k=0..n} ((-1)^(n-k))*a(n,k)*S(2*(n-k),x) with Chebyshev's S-polynomials. See A049310 for S(n,x). - Wolfdieter Lang, Mar 07 2007
From Wolfdieter Lang, Sep 18 2012: (Start)
The triangle T(n,k) appears also in the formula F(2*l+1)^(2*n+1) = (1/5^n)*Sum_{k=0..n} T(n,k)*F((2*(n-k)+1)*(2*l+1)), l >= 0, n >= 0, with F=A000045 (Fibonacci).
The signed triangle Ts(n,k):=T(n,k)*(-1)^k appears also in the formula
F(2*l)^(2*n+1) = (1/5^n)*Sum_{k=0..n} Ts(n,k)*F((2(n-k)+1)*2*l), l >= 0, n >= 0, with F=A000045 (Fibonacci).
This is Lemma 2 of the K. Ozeki reference, p. 108, written for odd and even indices separately.
(End)

Examples

			.......... / 1 \ .......... =A062344(0,0)=A034868(0,0),
......... / 1 . \ ......... =T(0,0)=A034868(1,0),
........ / 1 2 . \ ........ =A062344(1,0..1)=A034868(2,0..1),
....... / 1 3 ... \ ....... =T(1,0..1)=A034868(3,0..1),
...... / 1 4 6 ... \ ...... =A062344(2,0..2)=A034868(4,0..2),
..... / 1 5 10 .... \ ..... =T(2,0..2)=A034868(5,0..2),
.... / 1 6 15 20 ... \ .... =A062344(3,0..3)=A034868(6,0..3),
... / 1 7 21 35 ..... \ ... =T(3,0..3)=A034868(7,0..3),
.. / 1 8 28 56 70 .... \ .. =A062344(4,0..4)=A034868(8,0..4),
. / 1 9 36 84 126 ..... \ . =T(4,0..4)=A034868(9,0..4).
Row n=2:[1,5,10] appears in the expansion ((2*x)^5)/2 = T(5,x)+5*T(3,x)+10*T(1,x).
Row n=2:[1,5,10] appears in the expansion ((2*cos(phi))^5)/2 = cos(5*phi)+5*cos(3*phi)+10*cos(1*phi).
The signed row n=2:[1,-5,10] appears in the expansion ((2*sin(phi))^5)/2 = sin(5*phi)-5*sin(3*phi)+10*sin(phi).
The signed row n=2:[1,-5,10] appears therefore in the expansion (4-x^2)^2 = S(4,x)-5*S(2,x)+10*S(0,x).
Triangle T(n,k) starts:
  n\k 0  1   2   3    4     5     6     7     8     9  ...
  0   1
  1   1  3
  2   1  5  10
  3   1  7  21  35
  4   1  9  36  84  126
  5   1 11  55 165  330   462
  6   1 13  78 286  715  1287  1716
  7   1 15 105 455 1365  3003  5005  6435
  8   1 17 136 680 2380  6188 12376 19448 24310
  9   1 19 171 969 3876 11628 27132 50388 75582 92378
  ...  - _Wolfdieter Lang_, Sep 18 2012
Row n=2, with F(n)=A000045(n) (Fibonacci number), l >= 0, see a comment above:
F(2*l)^5   = (1*F(10*l) - 5*F(6*l) + 10*F(2*l))/25,
F(2*l+1)^5 = (1*F(10*l+5) + 5*F(6*l+3) + 10*F(2*l+1))/25.
- _Wolfdieter Lang_, Sep 19 2012
		

References

  • T. J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2nd ed., Wiley, New York, 1990, pp. 54-55, Ex. 1.5.31.

Crossrefs

Cf. A062344.
Odd numbered rows of A008314. Even numbered rows of A008314 are A127673.

Programs

  • Haskell
    a122366 n k = a122366_tabl !! n !! k
    a122366_row n = a122366_tabl !! n
    a122366_tabl = f 1 a007318_tabl where
       f x (_:bs:pss) = (take x bs) : f (x + 1) pss
    -- Reinhard Zumkeller, Mar 14 2014
  • Mathematica
    T[_, 0] = 1;
    T[n_, k_] := T[n, k] = T[n-1, k-1] 2n(2n+1)/(k(2n-k+1));
    Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Dec 01 2021 *)

Formula

T(n,0)=1; T(n,k) = T(n-1,k-1)*2*n*(2*n+1)/(k*(2*n-k+1)) for k > 0.
T(n,0)=1; for n > 0: T(n,1)=n+2; for n > 1: T(n,n) = T(n-1,n-2) + 3*T(n-1,n-1), T(n,k) = T(n-1,k-2) + 2*T(n-1,k-1) + T(n-1,k), 1 < k < n.
T(n,n) = A001700(n).
T(n,k) = A034868(2*n+1,k) = A007318(2*n+1,k), 0 <= k <= n;
G.f.: (2*y)/((y-1)*sqrt(1-4*x*y)-4*x*y^2+(1-4*x)*y+1). - Vladimir Kruchinin, Oct 30 2020

Extensions

Chebyshev and trigonometric comments from Wolfdieter Lang, Mar 07 2007.
Typo in comments fixed, thanks to Philippe Deléham, who indicated this.

A029760 A sum with next-to-central binomial coefficients of even order, Catalan related.

Original entry on oeis.org

1, 8, 47, 244, 1186, 5536, 25147, 112028, 491870, 2135440, 9188406, 39249768, 166656772, 704069248, 2961699667, 12412521388, 51854046982, 216013684528, 897632738722, 3721813363288, 15401045060572, 63616796642368, 262357557683422, 1080387930269464
Offset: 0

Views

Author

Keywords

Comments

Proof by induction.
a(n) = total area below paths consisting of steps east (1,0) and north (0,1) from (0,0) to (n+2,n+2) that stay weakly below y=x. For example, the two paths with n=0 are
. _|.....|
The first has area 1 below it, the second area 0 and so a(0)=1. - David Callan, Dec 09 2004
Convolution of A000346 with A001700. - Philippe Deléham, May 19 2009

Crossrefs

Programs

  • Mathematica
    a[n_] := (n+3)^2 CatalanNumber[n+2]/2 - 2^(2n+3);
    Table[a[n], {n, 0, 23}] (* Jean-François Alcover, Sep 25 2018 *)

Formula

a(n) = 4^(n+1)*Sum_{k=1..n+1} binomial(2k, k-1)/4^k = ((n+3)^2)*C(n+2)/2-2^(2*n+3), C = Catalan. Also a(n+1)=4*a(n)+binomial(2(n+2), n+1).
G.f.: (d/dx)c(x)/(1-4*x), where c(x) = g.f. for Catalan numbers; convolution of A001791 and powers of 4. G.f. also c(x)^2/(1-4*x)^(3/2); convolution of Catalan numbers A000108 C(n), n >= 1, with A002457; convolution of A008549(n), n >= 1, with A000984 (central binomial coefficients).
a(n) = Sum_{k=0..n+1} A039598(n+1,k)*k^2. - Philippe Deléham, Dec 16 2007

A068763 Irregular triangle of the Fibonacci polynomials of A011973 multiplied diagonally by the Catalan numbers.

Original entry on oeis.org

1, 1, 1, 2, 2, 5, 6, 1, 14, 20, 6, 42, 70, 30, 2, 132, 252, 140, 20, 429, 924, 630, 140, 5, 1430, 3432, 2772, 840, 70, 4862, 12870, 12012, 4620, 630, 14, 16796, 48620, 51480, 24024, 4620, 252, 58786, 184756, 218790
Offset: 0

Views

Author

Wolfdieter Lang, Mar 04 2002

Keywords

Comments

The row length sequence of this array is [1,2,2,3,3,4,4,5,5,...] = A008619(n+2), n>=0.
The row polynomials p(n,x) := Sum_{m=0..floor((n+1)/2)} a(n,m)*x^m produce, for x = (b-a^2)/a^2 (not 0), the two parameter family of sequences K(a,b; n) := (a^(n+1))*p(n,(b-a^2)/a^2) with g.f. K(a,b; x) := (1-sqrt(1-4*x*(a+x*(b-a^2))))/(2*x).
Some members are: K(1,1; n)=A000108(n) (Catalan), K(1,2; n)=A025227(n-1), K(2,1; n)=A025228(n-1), K(1,3; n)=A025229(n-1), K(3,1; n)=A025230(n-1). For a=b=2..10 the sequences K(a,a; n)/a are A068764-A068772.
The column sequences (without leading 0's) are: A000108 (Catalan), A000984 (central binomial), A002457, 2*A002802, 5*A020918, 14*A020920, 42*A020922, ...
a(n,m) is the number of ways to designate exactly m cherries over all binary trees with n internal nodes. A cherry is an internal node whose descendants are both external nodes. Cf. A091894 which gives the number of binary trees with m cherries. - Geoffrey Critzer, Jul 24 2020
This irregular triangle is essentially that of A011973 with its diagonals multiplied by the Catalan numbers of A000108. The diagonals of this triangle are then rows of the Pascal matrix A007318 multiplied by the Catalan numbers. - Tom Copeland, Dec 23 2023

Examples

			The irregular triangle begins:
   n\m    0     1     2     3    4   5
   0:     1
   1:     1     1
   2:     2     2
   3:     5     6     1
   4:    14    20     6
   5:    42    70    30     2
   6:   132   252   140    20
   7:   429   924   630   140    5
   8:  1430  3432  2772   840   70
   9:  4862 12870 12012  4620  630  14
  10: 16796 48620 51480 24024 4620 252
  ...
p(3,x) = 5 + 6*x + x^2.
		

Crossrefs

Cf. A025227(n-1) (row sums).
Cf. A000007(n) (alternating row sums).

Programs

  • Mathematica
    nn = 10; b[z_] := (1 - Sqrt[1 - 4 z])/(2 z);Map[Select[#, # > 0 &] &,
    CoefficientList[Series[v b[v z] /. v -> (1 + u z ), {z, 0, nn}], {z, u}]] // Grid (* Geoffrey Critzer, Jul 24 2020 *)

Formula

a(n, m) = binomial(n+1-m, m)*C(n-m) if 0 <= m <= floor((n+1)/2), otherwise 0, with C(n) := A000108(n) (Catalan).
G.f. for column m=1, 2, ...: (x^(2*m-1))*C(m-1)/(1-4*x)^((2*m-1)/2); m=0: c(x), g.f. for A000108 (Catalan).
G.f. for row polynomials p(n, x): c(z) + x*z*c(x*(z^2)/(1-4*z))/sqrt(1-4*z) = (1-sqrt(1-4*z*(1+x*z)))/(2*z), where c(x) is the g.f. of A000108 (Catalan).
G.f. for triangle: (1 - sqrt(1 - 4*x (1 + y*x)))/(2*x). - Geoffrey Critzer, Jul 24 2020
The series expansion of f(x) = (1 + 2sx - sqrt(1 + 4sx + 4d^2x^2))/(2x) at x = 0 is (s^2 - d^2) x + (2 d^2s - 2 s^3) x^2 + (d^4 - 6 d^2 s^2 + 5 s^4) x^3 + (-6 d^4 s + 20 d^2 s^3 - 14 s^5) x^4 + ..., containing the coefficients of this array. With s = (a+b)/2 and d = (a-b)/2, then f(x)/ab = g(x) = (1 + (a+b)x - sqrt((1+(a+b)x)^2 - 4abx^2))/(2abx) = x - (a + b) x^2 + (a^2 + 3 a b + b^2) x^3 - (a^3 + 6 a^2 b + 6 a b^2 + b^3) x^4 + ..., containing the Narayana polynomials of A001263, which can be simply transformed into A033282. The compositional inverse about the origin of g(x) is g^(-1)(x) = x/((1-ax)(1-bx)) = x/((1-(s+d)x)(1-(s-d)x)) = x + (a + b) x^2 + (a^2 + a b + b^2) x^3 + (a^3 + a^2 b + a b^2 + b^3) x^4 + ..., containing the complete homogeneous symmetric polynomials h_n(a,b) = (a^n - b^n)/(a-b), which are the polynomials of A034867 when expressed in s and d, e.g., ((s + d)^7 - (s - d)^7)/(2 d) = d^6 + 21 d^4 s^2 + 35 d^2 s^4 + 7 s^6. A133437 and A134264 for compositional inversion of o.g.f.s can be used to relate the sets of polynomials above. - Tom Copeland, Nov 28 2023

Extensions

Title changed by Tom Copeland, Dec 23 2023

A144635 a(n) = 5^n*Sum_{ k=0..n } binomial(2*k,k)/5^k.

Original entry on oeis.org

1, 7, 41, 225, 1195, 6227, 32059, 163727, 831505, 4206145, 21215481, 106782837, 536618341, 2693492305, 13507578125, 67693008145, 339066121115, 1697664211795, 8497396194275, 42522326235175, 212749477704695, 1064285646397915, 5323532330953295, 26625895085494075
Offset: 0

Views

Author

N. J. A. Sloane, Jan 21 2009

Keywords

Crossrefs

Programs

  • Mathematica
    Table[5^n Sum[Binomial[2k,k]/5^k,{k,0,n}],{n,0,30}] (* Harvey P. Dale, Aug 08 2011 *)
    Round@Table[5^(n + 1/2) - 2^(n + 1) (2 n + 1)!! Hypergeometric2F1[1, n + 3/2, n + 2, 4/5]/(5 (n + 1)!), {n, 0, 20}] (* Round is equivalent to FullSimplify here, but is much faster - Vladimir Reshetnikov, Oct 14 2016 *)
  • PARI
    a(n) = 5^n*sum(k=0, n, binomial(2*k,k)/5^k); \\ Michel Marcus, Oct 14 2016

Formula

From Vaclav Kotesovec, Jun 12 2013: (Start)
G.f.: 1/((1-5*x)*sqrt(1-4*x)).
Recurrence: n*a(n) = (9*n-2)*a(n-1) - 10*(2*n-1)*a(n-2).
a(n) ~ 5^(n+1/2). (End)
a(n) = 5^(n+1/2) - 2^(n+1)*(2*n+1)!!*hypergeom([1,n+3/2], [n+2], 4/5)/(5*(n+1)!). - Vladimir Reshetnikov, Oct 14 2016
a(n) = Sum_{k=0..n} C(2*n+1,n-k)*A000032(2*k+1). - Vladimir Kruchinin Jan 14 2025

A331431 Triangle read by rows: T(n,k) = (-1)^(n+k)*(n+k+1)*binomial(n,k)*binomial(n+k,k) for n >= k >= 0.

Original entry on oeis.org

1, -2, 6, 3, -24, 30, -4, 60, -180, 140, 5, -120, 630, -1120, 630, -6, 210, -1680, 5040, -6300, 2772, 7, -336, 3780, -16800, 34650, -33264, 12012, -8, 504, -7560, 46200, -138600, 216216, -168168, 51480, 9, -720, 13860, -110880, 450450, -1009008, 1261260, -823680, 218790
Offset: 0

Views

Author

N. J. A. Sloane, Jan 17 2020

Keywords

Comments

Tables I, III, IV on pages 92 and 93 of Ser have integer entries and are A331430, A331431 (the present sequence), and A331432.
Given the system of equations 1 = Sum_{j=0..n} H(i, j) * x(j) for i = 2..n+2 where H(i,j) = 1/(i+j-1) for 1 <= i,j <= n is the n X n Hilbert matrix, then the solutions are x(j) = T(n, j). - Michael Somos, Mar 20 2020 [Corrected by Petros Hadjicostas, Jul 09 2020]

Examples

			Triangle begins:
   1;
  -2,    6;
   3,  -24,    30;
  -4,   60,  -180,     140;
   5, -120,   630,   -1120,     630;
  -6,  210, -1680,    5040,   -6300,     2772;
   7, -336,  3780,  -16800,   34650,   -33264,   12012;
  -8,  504, -7560,   46200, -138600,   216216, -168168,   51480;
   9, -720, 13860, -110880,  450450, -1009008, 1261260, -823680, 218790;
  ...
		

References

  • J. Ser, Les Calculs Formels des Séries de Factorielles. Gauthier-Villars, Paris, 1933, p. 93. See Table III.

Crossrefs

Columns 1 is A331433 or equally A007531, column 2 is A331434 or equally A054559; the last three diagonals are A002738, A002736, A002457.
Cf. A000290 (row sums), A002457,, A100071, A108666 (alternating row sums), A109188 (diagonal sums), A331322, A331323, A331430, A331432.

Programs

  • Magma
    [(-1)^(n+k)*(k+1)*(2*k+1)*Binomial(n+k+1,n-k)*Catalan(k): k in [0..n], n in [0..15]]; // G. C. Greubel, Mar 22 2022
    
  • Maple
    gf := k -> (1+x)^(-2*(k+1)): ser := k -> series(gf(k), x, 32):
    T := (n, k) -> ((2*k+1)!/(k!)^2)*coeff(ser(k), x, n-k):
    seq(seq(T(n,k), k=0..n),n=0..7); # Peter Luschny, Jan 18 2020
    S:=(n,k)->(-1)^(n+k)*(n+k+1)!/((k!)^2*(n-k)!);
    rho:=n->[seq(S(n,k),k=0..n)];
    for n from 0 to 14 do lprint(rho(n)); od: # N. J. A. Sloane, Jan 18 2020
  • Mathematica
    Table[(-1)^(n+k)*(n+k+1)*Binomial[2*k,k]*Binomial[n+k,n-k], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Mar 22 2022 *)
  • Sage
    flatten([[(-1)^(n+k)*(2*k+1)*binomial(2*k,k)*binomial(n+k+1,n-k) for k in (0..n)] for n in (0..15)]) # G. C. Greubel, Mar 22 2022

Formula

T(n, 0) = (-1)^n*A000027(n+1).
T(n, 1) = A331433(n-1) = (-1)^(n+1)*A007531(n+2).
T(n, 2) = A331434(n-2) = (-1)^n*A054559(n+3).
T(n, n-2) = A002738(n-2).
T(n, n-1) = (-1)*A002736(n).
T(n, n) = A002457(n).
T(2*n, n) = (-1)^n*(3*n+1)!/(n!)^3 = (-1)^n*A331322(n).
Sum_{k=0..n} T(n, k) = A000290(n+1) (row sums).
Sum_{k=0..n} (-1)^k*T(n, k) = (-1)^n*A108666(n+1) (alternating row sums).
Sum_{k=0..n} T(n-k, k) = (-1)^n*A109188(n+1) (diagonal sums).
2^n*Sum_{k=0..n} T(n, k)/2^k = (-1)^floor(n/2)*A100071(n+1) (positive half sums).
(-2)^n*Sum_{k=0..n} T(n, k)/(-2)^k = A331323(n) (negative half sums).
T(n, k) = ((2*k+1)!/(k!)^2)*[x^(n-k)] (1+x)^(-2*(k+1)). - Georg Fischer and Peter Luschny, Jan 18 2020
T(n,k) = (-1)^(n+k)*(n+k+1)!/((k!)^2*(n-k)!), for n >= k >= 0. - N. J. A. Sloane, Jan 18 2020
From Petros Hadjicostas, Jul 09 2020: (Start)
Michael Somos's formulas above can be restated as
Sum_{k=0..n} T(n,k)/(i+k) = 1 for i = 1..n+1.
These are special cases of the following formula that is alluded to (in some way) in Ser's book:
1 - Sum_{k=0..n} T(n,k)/(x + k) = (x-1)*...*(x-(n + 1))/(x*(x+1)*...*(x+n)).
Because T(n,k) = (-1)^(n+1)*(n + k + 1)*A331430(n,k) and Sum_{k=0..n} A331430(n,k) = (-1)^(n+1), one may derive this formula from Ser's second formula stated in A331430. (End)
T(2*n+1, n) = (-2)*(-27)^n*Pochhammer(4/3, n)*Pochhammer(5/3, n)/(n!*(n+1)!). - G. C. Greubel, Mar 22 2022

Extensions

Several typos in the data corrected by Georg Fischer and Peter Luschny, Jan 18 2020
Definition changed by N. J. A. Sloane, Jan 18 2020

A002699 a(n) = n*2^(2*n-1).

Original entry on oeis.org

0, 2, 16, 96, 512, 2560, 12288, 57344, 262144, 1179648, 5242880, 23068672, 100663296, 436207616, 1879048192, 8053063680, 34359738368, 146028888064, 618475290624, 2611340115968, 10995116277760, 46179488366592, 193514046488576
Offset: 0

Views

Author

Keywords

Comments

Right side of binomial sum Sum(i * binomial(2*n, i), i=1..n) - Yong Kong (ykong(AT)curagen.com), Dec 26 2000
Coefficients of shifted Chebyshev polynomials.
Starting with offset 1 = 4th binomial transform of [2, 8, 0, 0, 0, ...]. - Gary W. Adamson, Jul 21 2009
Let P(A) be the power set of an n-element set A and B be the Cartesian product of P(A) with itself. Then a(n) = the sum of the size of the symmetric difference of x and y for every (x,y) of B. - Ross La Haye, Jan 04 2013
It's the relation [27] with T(n) in the document of Ross. Following the last comment of Ross, A002697 is the similar sequence when replacing "symmetric difference" by "intersection" and A212698 is the similar sequence when replacing "symmetric difference" by union. - Bernard Schott, Jan 04 2013
If Delta = Symmetric difference, here, X Delta Y and Y Delta X are considered as two distinct Cartesian products, if we want to consider that X Delta Y = X Delta Y is the same Cartesian product, see A002697. - Bernard Schott, Jan 15 2013

References

  • C. Lanczos, Applied Analysis. Prentice-Hall, Englewood Cliffs, NJ, 1956, p. 518.
  • A. P. Prudnikov, Yu. A. Brychkov and O.I. Marichev, "Integrals and Series", Volume 1: "Elementary Functions", Chapter 4: "Finite Sums", New York, Gordon and Breach Science Publishers, 1986-1992.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

Formula

a(n) = 2 * A002697(n). - Bernard Schott, Jan 04 2013
a(n) = A212698(n) - A002697(n)
a(n) = 8*a(n-1)-16*a(n-2) with n>1, a(0)=0, a(1)=2. - Vincenzo Librandi, Mar 20 2013
G.f.: (2*x)/(1 - 4*x)^2. - Harvey P. Dale, Jul 28 2021
E.g.f.: (exp(4*x) - 1)/2. - Stefano Spezia, Aug 04 2022

A302971 Triangle read by rows: T(n,k) is the numerator of R(n,k) defined implicitly by the identity Sum_{i=0..l-1} Sum_{j=0..m} R(m,j)*(l-i)^j*i^j = l^(2*m+1) holding for all l,m >= 0.

Original entry on oeis.org

1, 1, 6, 1, 0, 30, 1, -14, 0, 140, 1, -120, 0, 0, 630, 1, -1386, 660, 0, 0, 2772, 1, -21840, 18018, 0, 0, 0, 12012, 1, -450054, 491400, -60060, 0, 0, 0, 51480, 1, -11880960, 15506040, -3712800, 0, 0, 0, 0, 218790, 1, -394788954, 581981400, -196409840, 8817900, 0, 0, 0, 0, 923780, 1, -16172552880, 26003271294, -10863652800, 1031151660, 0, 0, 0, 0, 0, 3879876
Offset: 0

Views

Author

Kolosov Petro, Apr 16 2018

Keywords

Examples

			Triangle begins:
------------------------------------------------------------------------
k=   0          1         2         3    4     5      6      7       8
------------------------------------------------------------------------
n=0: 1;
n=1: 1,         6;
n=2: 1,         0,       30;
n=3: 1,       -14,        0,      140;
n=4: 1,      -120,        0,        0, 630;
n=5: 1,     -1386,      660,        0,   0, 2772;
n=6: 1,    -21840,    18018,        0,   0,    0, 12012;
n=7: 1,   -450054,   491400,   -60060,   0,    0,     0, 51480;
n=8: 1, -11880960, 15506040, -3712800,   0,    0,     0,     0, 218790;
		

Crossrefs

Items of second row are the coefficients in the definition of A287326.
Items of third row are the coefficients in the definition of A300656.
Items of fourth row are the coefficients in the definition of A300785.
T(n,n) gives A002457(n).
Denominators of R(n,k) are shown in A304042.
Row sums return A000079(2n+1) - 1.

Programs

  • Maple
    R := proc(n, k) if k < 0 or k > n then return 0 fi; (2*k+1)*binomial(2*k, k);
    if n = k then % else -%*add((-1)^j*R(n, j)*binomial(j, 2*k+1)*
    bernoulli(2*j-2*k)/(j-k), j=2*k+1..n) fi end: T := (n, k) -> numer(R(n, k)):
    seq(print(seq(T(n, k), k=0..n)), n=0..12);
    # Numerical check that S(m, n) = n^(2*m+1):
    S := (m, n) -> add(add(R(m, j)*(n-k)^j*k^j, j=0..m), k=0..n-1):
    seq(seq(S(m, n) - n^(2*m+1), n=0..12), m=0..12); # Peter Luschny, Apr 30 2018
  • Mathematica
    R[n_, k_] := 0
    R[n_, k_] := (2 k + 1)*Binomial[2 k, k]*
       Sum[R[n, j]*Binomial[j, 2 k + 1]*(-1)^(j - 1)/(j - k)*
       BernoulliB[2 j - 2 k], {j, 2 k + 1, n}] /; 2 k + 1 <= n
    R[n_, k_] := (2 n + 1)*Binomial[2 n, n] /; k == n;
    T[n_, k_] := Numerator[R[n, k]];
    (* Print Fifteen Initial rows of Triangle A302971 *)
    Column[ Table[ T[n, k], {n, 0, 15}, {k, 0, n}], Center]
  • PARI
    T(n, k) = if ((n>k) || (n<0), 0, if (k==n, (2*n+1)*binomial(2*n, n), if (2*n+1>k, 0, if (n==0, 1, (2*n+1)*binomial(2*n, n)*sum(j=2*n+1, k+1, T(j, k)*binomial(j, 2*n+1)*(-1)^(j-1)/(j-n)*bernfrac(2*j-2*n))))));
    tabl(nn) = for (n=0, nn, for (k=0, n, print1(numerator(T(k,n)), ", ")); print); \\ Michel Marcus, Apr 27 2018

Formula

Recurrence given by Max Alekseyev (see the MathOverflow link):
R(n, k) = 0 if k < 0 or k > n.
R(n, k) = (2k+1)*binomial(2k, k) if k = n.
R(n, k) = (2k+1)*binomial(2k, k)*Sum_{j=2k+1..n} R(n, j)*binomial(j, 2k+1)*(-1)^(j-1)/(j-k)*Bernoulli(2j-2k), otherwise.
T(n, k) = numerator(R(n, k)).

A304042 Triangle read by rows: T(n,k) is the denominator of R(n,k) defined implicitly by the identity Sum_{i=0..l-1} Sum_{j=0..m} R(m,j)*(l-i)^j*i^j = l^(2*m+1) holding for all l,m >= 0.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 0

Views

Author

Kolosov Petro, May 05 2018

Keywords

Examples

			Triangle begins:
-----------------------------------------------------
k=    0  1  2  3  4  5  6  7  8  9  10 11 12 13 14 15
-----------------------------------------------------
n=0:  1;
n=1:  1, 1;
n=2:  1, 1, 1;
n=3:  1, 1, 1, 1;
n=4:  1, 1, 1, 1, 1;
n=5:  1, 1, 1, 1, 1, 1;
n=6:  1, 1, 1, 1, 1, 1, 1;
n=7:  1, 1, 1, 1, 1, 1, 1, 1;
n=8:  1, 1, 1, 1, 1, 1, 1, 1, 1;
n=9:  1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=10: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=11: 1, 5, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1;
n=12: 1, 3, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=13: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=14: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
n=15: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1;
		

Crossrefs

Programs

  • Mathematica
    R[n_, k_] := 0
    R[n_, k_] := (2 k + 1)*Binomial[2 k, k]*
       Sum[R[n, j]*Binomial[j, 2 k + 1]*(-1)^(j - 1)/(j - k)*
       BernoulliB[2 j - 2 k], {j, 2 k + 1, n}] /; 2 k + 1 <= n
    R[n_, k_] := (2 n + 1)*Binomial[2 n, n] /; k == n;
    T[n_, k_] := Denominator[R[n, k]];
    (* Print Fifteen Initial rows of Triangle A304042 *)
    Column[ Table[ T[n, k], {n, 0, 15}, {k, 0, n}], Center]
  • PARI
    up_to = 1274; \\ = binomial(50+1,2)-1
    A304042aux(n, k) = if((k<0)||(k>n),0,(k+k+1)*binomial(2*k, k)*if(k==n,1,sum(j=k+k+1,n, A304042aux(n, j)*binomial(j, k+k+1)*((-1)^(j-1))/(j-k)*bernfrac(2*(j-k)))));
    A304042tr(n, k) = denominator(A304042aux(n, k));
    A304042list(up_to) = { my(v = vector(up_to), i=0); for(n=0,oo, for(k=0,n, if(i++ > up_to, return(v)); v[i] = A304042tr(n,k))); (v); };
    v304042 = A304042list(1+up_to);
    A304042(n) = v304042[1+n]; \\ Antti Karttunen, Nov 07 2018

Formula

Recurrence given by Max Alekseyev (see the MathOverflow link):
R(n, k) = 0 if k < 0 or k > n.
R(n, k) = (2k+1)*binomial(2k, k) if k = n.
R(n, k) = (2k+1)*binomial(2k, k)*Sum_{j=2k+1..n} R(n, j)*binomial(j, 2k+1)*(-1)^(j-1)/(j-k)*Bernoulli(2j-2k), otherwise.
T(n, k) = denominator(R(n, k)).

A020922 Expansion of 1/(1-4*x)^(11/2).

Original entry on oeis.org

1, 22, 286, 2860, 24310, 184756, 1293292, 8498776, 53117350, 318704100, 1848483780, 10418726760, 57302997180, 308554600200, 1630931458200, 8480843582640, 43464323361030, 219878341708740, 1099391708543700, 5439095821216200, 26651569523959380, 129450480544945560
Offset: 0

Views

Author

Keywords

Comments

Also convolution of A000984 with A040075, also convolution of A000302 with A020920, also convolution of A002457 with A038846, also convolution of A002697 with A020918, also convolution of A002802 with A038845. - Rui Duarte, Oct 08 2011

Crossrefs

Programs

  • GAP
    List([0..30], n-> Binomial(n+5, 5)*Binomial(2*n+10, n+5)/252); # G. C. Greubel, Jul 20 2019
  • Magma
    [(2*n+9)*(2*n+7)*(2*n+5)*(2*n+3)*(2*n+1)*Binomial(2*n, n)/945: n in [0..30]] // Vincenzo Librandi, Jul 05 2013
    
  • Mathematica
    CoefficientList[Series[1/(1-4x)^(11/2), {x,0,30}], x] (* Vincenzo Librandi, Jul 05 2013 *)
  • PARI
    vector(30, n, n--; m=n+5; binomial(m, 5)*binomial(2*m, m)/252) \\ G. C. Greubel, Jul 20 2019
    
  • Sage
    [binomial(n+5, 5)*binomial(2*n+10, n+5)/252 for n in (0..30)] # G. C. Greubel, Jul 20 2019
    

Formula

a(n) = binomial(n+5, 5)*A000984(n+5)/A000984(5), where A000984 are central binomial coefficients. - Wolfdieter Lang
From Rui Duarte, Oct 08 2011: (Start)
a(n) = ((2n+9)(2n+7)(2n+5)(2n+3)(2n+1)/(9*7*5*3*1)) * binomial(2n, n).
a(n) = binomial(2n+10, 10) * binomial(2n, n) / binomial(n+5, 5).
a(n) = binomial(n+5, 5) * binomial(2n+10, n+5) / binomial(10, 5).
a(n) = Sum_{ i_1+i_2+i_3+i_4+i_5+i_6+i_7+i_8+i_9+i_10+i_11 = n } f(i_1)* f(i_2)*f(i_3)*f(i_4)*f(i_5)*f(i_6)*f(i_7)*f(i_8)*f(i_9)*f(i_10)*f(i_11) with f(k)=A000984(k). (End)
Boas-Buck recurrence: a(n) = (22/n)*Sum_{k=0..n-1} 4^(n-k-1)*a(k), n >= 1, a(0) = 1. Proof from a(n) = A046521(n+5, 5). See a comment there. - Wolfdieter Lang, Aug 10 2017
From Amiram Eldar, Mar 25 2022: (Start)
Sum_{n>=0} 1/a(n) = 162*sqrt(3)*Pi - 30816/35.
Sum_{n>=0} (-1)^n/a(n) = 4500*sqrt(5)*log(phi) - 33888/7, where phi is the golden ratio (A001622). (End)

A021010 Triangle of coefficients of Laguerre polynomials L_n(x) (powers of x in decreasing order).

Original entry on oeis.org

1, -1, 1, 1, -4, 2, -1, 9, -18, 6, 1, -16, 72, -96, 24, -1, 25, -200, 600, -600, 120, 1, -36, 450, -2400, 5400, -4320, 720, -1, 49, -882, 7350, -29400, 52920, -35280, 5040, 1, -64, 1568, -18816, 117600, -376320, 564480, -322560, 40320, -1, 81, -2592, 42336, -381024, 1905120, -5080320, 6531840, -3265920, 362880
Offset: 0

Views

Author

Keywords

Comments

abs(T(n,k)) = k!*binomial(n,k)^2 = number of k-matchings of the complete bipartite graph K_{n,n}. Example: abs(T(2,2))=2 because in the bipartite graph K_{2,2} with vertex sets {A,B},{A',B'} we have the 2-matchings {AA',BB'} and {AB',BA'}. Row sums of the absolute values yield A002720. - Emeric Deutsch, Dec 25 2004

Examples

			   1;
  -1,   1;
   1,  -4,   2;
  -1,   9, -18,   6;
   1, -16,  72, -96,  24;
  ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 799.

Crossrefs

Cf. A002720, A021009, A009940 (row sums).
Central terms: A295383.

Programs

  • Magma
    [[(-1)^(n-k)*Factorial(k)*Binomial(n, k)^2: k in [0..n]]: n in [0..10]]; // G. C. Greubel, Feb 06 2018
  • Maple
    T:=(n,k)->(-1)^(n-k)*k!*binomial(n,k)^2: for n from 0 to 9 do seq(T(n,k),k=0..n) od; # yields sequence in triangular form # Emeric Deutsch, Dec 25 2004
  • Mathematica
    Flatten[ Table[ Reverse[ CoefficientList[n!*LaguerreL[n, x], x]], {n, 0, 9}]] (* Jean-François Alcover, Nov 24 2011 *)
  • PARI
    LaguerreL(n,v='x) = {
      my(x='x+O('x^(n+1)), t='t);
      subst(polcoeff(exp(-x*t/(1-x))/(1-x), n), 't, v);
    };
    concat(apply(n->Vec(n!*LaguerreL(n)), [0..9])) \\ Gheorghe Coserea, Oct 26 2017
    
  • PARI
    row(n) = Vec(n!*pollaguerre(n)); \\ Michel Marcus, Feb 06 2021
    

Formula

T(n, k) = (-1)^(n-k)*k!*binomial(n, k)^2. - Emeric Deutsch, Dec 25 2004
Previous Showing 41-50 of 146 results. Next