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.

Showing 1-10 of 12 results. Next

A002457 a(n) = (2n+1)!/n!^2.

Original entry on oeis.org

1, 6, 30, 140, 630, 2772, 12012, 51480, 218790, 923780, 3879876, 16224936, 67603900, 280816200, 1163381400, 4808643120, 19835652870, 81676217700, 335780006100, 1378465288200, 5651707681620, 23145088600920, 94684453367400, 386971244197200, 1580132580471900
Offset: 0

Views

Author

Keywords

Comments

Expected number of matches remaining in Banach's modified matchbox problem (counted when last match is drawn from one of the two boxes), multiplied by 4^(n-1). - Michael Steyer, Apr 13 2001
Hankel transform is (-1)^n*A014480(n). - Paul Barry, Apr 26 2009
Convolved with A000108: (1, 1, 1, 5, 14, 42, ...) = A000531: (1, 7, 38, 187, 874, ...). - Gary W. Adamson, May 14 2009
Convolution of A000302 and A000984. - Philippe Deléham, May 18 2009
1/a(n) is the integral of (x(1-x))^n on interval [0,1]. Apparently John Wallis computed these integrals for n=0,1,2,3,.... A004731, shifted left by one, gives numerators/denominators of related integrals (1-x^2)^n on interval [0,1]. - Marc van Leeuwen, Apr 14 2010
Extend the triangular peaks of Dyck paths of semilength n down to the baseline forming (possibly) larger and overlapping triangles. a(n) = sum of areas of these triangles. Also a(n) = triangular(n) * Catalan(n). - David Scambler, Nov 25 2010
Let H be the n X n Hilbert matrix H(i,j) = 1/(i+j-1) for 1 <= i,j <= n. Let B be the inverse matrix of H. The sum of the elements in row n of B equals a(n-1). - T. D. Noe, May 01 2011
Apparently the number of peaks in all symmetric Dyck paths with semilength 2n+1. - David Scambler, Apr 29 2013
Denominator of central elements of Leibniz's Harmonic Triangle A003506.
Central terms of triangle A116666. - Reinhard Zumkeller, Nov 02 2013
Number of distinct strings of length 2n+1 using n letters A, n letters B, and 1 letter C. - Hans Havermann, May 06 2014
Number of edges in the Hasse diagram of the poset of partitions in the n X n box ordered by containment (from Havermann's comment above, C represents the square added in the edge). - William J. Keith, Aug 18 2015
Let V(n, r) denote the volume of an n-dimensional sphere with radius r then V(n, 1/2^n) = V(n-1, 1/2^n) / a((n-1)/2) for all odd n. - Peter Luschny, Oct 12 2015
a(n) is the result of processing the n+1 row of Pascal's triangle A007318 with the method of A067056. Example: Let n=3. Given the 4th row of Pascal's triangle 1,4,6,4,1, we get 1*(4+6+4+1) + (1+4)*(6+4+1) + (1+4+6)*(4+1) + (1+4+6+4)*1 = 15+55+55+15 = 140 = a(3). - J. M. Bergot, May 26 2017
a(n) is the number of (n+1) X 2 Young tableaux with a two horizontal walls between the first and second column. If there is a wall between two cells, the entries may be decreasing; see [Banderier, Wallner 2021] and A000984 for one horizontal wall. - Michael Wallner, Jan 31 2022
a(n) is the number of facets of the symmetric edge polytope of the cycle graph on 2n+1 vertices. - Mariel Supina, May 12 2022
Diagonal of the rational function 1 / (1 - x - y)^2. - Ilya Gutkovskiy, Apr 23 2025

Examples

			G.f. = 1 + 6*x + 30*x^2 + 140*x^3 + 630*x^4 + 2772*x^5 + 12012*x^6 + 51480*x^7 + ...
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 159.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 83, Problem 25; p. 168, #30.
  • W. Feller, An Introduction to Probability Theory and Its Applications, Vol. I.
  • C. Jordan, Calculus of Finite Differences. Röttig and Romwalter, Budapest, 1939; Chelsea, NY, 1965, p. 449.
  • M. Klamkin, ed., Problems in Applied Mathematics: Selections from SIAM Review, SIAM, 1990; see pp. 127-129.
  • C. Lanczos, Applied Analysis. Prentice-Hall, Englewood Cliffs, NJ, 1956, p. 514.
  • 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.
  • J. Ser, Les Calculs Formels des Séries de Factorielles. Gauthier-Villars, Paris, 1933, p. 92.
  • 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).
  • J. Wallis, Operum Mathematicorum, pars altera, Oxford, 1656, pp 31,34 [Marc van Leeuwen, Apr 14 2010]

