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-7 of 7 results.

A126932 Binomial transform of A127358.

Original entry on oeis.org

1, 4, 15, 55, 199, 714, 2547, 9048, 32043, 113212, 399265, 1406079, 4946137, 17383162, 61048359, 214270215, 751691811, 2636004228, 9240836733, 32386215981, 113478349989, 397544907486, 1392493797765, 4876916883090, 17078574481941, 59802541979964
Offset: 0

Views

Author

Philippe Deléham, Mar 17 2007

Keywords

Comments

Hankel transform is (-1)^n.
Row sums of the Riordan array ((1-2*x)/(1+x+x^2), x/(1+x+x^2))^(-1). - Paul Barry, Nov 06 2008

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( (Sqrt(1-2*x-3*x^2) + 3*(1-3*x))/(2*(2-13*x+21*x^2)) )); // G. C. Greubel, Jan 29 2020
    
  • Maple
    seq(coeff(series( (sqrt(1-2*x-3*x^2) + 3*(1-3*x))/(2*(2-13*x+21*x^2)), x, n+1), x, n), n = 0..30); # G. C. Greubel, Jan 29 2020
  • Mathematica
    CoefficientList[Series[(Sqrt[-3*x^2-2*x+1]-3*(3*x-1))/(2*(21*x^2-13*x+2)), {x, 0, 30}], x] (* Vaclav Kotesovec, Feb 12 2014 *)
  • PARI
    my(x='x+O('x^30)); Vec( (sqrt(1-2*x-3*x^2) + 3*(1-3*x))/(2*(2-13*x+21*x^2)) ) \\ G. C. Greubel, Jan 29 2020
    
  • Sage
    def A126932_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( (sqrt(1-2*x-3*x^2) + 3*(1-3*x))/(2*(2-13*x+21*x^2)) ).list()
    A126932_list(30) # G. C. Greubel, Jan 29 2020

Formula

a(n+1) = 3*a(n) + A059738(n) with a(0)=1.
G.f: (sqrt(1-2*x-3*x^2) + 3*(1-3*x))/(2*(2-13*x+21*x^2)). - Paul Barry, Nov 06 2008
Conjecture: +2*n*a(n) -11*n*a(n-1) +4*(2*n+3)*a(n-2) +21*(n-2)*a(n-3)=0. - R. J. Mathar, Nov 24 2012
a(n) ~ 3 * 7^n / 2^(n+1). - Vaclav Kotesovec, Feb 12 2014

Extensions

Corrected and extended by Vincenzo Librandi, Feb 13 2014

A036987 Fredholm-Rueppel sequence.

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Keywords

Comments

Binary representation of the Kempner-Mahler number Sum_{k>=0} 1/2^(2^k) = A007404.
a(n) = (product of digits of n; n in binary notation) mod 2. This sequence is a transformation of the Thue-Morse sequence (A010060), since there exists a function f such that f(sum of digits of n) = (product of digits of n). - Ctibor O. Zizka, Feb 12 2008
a(n-1), n >= 1, the characteristic sequence for powers of 2, A000079, is the unique solution of the following formal product and formal power series identity: Product_{j>=1} (1 + a(j-1)*x^j) = 1 + Sum_{k>=1} x^k = 1/(1-x). The product is therefore Product_{l>=1} (1 + x^(2^l)). Proof. Compare coefficients of x^n and use the binary representation of n. Uniqueness follows from the recurrence relation given for the general case under A147542. - Wolfdieter Lang, Mar 05 2009
a(n) is also the number of orbits of length n for the map x -> 1-cx^2 on [-1,1] at the Feigenbaum critical value c=1.401155... . - Thomas Ward, Apr 08 2009
A054525 (Mobius transform) * A001511 = A036987 = A047999^(-1) * A001511 = the inverse of Sierpiński's gasket * the ruler sequence. - Gary W. Adamson, Oct 26 2009 [Of course this is only vaguely correct depending on how the fuzzy indexing in these formulas is made concrete. - R. J. Mathar, Jun 20 2014]
Characteristic function of A000225. - Reinhard Zumkeller, Mar 06 2012
Also parity of the Catalan numbers A000108. - Omar E. Pol, Jan 17 2012
For n >= 2, also the largest exponent k >= 0 such that n^k in binary notation does not contain both 0 and 1. Unlike for the decimal version of this sequence, A062518, where the terms are only conjectural, for this sequence the values of a(n) can be proved to be the characteristic function of A000225, as follows: n^k will contain both 0 and 1 unless n^k = 2^r-1 for some r. But this is a special case of Catalan's equation x^p = y^q-1, which was proved by Preda Mihăilescu to have no nontrivial solution except 2^3 = 3^2 - 1. - Christopher J. Smyth, Aug 22 2014
Image, under the coding a,b -> 1; c -> 0, of the fixed point, starting with a, of the morphism a -> ab, b -> cb, c -> cc. - Jeffrey Shallit, May 14 2016
Number of nonisomorphic Boolean algebras of order n+1. - Jianing Song, Jan 23 2020

