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

A106566 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 1, 1, 1, 1, 1, 1, 1, ... ] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ... ] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 5, 5, 3, 1, 0, 14, 14, 9, 4, 1, 0, 42, 42, 28, 14, 5, 1, 0, 132, 132, 90, 48, 20, 6, 1, 0, 429, 429, 297, 165, 75, 27, 7, 1, 0, 1430, 1430, 1001, 572, 275, 110, 35, 8, 1, 0, 4862, 4862, 3432, 2002, 1001, 429, 154, 44, 9, 1
Offset: 0

Views

Author

Philippe Deléham, May 30 2005

Keywords

Comments

Catalan convolution triangle; g.f. for column k: (x*c(x))^k with c(x) g.f. for A000108 (Catalan numbers).
Riordan array (1, xc(x)), where c(x) the g.f. of A000108; inverse of Riordan array (1, x*(1-x)) (see A109466).
Diagonal sums give A132364. - Philippe Deléham, Nov 11 2007

Examples

			Triangle begins:
  1;
  0,   1;
  0,   1,   1;
  0,   2,   2,  1;
  0,   5,   5,  3,  1;
  0,  14,  14,  9,  4,  1;
  0,  42,  42, 28, 14,  5, 1;
  0, 132, 132, 90, 48, 20, 6, 1;
From _Paul Barry_, Sep 28 2009: (Start)
Production array is
  0, 1,
  0, 1, 1,
  0, 1, 1, 1,
  0, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 1, 1 (End)
		

Crossrefs