Crossrefs

Cf. A000531 (Banach's original match problem).
Cf. A033876, A000984, A001803, A132818, A046521 (second column).
A diagonal of A331430.
The rightmost diagonal of the triangle A331431.

Programs

Formula

G.f.: (1-4x)^(-3/2) = 1F0(3/2;;4x).
a(n-1) = binomial(2*n, n)*n/2 = binomial(2*n-1, n)*n.
a(n-1) = 4^(n-1)*Sum_{i=0..n-1} binomial(n-1+i, i)*(n-i)/2^(n-1+i).
a(n) ~ 2*Pi^(-1/2)*n^(1/2)*2^(2*n)*{1 + 3/8*n^-1 + ...}. - Joe Keane (jgk(AT)jgk.org), Nov 21 2001
(2*n+2)!/(2*n!*(n+1)!) = (n+n+1)!/(n!*n!) = 1/beta(n+1, n+1) in A061928.
Sum_{i=0..n} i * binomial(n, i)^2 = n*binomial(2*n, n)/2. - Yong Kong (ykong(AT)curagen.com), Dec 26 2000
a(n) ~ 2*Pi^(-1/2)*n^(1/2)*2^(2*n). - Joe Keane (jgk(AT)jgk.org), Jun 07 2002
a(n) = 1/Integral_{x=0..1} x^n (1-x)^n dx. - Fred W. Helenius (fredh(AT)ix.netcom.com), Jun 10 2003
E.g.f.: exp(2*x)*((1+4*x)*BesselI(0, 2*x) + 4*x*BesselI(1, 2*x)). - Vladeta Jovovic, Sep 22 2003
a(n) = Sum_{i+j+k=n} binomial(2i, i)*binomial(2j, j)*binomial(2k, k). - Benoit Cloitre, Nov 09 2003
a(n) = (2*n+1)*A000984(n) = A005408(n)*A000984(n). - Zerinvary Lajos, Dec 12 2010
a(n-1) = Sum_{k=0..n} A039599(n,k)*A000217(k), for n >= 1. - Philippe Deléham, Jun 10 2007
Sum of (n+1)-th row terms of triangle A132818. - Gary W. Adamson, Sep 02 2007
Sum_{n>=0} 1/a(n) = 2*Pi/3^(3/2). - Jaume Oliver Lafont, Mar 07 2009
a(n) = Sum_{k=0..n} binomial(2k,k)*4^(n-k). - Paul Barry, Apr 26 2009
a(n) = A000217(n) * A000108(n). - David Scambler, Nov 25 2010
a(n) = f(n, n-3) where f is given in A034261.
a(n) = A005430(n+1)/2 = A002011(n)/4.
a(n) = binomial(2n+2, 2) * binomial(2n, n) / binomial(n+1, 1), a(n) = binomial(n+1, 1) * binomial(2n+2, n+1) / binomial(2, 1) = binomial(2n+2, n+1) * (n+1)/2. - Rui Duarte, Oct 08 2011
G.f.: (G(0) - 1)/(4*x) where G(k) = 1 + 2*x*((2*k + 3)*G(k+1) - 1)/(k + 1). - Sergei N. Gladkovskii, Dec 03 2011 [Edited by Michael Somos, Dec 06 2013]
G.f.: 1 - 6*x/(G(0)+6*x) where G(k) = 1 + (4*x+1)*k - 6*x - (k+1)*(4*k-2)/G(k+1); (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Aug 13 2012
G.f.: Q(0), where Q(k) = 1 + 4*(2*k + 1)*x*(2*k + 2 + Q(k+1))/(k+1). - Sergei N. Gladkovskii, May 10 2013 [Edited by Michael Somos, Dec 06 2013]
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - 4*x*(2*k+3)/(4*x*(2*k+3) + 2*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 06 2013
a(n) = 2^(4n)/Sum_{k=0..n} (-1)^k*C(2n+1,n-k)/(2k+1). - Mircea Merca, Nov 12 2013
a(n) = (2*n)!*[x^(2*n)] HeunC(0,0,-2,-1/4,7/4,4*x^2) where [x^n] f(x) is the coefficient of x^n in f(x) and HeunC is the Heun confluent function. - Peter Luschny, Nov 22 2013
0 = a(n) * (16*a(n+1) - 2*a(n+2)) + a(n+1) * (a(n+2) - 6*a(n+1)) for all n in Z. - Michael Somos, Dec 06 2013
a(n) = 4^n*binomial(n+1/2, 1/2). - Peter Luschny, Apr 24 2014
a(n) = 4^n*hypergeom([-2*n,-2*n-1,1/2],[-2*n-2,1],2)*(n+1)*(2*n+1). - Peter Luschny, Sep 22 2014
a(n) = 4^n*hypergeom([-n,-1/2],[1],1). - Peter Luschny, May 19 2015
a(n) = 2*4^n*Gamma(3/2+n)/(sqrt(Pi)*Gamma(1+n)). - Peter Luschny, Dec 14 2015
Sum_{n >= 0} 2^(n+1)/a(n) = Pi, related to Newton/Euler's Pi convergence transformation series. - Tony Foster III, Jul 28 2016. See the Weisstein Pi link, eq. (23). - Wolfdieter Lang, Aug 26 2016
Boas-Buck recurrence: a(n) = (6/n)*Sum_{k=0..n-1} 4^(n-k-1)*a(k), n >= 1, and a(0) = 1. Proof from a(n) = A046521(n+1,1). See comment in A046521. - Wolfdieter Lang, Aug 10 2017
a(n) = (1/3)*Sum_{i = 0..n+1} C(n+1,i)*C(n+1,2*n+1-i)*C(3*n+2-i,n+1) = (1/3)*Sum_{i = 0..2*n+1} (-1)^(i+1)*C(2*n+1,i)*C(n+i+1,i)^2. - Peter Bala, Feb 07 2018
a(n) = (2*n+1)*binomial(2*n, n). - Kolosov Petro, Apr 16 2018
a(n) = (-4)^n*binomial(-3/2, n). - Peter Luschny, Oct 23 2018
a(n) = 1 / Sum_{s=0..n} (-1)^s * binomial(n, s) / (n+s+1). - Kolosov Petro, Jan 22 2019
a(n) = Sum_{k = 0..n} (2*k + 1)*binomial(2*n + 1, n - k). - Peter Bala, Feb 25 2019
4^n/a(n) = Integral_{x=0..1} (1 - x^2)^n. - Michael Somos, Jun 13 2019
D-finite with recurrence: 0 = a(n)*(6 + 4*n) - a(n+1)*(n + 1) for all n in Z. - Michael Somos, Jun 13 2019
Sum_{n>=0} (-1)^n/a(n) = 4*arcsinh(1/2)/sqrt(5). - Amiram Eldar, Sep 10 2020
From Jianing Song, Apr 10 2022: (Start)
G.f. for {1/a(n)}: 4*arcsin(sqrt(x)/2) / sqrt(x*(4-x)).
E.g.f. for {1/a(n)}: exp(x/4)*sqrt(Pi/x)*erf(sqrt(x)/2). (End)
G.f. for {1/a(n)}: 4*arctan(sqrt(x/(4-x))) / sqrt(x*(4-x)). - Michael Somos, Jun 17 2023
a(n) = Sum_{k = 0..n} (-1)^(n+k) * (n + 2*k + 1)*binomial(n+k, k). This is the particular case m = 1 of the identity Sum_{k = 0..m*n} (-1)^k * (n + 2*k + 1) * binomial(n+k, k) = (-1)^(m*n) * (m*n + 1) * binomial((m+1)*n+1, n). Cf. A090816 and A306290. - Peter Bala, Nov 02 2024
a(n) = (1/Pi)*(2*n + 1)*(2^(2*n + 1))*Integral_{x=0..oo} 1/(x^2 + 1)^(n + 1) dx. - Velin Yanev, Jan 28 2025

A001803 Numerators in expansion of (1 - x)^(-3/2).

Original entry on oeis.org

1, 3, 15, 35, 315, 693, 3003, 6435, 109395, 230945, 969969, 2028117, 16900975, 35102025, 145422675, 300540195, 9917826435, 20419054425, 83945001525, 172308161025, 1412926920405, 2893136075115, 11835556670925
Offset: 0

Views

Author

Keywords

Comments

a(n) is the denominator of the integral from 0 to Pi of (sin(x))^(2*n+1). - James R. Buddenhagen, Aug 17 2008
a(n) is the denominator of (2n)!!/(2*n + 1)!! = 2^(2*n)*n!*n!/(2*n + 1)! (see Andersson). - N. J. A. Sloane, Jun 27 2011
a(n) = (2*n + 1)*A001790(n). A046161(n)/a(n) = 1, 2/3, 8/15, 16/35, 128/315, 256/693, ... is binomial transform of Madhava-Gregory-Leibniz series for Pi/4 (i.e., 1 - 1/3 + 1/5 - 1/7 + ... ). See A173384 and A173396. - Paul Curtz, Feb 21 2010
a(n) is the denominator of Integral_{x=-oo..oo} sech(x)^(2*n+2) dx. The corresponding numerator is A101926(n). - Mohammed Yaseen, Jul 25 2023

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. 798.
  • G. Prévost, Tables de Fonctions Sphériques. Gauthier-Villars, Paris, 1933, pp. 156-157.
  • 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).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 6, equation 6:14:9 at page 51.

Crossrefs

The denominator is given in A046161.
Largest odd divisors of A001800, A002011, A002457, A005430, A033876, A086228.
Bisection of A004731, A004735, A086116.
Second column of triangle A100258.
Cf. A002596 (numerators in expansion of (1-x)^(1/2)).
Cf. A161198 (triangle related to the series expansions of (1-x)^((-1-2*n)/2)).
A163590 is the odd part of the swinging factorial, A001790 at even indices. - Peter Luschny, Aug 01 2009

Programs

  • Julia
    A001803(n) = sum(<<(A001790(k), A005187(n) - A005187(k)) for k in 0:n) # Peter Luschny, Oct 03 2019
    
  • Magma
    A001803:= func< n | Numerator(Binomial(n+2,2)*Catalan(n+1)/4^n) >;
    [A001803(n): n in [0..30]]; // G. C. Greubel, Apr 27 2025
    
  • Maple
    swing := proc(n) option remember; if n = 0 then 1 elif irem(n, 2) = 1 then swing(n-1)*n else 4*swing(n-1)/n fi end:
    sigma := n -> 2^(add(i,i= convert(iquo(n,2),base,2))):
    a := n -> swing(2*n+1)/sigma(2*n+1); # Peter Luschny, Aug 01 2009
    A001803 := proc(n) (2*n+1)*binomial(2*n,n)/4^n ; numer(%) ; end proc: # R. J. Mathar, Jul 06 2011
    a := n -> denom(Pi*binomial(n, -1/2)): seq(a(n), n = 0..22); # Peter Luschny, Dec 06 2024
  • Mathematica
    Numerator/@CoefficientList[Series[(1-x)^(-3/2),{x,0,25}],x]  (* Harvey P. Dale, Feb 19 2011 *)
    Table[Denominator[Beta[1, n + 1, 1/2]], {n, 0, 22}] (* Gerry Martens, Nov 13 2016 *)
  • PARI
    a(n) = numerator((2*n+1)*binomial(2*n,n)/(4^n)); \\ Altug Alkan, Sep 06 2018
    
  • SageMath
    def A001803(n): return numerator((n+1)*binomial(2*n+2,n+1)/2^(2*n+1))
    print([A001803(n) for n in range(31)]) # G. C. Greubel, Apr 27 2025

Formula

a(n) = (2*n + 1)! /(n!^2*2^A000120(n)) = (n + 1)*binomial(2*n+2,n+1)/2^(A000120(n)+1). - Ralf Stephan, Mar 10 2004
From Johannes W. Meijer, Jun 08 2009: (Start)
a(n) = numerator( (2*n+1)*binomial(2*n,n)/(4^n) ).
(1 - x)^(-3/2) = Sum_{n>=0} ((2*n+1)*binomial(2*n,n)/4^n)*x^n. (End)
Truncations of rational expressions like those given by the numerator or denominator operators are artifacts in integer formulas and have many disadvantages. A pure integer formula follows. Let n$ denote the swinging factorial and sigma(n) = number of '1's in the base-2 representation of floor(n/2). Then a(n) = (2*n+1)$ / sigma(2*n+1) = A056040(2*n+1) / A060632(2*n+2). Simply said: This sequence gives the odd part of the swinging factorial at odd indices. - Peter Luschny, Aug 01 2009
a(n) = denominator(Pi*binomial(n, -1/2)). - Peter Luschny, Dec 06 2024

A088961 Zigzag matrices listed entry by entry.

Original entry on oeis.org

3, 5, 5, 5, 10, 14, 14, 7, 14, 21, 21, 7, 21, 35, 42, 48, 27, 9, 48, 69, 57, 36, 27, 57, 78, 84, 9, 36, 84, 126, 132, 165, 110, 44, 11, 165, 242, 209, 121, 55, 110, 209, 253, 220, 165, 44, 121, 220, 297, 330, 11, 55, 165, 330, 462
Offset: 1

Views

Author

Paul Boddington, Oct 28 2003

Keywords

Comments

For each n >= 1 the n X n matrix Z(n) is constructed as follows. The i-th row of Z(n) is obtained by generating a hexagonal array of numbers with 2*n+1 rows, 2*n numbers in the odd numbered rows and 2*n+1 numbers in the even numbered rows. The first row is all 0's except for two 1's in the i-th and the (2*n+1-i)th positions. The remaining rows are generated using the same rule for generating Pascal's triangle. The i-th row of Z(n) then consists of the first n numbers in the bottom row of our array.
For example the top row of Z(2) is [5,5], found from the array:
. 1 0 0 1
1 1 0 1 1
. 2 1 1 2
2 3 2 3 2
. 5 5 5 5
Zigzag matrices have remarkable properties. Here is a selection:
1) Z(n) is symmetric.
2) det(Z(n)) = A085527(n).
3) tr(Z(n)) = A033876(n-1).
4) If 2*n+1 is a power of a prime p then all entries of Z(n) are multiples of p.
5) If 4*n+1 is a power of a prime p then the dot product of any two distinct rows of Z(n) is a multiple of p.
6) It is always possible to move from the bottom left entry of Z(n) to the top right entry using only rightward and upward moves and visiting only odd numbers.
A001700(n) = last term of last row of Z(n): a(A000330(n-1)) = A001700(n); A230585(n) = first term of first row of Z(n): a(A056520(n-1)) = A230585(n); A051417(n) = greatest common divisor of entries of Z(n). - Reinhard Zumkeller, Oct 25 2013