Examples

			G.f. = 1 + x + x^3 + x^7 + x^15 + x^31 + x^63 + x^127 + x^255 + x^511 + ...
a(7) = 1 since 7 = 2^3 - 1, while a(10) = 0 since 10 is not of the form 2^k - 1 for any integer k.
		

Crossrefs

The first row of A073346. Occurs for first time in A073202 as row 6 (and again as row 8).
Congruent to any of the sequences A000108, A007460, A007461, A007463, A007464, A061922, A068068 reduced modulo 2. Characteristic function of A000225.
If interpreted with offset=1 instead of 0 (i.e., a(1)=1, a(2)=1, a(3)=0, a(4)=1, ...) then this is the characteristic function of 2^n (A000079) and as such occurs as the first row of A073265. Also, in that case the INVERT transform will produce A023359.
This is Guy Steele's sequence GS(1, 3), also GS(3, 1) (see A135416).
Cf. A054525, A047999. - Gary W. Adamson, Oct 26 2009

Programs

  • Haskell
    a036987 n = ibp (n+1) where
       ibp 1 = 1
       ibp n = if r > 0 then 0 else ibp n' where (n',r) = divMod n 2
    a036987_list = 1 : f [0,1] where f (x:y:xs) = y : f (x:xs ++ [x,x+y])
    -- Same list generator function as for a091090_list, cf. A091090.
    -- Reinhard Zumkeller, May 19 2015, Apr 13 2013, Mar 13 2013
    
  • Maple
    A036987:= n-> `if`(2^ilog2(n+1) = n+1, 1, 0):
    seq(A036987(n), n=0..128);
  • Mathematica
    RealDigits[ N[ Sum[1/10^(2^n), {n, 0, Infinity}], 110]][[1]]
    (* Recurrence: *)
    t[n_, 1] = 1; t[1, k_] = 1;
    t[n_, k_] := t[n, k] =
      If[n < k, If[n > 1 && k > 1, -Sum[t[k - i, n], {i, 1, n - 1}], 0],
       If[n > 1 && k > 1, Sum[t[n - i, k], {i, 1, k - 1}], 0]];
    Table[t[n, k], {k, n, n}, {n, 104}]
    (* Mats Granvik, Jun 03 2011 *)
    mb2d[n_]:=1 - Module[{n2 = IntegerDigits[n, 2]}, Max[n2] - Min[n2]]; Array[mb2d, 120, 0] (* Vincenzo Librandi, Jul 19 2019 *)
    Table[PadRight[{1},2^k,0],{k,0,7}]//Flatten (* Harvey P. Dale, Apr 23 2022 *)
  • PARI
    {a(n) =( n++) == 2^valuation(n, 2)}; /* Michael Somos, Aug 25 2003 */
    
  • PARI
    a(n) = !bitand(n, n+1); \\ Ruud H.G. van Tol, Apr 05 2023
    
  • Python
    from sympy import catalan
    def a(n): return catalan(n)%2 # Indranil Ghosh, May 25 2017
    
  • Python
    def A036987(n): return int(not(n&(n+1))) # Chai Wah Wu, Jul 06 2022

Formula

