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 16 results. Next

A002450 a(n) = (4^n - 1)/3.

Original entry on oeis.org

0, 1, 5, 21, 85, 341, 1365, 5461, 21845, 87381, 349525, 1398101, 5592405, 22369621, 89478485, 357913941, 1431655765, 5726623061, 22906492245, 91625968981, 366503875925, 1466015503701, 5864062014805, 23456248059221, 93824992236885, 375299968947541
Offset: 0

Views

Author

Keywords

Comments

For n > 0, a(n) is the degree (n-1) "numbral" power of 5 (see A048888 for the definition of numbral arithmetic). Example: a(3) = 21, since the numbral square of 5 is 5(*)5 = 101(*)101(base 2) = 101 OR 10100 = 10101(base 2) = 21, where the OR is taken bitwise. - John W. Layman, Dec 18 2001
a(n) is composite for all n > 2 and has factors x, (3*x + 2*(-1)^n) where x belongs to A001045. In binary the terms greater than 0 are 1, 101, 10101, 1010101, etc. - John McNamara, Jan 16 2002
Number of n X 2 binary arrays with path of adjacent 1's from upper left corner to right column. - R. H. Hardin, Mar 16 2002
The Collatz-function iteration started at a(n), for n >= 1, will end at 1 after 2*n+1 steps. - Labos Elemer, Sep 30 2002 [corrected by Wolfdieter Lang, Aug 16 2021]
Second binomial transform of A001045. - Paul Barry, Mar 28 2003
All members of sequence are also generalized octagonal numbers (A001082). - Matthew Vandermast, Apr 10 2003
Also sum of squares of divisors of 2^(n-1): a(n) = A001157(A000079(n-1)), for n > 0. - Paul Barry, Apr 11 2003
Binomial transform of A000244 (with leading zero). - Paul Barry, Apr 11 2003
Number of walks of length 2n between two vertices at distance 2 in the cycle graph C_6. For n = 2 we have for example 5 walks of length 4 from vertex A to C: ABABC, ABCBC, ABCDC, AFABC and AFEDC. - Herbert Kociemba, May 31 2004
Also number of walks of length 2n + 1 between two vertices at distance 3 in the cycle graph C_12. - Herbert Kociemba, Jul 05 2004
a(n+1) is the number of steps that are made when generating all n-step random walks that begin in a given point P on a two-dimensional square lattice. To make one step means to mark one vertex on the lattice (compare A080674). - Pawel P. Mazur (Pawel.Mazur(AT)pwr.wroc.pl), Mar 13 2005
a(n+1) is the sum of square divisors of 4^n. - Paul Barry, Oct 13 2005
a(n+1) is the decimal number generated by the binary bits in the n-th generation of the Rule 250 elementary cellular automaton. - Eric W. Weisstein, Apr 08 2006
a(n-1) / a(n) = percentage of wasted storage if a single image is stored as a pyramid with a each subsequent higher resolution layer containing four times as many pixels as the previous layer. n is the number of layers. - Victor Brodsky (victorbrodsky(AT)gmail.com), Jun 15 2006
k is in the sequence if and only if C(4k + 1, k) (A052203) is odd. - Paul Barry, Mar 26 2007
This sequence also gives the number of distinct 3-colorings of the odd cycle C(2*n - 1). - Keith Briggs, Jun 19 2007
All numbers of the form m*4^m + (4^m-1)/3 have the property that they are sums of two squares and also their indices are the sum of two squares. This follows from the identity m*4^m + (4^m-1)/3 = 4(4(..4(4m + 1) + 1) + 1) + 1 ..) + 1. - Artur Jasinski, Nov 12 2007
For n > 0, terms are the numbers that, in base 4, are repunits: 1_4, 11_4, 111_4, 1111_4, etc. - Artur Jasinski, Sep 30 2008
Let A be the Hessenberg matrix of order n, defined by: A[1, j] = 1, A[i, i] := 5, (i > 1), A[i, i - 1] = -1, and A[i, j] = 0 otherwise. Then, for n >= 1, a(n) = charpoly(A,1). - Milan Janjic, Jan 27 2010
This is the sequence A(0, 1; 3, 4; 2) = A(0, 1; 4, 0; 1) of the family of sequences [a, b : c, d : k] considered by G. Detlefs, and treated as A(a, b; c, d; k) in the W. Lang link given below. - Wolfdieter Lang, Oct 18 2010
6*a(n) + 1 is every second Mersenne number greater than or equal to M3, hence all Mersenne primes greater than M2 must be a 6*a(n) + 1 of this sequence. - Roderick MacPhee, Nov 01 2010
Smallest number having alternating bit sum n. Cf. A065359.
For n = 1, 2, ..., the last digit of a(n) is 1, 5, 1, 5, ... . - Washington Bomfim, Jan 21 2011
Rule 50 elementary cellular automaton generates this sequence. This sequence also appears in the second column of array in A173588. - Paul Muljadi, Jan 27 2011
Sequence found by reading the line from 0, in the direction 0, 5, ... and the line from 1, in the direction 1, 21, ..., in the square spiral whose edges are the Jacobsthal numbers A001045 and whose vertices are the numbers A000975. These parallel lines are two semi-diagonals in the spiral. - Omar E. Pol, Sep 10 2011
a(n), n >= 1, is also the inverse of 3, denoted by 3^(-1), Modd(2^(2*n - 1)). For Modd n see a comment on A203571. E.g., a(2) = 5, 3 * 5 = 15 == 1 (Modd 8), because floor(15/8) = 1 is odd and -15 == 1 (mod 8). For n = 1 note that 3 * 1 = 3 == 1 (Modd 2) because floor(3/2) = 1 and -3 == 1 (mod 2). The inverse of 3 taken Modd 2^(2*n) coincides with 3^(-1) (mod 2^(2*n)) given in A007583(n), n >= 1. - Wolfdieter Lang, Mar 12 2012
If an AVL tree has a leaf at depth n, then the tree can contain no more than a(n+1) nodes total. - Mike Rosulek, Nov 20 2012
Also, this is the Lucas sequence V(5, 4). - Bruno Berselli, Jan 10 2013
Also, for n > 0, a(n) is an odd number whose Collatz trajectory contains no odd number other than n and 1. - Jayanta Basu, Mar 24 2013
Sum_{n >= 1} 1/a(n) converges to (3*(log(4/3) - QPolyGamma[0, 1, 1/4]))/log(4) = 1.263293058100271... = A321873. - K. G. Stier, Jun 23 2014
Consider n spheres in R^n: the i-th one (i=1, ..., n) has radius r(i) = 2^(1-i) and the coordinates of its center are (0, 0, ..., 0, r(i), 0, ..., 0) where r(i) is in position i. The coordinates of the intersection point in the positive orthant of these spheres are (2/a(n), 4/a(n), 8/a(n), 16/a(n), ...). For example in R^2, circles centered at (1, 0) and (0, 1/2), and with radii 1 and 1/2, meet at (2/5, 4/5). - Jean M. Morales, May 19 2015
From Peter Bala, Oct 11 2015: (Start)
a(n) gives the values of m such that binomial(4*m + 1,m) is odd. Cf. A003714, A048716, A263132.
2*a(n) = A020988(n) gives the values of m such that binomial(4*m + 2, m) is odd.
4*a(n) = A080674(n) gives the values of m such that binomial(4*m + 4, m) is odd. (End)
Collatz Conjecture Corollary: Except for powers of 2, the Collatz iteration of any positive integer must eventually reach a(n) and hence terminate at 1. - Gregory L. Simay, May 09 2016
Number of active (ON, black) cells at stage 2^n - 1 of the two-dimensional cellular automaton defined by "Rule 598", based on the 5-celled von Neumann neighborhood. - Robert Price, May 16 2016
From Luca Mariot and Enrico Formenti, Sep 26 2016: (Start)
a(n) is also the number of coprime pairs of polynomials (f, g) over GF(2) where both f and g have degree n + 1 and nonzero constant term.
a(n) is also the number of pairs of one-dimensional binary cellular automata with linear and bipermutive local rule of neighborhood size n+1 giving rise to orthogonal Latin squares of order 2^m, where m is a multiple of n. (End)
Except for 0, 1 and 5, all terms are Brazilian repunits numbers in base 4, and so belong to A125134. For n >= 3, all these terms are composite because a(n) = {(2^n-1) * (2^n + 1)}/3 and either (2^n - 1) or (2^n + 1) is a multiple of 3. - Bernard Schott, Apr 29 2017
Given the 3 X 3 matrix A = [2, 1, 1; 1, 2, 1; 1, 1, 2] and the 3 X 3 unit matrix I_3, A^n = a(n)(A - I_3) + I_3. - Nicolas Patrois, Jul 05 2017
The binary expansion of a(n) (n >= 1) consists of n 1's alternating with n - 1 0's. Example: a(4) = 85 = 1010101_2. - Emeric Deutsch, Aug 30 2017
a(n) (n >= 1) is the viabin number of the integer partition [n, n - 1, n - 2, ..., 2, 1] (for the definition of viabin number see comment in A290253). Example: a(4) = 85 = 1010101_2; consequently, the southeast border of the Ferrers board of the corresponding integer partition is ENENENEN, where E = (1, 0), N = (0, 1); this leads to the integer partition [4, 3, 2, 1]. - Emeric Deutsch, Aug 30 2017
Numbers whose binary and Gray-code representations are both palindromes (i.e., intersection of A006995 and A281379). - Amiram Eldar, May 17 2021
Starting with n = 1 the sequence satisfies {a(n) mod 6} = repeat{1, 5, 3}. - Wolfdieter Lang, Jan 14 2022
Terms >= 5 are those q for which the multiplicative order of 2 mod q is floor(log_2(q)) + 2 (and which is 1 more than the smallest possible order for any q). - Tim Seuré, Mar 09 2024
The order of 2 modulo a(n) is 2*n for n >= 2. - Joerg Arndt, Mar 09 2024