The three triangles A059365, A106566 and A099039 are the same except for signs and the leading term.
See also A009766, A033184, A059365 for other versions.
The following are all versions of (essentially) the same Catalan triangle: A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    A106566:= func< n,k | n eq 0 select 1 else (k/n)*Binomial(2*n-k-1, n-k) >;
    [A106566(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Sep 06 2021
    
  • Maple
    A106566 := proc(n,k)
        if n = 0 then
            1;
        elif k < 0 or k > n then
            0;
        else
            binomial(2*n-k-1,n-k)*k/n ;
        end if;
    end proc: # R. J. Mathar, Mar 01 2015
  • Mathematica
    T[n_, k_] := Binomial[2n-k-1, n-k]*k/n; T[0, 0] = 1; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 18 2017 *)
    (* The function RiordanArray is defined in A256893. *)
    RiordanArray[1&, #(1-Sqrt[1-4#])/(2#)&, 11] // Flatten (* Jean-François Alcover, Jul 16 2019 *)
  • PARI
    {T(n, k) = if( k<=0 || k>n, n==0 && k==0, binomial(2*n - k, n) * k/(2*n - k))}; /* Michael Somos, Oct 01 2022 */
  • Sage
    def A106566(n, k): return 1 if (n==0) else (k/n)*binomial(2*n-k-1, n-k)
    flatten([[A106566(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Sep 06 2021
    

Formula

T(n, k) = binomial(2n-k-1, n-k)*k/n for 0 <= k <= n with n > 0; T(0, 0) = 1; T(0, k) = 0 if k > 0.
T(0, 0) = 1; T(n, 0) = 0 if n > 0; T(0, k) = 0 if k > 0; for k > 0 and n > 0: T(n, k) = Sum_{j>=0} T(n-1, k-1+j).
Sum_{j>=0} T(n+j, 2j) = binomial(2n-1, n), n > 0.
Sum_{j>=0} T(n+j, 2j+1) = binomial(2n-2, n-1), n > 0.
Sum_{k>=0} (-1)^(n+k)*T(n, k) = A064310(n). T(n, k) = (-1)^(n+k)*A099039(n, k).
Sum_{k=0..n} T(n, k)*x^k = A000007(n), A000108(n), A000984(n), A007854(n), A076035(n), A076036(n), A127628(n), A126694(n), A115970(n) for x = 0,1,2,3,4,5,6,7,8 respectively.
Sum_{k>=0} T(n, k)*x^(n-k) = C(x, n); C(x, n) are the generalized Catalan numbers.
Sum_{j=0..n-k} T(n+k,2*k+j) = A039599(n,k).
Sum_{j>=0} T(n,j)*binomial(j,k) = A039599(n,k).
Sum_{k=0..n} T(n,k)*A000108(k) = A127632(n).
Sum_{k=0..n} T(n,k)*(x+1)^k*x^(n-k) = A000012(n), A000984(n), A089022(n), A035610(n), A130976(n), A130977(n), A130978(n), A130979(n), A130980(n), A131521(n) for x= 0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Aug 25 2007
Sum_{k=0..n} T(n,k)*A000108(k-1) = A121988(n), with A000108(-1)=0. - Philippe Deléham, Aug 27 2007
Sum_{k=0..n} T(n,k)*(-x)^k = A000007(n), A126983(n), A126984(n), A126982(n), A126986(n), A126987(n), A127017(n), A127016(n), A126985(n), A127053(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Oct 27 2007
T(n,k)*2^(n-k) = A110510(n,k); T(n,k)*3^(n-k) = A110518(n,k). - Philippe Deléham, Nov 11 2007
Sum_{k=0..n} T(n,k)*A000045(k) = A109262(n), A000045: Fibonacci numbers. - Philippe Deléham, Oct 28 2008
Sum_{k=0..n} T(n,k)*A000129(k) = A143464(n), A000129: Pell numbers. - Philippe Deléham, Oct 28 2008
Sum_{k=0..n} T(n,k)*A100335(k) = A002450(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A100334(k) = A001906(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A099322(k) = A015565(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A106233(k) = A003462(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A151821(k+1) = A100320(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A082505(k+1) = A144706(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A000045(2k+2) = A026671(n). - Philippe Deléham, Feb 11 2009
Sum_{k=0..n} T(n,k)*A122367(k) = A026726(n). - Philippe Deléham, Feb 11 2009
Sum_{k=0..n} T(n,k)*A008619(k) = A000958(n+1). - Philippe Deléham, Nov 15 2009
Sum_{k=0..n} T(n,k)*A027941(k+1) = A026674(n+1). - Philippe Deléham, Feb 01 2014
G.f.: Sum_{n>=0, k>=0} T(n, k)*x^k*z^n = 1/(1 - x*z*c(z)) where c(z) the g.f. of A000108. - Michael Somos, Oct 01 2022

Extensions

Formula corrected by Philippe Deléham, Oct 31 2008
Corrected by Philippe Deléham, Sep 17 2009
Corrected by Alois P. Heinz, Aug 02 2012

A147703 Triangle [1,1,1,0,0,0,...] DELTA [1,0,0,0,...] with Deléham DELTA defined in A084938.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 5, 9, 5, 1, 13, 27, 20, 7, 1, 34, 80, 73, 35, 9, 1, 89, 234, 252, 151, 54, 11, 1, 233, 677, 837, 597, 269, 77, 13, 1, 610, 1941, 2702, 2225, 1199, 435, 104, 15, 1, 1597, 5523, 8533, 7943, 4956, 2158, 657, 135, 17, 1
Offset: 0

Views

Author

Paul Barry, Nov 10 2008

Keywords

Comments

Equal to A062110*A007318 when A062110 is regarded as a triangle read by rows.

Examples

			Triangle begins
   1;
   1,   1;
   2,   3,   1;
   5,   9,   5,   1;
  13,  27,  20,   7,  1;
  34,  80,  73,  35,  9,  1;
  89, 234, 252, 151, 54, 11, 1;
		

Crossrefs

Row sums are A006012. Diagonal sums are A147704.

Programs

  • Maple
    # The function RiordanSquare is defined in A321620:
    RiordanSquare(1 / (1 - x / (1 - x / (1 - x))), 10); # Peter Luschny, Jan 26 2020
  • Mathematica
    nmax=9; Flatten[CoefficientList[Series[CoefficientList[Series[(1 - 2*x)/(1 - (3 + y)*x + (1 + y)*x^2), {x, 0, nmax}], x], {y, 0, nmax}], y]] (* Indranil Ghosh, Mar 11 2017 *)

Formula

Riordan array ((1-2x)/(1-3x+x^2), x(1-x)/(1-3x+x^2)).
Sum_{k=0..n} T(n,k)*x^k = A152239(n), A152223(n), A152185(n), A152174(n), A152167(n), A152166(n), A152163(n), A000007(n), A001519(n), A006012(n), A081704(n), A082761(n), A147837(n), A147838(n), A147839(n), A147840(n), A147841(n), for x = -8, -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively. - Philippe Deléham, Dec 01 2008
G.f.: (1-2*x)/(1-(3+y)*x+(1+y)*x^2). - Philippe Deléham, Nov 26 2011
T(n,k) = 3*T(n-1,k) + T(n-1,k-1) - T(n-2,k) - T(n-2,k-1), for n > 1. - Philippe Deléham, Feb 12 2012
The Riordan square of the odd indexed Fibonacci numbers A001519. - Peter Luschny, Jan 26 2020

A122542 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 2, -1, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 2, 4, 1, 0, 2, 8, 6, 1, 0, 2, 12, 18, 8, 1, 0, 2, 16, 38, 32, 10, 1, 0, 2, 20, 66, 88, 50, 12, 1, 0, 2, 24, 102, 192, 170, 72, 14, 1, 0, 2, 28, 146, 360, 450, 292, 98, 16, 1, 0, 2, 32, 198, 608, 1002, 912, 462, 128, 18, 1
Offset: 0

Views

Author

Philippe Deléham, Sep 19 2006, May 28 2007

Keywords

Comments

Riordan array (1, x*(1+x)/(1-x)). Rising and falling diagonals are the tribonacci numbers A000213, A001590.

Examples

			Triangle begins:
  1;
  0, 1;
  0, 2,  1;
  0, 2,  4,   1;
  0, 2,  8,   6,   1;
  0, 2, 12,  18,   8,    1;
  0, 2, 16,  38,  32,   10,   1;
  0, 2, 20,  66,  88,   50,  12,   1;
  0, 2, 24, 102, 192,  170,  72,  14,   1;
  0, 2, 28, 146, 360,  450, 292,  98,  16,  1;
  0, 2, 32, 198, 608, 1002, 912, 462, 128, 18, 1;
		

Crossrefs

Other versions: A035607, A113413, A119800, A266213.
Sums include: A000007, A001333 (row), A001590 (diagonal), A007483, A057077 (signed row), A078016 (signed diagonal), A086901, A091928, A104934, A122558, A122690.

Programs

  • Haskell
    a122542 n k = a122542_tabl !! n !! k
    a122542_row n = a122542_tabl !! n
    a122542_tabl = map fst $ iterate
       (\(us, vs) -> (vs, zipWith (+) ([0] ++ us ++ [0]) $
                          zipWith (+) ([0] ++ vs) (vs ++ [0]))) ([1], [0, 1])
    -- Reinhard Zumkeller, Jul 20 2013, Apr 17 2013
    
  • Magma
    function T(n, k) // T = A122542
      if k eq 0 then return 0^n;
      elif k eq n then return 1;
      else return T(n-1,k) + T(n-1,k-1) + T(n-2,k-1);
      end if;
    end function;
    [T(n, k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 27 2024
  • Mathematica
    CoefficientList[#, y]& /@ CoefficientList[(1-x)/(1 - (1+y)x - y x^2) + O[x]^11, x] // Flatten (* Jean-François Alcover, Sep 09 2018 *)
    (* Second program *)
    T[n_, k_]:= T[n, k]= If[k==n, 1, If[k==0, 0, T[n-1,k-1] +T[n-1,k] +T[n-2,k- 1] ]]; (* T = A122542 *)
    Table[T[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Oct 27 2024 *)
  • Sage
    def A122542_row(n):
        @cached_function
        def prec(n, k):
            if k==n: return 1
            if k==0: return 0
            return prec(n-1,k-1)+2*sum(prec(n-i,k-1) for i in (2..n-k+1))
        return [prec(n, k) for k in (0..n)]
    for n in (0..10): print(A122542_row(n)) # Peter Luschny, Mar 16 2016
    

Formula

Sum_{k=0..n} x^k*T(n,k) = A000007(n), A001333(n), A104934(n), A122558(n), A122690(n), A091928(n) for x = 0, 1, 2, 3, 4, 5. - Philippe Deléham, Jan 25 2012
Sum_{k=0..n} 3^(n-k)*T(n,k) = A086901(n).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A007483(n-1), n >= 1. - Philippe Deléham, Oct 08 2006
T(2*n, n) = A123164(n).
T(n, k) = T(n-1,k) + T(n-1,k-1) + T(n-2,k-1), n > 1. - Philippe Deléham, Jan 25 2012
G.f.: (1-x)/(1-(1+y)*x-y*x^2). - Philippe Deléham, Mar 02 2012
From G. C. Greubel, Oct 27 2024: (Start)
Sum_{k=0..n} (-1)^k*T(n, k) = A057077(n).
Sum_{k=0..floor(n/2)} T(n-k, k) = A001590(n+1).
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = A078016(n). (End)

A122950 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 0, 1, -1, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 1, -1, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 0, 2, 0, 0, 1, 3, 0, 0, 0, 3, 5, 0, 0, 0, 1, 7, 8, 0, 0, 0, 0, 4, 15, 13, 0, 0, 0, 0, 1, 12, 30, 21, 0, 0, 0, 0, 0, 5, 31, 58, 34, 0, 0, 0, 0, 0, 1, 18, 73, 109, 55, 0, 0, 0, 0, 0, 0, 6, 54, 162, 201, 89, 0, 0, 0, 0, 0, 0, 1, 25, 145, 344, 365, 144, 0, 0, 0, 0, 0, 0, 0, 7, 85, 361, 707
Offset: 0

Views

Author

Philippe Deléham, Oct 25 2006

Keywords

Comments

Skew triangle associated with the Fibonacci numbers.

Examples

			Triangle begins:
  1;
  0, 1;
  0, 0, 2;
  0, 0, 1, 3;
  0, 0, 0, 3, 5;
  0, 0, 0, 1, 7, 8;
  0, 0, 0, 0, 4, 15, 13;
  0, 0, 0, 0, 1, 12, 30, 21;
  0, 0, 0, 0, 0,  5, 31, 58,  34;
  0, 0, 0, 0, 0,  1, 18, 73, 109,  55;
  0, 0, 0, 0, 0,  0,  6, 54, 162, 201,  89;
  0, 0, 0, 0, 0,  0,  1, 25, 145, 344, 365, 144;
  0, 0, 0, 0, 0,  0,  0,  7,  85, 361, 707, 655, 233;
		

Crossrefs

Cf. A055830 (another version).

Programs

  • Mathematica
    T[0, 0] = T[1, 1] = 1; T[, 0] = T[, 1] = 0; T[n_, n_] := Fibonacci[n+1]; T[n_, k_] /; 0 <= k <= n := T[n, k] = T[n-1, k-1] + T[n-2, k-1] + T[n-2, k-2]; T[, ] = 0;
    Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, Nov 29 2018 *)

Formula

Sum_{k=0..n} T(n,k) = A011782(n).
Sum_{n>=k} T(n,k) = A001333(k).
T(n,k) = 0 if k < 0 or if k > n, T(0,0) = 1, T(2,1) = 0, T(n,k) = T(n-1,k-1) + T(n-2,k-1) + T(n-2,k-2).
T(n,n) = Fibonacci(n+1) = A000045(n+1).
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A011782(n), A133592(n), A133594(n), A133642(n), A133646(n), A133678(n), A133679(n), A133680(n), A133681(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Jan 03 2008
G.f.: (1-y*x^2)/(1-y*x-y*(y+1)*x^2). - Philippe Deléham, Nov 26 2011

A130020 Triangle T(n,k), 0<=k<=n, read by rows given by [1,0,0,0,0,0,0,...] DELTA [0,1,1,1,1,1,1,...] where DELTA is the operator defined in A084938 .

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 5, 5, 0, 1, 4, 9, 14, 14, 0, 1, 5, 14, 28, 42, 42, 0, 1, 6, 20, 48, 90, 132, 132, 0, 1, 7, 27, 75, 165, 297, 429, 429, 0, 1, 8, 35, 110, 275, 572, 1001, 1430, 1430, 0, 1, 9, 44, 154, 429, 1001, 2002, 3432, 4862, 4862, 0
Offset: 0

Views

Author

Philippe Deléham, Jun 16 2007

Keywords

Comments

Reflected version of A106566.

Examples

			Triangle begins:
  1;
  1, 0;
  1, 1,  0;
  1, 2,  2,   0;
  1, 3,  5,   5,   0;
  1, 4,  9,  14,  14,    0;
  1, 5, 14,  28,  42,   42,    0;
  1, 6, 20,  48,  90,  132,  132,    0;
  1, 7, 27,  75, 165,  297,  429,  429,    0;
  1, 8, 35, 110, 275,  572, 1001, 1430, 1430,    0;
  1, 9, 44, 154, 429, 1001, 2002, 3432, 4862, 4862,  0;
  ...
		

Crossrefs

The following are all versions of (essentially) the same Catalan triangle: A009766, A030237, A033184, A047072, A059365, A099039, A106566, this sequence.
Cf. A000108 (Catalan numbers), A106566 (row reversal), A210736.

Programs

  • Magma
    A130020:= func< n,k | n eq 0 select 1 else (n-k)*Binomial(n+k-1, k)/n >;
    [A130020(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Jun 14 2022
    
  • Mathematica
    T[n_, k_]:= (n-k)Binomial[n+k-1, k]/n; T[0, 0] = 1;
    Table[T[n, k], {n, 0, 10}, {k, 0, n}]//Flatten (* Jean-François Alcover, Jun 14 2019 *)
  • PARI
    {T(n, k) = if( k<0 || k>=n, n==0 && k==0, binomial(n+k, n) * (n-k)/(n+k))}; /* Michael Somos, Oct 01 2022 */
  • Sage
    @CachedFunction
    def A130020(n, k):
        if n==k: return add((-1)^j*binomial(n, j) for j in (0..n))
        return add(A130020(n-1, j) for j in (0..k))
    for n in (0..10) :
        [A130020(n, k) for k in (0..n)]  # Peter Luschny, Nov 14 2012
    

Formula

T(n, k) = A106566(n, n-k).
Sum_{k=0..n} T(n,k) = A000108(n).
T(n, k) = (n-k)*binomial(n+k-1, k)/n with T(0, 0) = 1. - Jean-François Alcover, Jun 14 2019
Sum_{k=0..floor(n/2)} T(n-k, k) = A210736(n). - G. C. Greubel, Jun 14 2022
G.f.: Sum_{n>=0, k>=0} T(n, k)*x^k*z^n = 1/(1 - z*c(x*z)) where c(z) = g.f. of A000108.

A094665 Another version of triangular array in A083061: triangle T(n,k), 0<=k<=n, read by rows; given by [0, 1, 3, 6, 10, 15, 21, 28, ...] DELTA [1, 2, 3, 4, 5, 6, 7, 8, ...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 1, 3, 0, 4, 15, 15, 0, 34, 147, 210, 105, 0, 496, 2370, 4095, 3150, 945, 0, 11056, 56958, 111705, 107415, 51975, 10395, 0, 349504, 1911000, 4114110, 4579575, 2837835, 945945, 135135, 0, 14873104, 85389132, 197722980, 244909665, 178378200, 77567490, 18918900, 2027025
Offset: 0

Views

Author

Philippe Deléham, Jun 07 2004, Jun 12 2007

Keywords

Comments

Define polynomials P(n,x) = x(2x+1)P(n-1,x+1) - 2x^2P(n-1,x), P(0,x) = 1. Sequence gives triangle read by rows, defined by P(n,x) = Sum_{k = 0..n} T(n,k)*x^k. - Philippe Deléham, Jun 20 2004
From Johannes W. Meijer, May 24 2009: (Start)
In A160464 we defined the coefficients of the ES1 matrix by ES1[2*m-1,n=1] = 2*eta(2*m-1) and the recurrence relation ES1[2*m-1,n] = ((2*n-2)/(2*n-1))*(ES1[2*m-1,n-1] - ES1[2*m-3,n-1]/(n-1)^2) for m the positive and negative integers and n >= 1. As usual eta(m) = (1-2^(1-m))*zeta(m) with eta(m) the Dirichlet eta function and zeta(m) the Riemann zeta function. It is well-known that ES1[1-2*m,n=1] = (4^m-1)*(-bernoulli(2*m))/m for m >= 1. and together with the recurrence relation this leads to ES1[-1,n] = 0.5 for n >= 1.
We discovered that the n-th term of the row coefficients ES1[1-2*m,n] for m >= 1, can be generated with the rather simple polynomials RES1(1-2*m,n) = (-1)^(m+1)*ECGP(1-2*m, n)/2^m. This discovery was enabled by the recurrence relation for the RES1(1-2*m,n) which we derived from the recurrence relation for the ES1[2*m-1,n] coefficients and the fact that RES1(-1,n) = 0.5. The coefficients of the ECGP(1-2*m,n) polynomials led to this triangle and subsequently to triangle A083061. (End)
From David Callan, Jan 03 2011: (Start)
T(n,k) is the number of increasing 0-2 trees (A002105) on 2n edges in which the minimal path from the root has length k.
Proof. The number a(n,k) of such trees satisfies the recurrence a(0,0)=1, a(1,1)=1 and, counting by size of the subtree rooted at the smaller child of the root,
a(n,k) = Sum_{j=1..n-1} C(2n-1,j)*a(j,k-1)*a(n-1-j)
for 2<=k<=n, where a(n) = Sum_{k>=0} a(n,k) is the reduced tangent number A002105 (indexed from 0). The recurrence translates into the differential equation
F_x(x,y) = y*F(x,y)*G(x)
for the GF F(x,y) = Sum_{n,k>=0} a(n,k)x^(2n)/(2n)!*y^k, where G(x):=Sum_{n>=0} a(n)x^(2n+1)/(2n+1)! is known to be sqrt(2)*tan(x/sqrt(2)). The differential equation has solution F(x,y) = sec(x/sqrt(2))^(2y). (End)

Examples

			Triangle begins:
.1;
.0, 1;
.0, 1, 3;
.0, 4, 15, 15;
.0, 34, 147, 210, 105;
.0, 496, 2370, 4095, 3150, 945;
.0, 11056, 56958, 111705, 107415, 51975, 10395;
.0, 349504, 1911000, 4114110, 4579575, 2837835, 945945, 135135;
From _Johannes W. Meijer_, May 24 2009: (Start)
The first few ECGP(1-2*m,n) polynomials are: ECGP(-1,n) = 1; ECGP(-3,n) = n; ECGP(-5,n) = n + 3*n^2; ECGP(-7,n) = 4*n + 15*n^2+ 15*n^3 .
The first few RES1(1-2*m,n) are: RES1(-1,n) = (1/2)*(1); RES1(-3,n) = (-1/4)*(n); RES1(-5,n) = (1/8)*(n+3*n^2); RES1(-7,n) = (-1/16)*(4*n+15*n^2+15*n^3).
(End)
		

Crossrefs

From Johannes W. Meijer, May 24 2009 and Jun 27 2009: (Start)
A001147, A001880, A160470, A160471 and A160472 are the first five right hand columns.
Appears in A162005, A162006 and A162007.
(End)

Programs

  • Maple
    nmax:=7; imax := nmax: T1(0, x) := 1: T1(0, x+1) := 1: for i from 1 to imax do T1(i, x) := expand((2*x+1) * (x+1) * T1(i-1, x+1) - 2 * x^2 * T1(i-1, x)): dx:=degree(T1(i, x)): for k from 0 to dx do c(k) := coeff(T1(i, x), x, k) od: T1(i, x+1) := sum(c(j1)*(x+1)^(j1), j1=0..dx) od: for i from 0 to imax do for j from 0 to i do A083061(i, j) := coeff(T1(i, x), x, j) od: od: for n from 0 to nmax do for k from 0 to n do T(n+1, k+1) := A083061(n, k) od: od: T(0, 0):=1: for n from 1 to nmax do T(n, 0):=0 od: seq(seq(T(n, k), k=0..n), n=0..nmax);
    # Johannes W. Meijer, Jun 27 2009, revised Sep 23 2012
  • Mathematica
    nmax = 8;
    T[n_, k_] := SeriesCoefficient[Sec[x/Sqrt[2]]^(2y), {x, 0, 2n}, {y, 0, k}]* (2n)!;
    Table[T[n, k], {n, 0, nmax}, {k, 0, n}] // Flatten (* Jean-François Alcover, Aug 10 2018 *)

Formula

Sum_{k = 0..n} T(n, k) = A002105(n+1).
Sum_{k = 0..n} T(n, k)*2^(n-k) = A000364(n); Euler numbers.
Sum_{k = 0..n} T(n, k)*(-2)^(n-k) = 1.
RES1(1-2*m,n) = n^2*RES1(3-2*m,n)-n*(2*n+1)*RES1(3-2*m,n+1)/2 for m >= 2, with RES1(-1,n) = 0.5 for n >= 1. - Johannes W. Meijer, May 24 2009
G.f.: Sum_{n,k>=0} T(n,k)x^n/n!*y^k = sec(x/sqrt(2))^(2y).

Extensions

Term corrected by Johannes W. Meijer, Sep 23 2012

A167374 Triangle, read by rows, given by [ -1,1,0,0,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

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

Views

Author

Philippe Deléham, Nov 02 2009

Keywords

Comments

Riordan array (1-x,1) read by rows; Riordan inverse is (1/(1-x),1). Columns have g.f. (1-x)x^k. Diagonal sums are A033999. Unsigned version in A097806.
Table T(n,k) read by antidiagonals. T(n,1) = 1, T(n,2) = -1, T(n,k) = 0, k > 2. - Boris Putievskiy, Jan 17 2013
Finite difference operator (pair difference): left multiplication by T of a sequence arranged as a column vector gives a running forward difference, a(k+1)-a(k), or first finite difference (modulo sign), of the elements of the sequence. T^n gives the n-th finite difference (mod sign). T is the inverse of the summation matrix A000012 (regarded as lower triangular matrices). - Tom Copeland, Mar 26 2014

Examples

			Triangle begins:
   1;
  -1,  1;
   0, -1,  1;
   0,  0, -1,  1;
   0,  0,  0, -1,  1;
   0,  0,  0,  0, -1,  1; ...
Row number r (r>4) contains (r-2) times '0', then '-1' and '1'.
From _Boris Putievskiy_, Jan 17 2013: (Start)
The start of the sequence as a table:
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  ...
(End)
		

Crossrefs

Programs

  • Maple
    A167374 := proc(n,k)
        if k> n or k < n-1 then
            0;
        elif k = n then
            1;
        else
            -1 ;
        end if;
    end proc: # R. J. Mathar, Sep 07 2016
  • Mathematica
    Table[PadLeft[{-1, 1}, n], {n, 13}] // Flatten (* or *)
    MapIndexed[Take[#1, First@ #2] &, CoefficientList[Series[(1 - x)/(1 - x y), {x, 0, 12}], {x, y}]] // Flatten (* Michael De Vlieger, Nov 16 2016 *)
    T[n_, k_] := If[ k<0 || k>n, 0, Boole[n==k] - Boole[n==k+1]]; (* Michael Somos, Oct 01 2022 *)
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, (n==k) - (n==k+1))}; /* Michael Somos, Oct 01 2022 */

Formula

Sum_{k, 0<=k<=n} T(n,k)*x^k = A000007(n), A011782(n), A025192(n), A002001(n), A005054(n), A052934(n), A055272(n), A055274(n), A055275(n), A055268(n), A055276(n) for x = 1,2,3,4,5,6,7,8,9,10,11 respectively .
From Boris Putievskiy, Jan 17 2013: (Start)
a(n) = floor((A002260(n)+2)/(A003056(n)+2))*(-1)^(A002260(n)+A003056(n)+1), n>0.
a(n) = floor((i+2)/(t+2))*(-1)^(i+t+1), n > 0, where
i = n - t*(t+1)/2,
t = floor((-1 + sqrt(8*n-7))/2). (End)
T*A000012 = Identity matrix. T*A007318 = A097805. T*(A007318)^(-1)= signed A029653. - Tom Copeland, Mar 26 2014
G.f.: (1-x)/(1-x*y). - R. J. Mathar, Aug 11 2015
T = A130595*A156644 = M*T^(-1)*M = M*A000012*M, where M(n,k) = (-1)^n A130595(n,k). Note that M = M^(-1). Cf. A118800 and A097805. - Tom Copeland, Nov 15 2016

A152815 Triangle T(n,k), read by rows given by [1,0,-1,0,0,0,0,0,0,...] DELTA [0,1,-1,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 2, 1, 0, 0, 1, 2, 1, 0, 0, 0, 1, 3, 3, 1, 0, 0, 0, 1, 3, 3, 1, 0, 0, 0, 0, 1, 4, 6, 4, 1, 0, 0, 0, 0, 1, 4, 6, 4, 1, 0, 0, 0, 0, 0, 1, 5, 10, 10, 5, 1, 0, 0, 0, 0, 0, 1, 5, 10, 10, 5, 1, 0, 0, 0, 0, 0, 0, 1, 6, 15, 20, 15, 6, 1, 0, 0, 0, 0, 0, 0, 1, 6, 15, 20, 15, 6, 1, 0, 0, 0
Offset: 0

Views

Author

Philippe Deléham, Dec 13 2008

Keywords

Comments

Triangle read by rows, Pascal's triangle (A007318) rows repeated.
Riordan array (1/(1-x), x^2/(1-x^2)). - Philippe Deléham, Feb 27 2012

Examples

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

Crossrefs

Cf. A007318, A064861, A152198 (another version), A000931 (diagonal sums), A016116 (row sums).

Programs

  • Haskell
    a152815 n k = a152815_tabl !! n !! k
    a152815_row n = a152815_tabl !! n
    a152815_tabl = [1] : [1,0] : t [1,0] where
       t ys = zs : zs' : t zs' where
         zs' = zs ++ [0]; zs = zipWith (+) ([0] ++ ys) (ys ++ [0])
    -- Reinhard Zumkeller, Feb 28 2012
    
  • Mathematica
    m = 13;
    (* DELTA is defined in A084938 *)
    DELTA[Join[{1, 0, -1}, Table[0, {m}]], Join[{0, 1, -1}, Table[0, {m}]], m] // Flatten (* Jean-François Alcover, Feb 19 2020 *)
    T[n_, k_] := If[n<0, 0, Binomial[Floor[n/2], k]]; (* Michael Somos, Oct 01 2022 *)
  • PARI
    {T(n, k) = if(n<0, 0, binomial(n\2, k))}; /* Michael Somos, Oct 01 2022 */

Formula

T(n,k) = T(n-1,k) + ((1+(-1)^n)/2)*T(n-1,k-1).
G.f.: (1+x)/(1-(1+y)*x^2).
Sum_{k=0..n} T(n,k)*x^k = A000012(n), A016116(n), A108411(n), A213173(n), A074872(n+1) for x = 0,1,2,3,4 respectively. - Philippe Deléham, Nov 26 2011, Apr 22 2013

Extensions

Example corrected by Philippe Deléham, Dec 13 2008

A111146 Triangle T(n,k), read by rows, given by [0, 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, ...] DELTA [1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 0, 2, 0, 0, 1, 4, 0, 0, 2, 5, 8, 0, 0, 6, 15, 17, 16, 0, 0, 24, 62, 68, 49, 32, 0, 0, 120, 322, 359, 243, 129, 64, 0, 0, 720, 2004, 2308, 1553, 756, 321, 128, 0, 0, 5040, 14508, 17332, 11903, 5622, 2151, 769, 256, 0, 0, 40320, 119664
Offset: 0

Views

Author

Philippe Deléham, Oct 19 2005

Keywords

Comments

Let R(m,n,k), 0<=k<=n, the Riordan array (1, x*g(x)) where g(x) is g.f. of the m-fold factorials . Then Sum_{k, 0<=k<=n} = R(m,n,k) = Sum_{k, 0<=k<=n} T(n,k)*m^(n-k).
For m = -1, R(-1,n,k) is A026729(n,k).
For m = 0, R(0,n,k) is A097805(n,k).
For m = 1, R(1,n,k) is A084938(n,k).
For m = 2, R(2,n,k) is A111106(n,k).

Examples

			Triangle begins:
.1;
.0, 1;
.0, 0, 2;
.0, 0, 1, 4;
.0, 0, 2, 5, 8;
.0, 0, 6, 15, 17, 16;
.0, 0, 24, 62, 68, 49, 32;
.0, 0, 120, 322, 359, 243, 129, 64;
.0, 0, 720, 2004, 2308, 1553, 756, 321, 128;
.0, 0, 5040, 14508, 17332, 11903, 5622, 2151, 769, 256;
.0, 0, 40320, 119664, 148232, 105048, 49840, 18066, 5756, 1793, 512;
....................................................................
At y=2: Sum_{k=0..n} 2^k*T(n,k) = A113327(n) where (1 + 2*x + 8*x^2 + 36*x^3 +...+ A113327(n)*x^n +..) = 1/(1 - 2/1!*x*(1! + 2!*x + 3!*x^2 + 4!*x^3 +..) ).
At y=3: Sum_{k=0..n} 3^k*T(n,k) = A113328(n) where (1 + 3*x + 18*x^2 + 117*x^3 +...+ A113328(n)*x^n +..) = 1/(1 - 3/2!*x*(2! + 3!*x + 4!*x^2 + 5!*x^3 +..) ).
At y=4: Sum_{k=0..n} 4^k*T(n,k) = A113329(n) where (1 + 4*x + 32*x^2 + 272*x^3 +...+ A113329(n)*x^n +..) = 1/(1 - 4/3!*x*(3! + 4!*x + 5!*x^2 + 6!*x^3 +..) ).
		

Crossrefs

Cf. m-fold factorials : A000142, A001147, A007559, A007696, A008548, A008542.
Cf. A113326, A113327 (y=2), A113328 (y=3), A113329 (y=4), A113330 (y=5), A113331 (y=6).

Programs

  • Mathematica
    T[n_, k_] := Module[{x = X + X*O[X]^n, y = Y + Y*O[Y]^k}, A = 1/(1 - x*y*Sum[x^j*Product[y + i, {i, 0, j - 1}], {j, 0, n}]); Coefficient[ Coefficient[A, X, n], Y, k]];
    Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 26 2019, from PARI *)
  • PARI
    {T(n,k)=local(x=X+X*O(X^n),y=Y+Y*O(Y^k)); A=1/(1-x*y*sum(j=0,n,x^j*prod(i=0,j-1,y+i))); return(polcoeff(polcoeff(A,n,X),k,Y))} (Hanna)

Formula

Sum_{k, 0<=k<=n} (-1)^(n-k)*T(n, k) = A000045(n+1), Fibonacci numbers.
Sum_{k, 0<=k<=n} T(n, k) = A051295(n).
Sum_{k, 0<=k<=n} 2^(n-k)*T(n, k) = A112934(n).
T(0, 0) = 1, T(n, n) = 2^(n-1).
G.f.: A(x, y) = 1/(1 - x*y*Sum_{j>=0} (y-1+j)!/(y-1)!*x^j ). - Paul D. Hanna, Oct 26 2005

A165253 Triangle T(n,k), read by rows given by [1,0,1,0,0,0,0,0,0,...] DELTA [0,1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 6, 5, 1, 0, 1, 10, 15, 7, 1, 0, 1, 15, 35, 28, 9, 1, 0, 1, 21, 70, 84, 45, 11, 1, 0, 1, 28, 126, 210, 165, 66, 13, 1, 0, 1, 36, 210, 462, 495, 286, 91, 15, 1, 0, 1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1, 0, 1, 55, 495, 1716, 3003, 3003, 1820, 680
Offset: 0

Views

Author

Philippe Deléham, Sep 10 2009

Keywords

Comments

Mirror image of triangle in A121314.

Examples

			Triangle begins:
  1;
  1,    0;
  1,    1,    0;
  1,    3,    1,    0;
  1,    6,    5,    1,    0;
  1,   10,   15,    7,    1,    0;
  1,   15,   35,   28,    9,    1,    0;
  1,   21,   70,   84,   45,   11,    1,    0;
  1,   28,  126,  210,  165,   66,   13,    1,    0;
  1,   36,  210,  462,  495,  286,   91,   15,    1,    0,
  1,   45,  330,  924, 1287, 1001,  455,  120,   17,    1,    0;
		

Crossrefs

Programs

  • Mathematica
    m = 13;
    (* DELTA is defined in A084938 *)
    DELTA[Join[{1, 0, 1}, Table[0, {m}]], Join[{0, 1}, Table[0, {m}]], m] // Flatten (* Jean-François Alcover, Feb 19 2020 *)

Formula

T(0,0)=1, T(n,k) = binomial(n-1+k,2k) for n >= 1.
Sum {k=0..n} T(n,k)*x^k = A000012(n), A001519(n), A001835(n), A004253(n), A001653(n), A049685(n-1), A070997(n-1), A070998(n-1), A072256(n), A078922(n), A077417(n-1), A085260(n), A001570(n) for x = 0,1,2,3,4,5,6,7,8,9,10,11,12 respectively.
Sum_{k=0..n} T(n,k)*x^(n-k) = A000007(n), A001519(n), A047849(n), A165310(n), A165311(n), A165312(n), A165314(n), A165322(n), A165323(n), A165324(n) for x= 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Sep 26 2009
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k), T(0,0)=T(1,0)=1, T(1,1)=0. - Philippe Deléham, Feb 18 2012
G.f.: (1-x-y*x)/((1-x)^2-y*x). - Philippe Deléham, Feb 19 2012
Showing 1-10 of 637 results. Next