1 followed by a string of 2^k - 1 0's. Also a(n)=1 iff n = 2^m - 1.
a(n) = a(floor(n/2)) * (n mod 2) for n>0 with a(0)=1. - Reinhard Zumkeller, Aug 02 2002 [Corrected by Mikhail Kurkov, Jul 16 2019]
Sum_{n>=0} 1/10^(2^n) = 0.110100010000000100000000000000010...
1 if n=0, floor(log_2(n+1)) - floor(log_2(n)) otherwise. G.f.: (1/x) * Sum_{k>=0} x^(2^k) = Sum_{k>=0} x^(2^k-1). - Ralf Stephan, Apr 28 2003
a(n) = 1 - A043545(n). - Michael Somos, Aug 25 2003
a(n) = -Sum_{d|n+1} mu(2*d). - Benoit Cloitre, Oct 24 2003
Dirichlet g.f. for right-shifted sequence: 2^(-s)/(1-2^(-s)).
a(n) = A000108(n) mod 2 = A001405(n) mod 2. - Paul Barry, Nov 22 2004
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(n,k)*Sum_{j=0..k} binomial(k, 2^j-1). - Paul Barry, Jun 01 2006
A000523(n+1) = Sum_{k=1..n} a(k). - Mitch Harris, Jul 22 2011
a(n) = A209229(n+1). - Reinhard Zumkeller, Mar 07 2012
a(n) = Sum_{k=1..n} A191898(n,k)*cos(Pi*(n-1)*(k-1))/n; (conjecture). - Mats Granvik, Mar 04 2013
a(n) = A000035(A000108(n)). - Omar E. Pol, Aug 06 2013
a(n) = 1 iff n=2^k-1 for some k, 0 otherwise. - M. F. Hasler, Jun 20 2014
a(n) = ceiling(log_2(n+2)) - ceiling(log_2(n+1)). - Gionata Neri, Sep 06 2015
From John M. Campbell, Jul 21 2016: (Start)
a(n) = (A000168(n-1) mod 2).
a(n) = (A000531(n+1) mod 2).
a(n) = (A000699(n+1) mod 2).
a(n) = (A000891(n) mod 2).
a(n) = (A000913(n-1) mod 2), for n>1.
a(n) = (A000917(n-1) mod 2), for n>0.
a(n) = (A001142(n) mod 2).
a(n) = (A001246(n) mod 2).
a(n) = (A001246(n) mod 4).
a(n) = (A002057(n-2) mod 2), for n>1.
a(n) = (A002430(n+1) mod 2). (End)
a(n) = 2 - A043529(n). - Antti Karttunen, Nov 19 2017
a(n) = floor(1+log(n+1)/log(2)) - floor(log(2n+1)/log(2)). - Adriano Caroli, Sep 22 2019
This is also the decimal expansion of -Sum_{k>=1} mu(2*k)/(10^k - 1), where mu is the Möbius function (A008683). - Amiram Eldar, Jul 12 2020

Extensions

Edited by M. F. Hasler, Jun 20 2014

A061554 Square table read by antidiagonals: a(n,k) = binomial(n+k, floor(k/2)).

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 6, 4, 4, 1, 1, 10, 10, 5, 5, 1, 1, 20, 15, 15, 6, 6, 1, 1, 35, 35, 21, 21, 7, 7, 1, 1, 70, 56, 56, 28, 28, 8, 8, 1, 1, 126, 126, 84, 84, 36, 36, 9, 9, 1, 1, 252, 210, 210, 120, 120, 45, 45, 10, 10, 1, 1, 462, 462, 330, 330, 165, 165, 55, 55, 11, 11, 1, 1
Offset: 0

Views

Author

Henry Bottomley, May 17 2001

Keywords

Comments