Examples

			Apply Collatz iteration to 9: 9, 28, 14, 7, 22, 11, 34, 17, 52, 26, 13, 40, 20, 10, 5 and hence 16, 8, 4, 2, 1.
Apply Collatz iteration to 27: 27, 82, 41, 124, 62, 31, 94, 47, 142, 71, 214, 107, 322, 161, 484, 242, 121, 364, 182, 91, 274, 137, 412, 206, 103, 310, 155, 466, 233, 700, 350, 175, 526, 263, 790, 395, 1186, 593, 1780, 890, 445, 1336, 668, 334, 167, 502, 251, 754, 377, 1132, 566, 283, 850, 425, 1276, 638, 319, 958, 479, 1438, 719, 2158, 1079, 3238, 1619, 4858, 2429, 7288, 3644, 1822, 911, 2734, 1367, 4102, 2051, 6154, 3077, 9232, 4616, 2308, 1154, 577, 1732, 866, 433, 1300, 650, 325, 976, 488, 244, 122, 61, 184, 92, 46, 23, 70, 35, 106, 53, 160, 80, 40, 20, 10, 5 and hence 16, 8, 4, 2, 1. [Corrected by _Sean A. Irvine_ at the suggestion of Stephen Cornelius, Mar 04 2024]
a(5) = (4^5 - 1)/3 = 341 = 11111_4 = {(2^5 - 1) * (2^5 + 1)}/3 = 31 * 33/3 = 31 * 11. - _Bernard Schott_, Apr 29 2017
		

