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

A181474 Sequence related to the Hankel transform of A105523(n+5).

Original entry on oeis.org

1, 2, 10, 15, 42, 56, 120, 150, 275, 330, 546, 637, 980, 1120, 1632, 1836, 2565, 2850, 3850, 4235, 5566, 6072, 7800, 8450, 10647, 11466, 14210, 15225, 18600, 19840, 23936, 25432, 30345, 32130, 37962, 40071, 46930
Offset: 0

Views

Author

Paul Barry, Oct 22 2010

Keywords

Comments

The Hankel transform of A105523(n+5) is (-1)^C(n+2,2)a(n+1).

Programs

  • Mathematica
    CoefficientList[Series[(1 + x + 4 x^2 + x^3 + x^4)/((1 - x)^5 (1 + x)^4), {x, 0, 36}], x] (* Michael De Vlieger, Jul 25 2023 *)

Formula

G.f.: (1+x+4x^2+x^3+x^4)/((1-x)^5(1+x)^4);
a(n) = 2^n*(n+2)*(n+3)*Gamma(floor(n/2)+3)*Gamma(floor((n+1)/2)+1/2)/(12n!*sqrt(Pi)) (suggested by WolframAlpha).
a(n) = +a(n-1) +4*a(n-2) -4*a(n-3) -6*a(n-4) +6*a(n-5) +4*a(n-6) -4*a(n-7) -a(n-8) +a(n-9). a(n) = (n+3)*(n+2)*(2*n^2+2*(-1)^n*n+10*n+5*(-1)^n+11)/96. [R. J. Mathar, Oct 23 2010]

A011782 Coefficients of expansion of (1-x)/(1-2*x) in powers of x.

Original entry on oeis.org

1, 1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, 16384, 32768, 65536, 131072, 262144, 524288, 1048576, 2097152, 4194304, 8388608, 16777216, 33554432, 67108864, 134217728, 268435456, 536870912, 1073741824, 2147483648, 4294967296, 8589934592
Offset: 0

Views

Author

Lee D. Killough (killough(AT)wagner.convex.com)

Keywords

Comments

Apart from initial term, same as A000079 (powers of 2).
Number of compositions (ordered partitions) of n. - Toby Bartels, Aug 27 2003
Number of ways of putting n unlabeled items into (any number of) labeled boxes where every box contains at least one item. Also "unimodal permutations of n items", i.e., those which rise then fall. (E.g., for three items: ABC, ACB, BCA and CBA are unimodal.) - Henry Bottomley, Jan 17 2001
Number of permutations in S_n avoiding the patterns 213 and 312. - Tuwani Albert Tshifhumulo, Apr 20 2001. More generally (see Simion and Schmidt), the number of permutations in S_n avoiding (i) the 123 and 132 patterns; (ii) the 123 and 213 patterns; (iii) the 132 and 213 patterns; (iv) the 132 and 231 patterns; (v) the 132 and 312 patterns; (vi) the 213 and 231 patterns; (vii) the 213 and 312 patterns; (viii) the 231 and 312 patterns; (ix) the 231 and 321 patterns; (x) the 312 and 321 patterns.
a(n+2) is the number of distinct Boolean functions of n variables under action of symmetric group.
Number of unlabeled (1+2)-free posets. - Detlef Pauly, May 25 2003
Image of the central binomial coefficients A000984 under the Riordan array ((1-x), x*(1-x)). - Paul Barry, Mar 18 2005
Binomial transform of (1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...); inverse binomial transform of A007051. - Philippe Deléham, Jul 04 2005
Also, number of rationals in [0, 1) whose binary expansions terminate after n bits. - Brad Chalfan, May 29 2006
Equals row sums of triangle A144157. - Gary W. Adamson, Sep 12 2008
Prepend A089067 with a 1, getting (1, 1, 3, 5, 13, 23, 51, ...) as polcoeff A(x); then (1, 1, 2, 4, 8, 16, ...) = A(x)/A(x^2). - Gary W. Adamson, Feb 18 2010
An elephant sequence, see A175655. For the central square four A[5] vectors, with decimal values 2, 8, 32 and 128, lead to this sequence. For the corner squares these vectors lead to the companion sequence A094373. - Johannes W. Meijer, Aug 15 2010
From Paul Curtz, Jul 20 2011: (Start)
Array T(m,n) = 2*T(m,n-1) + T(m-1,n):
1, 1, 2, 4, 8, 16, ... = a(n)
1, 3, 8, 20, 48, 112, ... = A001792,
1, 5, 18, 56, 160, 432, ... = A001793,
1, 7, 32, 120, 400, 1232, ... = A001794,
1, 9, 50, 220, 840, 2912, ... = A006974, followed with A006975, A006976, gives nonzero coefficients of Chebyshev polynomials of first kind A039991 =
1,
1, 0,
2, 0, -1,
4, 0, -3, 0,
8, 0, -8, 0, 1.
T(m,n) third vertical: 2*n^2, n positive (A001105).
Fourth vertical appears in Janet table even rows, last vertical (A168342 array, A138509, rank 3, 13, = A166911)). (End)
A131577(n) and differences are:
0, 1, 2, 4, 8, 16,
1, 1, 2, 4, 8, 16, = a(n),
0, 1, 2, 4, 8, 16,
1, 1, 2, 4, 8, 16.
Number of 2-color necklaces of length 2n equal to their complemented reversal. For length 2n+1, the number is 0. - David W. Wilson, Jan 01 2012
Edges and also central terms of triangle A198069: a(0) = A198069(0,0) and for n > 0: a(n) = A198069(n,0) = A198069(n,2^n) = A198069(n,2^(n-1)). - Reinhard Zumkeller, May 26 2013
These could be called the composition numbers (see the second comment) since the equivalent sequence for partitions is A000041, the partition numbers. - Omar E. Pol, Aug 28 2013
Number of self conjugate integer partitions with exactly n parts for n>=1. - David Christopher, Aug 18 2014
The sequence is the INVERT transform of (1, 0, 1, 1, 2, 3, 5, 8, 13, 21, 34, ...). - Gary W. Adamson, Jul 16 2015
Number of threshold graphs on n nodes [Hougardy]. - Falk Hüffner, Dec 03 2015
Number of ternary words of length n in which binary subwords appear in the form 10...0. - Milan Janjic, Jan 25 2017
a(n) is the number of words of length n over an alphabet of two letters, of which one letter appears an even number of times (the empty word of length 0 is included). See the analogous odd number case in A131577, and the Balakrishnan reference in A006516 (the 4-letter odd case), pp. 68-69, problems 2.66, 2.67 and 2.68. - Wolfdieter Lang, Jul 17 2017
Number of D-equivalence classes of Łukasiewicz paths. Łukasiewicz paths are D-equivalent iff the positions of pattern D are identical in these paths. - Sergey Kirgizov, Apr 08 2018
Number of color patterns (set partitions) for an oriented row of length n using two or fewer colors (subsets). Two color patterns are equivalent if we permute the colors. For a(4)=8, the 4 achiral patterns are AAAA, AABB, ABAB, and ABBA; the 4 chiral patterns are the 2 pairs AAAB-ABBB and AABA-ABAA. - Robert A. Russell, Oct 30 2018
The determinant of the symmetric n X n matrix M defined by M(i,j) = (-1)^max(i,j) for 1 <= i,j <= n is equal to a(n) * (-1)^(n*(n+1)/2). - Bernard Schott, Dec 29 2018
For n>=1, a(n) is the number of permutations of length n whose cyclic representations can be written in such a way that when the cycle parentheses are removed what remains is 1 through n in natural order. For example, a(4)=8 since there are exactly 8 permutations of this form, namely, (1 2 3 4), (1)(2 3 4), (1 2)(3 4), (1 2 3)(4), (1)(2)(3 4), (1)(2 3)(4), (1 2)(3)(4), and (1)(2)(3)(4). Our result follows readily by conditioning on k, the number of parentheses pairs of the form ")(" in the cyclic representation. Since there are C(n-1,k) ways to insert these in the cyclic representation and since k runs from 0 to n-1, we obtain a(n) = Sum_{k=0..n-1} C(n-1,k) = 2^(n-1). - Dennis P. Walsh, May 23 2020
Maximum number of preimages that a permutation of length n + 1 can have under the consecutive-231-avoiding stack-sorting map. - Colin Defant, Aug 28 2020
a(n) is the number of occurrences of the empty set {} in the von Neumann ordinals from 0 up to n. Each ordinal k is defined as the set of all smaller ordinals: 0 = {}, 1 = {0}, 2 = {0,1}, etc. Since {} is the foundational element of all ordinals, the total number of times it appears grows as powers of 2. - Kyle Wyonch, Mar 30 2025