Equivalently, a triangle read by rows, where the rows are obtained by sorting the elements of the rows of Pascal's triangle (A007318) into descending order. - Philippe Deléham, May 21 2005
Equivalently, as a triangle read by rows, this is T(n,k)=binomial(n,floor((n-k)/2)); column k then has e.g.f. Bessel_I(k,2x)+Bessel_I(k+1,2x). - Paul Barry, Feb 28 2006
Antidiagonal sums are A037952(n+1) = C(n+1,[n/2]). Matrix inverse is the row reversal of triangle A066170. Eigensequence is A125094(n) = Sum_{k=0..n-1} A125093(n-1,k)*A125094(k). - Paul D. Hanna, Nov 20 2006
Riordan array (1/(1-x-x^2*c(x^2)),x*c(x^2)); where c(x)=g.f.for Catalan numbers A000108. - Philippe Deléham, Mar 17 2007
Triangle T(n,k), 0<=k<=n, read by rows given by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+T(n-1,k+1) for k>=1. - Philippe Deléham, Mar 27 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
T(n,k) is the number of paths from (0,k) to some (n,m) which never dip below y=0, touch y=0 at least once and are made up only of the steps (1,1) and (1,-1). This can be proved using the recurrence supplied by Deléham. - Gerald McGarvey, Oct 15 2008
Triangle read by rows = partial sums of A053121 terms starting from the right. - Gary W. Adamson, Oct 24 2008
As a subset of the "family of triangles" (Deleham comment of Sep 25 2007), beginning with a signed variant of A061554, M = (-1,0) = (1; -1, 1; 2, -1, 1; -3, 3, -1, 1; ...) successive binomial transforms of M yield (0,1) - A089942; (1,2) - A039599; (2,3) - A124733; (3,4) - A124574; (4,5) - A126331; ... such that the binomial transform of the triangle generated from (n,n+1) = the triangle generated from (n+1,n+2). Similarly, another subset beginning with A053121 - (0,0), and taking successive binomial transforms yields (1,1) - A064189; (2,2) - A039598; (3,3) - A091965, ... By rows, the triangle generated from (n,n) can be obtained by taking pairwise sums from the (n-1,n) triangle starting from the right. For example, row 2 of (1,2) - A039599 = (2, 3, 1); and taking pairwise sums from the right we obtain (5, 4, 1) = row 2 of (2,2) - A039598. - Gary W. Adamson, Aug 04 2011
The triangle by rows (n) with alternating signs (+-+...) from the top as a set of simultaneous equations solves for diagonal lengths of odd N (N = 2n+1) regular polygons. The constants in each case are powers of c = 2*cos(2*Pi/N). By way of example, the first 3 rows relate to the heptagon and the simultaneous equations are (1,0,0) = 1; (-1,1,0) = c = 1.24697...; and (2,-1,1) = c^2. The answers are 1, 2.24697..., and 1.801...; the 3 distinct diagonal lengths of the heptagon with edge = 1. - Gary W. Adamson, Sep 07 2011

Examples

			The array starts:
   1, 1,  2,  3,  6, 10,  20,  35,   70,  126, ...
   1, 1,  3,  4, 10, 15,  35,  56,  126,  210, ...
   1, 1,  4,  5, 15, 21,  56,  84,  210,  330, ...
   1, 1,  5,  6, 21, 28,  84, 120,  330,  495, ...
   1, 1,  6,  7, 28, 36, 120, 165,  495,  715, ...
   1, 1,  7,  8, 36, 45, 165, 220,  715, 1001, ...
   1, 1,  8,  9, 45, 55, 220, 286, 1001, 1365, ...
   1, 1,  9, 10, 55, 66, 286, 364, 1365, 1820, ...
   1, 1, 10, 11, 66, 78, 364, 455, 1820, 2380, ...
   1, 1, 11, 12, 78, 91, 455, 560, 2380, 3060, ...
Triangle (antidiagonal) version begins:
    1;
    1,   1;
    2,   1,   1;
    3,   3,   1,   1;
    6,   4,   4,   1,   1;
   10,  10,   5,   5,   1,   1;
   20,  15,  15,   6,   6,   1,  1;
   35,  35,  21,  21,   7,   7,  1,  1;
   70,  56,  56,  28,  28,   8,  8,  1,  1;
  126, 126,  84,  84,  36,  36,  9,  9,  1,  1;
  252, 210, 210, 120, 120,  45, 45, 10, 10,  1, 1;
  462, 462, 330, 330, 165, 165, 55, 55, 11, 11, 1, 1; ...
Matrix inverse begins:
   1;
  -1,  1;
  -1, -1,   1;
   1, -2,  -1,   1;
   1,  2,  -3,  -1,  1;
  -1,  3,   3,  -4, -1,  1;
  -1, -3,   6,   4, -5, -1,  1;
   1, -4,  -6,  10,  5, -6, -1,  1;
   1,  4, -10, -10, 15,  6, -7, -1, 1; ...
From _Paul Barry_, May 21 2009: (Start)
Production matrix is
  1, 1,
  1, 0, 1,
  0, 1, 0, 1,
  0, 0, 1, 0, 1,
  0, 0, 0, 1, 0, 1,
  0, 0, 0, 0, 1, 0, 1,
  0, 0, 0, 0, 0, 1, 0, 1 (End)
		

Crossrefs

Rows are A001405, A037952, A037955, A037951, A037956, A037953, A037957 etc. Columns are truncated pairs of A000012, A000027, A000217, A000292, A000332, A000389, A000579, etc. Main diagonal is alternate values of A051036.