References

  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 112.
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 217.
  • 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

Partial sums of powers of 4, A000302.
When converted to binary, this gives A094028.
Subsequence of A003714.
Primitive factors: A129735.

Programs

  • GAP
    List([0..25], n -> (4^n-1)/3); # Muniru A Asiru, Feb 18 2018
    
  • Haskell
    a002450 = (`div` 3) . a024036
    a002450_list = iterate ((+ 1) . (* 4)) 0
    -- Reinhard Zumkeller, Oct 03 2012
    
  • Magma
    [ (4^n-1)/3: n in [0..25] ]; // Klaus Brockhaus, Oct 28 2008
    
  • Magma
    [n le 2 select n-1 else 5*Self(n-1)-4*Self(n-2): n in [1..70]]; // Vincenzo Librandi, Jun 13 2015
    
  • Maple
    [seq((4^n-1)/3,n=0..40)];
    A002450:=1/(4*z-1)/(z-1); # Simon Plouffe in his 1992 dissertation, dropping the initial zero
  • Mathematica
    Table[(4^n - 1)/3, {n, 0, 127}] (* Vladimir Joseph Stephan Orlovsky, Sep 29 2008 *)
    LinearRecurrence[{5, -4}, {0, 1}, 30] (* Harvey P. Dale, Jun 23 2013 *)
  • Maxima
    makelist((4^n-1)/3, n, 0, 30); /* Martin Ettl, Nov 05 2012 */
    
  • PARI
    a(n) = (4^n-1)/3;
    
  • PARI
    my(z='z+O('z^40)); Vec(z/((1-z)*(1-4*z))) \\ Altug Alkan, Oct 11 2015
    
  • Python
    def A002450(n): return ((1<<(n<<1))-1)//3 # Chai Wah Wu, Jan 29 2023
  • Scala
    ((List.fill(20)(4: BigInt)).scanLeft(1: BigInt)( * )).scanLeft(0: BigInt)( + ) // Alonso del Arte, Sep 17 2019
    

Formula