Examples

			G.f. = 1 + x + 2*x^2 + 4*x^3 + 8*x^4 + 16*x^5 + 32*x^6 + 64*x^7 + 128*x^8 + ...
    ( -1   1  -1)
det (  1   1   1)  = 4
    ( -1  -1  -1)
		

References

  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem, Mathematics and Computer Education Journal, Vol. 31, No. 1, pp. 24-28, Winter 1997.
  • S. Kitaev, Patterns in Permutations and Words, Springer-Verlag, 2011. see p. 399 Table A.7
  • Xavier Merlin, Methodix Algèbre, Ellipses, 1995, p. 153.

Crossrefs

Sequences with g.f.'s of the form ((1-x)/(1-2*x))^k: this sequence (k=1), A045623 (k=2), A058396 (k=3), A062109 (k=4), A169792 (k=5), A169793 (k=6), A169794 (k=7), A169795 (k=8), A169796 (k=9), A169797 (k=10).
Cf. A005418 (unoriented), A122746(n-3) (chiral), A016116 (achiral).
Row sums of triangle A100257.
A row of A160232.
Row 2 of A278984.

Programs

  • Haskell
    a011782 n = a011782_list !! n
    a011782_list = 1 : scanl1 (+) a011782_list
    -- Reinhard Zumkeller, Jul 21 2013
    
  • Magma
    [Floor((1+2^n)/2): n in [0..35]]; // Vincenzo Librandi, Aug 21 2011
    
  • Maple
    A011782:= n-> ceil(2^(n-1)): seq(A011782(n), n=0..50); # Wesley Ivan Hurt, Feb 21 2015
    with(PolynomialTools):  A011782:=seq(coeftayl((1-x)/(1-2*x), x = 0, k),k=0..10^2); # Muniru A Asiru, Sep 26 2017
  • Mathematica
    f[s_] := Append[s, Ceiling[Plus @@ s]]; Nest[f, {1}, 32] (* Robert G. Wilson v, Jul 07 2006 *)
    CoefficientList[ Series[(1-x)/(1-2x), {x, 0, 32}], x] (* Robert G. Wilson v, Jul 07 2006 *)
    Table[Sum[StirlingS2[n, k], {k,0,2}], {n, 0, 30}] (* Robert A. Russell, Apr 25 2018 *)
    Join[{1},NestList[2#&,1,40]] (* Harvey P. Dale, Dec 06 2018 *)
  • PARI
    {a(n) = if( n<1, n==0, 2^(n-1))};
    
  • PARI
    Vec((1-x)/(1-2*x) + O(x^30)) \\ Altug Alkan, Oct 31 2015
    
  • Python
    def A011782(n): return 1 if n == 0 else 2**(n-1) # Chai Wah Wu, May 11 2022
  • Sage
    [sum(stirling_number2(n,j) for j in (0..2)) for n in (0..35)] # G. C. Greubel, Jun 02 2020
    

Formula

a(0) = 1, a(n) = 2^(n-1).
G.f.: (1 - x) / (1 - 2*x) = 1 / (1 - x / (1 - x)). - Michael Somos, Apr 18 2012
E.g.f.: cosh(z)*exp(z) = (exp(2*z) + 1)/2.
a(0) = 1 and for n>0, a(n) = sum of all previous terms.
a(n) = Sum_{k=0..n} binomial(n, 2*k). - Paul Barry, Feb 25 2003
a(n) = Sum_{k=0..n} binomial(n,k)*(1+(-1)^k)/2. - Paul Barry, May 27 2003
a(n) = floor((1+2^n)/2). - Toby Bartels (toby+sloane(AT)math.ucr.edu), Aug 27 2003
G.f.: Sum_{i>=0} x^i/(1-x)^i. - Jon Perry, Jul 10 2004
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(k+1, n-k)*binomial(2*k, k). - Paul Barry, Mar 18 2005
a(n) = Sum_{k=0..floor(n/2)} A055830(n-k, k). - Philippe Deléham, Oct 22 2006
a(n) = Sum_{k=0..n} A098158(n,k). - Philippe Deléham, Dec 04 2006
G.f.: 1/(1 - (x + x^2 + x^3 + ...)). - Geoffrey Critzer, Aug 30 2008
a(n) = A000079(n) - A131577(n).
a(n) = A173921(A000079(n)). - Reinhard Zumkeller, Mar 04 2010
a(n) = Sum_{k=2^n..2^(n+1)-1} A093873(k)/A093875(k), sums of rows of the full tree of Kepler's harmonic fractions. - Reinhard Zumkeller, Oct 17 2010
E.g.f.: (exp(2*x)+1)/2 = (G(0) + 1)/2; G(k) = 1 + 2*x/(2*k+1 - x*(2*k+1)/(x + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Dec 03 2011
A051049(n) = p(n+1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, Apr 18 2012
A008619(n) = p(-1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, Apr 18 2012
INVERT transform is A122367. MOBIUS transform is A123707. EULER transform of A059966. PSUM transform is A000079. PSUMSIGN transform is A078008. BINOMIAL transform is A007051. REVERT transform is A105523. A002866(n) = a(n)*n!. - Michael Somos, Apr 18 2012
G.f.: U(0), where U(k) = 1 + x*(k+3) - x*(k+2)/U(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Oct 10 2012
a(n) = A000041(n) + A056823(n). - Omar E. Pol, Aug 31 2013
E.g.f.: E(0), where E(k) = 1 + x/( 2*k+1 - x/E(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Dec 25 2013
G.f.: 1 + x/(1 + x)*( 1 + 3*x/(1 + 3*x)*( 1 + 5*x/(1 + 5*x)*( 1 + 7*x/(1 + 7*x)*( 1 + ... )))). - Peter Bala, May 27 2017
a(n) = Sum_{k=0..2} stirling2(n, k).
G.f.: Sum_{j=0..k} A248925(k,j)*x^j / Product_{j=1..k} 1-j*x with k=2. - Robert A. Russell, Apr 25 2018
a(n) = A053120(n, n), n >= 0, (main diagonal of triangle of Chebyshev's T polynomials). - Wolfdieter Lang, Nov 26 2019

Extensions

Additional comments from Emeric Deutsch, May 14 2001
Typo corrected by Philippe Deléham, Oct 25 2008

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

A090181 Triangle of Narayana (A001263) with 0 <= k <= n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 6, 6, 1, 0, 1, 10, 20, 10, 1, 0, 1, 15, 50, 50, 15, 1, 0, 1, 21, 105, 175, 105, 21, 1, 0, 1, 28, 196, 490, 490, 196, 28, 1, 0, 1, 36, 336, 1176, 1764, 1176, 336, 36, 1, 0, 1, 45, 540, 2520, 5292, 5292, 2520, 540, 45, 1, 0, 1, 55, 825, 4950, 13860
Offset: 0

Views

Author

Philippe Deléham, Jan 19 2004

Keywords

Comments

Number of Dyck n-paths with exactly k peaks. - Peter Luschny, May 10 2014

Examples

			Triangle starts:
[0] 1;
[1] 0, 1;
[2] 0, 1,  1;
[3] 0, 1,  3,   1;
[4] 0, 1,  6,   6,    1;
[5] 0, 1, 10,  20,   10,    1;
[6] 0, 1, 15,  50,   50,   15,    1;
[7] 0, 1, 21, 105,  175,  105,   21,   1;
[8] 0, 1, 28, 196,  490,  490,  196,  28,  1;
[9] 0, 1, 36, 336, 1176, 1764, 1176, 336, 36, 1;
		

Crossrefs

Mirror image of triangle A131198. A000108 (row sums, Catalan).
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000108(n), A006318(n), A047891(n+1), A082298(n), A082301(n), A082302(n), A082305(n), A082366(n), A082367(n) for x=0,1,2,3,4,5,6,7,8,9. - Philippe Deléham, Aug 10 2006
Sum_{k=0..n} x^(n-k)*T(n,k) = A090192(n+1), A000012(n), A000108(n), A001003(n), A007564(n), A059231(n), A078009(n), A078018(n), A081178(n), A082147(n), A082181(n), A082148(n), A082173(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11. - Philippe Deléham, Oct 21 2006
Sum_{k=0..n} T(n,k)*x^k*(x-1)^(n-k) = A000012(n), A006318(n), A103210(n), A103211(n), A133305(n), A133306(n), A133307(n), A133308(n), A133309(n) for x = 1, 2, 3, 4, 5, 6, 7, 8, 9, respectively. - Philippe Deléham, Oct 20 2007

Programs

  • Magma
    [[(&+[(-1)^(j-k)*Binomial(2*n-j,j)*Binomial(j,k)*Binomial(2*n-2*j,n-j)/(n-j+1): j in [0..n]]): k in [0..n]]: n in [0..10]];
  • Maple
    A090181 := (n,k) -> binomial(n,n-k)*binomial(n-1,n-k)/(n-k+1):
    seq(print( seq(A090181(n,k),k=0..n)),n=0..5); # Peter Luschny, May 10 2014
    egf := 1+int((sqrt(t)*exp((1+t)*x)*BesselI(1,2*sqrt(t)*x))/x,x);
    s := n -> n!*coeff(series(egf,x,n+2),x,n);
    seq(print(seq(coeff(s(n),t,j),j=0..n)),n=0..9); # Peter Luschny, Oct 30 2014
    T := proc(n, k) option remember; if k = n or k = 1 then 1 elif k < 1 then 0 else (2*n/k - 1) * T(n-1, k-1) + T(n-1, k) fi end:
    for n from 0 to 8 do seq(T(n, k), k = 0..n) od;  # Peter Luschny, Dec 31 2024
  • Mathematica
    Flatten[Table[Sum[(-1)^(j-k) * Binomial[2n-j,j] * Binomial[j,k] * CatalanNumber[n-j], {j, 0, n}], {n,0,11},{k,0,n}]] (* Indranil Ghosh, Mar 05 2017 *)
    p[0, ] := 1; p[1, x] := x; p[n_, x_] := ((2 n - 1) (1 + x) p[n - 1, x] - (n - 2) (x - 1)^2 p[n - 2, x]) / (n + 1);
    Table[CoefficientList[p[n, x], x], {n, 0, 9}] // TableForm (* Peter Luschny, Apr 26 2022 *)
  • PARI
    c(n) = binomial(2*n,n)/ (n+1);
    tabl(nn) = {for(n=0, nn, for(k=0, n, print1(sum(j=0, n, (-1)^(j-k) * binomial(2*n-j,j) * binomial(j,k) * c(n-j)),", ");); print(););};
    tabl(11); \\ Indranil Ghosh, Mar 05 2017
    
  • Python
    from functools import cache
    @cache
    def Trow(n):
        if n == 0: return [1]
        if n == 1: return [0, 1]
        if n == 2: return [0, 1, 1]
        A = Trow(n - 2) + [0, 0]
        B = Trow(n - 1) + [1]
        for k in range(n - 1, 1, -1):
            B[k] = (((B[k] + B[k - 1]) * (2 * n - 1)
                   - (A[k] - 2 * A[k - 1] + A[k - 2]) * (n - 2)) // (n + 1))
        return B
    for n in range(10): print(Trow(n)) # Peter Luschny, May 02 2022
    
  • Sage
    def A090181_row(n):
        U = [0]*(n+1)
        for d in DyckWords(n):
            U[d.number_of_peaks()] +=1
        return U
    for n in range(8): A090181_row(n) # Peter Luschny, May 10 2014
    

Formula

Triangle T(n, k), read by rows, given by [0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] where DELTA is the operator defined in A084938. T(0, 0) = 1, T(n, 0) = 0 for n>0, T(n, k) = C(n-1, k-1)*C(n, k-1)/k for k>0.
Sum_{j>=0} T(n,j)*binomial(j,k) = A060693(n,k). - Philippe Deléham, May 04 2007
Sum_{k=0..n} T(n,k)*10^k = A143749(n+1). - Philippe Deléham, Oct 14 2008
From Paul Barry, Nov 10 2008: (Start)
Coefficient array of the polynomials P(n,x) = x^n*2F1(-n,-n+1;2;1/x).
T(n,k) = Sum_{j=0..n} (-1)^(j-k)*C(2n-j,j)*C(j,k)*A000108(n-j). (End)
Sum_{k=0..n} T(n,k)*5^k*3^(n-k) = A152601(n). - Philippe Deléham, Dec 10 2008
Sum_{k=0..n} T(n,k)*(-2)^k = A152681(n); Sum_{k=0..n} T(n,k)*(-1)^k = A105523(n). - Philippe Deléham, Feb 03 2009
Sum_{k=0..n} T(n,k)*2^(n+k) = A156017(n). - Philippe Deléham, Nov 27 2011
T(n, k) = C(n,n-k)*C(n-1,n-k)/(n-k+1). - Peter Luschny, May 10 2014
E.g.f.: 1+Integral((sqrt(t)*exp((1+t)*x)*BesselI(1,2*sqrt(t)*x))/x dx). - Peter Luschny, Oct 30 2014
G.f.: (1+x-x*y-sqrt((1-x*(1+y))^2-4*y*x^2))/(2*x). - Alois P. Heinz, Nov 28 2021, edited by Ron L.J. van den Burg, Dec 19 2021
T(n, k) = [x^k] (((2*n - 1)*(1 + x)*p(n-1, x) - (n - 2)*(x - 1)^2*p(n-2, x))/(n + 1)) with p(0, x) = 1 and p(1, x) = x. - Peter Luschny, Apr 26 2022
Recursion based on rows (see the Python program):
T(n, k) = (((B(k) + B(k-1))*(2*n - 1) - (A(k) - 2*A(k-1) + A(k-2))*(n-2))/(n+1)), where A(k) = T(n-2, k) and B(k) = T(n-1, k), for n >= 3. # Peter Luschny, May 02 2022

A120730 Another version of Catalan triangle A009766.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 2, 1, 0, 0, 2, 3, 1, 0, 0, 0, 5, 4, 1, 0, 0, 0, 5, 9, 5, 1, 0, 0, 0, 0, 14, 14, 6, 1, 0, 0, 0, 0, 14, 28, 20, 7, 1, 0, 0, 0, 0, 0, 42, 48, 27, 8, 1, 0, 0, 0, 0, 0, 42, 90, 75, 35, 9, 1, 0, 0, 0, 0, 0, 0, 132, 165, 110, 44, 10, 1
Offset: 0

Views

Author

Philippe Deléham, Aug 17 2006, corrected Sep 15 2006

Keywords

Comments

Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 1, -1, 0, 0, 1, -1, 0, 0, 1, -1, 0, 0, ...] DELTA [1, 0, 0, -1, 1, 0, 0, -1, 1, 0, 0, -1, 1, ...] where DELTA is the operator defined in A084938.
Aerated version gives A165408. - Philippe Deléham, Sep 22 2009
T(n,k) is the number of length n left factors of Dyck paths having k up steps. Example: T(5,4)=4 because we have UDUUU, UUDUU, UUUDU, and UUUUD, where U=(1,1) and D=(1,-1). - Emeric Deutsch, Jun 19 2011
With zeros omitted: 1,1,1,1,2,1,2,3,1,5,4,1,... = A008313. - Philippe Deléham, Nov 02 2011

Examples

			As a triangle, this begins:
  1;
  0,  1;
  0,  1,  1;
  0,  0,  2,  1;
  0,  0,  2,  3,  1;
  0,  0,  0,  5,  4,  1;
  0,  0,  0,  5,  9,  5,  1;
  0,  0,  0,  0, 14, 14,  6,  1;
  ...
		

Crossrefs

Programs

  • Magma
    A120730:= func< n,k | n gt 2*k select 0 else Binomial(n, k)*(2*k-n+1)/(k+1) >;
    [A120730(n,k): k in [0..n], n in [0..13]]; // G. C. Greubel, Nov 07 2022
    
  • Maple
    G := 4*z/((2*z-1+sqrt(1-4*z^2*t))*(1+sqrt(1-4*z^2*t))): Gser := simplify(series(G, z = 0, 13)): for n from 0 to 12 do P[n] := sort(coeff(Gser, z, n)) end do: for n from 0 to 12 do seq(coeff(P[n], t, k), k = 0 .. n) end do; # yields sequence in triangular form  # Emeric Deutsch, Jun 19 2011
    # second Maple program:
    b:= proc(x, y) option remember; `if`(y<0 or y>x, 0,
         `if`(x=0, 1, add(b(x-1, y+j), j=[-1, 1])))
        end:
    T:= (n, k)-> b(n, 2*k-n):
    seq(seq(T(n, k), k=0..n), n=0..14);  # Alois P. Heinz, Oct 13 2022
  • Mathematica
    b[x_, y_]:= b[x, y]= If[y<0 || y>x, 0, If[x==0, 1, Sum[b[x-1, y+j], {j, {-1, 1}}] ]];
    T[n_, k_] := b[n, 2 k - n];
    Table[Table[T[n, k], {k, 0, n}], {n, 0, 14}] // Flatten (* Jean-François Alcover, Oct 21 2022, after Alois P. Heinz *)
    T[n_, k_]:= If[n>2*k, 0, Binomial[n, k]*(2*k-n+1)/(k+1)];
    Table[T[n, k], {n,0,13}, {k,0,n}]//Flatten (* G. C. Greubel, Nov 07 2022 *)
  • SageMath
    def A120730(n,k): return 0 if (n>2*k) else binomial(n, k)*(2*k-n+1)/(k+1)
    flatten([[A120730(n,k) for k in range(n+1)] for n in range(14)]) # G. C. Greubel, Nov 07 2022

Formula

G.f.: G(t,z) = 4*z/((2*z-1+sqrt(1-4*t*z^2))*(1+sqrt(1-4*t*z^2))). - Emeric Deutsch, Jun 19 2011
Sum_{k=0..n} x^k*T(n,n-k) = A001405(n), A126087(n), A128386(n), A121724(n), A128387(n), A132373(n), A132374(n), A132375(n), A121725(n) for x=1,2,3,4,5,6,7,8,9 respectively. [corrected by Philippe Deléham, Oct 16 2008]
T(2*n,n) = A000108(n); A000108: Catalan numbers.
From Philippe Deléham, Oct 18 2008: (Start)
Sum_{k=0..n} T(n,k)^2 = A000108(n) and Sum_{n>=k} T(n,k) = A000108(k+1).
Sum_{k=0..n} T(n,k)^3 = A003161(n).
Sum_{k=0..n} T(n,k)^4 = A129123(n). (End)
Sum_{k=0..n}, T(n,k)*x^k = A000007(n), A001405(n), A151281(n), A151162(n), A151254(n), A156195(n), A156361(n), A156362(n), A156566(n), A156577(n) for x=0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Feb 10 2009
From G. C. Greubel, Nov 07 2022: (Start)
T(n, k) = 0 if n > 2*k, otherwise binomial(n, k)*(2*k-n+1)/(k+1).
Sum_{k=0..n} (-1)^k*T(n,k) = A105523(n).
Sum_{k=0..n} (-1)^k*T(n,k)^2 = -A132889(n), n >= 1.
Sum_{k=0..floor(n/2)} T(n-k, k) = A357654(n).
T(n, n-1) = A001477(n).
T(n, n-2) = [n=2] + A000096(n-3), n >= 2.
T(n, n-3) = 2*[n<5] + A005586(n-5), n >= 3.
T(n, n-4) = 5*[n<7] - 2*[n=4] + A005587(n-7), n >= 4.
T(2*n+1, n+1) = A000108(n+1), n >= 0.
T(2*n-1, n+1) = A099376(n-1), n >= 1. (End)

A097331 Expansion of 1 + 2x/(1 + sqrt(1 - 4x^2)).

Original entry on oeis.org

1, 1, 0, 1, 0, 2, 0, 5, 0, 14, 0, 42, 0, 132, 0, 429, 0, 1430, 0, 4862, 0, 16796, 0, 58786, 0, 208012, 0, 742900, 0, 2674440, 0, 9694845, 0, 35357670, 0, 129644790, 0, 477638700, 0, 1767263190, 0, 6564120420, 0, 24466267020, 0, 91482563640, 0, 343059613650, 0
Offset: 0

Views

Author

Paul Barry, Aug 05 2004

Keywords

Comments

Binomial transform is A097332. Second binomial transform is A014318.
Essentially the same as A126120. - R. J. Mathar, Jun 15 2008
Hankel transform is A087960(n) = (-1)^binomial(n+1,2). - Paul Barry, Aug 10 2009

Programs

  • Maple
    A097331_list := proc(n) local j, a, w; a := array(0..n); a[0] := 1;
    for w from 1 to n do a[w]:=a[w-1]-(-1)^w*add(a[j]*a[w-j-1],j=1..w-1) od; convert(a,list)end: A097331_list(48); # Peter Luschny, May 19 2011
  • Mathematica
    a[0] = 1; a[n_?OddQ] := CatalanNumber[(n-1)/2]; a[] = 0; Table[a[n], {n, 0, 48}] (* _Jean-François Alcover, Jul 24 2013 *)
  • Sage
    def A097331_list(n) :
        D = [0]*(n+2); D[1] = 1
        b = True; h = 1; R = []
        for i in range(2*n-1) :
            if b :
                for k in range(h,0,-1) : D[k] -= D[k-1]
                h += 1; R.append(abs(D[1]))
            else :
                for k in range(1,h, 1) : D[k] += D[k+1]
            b = not b
        return R
    A097331_list(49) # Peter Luschny, Jun 03 2012

Formula

a(n) = 0^n + Catalan((n-1)/2)(1-(-1)^n)/2.
Unsigned version of A090192, A105523. - Philippe Deléham, Sep 29 2006
From Paul Barry, Aug 10 2009: (Start)
G.f.: 1+xc(x^2), c(x) the g.f. of A000108;
G.f.: 1/(1-x/(1+x/(1+x/(1-x/(1-x/(1+x/(1+x/(1-x/(1-x/(1+... (continued fraction);
G.f.: 1+x/(1-x^2/(1-x^2/(1-x^2/(1-x^2/(1-... (continued fraction). (End)
G.f.: 1/(1-z/(1-z/(1-z/(...)))) where z=x/(1+2*x) (continued fraction); more generally g.f. C(x/(1+2*x)) where C(x) is the g.f. for the Catalan numbers (A000108). - Joerg Arndt, Mar 18 2011
Conjecture: (n+1)*a(n) + n*a(n-1) + 4*(-n+2)*a(n-2) + 4*(-n+3)*a(n-3)=0. - R. J. Mathar, Dec 02 2012
Recurrence: (n+3)*a(n+2) = 4*n*a(n), a(0)=a(1)=1. For nonzero terms, a(n) ~ 2^(n+1)/((n+1)^(3/2)*sqrt(2*Pi)). - Fung Lam, Mar 17 2014

A106195 Riordan array (1/(1-2*x), x*(1-x)/(1-2*x)).

Original entry on oeis.org

1, 2, 1, 4, 3, 1, 8, 8, 4, 1, 16, 20, 13, 5, 1, 32, 48, 38, 19, 6, 1, 64, 112, 104, 63, 26, 7, 1, 128, 256, 272, 192, 96, 34, 8, 1, 256, 576, 688, 552, 321, 138, 43, 9, 1, 512, 1280, 1696, 1520, 1002, 501, 190, 53, 10, 1, 1024, 2816, 4096, 4048, 2972, 1683, 743, 253, 64, 11
Offset: 0

Views

Author

Gary W. Adamson, Apr 24 2005; Paul Barry, May 21 2006

Keywords

Comments

Extract antidiagonals from the product P * A, where P = the infinite lower triangular Pascal's triangle matrix; and A = the Pascal's triangle array:
1, 1, 1, 1, ...
1, 2, 3, 4, ...
1, 3, 6, 10, ...
1, 4, 10, 20, ...
...
Row sums are Fibonacci(2n+2). Diagonal sums are A006054(n+2). Row sums of inverse are A105523. Product of Pascal triangle A007318 and A046854.
A106195 with an appended column of ones = A055587. Alternatively, k-th column (k=0, 1, 2) is the binomial transform of bin(n, k).
T(n,k) is the number of ideals in the fence Z(2n) having k elements of rank 1. - Emanuele Munarini, Mar 22 2011
Subtriangle of the triangle given by (0, 2, 0, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 22 2012

Examples

			Triangle begins
   1;
   2,   1;
   4,   3,   1;
   8,   8,   4,  1;
  16,  20,  13,  5,  1;
  32,  48,  38, 19,  6, 1;
  64, 112, 104, 63, 26, 7, 1;
(0, 2, 0, 0, 0, ...) DELTA (1, 0, -1/2, 1/2, 0, 0, ...) begins :
  1;
  0,  1;
  0,  2,   1;
  0,  4,   3,   1;
  0,  8,   8,   4,  1;
  0, 16,  20,  13,  5,  1;
  0, 32,  48,  38, 19,  6, 1;
  0, 64, 112, 104, 63, 26, 7, 1. - _Philippe Deléham_, Mar 22 2012
		

Crossrefs

Column 0 = 1, 2, 4...; (binomial transform of 1, 1, 1...); column 1 = 1, 3, 8, 20...(binomial transform of 1, 2, 3...); column 2: 1, 4, 13, 38...= binomial transform of bin(n, 2): 1, 3, 6...

Programs

  • Haskell
    a106195 n k = a106195_tabl !! n !! k
    a106195_row n = a106195_tabl !! n
    a106195_tabl = [1] : [2, 1] : f [1] [2, 1] where
       f us vs = ws : f vs ws where
         ws = zipWith (-) (zipWith (+) ([0] ++ vs) (map (* 2) vs ++ [0]))
                          ([0] ++ us ++ [0])
    -- Reinhard Zumkeller, Dec 16 2013
    
  • Magma
    [ (&+[Binomial(n-k, n-j)*Binomial(j, k): j in [0..n]]): k in [0..n], n in [0..10]]; // G. C. Greubel, Mar 15 2020
    
  • Maple
    T := (n, k) -> hypergeom([-n+k, k+1],[1],-1):
    seq(lprint(seq(simplify(T(n, k)), k=0..n)), n=0..7); # Peter Luschny, May 20 2015
  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + v[n - 1, x]
    v[n_, x_] := u[n - 1, x] + (x + 1) v[n - 1, x]
    Table[Factor[u[n, x]], {n, 1, z}]
    Table[Factor[v[n, x]], {n, 1, z}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]  (* A207605 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]  (* A106195 *)
    (* Clark Kimberling, Feb 19 2012 *)
    Table[Hypergeometric2F1[-n+k, k+1, 1, -1], {n, 0, 12}, {k, 0, n}]//Flatten (* G. C. Greubel, Mar 15 2020 *)
  • Maxima
    create_list(sum(binomial(i,k)*binomial(n-k,n-i),i,0,n),n,0,8,k,0,n); /* Emanuele Munarini, Mar 22 2011 */
    
  • Python
    from sympy import Poly, symbols
    x = symbols('x')
    def u(n, x): return 1 if n==1 else u(n - 1, x) + v(n - 1, x)
    def v(n, x): return 1 if n==1 else u(n - 1, x) + (x + 1)*v(n - 1, x)
    def a(n): return Poly(v(n, x), x).all_coeffs()[::-1]
    for n in range(1, 13): print(a(n)) # Indranil Ghosh, May 28 2017
    
  • Python
    from mpmath import hyp2f1, nprint
    def T(n, k): return hyp2f1(k - n, k + 1, 1, -1)
    for n in range(13): nprint([int(T(n, k)) for k in range(n + 1)]) # Indranil Ghosh, May 28 2017, after formula from Peter Luschny
    
  • Sage
    [[sum(binomial(n-k,n-j)*binomial(j,k) for j in (0..n)) for k in (0..n)] for n in (0..10)] # G. C. Greubel, Mar 15 2020

Formula

T(n,k) = Sum_{j=0..n} C(n-k,n-j)*C(j,k).
From Emanuele Munarini, Mar 22 2011: (Start)
T(n,k) = Sum_{i=0..n-k} C(k,i)*C(n-k,i)*2^(n-k-i).
T(n,k) = Sum_{i=0..n-k} C(k,i)*C(n-i,k)*(-1)^i*2^(n-k-i).
Recurrence: T(n+2,k+1) = 2*T(n+1,k+1)+T(n+1,k)-T(n,k). (End)
From Clark Kimberling, Feb 19 2012: (Start)
Define u(n,x) = u(n-1,x)+v(n-1,x), v(n,x) = u(n-1,x)+(x+1)*v(n-1,x),
where u(1,x)=1, v(1,x)=1. Then v matches A106195 and u matches A207605. (End)
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k-1). - Philippe Deléham, Mar 22 2012
T(n+k,k) is the coefficient of x^n y^k in 1/(1-2x-y+xy). - Ira M. Gessel, Oct 30 2012
T(n, k) = A208341(n+1,n-k+1), k = 0..n. - Reinhard Zumkeller, Dec 16 2013
T(n, k) = hypergeometric_2F1(-n+k, k+1, 1 , -1). - Peter Luschny, May 20 2015
G.f. 1/(1-2*x+x^2*y-x*y). - R. J. Mathar, Aug 11 2015
Sum_{k=0..n} T(n, k) = Fibonacci(2*n+2) = A088305(n+1). - G. C. Greubel, Mar 15 2020

Extensions

Edited by N. J. A. Sloane, Apr 09 2007, merging two sequences submitted independently by Gary W. Adamson and Paul Barry

A108623 G.f. satisfies x = A(x)*(1-A(x))/(1-A(x)-(A(x))^2).

Original entry on oeis.org

1, 0, -1, -1, 1, 4, 3, -8, -23, -10, 67, 153, 9, -586, -1081, 439, 5249, 7734, -7941, -47501, -53791, 105314, 430119, 343044, -1249799, -3866556, -1730017, 13996097, 34243897, 1947204, -150962373, -296101864, 121857185, 1582561870
Offset: 1

Views

Author

Christian G. Bower, Jun 12 2005

Keywords

Comments

Row sums of inverse of Riordan array (1/(1-x-x^2), x*(1-x)/(1-x-x^2)) (Cf. A053538). - Paul Barry, Nov 01 2006

Examples

			G.f. = x - x^3 - x^4 + x^5 + 4*x^6 + 3*x^7 - 8*x^8 - 23*x^9 - 10*x^10 + ...
		

Crossrefs

Except for signs, same as A108624.

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 41);
    Coefficients(R!( (1+x-Sqrt(1-2*x+5*x^2))/(2*(1-x)) )); // G. C. Greubel, Oct 20 2023
    
  • Maple
    # Using function CompInv from A357588.
    CompInv(34, n -> ifelse(n=-1, 1, combinat:-fibonacci(n-2))); # Peter Luschny, Oct 05 2022
  • Mathematica
    CoefficientList[Series[(1+x-Sqrt[1-2*x+5*x^2])/(2*x*(1-x)), {x, 0, 20}], x] (* Vaclav Kotesovec, Feb 08 2014 *)
    a[ n_] := SeriesCoefficient[ (1 + x - Sqrt[1 - 2 x + 5 x^2]) / (2 (1 - x)), {x, 0, n}]; (* Michael Somos, May 19 2014 *)
    a[ n_] := If[ n < 1, 0, SeriesCoefficient[ InverseSeries[ Series[ (x - x^2) / (1 - x - x^2), {x, 0, n}]], {x, 0, n}]]; (* Michael Somos, May 19 2014 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( (1 + x - sqrt(1 - 2*x + 5*x^2 + x^2 * O(x^n))) / (2 * (1 - x)), n))}; /* Michael Somos, May 19 2014 */
    
  • PARI
    {b(n) = if( n<1, 0, polcoeff( serreverse( (x - x^2) / (1 - x - x^2) + x * O(x^n)), n))}; /* Michael Somos, May 19 2014 */
    
  • SageMath
    def A108623_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( (1+x-sqrt(1-2*x+5*x^2))/(2*(1-x)) ).list()
    a=A108623_list(41); a[1:] # G. C. Greubel, Oct 20 2023

Formula

Binomial transform of A105523. - Paul Barry, Nov 01 2006
G.f.: (1+x-sqrt(1-2*x+5*x^2))/(2*(1-x)). - Paul Barry, Nov 01 2006
Conjecture: n*a(n) +3*(1-n)*a(n-1) +(7*n-18)*a(n-2) +5*(3-n)*a(n-3)=0. - R. J. Mathar, Nov 15 2011
Lim sup_{n->infinity} |a(n)|^(1/n) = sqrt(5). - Vaclav Kotesovec, Feb 08 2014
Series reversion of g.f. of A212804. - Michael Somos, May 19 2014
G.f.: x / (1 - x + x /(1 - x / (1 - x + x / (1 - x / ...)))). - Michael Somos, May 19 2014
0 = a(n)*(25*a(n+1) - 50*a(n+2) + 45*a(n+3) - 20*a(n+4)) + a(n+1)*(-20*a(n+1) + 34*a(n+2) - 44*a(n+3) + 25*a(n+4)) + a(n+2)*(12*a(n+2) - 2*a(n+3) - 6*a(n+4)) + a(n+3)*(a(n+4)) if n>=0. - Michael Somos, May 19 2014

A106180 Matrix inverse of number triangle A046854.

Original entry on oeis.org

1, -1, 1, 0, -1, 1, 1, -1, -1, 1, 0, 2, -2, -1, 1, -2, 2, 3, -3, -1, 1, 0, -5, 5, 4, -4, -1, 1, 5, -5, -9, 9, 5, -5, -1, 1, 0, 14, -14, -14, 14, 6, -6, -1, 1, -14, 14, 28, -28, -20, 20, 7, -7, -1, 1, 0, -42, 42, 48, -48, -27, 27
Offset: 0

Views

Author

Paul Barry, Apr 24 2005

Keywords

Comments

First column is A105523; second column is A106181.
Triangle T(n,k), 0 <= k <= n, read by rows given by [ -1, 1, -1, 1, -1, 1, -1, 1, -1, 1,...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Sep 29 2006
A124448*A007318 as infinite lower triangular matrices. - Philippe Deléham, Oct 16 2007

Examples

			Triangle begins
   1;
  -1,  1;
   0, -1,  1;
   1, -1, -1,  1;
   0,  2, -2, -1,  1;
  -2,  2,  3, -3, -1,  1;
   0, -5,  5,  4, -4, -1,  1;
		

Crossrefs

Cf. A000108.

Formula

Riordan array (1-y, y) where y=-(1-sqrt(1+4x^2))/(2x).
Sum_{k=0..n} abs(T(n,k)) = A063886(n). - Philippe Deléham, Oct 06 2006
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,k-1) - T(n-1,k+1) for k >= 1. - Philippe Deléham, Oct 27 2007
T(2n,0) = A000007(n); T(2n+2,2k+2) = -T(2n+2,2k+1) = (-1)^(n-k)*A039598(n,k); T(2n+1,2k+1) = -T(2n+1,2k) = (-1)^(n-k)*A039599(n,k). - Philippe Deléham, Oct 29 2007
Sum_{k>=0} T(m,k)*T(n,k)*(-1)^k = T(m+n,0) = A105523(m+n). - Philippe Deléham, Jan 24 2010

A165408 An aerated Catalan triangle.

Original entry on oeis.org

1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 2, 0, 1, 0, 0, 0, 3, 0, 1, 0, 0, 2, 0, 4, 0, 1, 0, 0, 0, 5, 0, 5, 0, 1, 0, 0, 0, 0, 9, 0, 6, 0, 1, 0, 0, 0, 5, 0, 14, 0, 7, 0, 1, 0, 0, 0, 0, 14, 0, 20, 0, 8, 0, 1, 0, 0, 0, 0, 0, 28, 0, 27, 0, 9, 0, 1, 0, 0, 0, 0, 14, 0, 48, 0, 35, 0, 10, 0, 1, 0, 0, 0, 0, 0, 42, 0, 75, 0, 44, 0, 11, 0, 1
Offset: 0

Views

Author

Paul Barry, Sep 17 2009

Keywords

Comments

Aeration of A120730. Row sums are A165407.
T(n,k) is the number of lattice paths from (0,0) to (k,(n-k)/2) that do not go above the diagonal x=y using steps in {(1,0), (0,1)}. - Alois P. Heinz, Sep 20 2022

Examples

			Triangle T(n,k) begins:
  1;
  0, 1;
  0, 0, 1;
  0, 1, 0, 1;
  0, 0, 2, 0,  1;
  0, 0, 0, 3,  0,  1;
  0, 0, 2, 0,  4,  0,  1;
  0, 0, 0, 5,  0,  5,  0,  1;
  0, 0, 0, 0,  9,  0,  6,  0,  1;
  0, 0, 0, 5,  0, 14,  0,  7,  0, 1;
  0, 0, 0, 0, 14,  0, 20,  0,  8, 0,  1;
  0, 0, 0, 0,  0, 28,  0, 27,  0, 9,  0, 1;
  0, 0, 0, 0, 14,  0, 48,  0, 35, 0, 10, 0, 1;
  ...
		

Crossrefs

Programs

  • Magma
    A165408:= func< n,k | n le 3*k select Binomial(Floor((n+k)/2), k)*((3*k-n)/2 +1)*(1+(-1)^(n-k))/(2*(k+1)) else 0 >;
    [A165408(n,k): k in [0..n], n in [0..15]]; // G. C. Greubel, Nov 09 2022
    
  • Maple
    b:= proc(x, y) option remember; `if`(y<=x, `if`(x=0, 1,
          b(x-1, y)+`if`(y>0, b(x, y-1), 0)), 0)
        end:
    T:= (n, k)-> `if`((n-k)::even, b(k, (n-k)/2), 0):
    seq(seq(T(n, k), k=0..n), n=0..14);  # Alois P. Heinz, Sep 20 2022
  • Mathematica
    b[x_, y_]:= b[x, y]= If[y<=x, If[x==0, 1, b[x-1, y] + If[y>0, b[x, y-1], 0]], 0];
    T[n_, k_]:= If[EvenQ[n-k], b[k, (n-k)/2], 0];
    Table[T[n, k], {n,0,14}, {k,0,n}]//Flatten (* Jean-François Alcover, Oct 08 2022, after Alois P. Heinz *)
  • SageMath
    def A165408(n,k): return 0 if (n>3*k) else binomial(int((n+k)/2), k)*((3*k-n+2)/2 )*(1+(-1)^(n-k))/(2*(k+1))
    flatten([[A165408(n,k) for k in range(n+1)] for n in range(16)]) # G. C. Greubel, Nov 09 2022

Formula

T(n,k) = if(n<=3k, C((n+k)/2, k)*((3*k-n)/2 + 1)*(1 + (-1)^(n-k))/(2*(k+1)), 0).
G.f.: 1/(1-x*y-x^3*y/(1-x^3*y/(1-x^3*y/(1-x^3*y/(1-... (continued fraction).
Sum_{k=0..n} T(n, k) = A165407(n).
From G. C. Greubel, Nov 09 2022: (Start)
Sum_{k=0..floor(n/2)} T(n-k, k) = (1+(-1)^n)*A001405(n/2)/2.
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = (1+(-1)^n)*A105523(n/2)/2.
Sum_{k=0..n} (-1)^k*T(n, k) = (-1)^n*A165407(n).
Sum_{k=0..n} 2^k*T(n, k) = A165409(n).
T(n, n-2) = A001477(n-2), n >= 2.
T(2*n, n) = (1+(-1)^n)*A174687(n/2)/2.
T(2*n, n+1) = (1-(-1)^n)*A262394(n/2)/2.
T(2*n, n-1) = (1+(-1)^n)*A236194(n/2)/2
T(3*n-2, n) = A000108(n), n >= 1. (End)
Showing 1-10 of 19 results. Next