Programs

  • Maple
    T := proc(n, k) option remember;
    if n = k then 1 elif k < 0 or n < 0 or k > n then 0
    elif k = 0 then T(n-1, 0) + T(n-1, 1) else T(n-1, k-1) + T(n-1, k+1) fi end:
    for n from 0 to 9 do seq(T(n, k), k = 0..n) od; # Peter Luschny, May 25 2021
  • Mathematica
    t[n_, k_] = Binomial[n, Floor[(n+1)/2 - (-1)^(n-k)*(k+1)/2]]; Flatten[Table[t[n, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, May 31 2011 *)
  • PARI
    T(n,k)=binomial(n,(n+1)\2-(-1)^(n-k)*((k+1)\2))

Formula

As a triangle: T(n,k) = binomial(n,m) where m = floor((n+1)/2 - (-1)^(n-k)*(k+1)/2).
a(0, k) = binomial(k, floor(k/2)) = A001405(k); for n>0 T(n, k) = T(n+1, k-2) + T(n-1, k).
n-th row = M^n * V, where M = the infinite tridiagonal matrix with all 1's in the super and subdiagonals and (1,0,0,0,...) in the main diagonal. V = the infinite vector [1,0,0,0,...]. Example: (3,3,1,1,0,0,0,...) = M^3 * V. - Gary W. Adamson, Nov 04 2006
Sum_{k=0..n} T(m,k)*T(n,k) = T(m+n,0) = A001405(m+n). - Philippe Deléham, Feb 26 2007
Sum_{k=0..n} T(n,k)=2^n. - Philippe Deléham, Mar 27 2007
Sum_{k=0..n} T(n,k)*x^k = A127361(n), A126869(n), A001405(n), A000079(n), A127358(n), A127359(n), A127360(n) for x = -2, -1, 0, 1, 2, 3, 4 respectively. - Philippe Deléham, Dec 04 2009

Extensions

Entry revised by N. J. A. Sloane, Nov 22 2006

A126075 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) = 2*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + T(n-1,k+1) for k >= 1.

Original entry on oeis.org

1, 2, 1, 5, 2, 1, 12, 6, 2, 1, 30, 14, 7, 2, 1, 74, 37, 16, 8, 2, 1, 185, 90, 45, 18, 9, 2, 1, 460, 230, 108, 54, 20, 10, 2, 1, 1150, 568, 284, 128, 64, 22, 11, 2, 1, 2868, 1434, 696, 348, 150, 75, 24, 12, 2, 1
Offset: 0

Views

Author

Philippe Deléham, Mar 02 2007

Keywords

Comments

Riordan array (c(x^2)/(1-2xc(x^2)),xc(x^2)) where c(x)=g.f. of Catalan numbers A000108. - Philippe Deléham, Mar 18 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

Examples

			Triangle begins:
     1;
     2,    1;
     5,    2,   1;
    12,    6,   2,   1;
    30,   14,   7,   2,   1;
    74,   37,  16,   8,   2,  1;
   185,   90,  45,  18,   9,  2,  1;
   460,  230, 108,  54,  20, 10,  2,  1;
  1150,  568, 284, 128,  64, 22, 11,  2, 1;
  2868, 1434, 696, 348, 150, 75, 24, 12, 2, 1;
		

Crossrefs

Programs

  • Maple
    A126075 := proc (n, k)
    add( 2^(n-k-2*j)*binomial(n, j), j = 0..floor((n-k)/2) ) - add( 2^(n-k-2-2*j)*binomial(n, j), j = 0..floor((n-k-2)/2) )
    end proc:
    # display sequence in triangular form
    for n from 0 to 10 do seq(A126075(n, k), k = 0..n) end do;
    # Peter Bala, Feb 20 2018
  • Mathematica
    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, 2, 0], {n, 0, 49}, {k, 0, n}] // Flatten  (* G. C. Greubel, Apr 21 2017 *)

Formula