From Wolfdieter Lang, Apr 24 2001: (Start)
a(n+1) = Sum_{m = 0..n} A060921(n, m).
G.f.: x/((1-x)*(1-4*x)). (End)
a(n) = Sum_{k = 0..n-1} 4^k; a(n) = A001045(2*n). - Paul Barry, Mar 17 2003
E.g.f.: (exp(4*x) - exp(x))/3. - Paul Barry, Mar 28 2003
a(n) = (A007583(n) - 1)/2. - N. J. A. Sloane, May 16 2003
a(n) = A000975(2*n)/2. - N. J. A. Sloane, Sep 13 2003
a(n) = A084160(n)/2. - N. J. A. Sloane, Sep 13 2003
a(n+1) = 4*a(n) + 1, with a(0) = 0. - Philippe Deléham, Feb 25 2004
a(n) = Sum_{i = 0..n-1} C(2*n - 1 - i, i)*2^i. - Mario Catalani (mario.catalani(AT)unito.it), Jul 23 2004
a(n+1) = Sum_{k = 0..n} binomial(n+1, k+1)*3^k. - Paul Barry, Aug 20 2004
a(n) = center term in M^n * [1 0 0], where M is the 3 X 3 matrix [1 1 1 / 1 3 1 / 1 1 1]. M^n * [1 0 0] = [A007583(n-1) a(n) A007583(n-1)]. E.g., a(4) = 85 since M^4 * [1 0 0] = [43 85 43] = [A007583(3) a(4) A007583(3)]. - Gary W. Adamson, Dec 18 2004
a(n) = Sum_{k = 0..n, j = 0..n} C(n, j)*C(j, k)*A001045(j - k). - Paul Barry, Feb 15 2005
a(n) = Sum_{k = 0..n} C(n, k)*A001045(n-k)*2^k = Sum_{k = 0..n} C(n, k)*A001045(k)*2^(n-k). - Paul Barry, Apr 22 2005
a(n) = A125118(n, 3) for n > 2. - Reinhard Zumkeller, Nov 21 2006
a(n) = Sum_{k = 0..n} 2^(n - k)*A128908(n, k), n >= 1. - Philippe Deléham, Oct 19 2008
a(n) = Sum_{k = 0..n} A106566(n, k)*A100335(k). - Philippe Deléham, Oct 30 2008
If we define f(m, j, x) = Sum_{k = j..m} binomial(m, k)*stirling2(k, j)*x^(m - k) then a(n-1) = f(2*n, 4, -2), n >= 2. - Milan Janjic, Apr 26 2009
a(n) = A014551(n) * A001045(n). - R. J. Mathar, Jul 08 2009
a(n) = 4*a(n-1) + a(n-2) - 4*a(n-3) = 5*a(n-1) - 4*a(n-2), a(0) = 0, a(1) = 1, a(2) = 5. - Wolfdieter Lang, Oct 18 2010
a(0) = 0, a(n+1) = a(n) + 2^(2*n). - Washington Bomfim, Jan 21 2011
A036555(a(n)) = 2*n. - Reinhard Zumkeller, Jan 28 2011
a(n) = Sum_{k = 1..floor((n+2)/3)} C(2*n + 1, n + 2 - 3*k). - Mircea Merca, Jun 25 2011
a(n) = Sum_{i = 1..n} binomial(2*n + 1, 2*i)/3. - Wesley Ivan Hurt, Mar 14 2015
a(n+1) = 2^(2*n) + a(n), a(0) = 0. - Ben Paul Thurston, Dec 27 2015
a(k*n)/a(n) = 1 + 4^n + ... + 4^((k-1)*n). - Gregory L. Simay, Jun 09 2016
Dirichlet g.f.: (PolyLog(s, 4) - zeta(s))/3. - Ilya Gutkovskiy, Jun 26 2016
A000120(a(n)) = n. - André Dalwigk, Mar 26 2018
a(m) divides a(m*n), in particular: a(2*n) == 0 (mod 5), a(3*n) == 0 (mod 3*7), a(5*n) == 0 (mod 11*31), etc. - M. F. Hasler, Oct 19 2018
a(n) = 4^(n-1) + a(n-1). - Bob Selcoe, Jan 01 2020
a(n) = A178415(1, n) = A347834(1, n-1), arrays, for n >= 1. - Wolfdieter Lang, Nov 29 2021
a(n) = A000225(2*n)/3. - John Keith, Jan 22 2022
a(n) = A080674(n) + 1 = A047849(n) - 1 = A163834(n) - 2 = A155701(n) - 3 = A163868(n) - 4 = A156605(n) - 7. - Ray Chandler, Jun 16 2023
From Peter Bala, Jul 23 2025: (Start)
The following are examples of telescoping products. Cf. A016153:
Product_{k = 1..2*n} 1 + 2^k/a(k+1) = a(n+1)/A007583(n) = (4^(n+1) - 1)/(2*4^n + 1).
Hence, Product_{k >= 1} 1 + 2^k/a(k+1) = 2.
Product_{k >= 1} 1 - 2^k/a(k+1) = 2/5, since 1 - 2^n/a(n+1) = b(n)/b(n-1), where b(n) = 2 - 3/(1 - 2^(n+1)).
Product_{k >= 1} 1 + (-2)^k/a(k+1) = 2/3, since 1 + (-2)^n/a(n+1) = c(n)/c(n-1), where c(n) = 2 - 1/(1 + (-2)^(n+1)).
Product_{k >= 1} 1 - (-2)^k/a(k+1) = 6/5, since 1 - (-2)^n/a(n+1) = d(n)/d(n-1), where d(n) = 2 - 1/(1 - (-2)^(n+1)). (End)

A111418 Right-hand side of odd-numbered rows of Pascal's triangle.

Original entry on oeis.org

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

Views

Author

Philippe Deléham, Nov 13 2005

Keywords

Comments

Riordan array (c(x)/sqrt(1-4*x),x*c(x)^2) where c(x) is g.f. of A000108. Unsigned version of A113187. Diagonal sums are A014301(n+1).
Triangle T(n,k),0<=k<=n, read by rows defined by :T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=3*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+2*T(n-1,k)+T(n-1,k+1) for k>=1. - Philippe Deléham, Mar 22 2007
Reversal of A122366. - Philippe Deléham, Mar 22 2007
Column k has e.g.f. exp(2x)(Bessel_I(k,2x)+Bessel_I(k+1,2x)). - Paul Barry, Jun 06 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1 . Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Diagonal sums are A014301(n+1). - Paul Barry, Mar 08 2011
This triangle T(n,k) appears in the expansion of odd powers of Fibonacci numbers F=A000045 in terms of F-numbers with multiples of odd numbers as indices. See the Ozeki reference, p. 108, Lemma 2. The formula is: F_l^(2*n+1) = sum(T(n,k)*(-1)^((n-k)*(l+1))* F_{(2*k+1)*l}, k=0..n)/5^n, n >= 0, l >= 0. - Wolfdieter Lang, Aug 24 2012
Central terms give A052203. - Reinhard Zumkeller, Mar 14 2014
This triangle appears in the expansion of (4*x)^n in terms of the polynomials Todd(n, x):= T(2*n+1, sqrt(x))/sqrt(x) = sum(A084930(n,m)*x^m), n >= 0. This follows from the inversion of the lower triangular Riordan matrix built from A084930 and comparing the g.f. of the row polynomials. - Wolfdieter Lang, Aug 05 2014
From Wolfdieter Lang, Aug 15 2014: (Start)
This triangle is the inverse of the signed Riordan triangle (-1)^(n-m)*A111125(n,m).
This triangle T(n,k) appears in the expansion of x^n in terms of the polynomials todd(k, x):= T(2*k+1, sqrt(x)/2)/(sqrt(x)/2) = S(k, x-2) - S(k-1, x-2) with the row polynomials T and S for the triangles A053120 and A049310, respectively: x^n = sum(T(n,k)*todd(k, x), k=0..n). Compare this with the preceding comment.
The A- and Z-sequences for this Riordan triangle are [1, 2, 1, repeated 0] and [3, 1, repeated 0]. For A- and Z-sequences for Riordan triangles see the W. Lang link under A006232. This corresponds to the recurrences given in the Philippe Deléham, Mar 22 2007 comment above. (End)