Examples

			The first five values are 3, 5, 5, 5, 10 because the first two zigzag matrices are [[3]] and [[5,5],[5,10]].
		

Crossrefs

Programs

  • Haskell
    a088961 n = a088961_list !! (n-1)
    a088961_list = concat $ concat $ map f [1..] where
       f x = take x $ g (take x (1 : [0,0..])) where
         g us = (take x $ g' us) : g (0 : init us)
         g' vs = last $ take (2 * x + 1) $
                        map snd $ iterate h (0, vs ++ reverse vs)
       h (p,ws) = (1 - p, drop p $ zipWith (+) ([0] ++ ws) (ws ++ [0]))
    -- Reinhard Zumkeller, Oct 25 2013
  • Mathematica
    Flatten[Table[Binomial[2n,n+j-i]-Binomial[2n,n+i+j]+ Binomial[2n, 3n+1-i-j], {n,5},{i,n},{j,n}]] (* Harvey P. Dale, Dec 15 2011 *)

Formula

The ij entry of Z(n) is binomial(2*n, n+j-i) - binomial(2*n, n+i+j) + binomial(2*n, 3*n+1-i-j).

A127952 Triangle read by rows, T(n,k) = (n+1)*binomial(n-1,k-1).

Original entry on oeis.org

1, 0, 2, 0, 3, 3, 0, 4, 8, 4, 0, 5, 15, 15, 5, 0, 6, 24, 36, 24, 6, 0, 7, 35, 70, 70, 35, 7, 0, 8, 48, 120, 160, 120, 48, 8, 0, 9, 63, 189, 315, 315, 189, 63, 9, 0, 10, 80, 280, 560, 700, 560, 280, 80, 10, 0, 11, 99, 396, 924, 1386, 1386, 924, 396, 99, 11
Offset: 0

Views

Author

Gary W. Adamson, Feb 09 2007

Keywords

Comments

Row sums = A057711, starting (1, 2, 6, 16, 40, 96, ...).
T(2n,n) gives A033876(n-1) for n > 0. - Alois P. Heinz, Sep 04 2014

Examples

			First few rows of the triangle:
  1;
  0, 2;
  0, 3,  3;
  0, 4,  8,  4;
  0, 5, 15, 15,  5;
  0, 6, 24, 36, 24,  6;
  0, 7, 35, 70, 70, 35,  7;
  ...
		

Crossrefs

Programs

  • Magma
    [[n le 0 select 1 else (n+1)*Binomial(n-1,k-1): k in [0..n]]: n in [0..10]]; // G. C. Greubel, May 05 2018
  • Maple
    T := (n,k) -> (n+1)*binomial(n-1, k-1);
    seq(print(seq(T(n,k),k= 0..n)),n=0..6); # Peter Luschny, Sep 02 2014
  • Mathematica
    Table[(n+1)*Binomial[n-1, k-1], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, May 05 2018 *)
  • PARI
    for(n=0,10, for(k=0,n, print1(if(n==0, 1,(n+1)*binomial(n-1,k-1)), ", "))) \\ G. C. Greubel, May 05 2018
    

Extensions

Name corrected after a suggestion of Joerg Arndt by Peter Luschny, Sep 02 2014

A275324 Expansion of (x*(1-4*x^2)^(-3/2) + (1-4*x^2)^(-1/2) + x + 1)/2.

Original entry on oeis.org

1, 1, 1, 3, 3, 15, 10, 70, 35, 315, 126, 1386, 462, 6006, 1716, 25740, 6435, 109395, 24310, 461890, 92378, 1939938, 352716, 8112468, 1352078, 33801950, 5200300, 140408100, 20058300, 581690700, 77558760, 2404321560, 300540195, 9917826435, 1166803110, 40838108850
Offset: 0

Views

Author

Peter Luschny, Aug 15 2016

Keywords

Crossrefs

Programs

  • Maple
    st := (x*(1-4*x^2)^(-3/2)+(1-4*x^2)^(-1/2)+x+1)/2: series(st,x,36):
    PolynomialTools:-CoefficientList(convert(%,polynom),x);
  • Mathematica
    Table[If[n<2,1,n!/Quotient[n,2]!^2/2], {n,0,30}]
    CoefficientList[Series[(x*(1 - 4*x^2)^(-3/2) + (1 - 4*x^2)^(-1/2) + x + 1)/2, {x, 0, 50}], x] (* G. C. Greubel, Aug 15 2016 *)
  • Sage
    def A275324():
        r, n = 2, 1
        yield 1
        yield 1
        while True:
            n += 1
            r *= 4/n if is_even(n) else n
            yield r // 4
    a = A275324(); print([next(a) for i in range(16)])

Formula

Interweaved from (1+(1-4*x)^(-1/2))/2 (compare A088218 & A001700) and (1+(1-4*x)^(-3/2))/2 (compare A033876).
E.g.f.: (1 + x)*(1 + BesselI(0, 2*x))/2.
For a recurrence see the Sage script.
a(n) = A056040(n)/2 for n>=2.
From Amiram Eldar, Mar 04 2023: (Start)
Sum_{n>=0} 1/a(n) = 2/3 + 16*Pi/(9*sqrt(3)).
Sum_{n>=0} (-1)^n/a(n) = 8/3 - 8*Pi/(9*sqrt(3)). (End)

A296769 Row sums of A296662.

Original entry on oeis.org

1, 7, 38, 187, 874, 3958, 17548, 76627, 330818, 1415650, 6015316, 25413342, 106853668, 447472972, 1867450648, 7770342787, 32248174258, 133530264682, 551793690628, 2276098026922, 9373521044908, 38546133661492, 158301250009768, 649328801880622
Offset: 0

Views

Author

Peter Luschny, Dec 20 2017

Keywords

Crossrefs

Cf. A000531, A296662, bisection of A296663.

Programs

  • Maple
    a := n -> 2^(1 + 2*n)*((2*GAMMA(5/2 + n))/(sqrt(Pi)*GAMMA(2 + n)) - 1):
    seq(a(n), n=0..23);
    # alternative
    A296769 := proc(n)
        2^(2*n)*(doublefactorial(2*n+3)/(1+n)/doublefactorial(2*n)-2) ;
    end proc:
    seq(A296769(n),n=0..10) ; # R. J. Mathar, Jan 03 2018

Formula

a(n) = 2^(1 + 2*n)*((2*Gamma(5/2 + n))/(sqrt(Pi)*Gamma(2 + n)) - 1).
a(n) ~ 4^(n+1)*(sqrt(n/Pi) - 1/2).
a(n) = A000531(n+1). - R. J. Mathar, Jan 03 2018
a(n) = A033876(n)-2^(2*n+1). - R. J. Mathar, Jan 03 2018

A091811 Array read by rows: T(n,k) = binomial(n+k-2,k-1)*binomial(2*n-1,n-k).

Original entry on oeis.org

1, 3, 2, 10, 15, 6, 35, 84, 70, 20, 126, 420, 540, 315, 70, 462, 1980, 3465, 3080, 1386, 252, 1716, 9009, 20020, 24024, 16380, 6006, 924, 6435, 40040, 108108, 163800, 150150, 83160, 25740, 3432, 24310, 175032, 556920, 1021020, 1178100, 875160
Offset: 1

Views

Author

Benoit Cloitre, Mar 18 2004

Keywords

Comments

Alternating sum of elements of n-th row = 1.
If a certain event has a probability p of occurring in any given trial, the probability of its occurring at least n times in 2n-1 trials is Sum_{k=1..n} T(n,k)*(-1)^(k-1)*p^(n+k-1). For example, the probability of its occurring at least 4 out of 7 times is 35p^4 - 84p^5 + 70p^6 - 20p^7. - Matthew Vandermast, Jun 05 2004
With the row polynomial defined as R(n,x) = Sum_{k = 1..n} T(n,k)*x^k, the row polynomial is related to the regularized incomplete Beta function I_x(a,b), through the relation R(n,x) = -(-x)^{-n+1}*I_{-x}(n,n). - Leo C. Stein, Jun 06 2019

Examples

			Triangle starts:
    1,
    3,   2,
   10,  15,   6,
   35,  84,  70,  20,
  126, 420, 540, 315, 70,
  ...
		

Crossrefs

Cf. A001700 (first column), A002740 (second column), A000984 (main diagonal), A033876 (second diagonal), A178792 (row sums).

Programs

  • Magma
    [[Binomial(n+k-2,k-1)*Binomial(2*n-1,n-k): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Jun 15 2015
  • Mathematica
    t[n_, k_] := Binomial[n+k-2, k-1]*Binomial[2n-1, n-k]; Table[t[n, k], {n, 1, 9}, {k, 1, n}] // Flatten (* Jean-François Alcover, Dec 06 2012 *)
  • PARI
    T(x,y)=binomial(x+y-2,y-1)*binomial(2*x-1,x-y)
    

Formula

From Peter Bala, Apr 10 2012: (Start)
O.g.f.: x*t*(1+2*x-sqrt(1-4*t*(x+1)))/(2*(x+t)*sqrt(1-4*t*(x+1))) = x*t + (3*x+2*x^2)*t^2 + (10*x+15*x^2+6*x^3)*t^3 + ....
Sum_{k = 1..n} (-1)^(k-1)*T(n,k)*2^(n-k) = 4^(n-1).
Row polynomial R(n+1,x) = ((2*n+1)!/n!^2)*x*Integral_{y = 0..1} (y*(1+x*y))^n dy. Row sums A178792. (End)

A376484 Array read by ascending antidiagonals: A(n,k)=4^k*Sum_{j=1..n} sin(2*j*Pi/(2*n+1))^(2*k).

Original entry on oeis.org

0, 1, 0, 2, 3, 0, 3, 5, 9, 0, 4, 7, 15, 27, 0, 5, 9, 21, 50, 81, 0, 6, 11, 27, 70, 175, 243, 0, 7, 13, 33, 90, 245, 625, 729, 0, 8, 15, 39, 110, 315, 882, 2250, 2187, 0, 9, 17, 45, 130, 385, 1134, 3234, 8125, 6561, 0, 10, 19, 51, 150, 455, 1386, 4158, 12005, 29375, 19683, 0
Offset: 0

Views

Author

Cheng-Jun Li, Sep 24 2024

Keywords

Comments

It is only a conjecture that A(n,k) is always an integer.
It appears that A(n,k) is divisible by 2*n+1 when n, k are positive integers.

Examples

			For n = 0 to 10 and k = 0 to 10, A(n, k) shows as below :
  0  0  0   0   0    0    0     0      0      0       0
  1  3  9  27  81  243  729  2187   6561  19683   59049
  2  5 15  50 175  625 2250  8125  29375 106250  384375
  3  7 21  70 245  882 3234 12005  44933 169099  638666
  4  9 27  90 315 1134 4158 15444  57915 218781  831222
  5 11 33 110 385 1386 5082 18876  70785 267410 1016158
  6 13 39 130 455 1638 6006 22308  83655 316030 1200914
  7 15 45 150 525 1890 6930 25740  96525 364650 1385670
  8 17 51 170 595 2142 7854 29172 109395 413270 1570426
  9 19 57 190 665 2394 8778 32604 122265 461890 1755182
 10 21 63 210 735 2646 9702 36036 135135 510510 1939938
		

Crossrefs

Conjectures: This array is related to existing sequences as follows: (Start)
Rows: A000004, A000244, A020876, A322459 (with alternate signs).
Columns: A001477, A005408, A016945.
Main Diagonals: A033876.
A(0,k) = A000004, A(1,k) = A000244, A(2,k) = A020876, A(3,k) = (-1)^k * A322459 (First 4 rows of the array).
A(n,0) = A001477(n > 0), A(n,1) = A005408(n > 0), A(n,2) = A016945(n > 0) (First 3 columns of the array).
A(n,n) = A033876(n > 0) (Main diagonal from top left corner). (End)

Programs

  • PARI
    A(n,k) = 4^k*sum(j=1,n,(sin(2*j*Pi/(2*n+1)))^(2*k))

A092370 Triangle read by rows: T(n,k)=(1/2)*C(n+k,k)*C(n,n-k).

Original entry on oeis.org

1, 3, 3, 6, 15, 10, 10, 45, 70, 35, 15, 105, 280, 315, 126, 21, 210, 840, 1575, 1386, 462, 28, 378, 2100, 5775, 8316, 6006, 1716, 36, 630, 4620, 17325, 36036, 42042, 25740, 6435, 45, 990, 9240, 45045, 126126, 210210, 205920, 109395, 24310, 55, 1485, 17160
Offset: 1

Views

Author

Benoit Cloitre, Mar 20 2004

Keywords

Crossrefs

First column = A000217, second column = A050534, main diagonal = A001700, second diagonal = A033876.
Cf. A063007.

Programs

  • Mathematica
    Table[(Binomial[n+k,k]Binomial[n,n-k])/2,{n,10},{k,n}]//Flatten (* Harvey P. Dale, Dec 29 2021 *)
  • PARI
    T(n,k)=(1/2)*binomial(n+k,k)*binomial(n,n-k)

Extensions

Definition corrected by Harvey P. Dale, Dec 29 2021

A242172 a(n) = 2^n*binomial((n + 1 + (n mod 2))/2, 1/2).

Original entry on oeis.org

1, 3, 6, 15, 30, 70, 140, 315, 630, 1386, 2772, 6006, 12012, 25740, 51480, 109395, 218790, 461890, 923780, 1939938, 3879876, 8112468, 16224936, 33801950, 67603900, 140408100, 280816200, 581690700, 1163381400, 2404321560, 4808643120, 9917826435, 19835652870
Offset: 0

Views

Author

Peter Luschny, May 06 2014

Keywords

Crossrefs

Programs

  • Maple
    a := n -> 2^n*binomial((n+1+(n mod 2))/2, 1/2); seq(a(n), n=0..29);
  • Mathematica
    a[n_] := 2^n*Binomial[(n + 1 + Mod[n, 2])/2, 1/2]; Array[a, 33, 0] (* Amiram Eldar, Mar 04 2023 *)

Formula

a(2*n) = A002457(n).
a(2*n+1) = A033876(n).
a(2*n+2)/2 = a(2*n+1).
Conjecture: (n+1)*a(n) -2*a(n-1) +4*(-n-1)*a(n-2)=0. - R. J. Mathar, May 11 2014
a(n) = A100071(n+2)/2. - Michel Marcus, Sep 14 2015
Sum_{n>=0} 1/a(n) = 2*Pi/sqrt(3) - 2. - Amiram Eldar, Mar 04 2023
a(n) = (n+2)*binomial(n+1,ceiling(n/2))/2. - Wesley Ivan Hurt, Nov 23 2023
Showing 1-10 of 12 results. Next