Sum_{k=0..n} T(n,k) = A127358(n). T(n,0)=A054341(n).
Sum_{k=0..n} T(n,k)*(-k+1) = 2^n. - Philippe Deléham, Mar 25 2007
From Peter Bala, Feb 20 2018: (Start)
T(n,k) = Sum_{j = 0..floor((n-k)/2)} 2^(n-k-2*j)*binomial(n, j) - Sum_{j = 0..floor((n-k-2)/2)} 2^(n-k-2-2*j)*binomial(n, j), 0 <= k <= n. - Peter Bala, Feb 20 2018
The n-th row polynomial in descending powers of x is the n-th Taylor polynomial of the rational function (1 - x^2)/(1 - 2*x) * (1 + x^2)^n about 0. For example, for n = 4, (1 - x^2)/(1 - 2*x) * (1 + x^2)^4 = (30*x^4 + 14*x*3 + 7*x^2 + 2*x + 1) + O(x^5). (End)

A054341 Row sums of triangle A054336 (central binomial convolutions).

Original entry on oeis.org

1, 2, 5, 12, 30, 74, 185, 460, 1150, 2868, 7170, 17904, 44760, 111834, 279585, 698748, 1746870, 4366460, 10916150, 27287944, 68219860, 170541252, 426353130, 1065853432, 2664633580, 6661479944, 16653699860, 41633878200, 104084695500, 260210401530, 650526003825
Offset: 0

Views

Author

Wolfdieter Lang, Mar 13 2000

Keywords

Comments

a(n) = # Dyck (n+1)-paths all of whose components are symmetric. A strict Dyck path is one with exactly one return to ground level (necessarily at the end). Every nonempty Dyck path is expressible uniquely as a concatenation of one or more strict Dyck paths, called its components. - David Callan, Mar 02 2005
a(n) = # 2-Motzkin paths (i.e., Motzkin paths with blue and red level steps) with no level steps at positive height. Example: a(2)=5 because, denoting U=(1,1), D=(1,-1), B=blue (1,0), R=red (1,0), we have BB, BR, RB, RR, and UD. - Emeric Deutsch, Jun 07 2011
Inverse Chebyshev transform of the second kind applied to 2^n. This maps g(x) -> c(x^2)g(xc(x^2)). - Paul Barry, Sep 14 2005
Hankel transform of this sequence gives A000012 = [1,1,1,1,1,1,1,...]. - Philippe Deléham, Oct 24 2007
Inverse binomial transform of A059738. - Philippe Deléham, Nov 24 2009

Examples

			a(4) = 30, the upper left term of M^4.
		

Crossrefs

Programs

  • Maple
    b:= proc(x, y) option remember; `if`(x=0, 1,
          b(x-1, 0)+`if`(y>0, b(x-1, y-1), 0)+b(x-1, y+1))
        end:
    a:= n-> b(n, 0):
    seq(a(n), n=0..31);  # Alois P. Heinz, Jan 23 2024
  • Mathematica
    CoefficientList[Series[2/(1-4*x+Sqrt[1-4*x^2]), {x, 0, 20}], x] (* Vaclav Kotesovec, Feb 13 2014 *)

Formula

a(n) = Sum_{m=0..n} A054336(n, m).
G.f.: 1/(1-2*x-x^2*c(x^2)), where c(x) = g.f. for Catalan numbers A000108.
From_Paul Barry_, Sep 14 2005: (Start)
G.f.: c(x^2)/(1-2*x*c(x^2));
a(n) = Sum_{k=0..n} binomial(n,(n-k)/2)*(1 + (-1)^(n+k))*2^k*(k+1)/(n+k+2). (End)
G.f.: 2/(1-4*x+sqrt(1-4*x^2)). - Ira M. Gessel, Oct 27 2013
a(n) = A127358(n+1) - 2*A127358(n). - Philippe Deléham, Mar 02 2007
a(n) = A126075(n,0). - Philippe Deléham, Nov 24 2009
a(n) = Sum_{k=0..n} A053121(n,k)*2^k. - Philippe Deléham, Nov 28 2009
From Gary W. Adamson, Sep 07 2011: (Start)
a(n) is the upper left term of M^n, M is an infinite square production matrix as follows:
2, 1, 0, 0, 0, ...
1, 0, 2, 0, 0, ...
0, 1, 0, 1, 0, ...
0, 0, 1, 0, 1, ...
0, 0, 0, 1, 0, ...
... (End)
Conjecture: 2*(n+1)*a(n) +5*(-n-1)*a(n-1) +8*(-n+2)*a(n-2) +20*(n-2)*a(n-3)=0. - R. J. Mathar, Nov 30 2012
a(n) ~ 3 * 5^n / 2^(n+2). - Vaclav Kotesovec, Feb 13 2014