Examples

			From _Wolfdieter Lang_, Aug 05 2014: (Start)
The triangle T(n,k) begins:
n\k      0      1      2      3     4     5    6    7   8  9  10 ...
0:       1
1:       3      1
2:      10      5      1
3:      35     21      7      1
4:     126     84     36      9     1
5:     462    330    165     55    11     1
6:    1716   1287    715    286    78    13    1
7:    6435   5005   3003   1365   455   105   15    1
8:   24310  19448  12376   6188  2380   680  136   17   1
9:   92378  75582  50388  27132 11628  3876  969  171  19  1
10: 352716 293930 203490 116280 54264 20349 5985 1330 210 21   1
...
Expansion examples (for the Todd polynomials see A084930 and a comment above):
(4*x)^2 = 10*Todd(n,  0) + 5*Todd(n, 1) + 1*Todd(n, 2) = 10*1 + 5*(-3 + 4*x) + 1*(5 - 20*x + 16*x^2).
(4*x)^3 =  35*1 + 21*(-3 + 4*x) + 7*(5 - 20*x + 16*x^2) + (-7 + 56*x - 112*x^2 +64*x^3)*1. (End)
---------------------------------------------------------------------
Production matrix is
3, 1,
1, 2, 1,
0, 1, 2, 1,
0, 0, 1, 2, 1,
0, 0, 0, 1, 2, 1,
0, 0, 0, 0, 1, 2, 1,
0, 0, 0, 0, 0, 1, 2, 1,
0, 0, 0, 0, 0, 0, 1, 2, 1,
0, 0, 0, 0, 0, 0, 0, 1, 2, 1
- _Paul Barry_, Mar 08 2011
Application to odd powers of Fibonacci numbers F, row n=2:
F_l^5 = (10*(-1)^(2*(l+1))*F_l + 5*(-1)^(1*(l+1))*F_{3*l} + 1*F_{5*l})/5^2, l >= 0. - _Wolfdieter Lang_, Aug 24 2012
		

Crossrefs

Programs

  • Haskell
    a111418 n k = a111418_tabl !! n !! k
    a111418_row n = a111418_tabl !! n
    a111418_tabl = map reverse a122366_tabl
    -- Reinhard Zumkeller, Mar 14 2014
  • Mathematica
    Table[Binomial[2*n+1, n-k], {n,0,10}, {k,0,n}] (* G. C. Greubel, May 22 2017 *)
    T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0,
    T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]];
    Table[T[n, k, 3, 2], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, May 22 2017 *)

Formula

T(n, k) = C(2*n+1, n-k).
Sum_{k=0..n} T(n, k) = 4^n.
Sum_{k, 0<=k<=n}(-1)^k *T(n,k) = binomial(2*n,n) = A000984(n). - Philippe Deléham, Mar 22 2007
T(n,k) = sum{j=k..n, C(n,j)*2^(n-j)*C(j,floor((j-k)/2))}. - Paul Barry, Jun 06 2007
Sum_{k, k>=0} T(m,k)*T(n,k) = T(m+n,0)= A001700(m+n). - Philippe Deléham, Nov 22 2009
G.f. row polynomials: ((1+x) - (1-x)/sqrt(1-4*z))/(2*(x - (1+x)^2*z))
(see the Riordan property mentioned in a comment above). - Wolfdieter Lang, Aug 05 2014

A262977 a(n) = binomial(4*n-1,n).

Original entry on oeis.org

1, 3, 21, 165, 1365, 11628, 100947, 888030, 7888725, 70607460, 635745396, 5752004349, 52251400851, 476260169700, 4353548972850, 39895566894540, 366395202809685, 3371363686069236, 31074067324187580, 286845713747883300, 2651487106659130740, 24539426037817994160
Offset: 0

Views

Author

Vladimir Kruchinin, Oct 06 2015

Keywords

Comments

From Gus Wiseman, Sep 28 2022: (Start)
Also the number of integer compositions of 4n with alternating sum 2n, where the alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. These compositions are ranked by A348614. The a(12) = 21 compositions are:
(6,2) (1,2,5) (1,1,5,1) (1,1,1,1,4)
(2,2,4) (2,1,4,1) (1,1,2,1,3)
(3,2,3) (3,1,3,1) (1,1,3,1,2)
(4,2,2) (4,1,2,1) (1,1,4,1,1)
(5,2,1) (5,1,1,1) (2,1,1,1,3)
(2,1,2,1,2)
(2,1,3,1,1)
(3,1,1,1,2)
(3,1,2,1,1)
(4,1,1,1,1)
The following pertain to this interpretation:
- The case of partitions is A000712, reverse A006330.
- Allowing any alternating sum gives A013777 (compositions of 4n).
- A011782 counts compositions of n.
- A034871 counts compositions of 2n with alternating sum 2k.
- A097805 counts compositions by alternating (or reverse-alternating) sum.
- A103919 counts partitions by sum and alternating sum (reverse: A344612).
- A345197 counts compositions by length and alternating sum.
(End)

Crossrefs

Programs

  • Magma
    [Binomial(4*n-1,n): n in [0..20]]; // Vincenzo Librandi, Oct 06 2015
    
  • Mathematica
    Table[Binomial[4 n - 1, n], {n, 0, 40}] (* Vincenzo Librandi, Oct 06 2015 *)
  • Maxima
    B(x):=sum(binomial(4*n-1,n-1)*3/(4*n-1)*x^n,n,1,30);
    taylor(x*diff(B(x),x,1)/B(x),x,0,20);
    
  • PARI
    a(n) = binomial(4*n-1,n); \\ Michel Marcus, Oct 06 2015