A107430 Triangle read by rows: row n is row n of Pascal's triangle (A007318) sorted into increasing order.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 4, 6, 1, 1, 5, 5, 10, 10, 1, 1, 6, 6, 15, 15, 20, 1, 1, 7, 7, 21, 21, 35, 35, 1, 1, 8, 8, 28, 28, 56, 56, 70, 1, 1, 9, 9, 36, 36, 84, 84, 126, 126, 1, 1, 10, 10, 45, 45, 120, 120, 210, 210, 252, 1, 1, 11, 11, 55, 55, 165, 165, 330, 330, 462, 462, 1
Offset: 0

Views

Author

Philippe Deléham, May 21 2005

Keywords

Comments

By rows, equals partial sums of A053121 reversed rows. Example: Row 4 of A053121 = (2, 0, 3, 0, 1) -> (1, 0, 3, 0, 2) -> (1, 1, 4, 4, 6). - Gary W. Adamson, Dec 28 2008, edited by Michel Marcus, Sep 22 2015

Examples

			Triangle begins:
1;
1,1;
1,1,2;
1,1,3,3;
1,1,4,4,6;
		

Crossrefs

A061554 is similar but with rows sorted into decreasing order.
Cf. A034868.
Cf. A053121. - Gary W. Adamson, Dec 28 2008
Cf. A103284.

Programs

  • Haskell
    import Data.List (sort)
    a107430 n k = a107430_tabl !! n !! k
    a107430_row n = a107430_tabl !! n
    a107430_tabl = map sort a007318_tabl
    -- Reinhard Zumkeller, May 26 2013
    
  • Magma
    /* As triangle */ [[Binomial(n,Floor(k/2)) : k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 22 2015
    
  • Maple
    for n from 0 to 10 do sort([seq(binomial(n,k),k=0..n)]) od; # yields sequence in triangular form. - Emeric Deutsch, May 28 2005
  • Mathematica
    Flatten[ Table[ Sort[ Table[ Binomial[n, k], {k, 0, n}]], {n, 0, 12}]] (* Robert G. Wilson v, May 28 2005 *)
  • PARI
    for(n=0,20, for(k=0,n, print1(binomial(n,floor(k/2)), ", "))) \\ G. C. Greubel, May 22 2017

Formula

T(n,k) = C(n,floor(k/2)). - Paul Barry, Dec 15 2006; corrected by Philippe Deléham, Mar 15 2007
Sum_{k=0..n} T(n,k)*x^(n-k) = A127363(n), A127362(n), A127361(n), A126869(n), A001405(n), A000079(n), A127358(n), A127359(n), A127360(n) for x=-4,-3,-2,-1,0,1,2,3,4 respectively. - Philippe Deléham, Mar 29 2007

Extensions

More terms from Emeric Deutsch and Robert G. Wilson v, May 28 2005

A126125 Triangle equal to the matrix square of the triangle binomial(n,floor((n+1-(-1)^(n+k)*(k+1))/2)).

Original entry on oeis.org

1, 2, 1, 5, 2, 1, 11, 7, 2, 1, 27, 15, 9, 2, 1, 61, 44, 19, 11, 2, 1, 149, 97, 65, 23, 13, 2, 1, 342, 267, 141, 90, 27, 15, 2, 1, 835, 599, 433, 193, 119, 31, 17, 2, 1, 1939, 1598, 956, 655, 253, 152, 35, 19, 2, 1, 4740, 3631
Offset: 0

Views

Author

Gary W. Adamson, Dec 17 2006

Keywords

Comments

The triangular view of A061554 is:
1;
1, 1;
2, 1, 1;
3, 3, 1, 1;
6, 4, 4, 1, 1;
The matrix square of this (infinite) triangle is this triangle here.

Examples

			First few rows of the triangle are:
1;
2, 1;
5, 2, 1;
11, 7, 2, 1;
27, 15, 9, 2, 1;
61, 44, 19, 11, 2, 1;
..,
		

Crossrefs

Cf. A061554, A127358 (row sums).

Programs

  • Maple
    A061554 := proc(n,k) local m ; m := floor((n+1)/2 - (-1)^(n-k)*(k+1)/2) ; binomial(n,m) ; end proc:
    A126125 := proc(n,k) add(A061554(n,j)*A061554(j,k),j=k..n) ; end proc: # R. J. Mathar, Sep 17 2011
Showing 1-7 of 7 results.