Formula

G.f.: A(x)=x*B'(x)/B(x), where B(x) if g.f. of A006632.
a(n) = Sum_{k=0..n}(binomial(n-1,n-k)*binomial(3*n,k)).
a(n) = 3*A224274(n), for n > 0. - Michel Marcus, Oct 12 2015
From Peter Bala, Nov 04 2015: (Start)
The o.g.f. equals f(x)/g(x), where f(x) is the o.g.f. for A005810 and g(x) is the o.g.f. for A002293. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(4*n + k,n). Cf. A005810 (k = 0), A052203 (k = 1), A257633 (k = 2), A224274 (k = 3) and A004331 (k = 4). (End)
a(n) = [x^n] 1/(1 - x)^(3*n). - Ilya Gutkovskiy, Oct 03 2017
a(n) = A071919(3n-1,n+1) = A097805(4n,n+1). - Gus Wiseman, Sep 28 2022
From Peter Bala, Feb 14 2024: (Start)
a(n) = (-1)^n * binomial(-3*n, n).
a(n) = hypergeom([1 - 3*n, -n], [1], 1).
The g.f. A(x) satisfies A(x/(1 + x)^4) = 1/(1 - 3*x). (End)
a(n) = Sum_{k = 0..n} binomial(2*n+k-1, k)*binomial(2*n-k-1, n-k). - Peter Bala, Sep 16 2024
G.f.: 1/(4-3*g) where g = 1+x*g^4 is the g.f. of A002293. - Seiichi Manyama, Aug 17 2025

A224274 a(n) = binomial(4*n,n)/4.

Original entry on oeis.org

1, 7, 55, 455, 3876, 33649, 296010, 2629575, 23535820, 211915132, 1917334783, 17417133617, 158753389900, 1451182990950, 13298522298180, 122131734269895, 1123787895356412, 10358022441395860, 95615237915961100, 883829035553043580, 8179808679272664720, 75788358475481302185
Offset: 1

Views

Author

Gary Detlefs, Apr 02 2013

Keywords

Comments

In general, binomial(k*n,n)/k = binomial(k*n-1,n-1).
Sequences in the OEIS related to this identity are:
. C(2n,n) = A000984, C(2n,n)/2 = A001700;
. C(3n,n) = A005809, C(3n,n)/3 = A025174;
. C(4n,n) = A005810, C(4n,n)/4 = a(n);
. C(5n,n) = A001449, C(5n,n)/5 = A163456;
. C(6n,n) = A004355, C(6n,n)/6 is not in the OEIS.
Conjecture: a(n) == 1 (mod n^3) iff n is an odd prime.
It is known that a(p) == 1(mod p^3) for prime p >= 3. See Mestrovic, Section 3. - Peter Bala, Oct 09 2015

Examples

			For n=2, binomial(4*n,n) = binomial(8,2) = 8*7/2 = 28, so a(2) = 28/4 = 7. - _Michael B. Porter_, Jul 12 2016
		

Crossrefs

Programs

  • Magma
    [Binomial(4*n,n) div 4: n in [1..25]]; // Vincenzo Librandi, Jun 03 2015
  • Maple
    seq(binomial(4*n,n)/4, n=1..17);
  • Mathematica
    Table[Binomial[4 n, n]/4, {n, 30}] (* Vincenzo Librandi, Jun 03 2015 *)
  • PARI
    a(n) = binomial(4*n,n)/4; /* Joerg Arndt, Apr 02 2013 */
    

Formula

a(n) = binomial(4*n,n)/4 = A005810(n)/4.
a(n) = binomial(4*n-1,n-1).
G.f.: A(x) = B'(x)/B(x), where B(x) = 1 + x*B(x)^4 is g.f. of A002293. - Vladimir Kruchinin, Aug 13 2015
From Peter Bala, Oct 08 2015: (Start)
a(n) = 1/2*[x^n] (C(x)^2)^n, where C(x) = (1 - sqrt(1 - 4*x))/(2*x) is the o.g.f. for the Catalan numbers A000108. Cf. A163456.
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + x + 4*x^2 + 22*x^3 + ... is the o.g.f. for A002293.
exp( 2*Sum_{n >= 1} a(n)*x^n/n ) = 1 + 2*x + 9*x^2 + 52*x^3 + ... is the o.g.f. for A069271. (End)
From Peter Bala, Nov 04 2015: (Start)
With an offset of 1, the o.g.f. equals f(x)*g(x)^3, where f(x) is the o.g.f. for A005810 and g(x) is the o.g.f. for A002293. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(4*n + k,n). Cf. A262977 (k = -1), A005810 (k = 0), A052203 (k = 1), A257633 (k = 2) and A004331 (k = 4). (End)
a(n) = 1/5*[x^n] (1 + x)/(1 - x)^(3*n + 1) = 1/5*[x^n]( 1/C(-x) )^(5*n), where C(x) = (1 - sqrt(1 - 4*x))/(2*x) is the o.g.f. for the Catalan numbers A000108. Cf. A227726. - Peter Bala, Jul 12 2016
a(n) ~ 2^(8*n-3/2)*3^(-3*n-1/2)*n^(-1/2)/sqrt(Pi). - Ilya Gutkovskiy, Jul 12 2016
O.g.f.: A(x) = f(x)/(1 - 3*f(x)), where f(x) = series reversion (x/(1 + x)^4) = x + 4*x^2 + 22*x^3 + 140*x^4 + 969*x^5 + ... is the o.g.f. of A002293 with the initial term omitted. Cf. A025174. - Peter Bala, Feb 03 2022
Right-hand side of the identities (1/3)*Sum_{k = 0..n} (-1)^(n+k)*C(x*n,n-k)*C((x+3)*n+k-1,k) = C(4*n,n)/4 and (1/4)*Sum_{k = 0..n} (-1)^k*C(x*n,n-k)*C((x-4)*n+k-1,k) = C(4*n,n)/4, both valid for n >= 1 and x arbitrary. - Peter Bala, Feb 28 2022
Right-hand side of the identity (1/3)*Sum_{k = 0..2*n} (-1)^k*binomial(5*n-k-1,2*n-k)*binomial(3*n+k-1,k) = binomial(4*n,n)/4. - Peter Bala, Mar 09 2022
a(n) = [x^n] G(x)^n, where G(x) = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + ... is the g.f. of A001764. - Peter Bala, Oct 17 2024

A004331 Binomial coefficient C(4n,n-1).

Original entry on oeis.org

1, 8, 66, 560, 4845, 42504, 376740, 3365856, 30260340, 273438880, 2481256778, 22595200368, 206379406870, 1889912732400, 17345898649800, 159518999862720, 1469568786235308, 13559593014190944, 125288932441604200
Offset: 1

Views

Author

Keywords

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. 828.

Crossrefs

Programs

  • Maple
    #A004331
    seq(binomial(4*n - 1,n), n = 0..20);
  • Mathematica
    a[n_] := Binomial[4*n, n - 1]; Array[a, 19] (* Amiram Eldar, May 09 2020 *)
  • PARI
    vector(30, n, binomial(4*n, n-1)) \\ Altug Alkan, Nov 05 2015

Formula

G.f.: (g^2-g)/(4-3*g) where g = 1+x*g^4 is the g.f. of A002293. - Mark van Hoeij, Nov 11 2011
With an offset of 0, the o.g.f. equals f(x)*g(x)^4, where f(x) is the o.g.f. for A005810 and g(x) is the o.g.f. for A002293. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(4*n + k,n). Cf. A262977 (k = -1), A005810 (k = 0), A052203 (k = 1), A257633 (k = 2) and A224274 (k = 3). - Peter Bala, Nov 04 2015
D-finite with recurrence 3*(n-1)*(3*n-1)*(3*n+1)*a(n) -8*(4*n-3)*(2*n-1)*(4*n-1)*a(n-1)=0. - R. J. Mathar, Mar 19 2025

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.

A079589 a(n) = C(5*n+1,n).

Original entry on oeis.org

1, 6, 55, 560, 5985, 65780, 736281, 8347680, 95548245, 1101716330, 12777711870, 148902215280, 1742058970275, 20448884000160, 240719591939480, 2840671544105280, 33594090947249085, 398039194165652550, 4724081931321677925, 56151322242892212960, 668324943343021950370
Offset: 0

Views

Author

Benoit Cloitre, Jan 26 2003

Keywords

Comments

a(n) is the number of paths from (0,0) to (5n,n) taking north and east steps while avoiding exactly 2 consecutive north steps. - Shanzhen Gao, Apr 15 2010

Crossrefs

Programs

  • Magma
    [Binomial(5*n+1, n): n in [0..20]]; // Vincenzo Librandi, Aug 07 2014
  • Maple
    seq(binomial(5*n+1,n),n=0..100); # Robert Israel, Aug 07 2014
  • Mathematica
    Table[Binomial[5n+1,n],{n,0,20}]  (* Harvey P. Dale, Jan 23 2011 *)

Formula

a(n) is asymptotic to c*(3125/256)^n/sqrt(n) with c=0.557.... [c = 5^(3/2)/(sqrt(Pi)*2^(7/2)) = 0.55753878629774... - Vaclav Kotesovec, Feb 14 2019 and Aug 20 2025]
8*n*(4*n+1)*(2*n-1)*(4*n-1)*a(n) -5*(5*n+1)*(5*n-3)*(5*n-2)*(5*n-1)*a(n-1)=0. - R. J. Mathar, Jul 17 2014
G.f.: hypergeom([2/5, 3/5, 4/5, 6/5], [1/2, 3/4, 5/4], (3125/256)*x). - Robert Israel, Aug 07 2014
a(n) = [x^n] 1/(1 - x)^(2*(2*n+1)). - Ilya Gutkovskiy, Oct 10 2017
From Seiichi Manyama, Aug 16 2025: (Start)
a(n) = Sum_{k=0..n} binomial(5*n-k,n-k).
G.f.: 1/(1 - x*g^3*(5+g)) where g = 1+x*g^5 is the g.f. of A002294.
G.f.: g^2/(5-4*g) where g = 1+x*g^5 is the g.f. of A002294.
G.f.: B(x)^2/(1 + 4*(B(x)-1)/5), where B(x) is the g.f. of A001449. (End)

A257633 a(n) = binomial(4*n + 2,n).

Original entry on oeis.org

1, 6, 45, 364, 3060, 26334, 230230, 2035800, 18156204, 163011640, 1471442973, 13340783196, 121399651100, 1108176102180, 10142940735900, 93052749919920, 855420636763836, 7877932561061640, 72667580816130436, 671262558647881200, 6208770443303347920
Offset: 0

Views

Author

Peter Bala, Nov 04 2015

Keywords

Crossrefs

Programs

  • Maple
    #A257633
    seq(binomial(4*n + 2,n), n = 0..20);
  • Mathematica
    Table[Binomial[4*n + 2, n], {n, 0, 120}] (* Michael De Vlieger, Apr 11 2025 *)
  • PARI
    vector(30, n, n--; binomial(4*n+2, n)) \\ Altug Alkan, Nov 05 2015

Formula

The o.g.f. equals f(x)*g(x)^2, where f(x) is the o.g.f. for A005810 and g(x) is the o.g.f. for A002293. More generally, f(x)*g(x)^k is the o.g.f. for the sequence binomial(4*n + k,n). Cf. A262977 (k = -1), A005810 (k = 0), A052203 (k = 1), A224274 (k = 3) and A004331 (k = 4).

A079590 a(n) = C(6*n+1,n).

Original entry on oeis.org

1, 7, 78, 969, 12650, 169911, 2324784, 32224114, 450978066, 6358402050, 90177170226, 1285063345176, 18385569737808, 263926640438545, 3799541229226200, 54834293825867556, 793067310934426018
Offset: 0

Views

Author

Benoit Cloitre, Jan 26 2003

Keywords

Comments

a(n) is the number of paths from (0,0) to (6n,n) taking north and east steps and avoiding north^{=2}. - Shanzhen Gao, Apr 15 2010

Crossrefs

Cf. A052203.

Programs

  • Magma
    [Binomial(6*n+1, n): n in [0..20]]; // Vincenzo Librandi, Aug 07 2014
  • Maple
    seq(binomial(6*n+1,n), n = 1 .. 100); # Robert Israel, Aug 07 2014
  • Mathematica
    Table[Binomial[6 n + 1, n], {n, 0, 20}] (* Vincenzo Librandi, Aug 07 2014 *)

Formula

a(n) is asymptotic to c*sqrt(n)*(46656/3125)^n with c= 0.5243... [c = 2*3^(3/2)/(5^(3/2)*sqrt(Pi)) = 0.52442324668419795386... - Vaclav Kotesovec, Feb 14 2019]
G.f.: hypergeom([1/3, 1/2, 2/3, 5/6, 7/6], [2/5, 3/5, 4/5, 6/5], (46656/3125)*x). - Robert Israel, Aug 07 2014

A117178 Riordan array (c(x^2)/sqrt(1-4*x^2), x*c(x^2)), c(x) the g.f. of A000108.

Original entry on oeis.org

1, 0, 1, 3, 0, 1, 0, 4, 0, 1, 10, 0, 5, 0, 1, 0, 15, 0, 6, 0, 1, 35, 0, 21, 0, 7, 0, 1, 0, 56, 0, 28, 0, 8, 0, 1, 126, 0, 84, 0, 36, 0, 9, 0, 1, 0, 210, 0, 120, 0, 45, 0, 10, 0, 1, 462, 0, 330, 0, 165, 0, 55, 0, 11, 0, 1, 0, 792, 0, 495, 0, 220, 0, 66, 0, 12, 0, 1
Offset: 0

Views

Author

Paul Barry, Mar 01 2006

Keywords

Comments

Row sums are A058622(n+1). Diagonal sums are A001791(n+1), with interpolated zeros. Inverse is A117179.
De-aerated and rows reversed, this matrix apparently becomes A014462. The nonzero antidiagonals are embedded in several entries and apparently contain partial sums of previous nonzero antidiagonals. - Tom Copeland, May 30 2017

Examples

			Triangle begins
    1;
    0,  1;
    3,  0,  1;
    0,  4,  0,  1;
   10,  0,  5,  0,  1;
    0, 15,  0,  6,  0,  1;
   35,  0, 21,  0,  7,  0,  1;
    0, 56,  0, 28,  0,  8,  0,  1;
  126,  0, 84,  0, 36,  0,  9,  0,  1;
		

Crossrefs

Programs

  • Magma
    [(1+(-1)^(n-k))*Binomial(n+1, Floor((n-k)/2))/2: k in [0..n], n in [0..15]]; // G. C. Greubel, Aug 08 2022
    
  • Mathematica
    T[n_, k_]:= Binomial[n+1, (n-k)/2]*(1+(-1)^(n-k))/2;
    Table[T[n, k], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Aug 08 2022 *)
  • SageMath
    def A117178(n,k): return (1 + (-1)^(n-k))*binomial(n+1, (n-k)//2)/2
    flatten([[A117178(n,k) for k in (0..n)] for n in (0..15)]) # G. C. Greubel, Aug 08 2022

Formula

T(n,k) = C(n+1, (n-k)/2) * (1 + (-1)^(n-k))/2.
Column k has e.g.f. Bessel_I(k,2x) + Bessel_I(k+2, 2x).
From G. C. Greubel, Aug 08 2022: (Start)
Sum_{k=0..n} T(n, k) = A058622(n+1).
Sum_{k=0..floor(n/2)} T(n-k, k) = ((1+(-1)^n)/2) * A001791((n+2)/2).
T(2*n, n) = ((1+(-1)^n)/2) * A052203(n/2).
T(2*n+1, n) = ((1-(-1)^n)/2) * A224274((n+1)/2).
T(2*n-1, n-1) = ((1+(-1)^n)/2) * A224274(n/2). (End)
Showing 1-10 of 16 results. Next