A112556 -id:A112556 - cp's OEIS frontend

A112555 Triangle T, read by rows, such that the m-th matrix power satisfies T^m = I + m*(T - I) and consequently the matrix logarithm satisfies log(T) = T - I, where I is the identity matrix.

1, 1, 1, -1, 0, 1, 1, 1, 1, 1, -1, -2, -2, 0, 1, 1, 3, 4, 2, 1, 1, -1, -4, -7, -6, -3, 0, 1, 1, 5, 11, 13, 9, 3, 1, 1, -1, -6, -16, -24, -22, -12, -4, 0, 1, 1, 7, 22, 40, 46, 34, 16, 4, 1, 1, -1, -8, -29, -62, -86, -80, -50, -20, -5, 0, 1, 1, 9, 37, 91, 148, 166, 130, 70, 25, 5, 1, 1, -1, -10, -46, -128, -239, -314, -296, -200, -95, -30, -6, 0

Offset: 0

Paul D. Hanna, Sep 21 2005

Signed version of A108561. Row sums equal A084247. The n-th unsigned row sum = A001045(n) + 1 (Jacobsthal numbers). Central terms of even-indexed rows are a signed version of A072547. Sums of squared terms in rows yields A112556, which equals the first differences of the unsigned central terms.
Equals row reversal of triangle A112468 up to sign, where A112468 is the Riordan array (1/(1-x),x/(1+x)). - Paul D. Hanna, Jan 20 2006
The elements here match A108561 in absolute value, but the signs are crucial to the properties that the matrix A112555 exhibits; the main property being T^m = I + m*(T - I). This property is not satisfied by A108561. - Paul D. Hanna, Nov 10 2009
Eigensequence of the triangle = A140165. - Gary W. Adamson, Jan 30 2009
Triangle T(n,k), read by rows, given by [1,-2,0,0,0,0,0,0,0,...] DELTA [1,0,-1,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Sep 17 2009

			Triangle T begins:
   1;
   1,   1;
  -1,   0,   1;
   1,   1,   1,   1;
  -1,  -2,  -2,   0,   1;
   1,   3,   4,   2,   1,   1;
  -1,  -4,  -7,  -6,  -3,   0,   1;
   1,   5,  11,  13,   9,   3,   1,   1;
  -1,  -6, -16, -24, -22, -12,  -4,   0,   1;
   1,   7,  22,  40,  46,  34,  16,   4,   1,   1;
  -1,  -8, -29, -62, -86, -80, -50, -20,  -5,   0,   1;
  ...
Matrix log, log(T) = T - I, begins:
   0;
   1,  0;
  -1,  0,  0;
   1,  1,  1,  0;
  -1, -2, -2,  0,  0;
   1,  3,  4,  2,  1,  0;
  -1, -4, -7, -6, -3,  0,  0;
  ...
Matrix inverse, T^-1 = 2*I - T, begins:
   1;
  -1,  1;
   1,  0,  1;
  -1, -1, -1,  1;
   1,  2,  2,  0,  1;
  -1, -3, -4, -2, -1,  1;
  ...
where adjacent sums in row n of T^-1 gives row n+1 of T.

Paul D. Hanna, Table of n, a(n) for n = 0..1080

Cf. A112468 (reversed rows), A108561, A084247, A001045, A072547, A112556, A279006, A140165.
From Philippe Deléham, Oct 07 2009: (Start)
Sum_{k=0..n} T(n, k)*x^(n-k) = A165760(n), A165759(n), A165758(n), A165755(n), A165752(n), A165746(n), A165751(n), A165747(n), A000007(n), A000012(n), A084247(n), A165553(n), A165622(n), A165625(n), A165638(n), A165639(n), A165748(n), A165749(n), A165750(n) for x= -9,-8,-7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9 respectively.
Sum_{k=0..n} T(n, k)*x^k = A166157(n), A166153(n), A166152(n), A166149(n), A166036(n), A166035(n), A091004(n+1), A077925(n), A000007(n), A165326(n), A084247(n), A165405(n), A165458(n), A165470(n), A165491(n), A165505(n), A165506(n), A165510(n), A165511(n) for x = -9,-8,-7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9 respectively.  (End)

Clear[t]; t[0, 0] = 1; t[n_, 0] = (-1)^(Mod[n, 2]+1); t[n_, n_] = 1; t[n_, k_] /; k == n-1 := t[n, k] = Mod[n, 2]; t[n_, k_] /; 0 < k < n-1 := t[n, k] = -t[n-1, k] - t[n-1, k-1]; Table[t[n, k], {n, 0, 13}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 06 2013 *)

{T(n,k)=local(x=X+X*O(X^n),y=Y+Y*O(Y^k)); polcoeff( polcoeff( (1+2*x+x*y)/((1-x*y)*(1+x+x*y)),n,X),k,Y)}
for(n=0,12, for(k=0,n, print1(T(n,k),", "));print(""))

{T(n,k)=local(m=1,x=X+X*O(X^n),y=Y+Y*O(Y^k)); polcoeff(polcoeff(1/(1-x*y) + m*x/((1-x*y)*(1+x+x*y)),n,X),k,Y)}
for(n=0,12, for(k=0,n, print1(T(n,k),", "));print(""))

def A112555_row(n):
    @cached_function
    def prec(n, k):
        if k==n: return 1
        if k==0: return 0
        return -prec(n-1,k-1)-sum(prec(n,k+i-1) for i in (2..n-k+1))
    return [(-1)^(n-k+1)*prec(n+1, k) for k in (1..n+1)]
for n in (0..12): print(A112555_row(n)) # Peter Luschny, Mar 16 2016

G.f.: 1/(1-x*y) + x/((1-x*y)*(1+x+x*y)).
The m-th matrix power T^m has the g.f.: 1/(1-x*y) + m*x/((1-x*y)*(1+x+x*y)).
Recurrence: T(n, k) = [T^-1](n-1, k) + [T^-1](n-1, k-1), where T^-1 is the matrix inverse of T.
From Peter Bala, Jun 23 2025: (Start)
T^z = exp(z*log(T)) = I + z*(T - I) for arbitrary complex z, where I is the identity array.
exp(T) = e*T. More generally, exp(z * T^u) = exp(z)*T^(u*z) = exp(z)*I + u*z*exp(z)*(T - I).
sin(z * T^u) =  sin(z)*I + u*z*cos(z)*(T - I).
cos(z * T^u) =  cos(z)*I - u*z*sin(z)*(T - I).
tan(z * T^u) =  tan(z)*I + u*z*sec(z)^2*(T - I).
Chebyshev_T(n, T^u) = I + (n^2)*u*(T - I) and
Legendre_P(n, T^u) = I + (n*(n+1)/2)*u*(T - I).
More generally, for n >= 1,
Chebyshev_T(n, z*T^u) = Chebyshev_T(n, z)*I + n*u*z*Chebyshev_U(n-1, z)*(T - I) and
Legendre_P(n, z*T^u) = Legendre_P(n, z)*I + u*Q(n, z)*(T - I), where Q(1, z) = z and Q(n, z) = n*Legendre_P(n, z) + Q(n-1, z)/z for n > 1.
All the above properties may also hold for the triangle A279006. (End)

A112468 Riordan array (1/(1-x), x/(1+x)).

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 1, 1, -1, 1, 1, 0, 2, -2, 1, 1, 1, -2, 4, -3, 1, 1, 0, 3, -6, 7, -4, 1, 1, 1, -3, 9, -13, 11, -5, 1, 1, 0, 4, -12, 22, -24, 16, -6, 1, 1, 1, -4, 16, -34, 46, -40, 22, -7, 1, 1, 0, 5, -20, 50, -80, 86, -62, 29, -8, 1, 1, 1, -5, 25, -70, 130, -166, 148, -91, 37, -9, 1, 1, 0, 6, -30, 95, -200, 296, -314, 239, -128, 46, -10, 1

Offset: 0

Paul Barry, Sep 06 2005

Row sums are A040000. Diagonal sums are A112469. Inverse is A112467. Row sums of k-th power are 1, k+1, k+1, k+1, .... Note that C(n,k) = Sum_{j=0..n-k} C(n-j-1, n-k-j).
Equals row reversal of triangle A112555 up to sign, where log(A112555) = A112555 - I. Unsigned row sums equals A052953 (Jacobsthal numbers + 1). Central terms of even-indexed rows are a signed version of A072547. Sums of squared terms in rows yields A112556, which equals the first differences of the unsigned central terms. - Paul D. Hanna, Jan 20 2006
Sum_{k=0..n} T(n,k)*x^k = A000012(n), A040000(n), A005408(n), A033484(n), A048473(n), A020989(n), A057651(n), A061801(n), A238275(n), A238276(n), A138894(n), A090843(n), A199023(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 respectively (see the square array in A112739). - Philippe Deléham, Feb 22 2014

			Triangle starts
  1;
  1,  1;
  1,  0,  1;
  1,  1, -1,  1;
  1,  0,  2, -2,  1;
  1,  1, -2,  4, -3,  1;
  1,  0,  3, -6,  7, -4,  1;
Matrix log begins:
  0;
  1,  0;
  1,  0,  0;
  1,  1, -1,  0;
  1,  1,  1, -2,  0;
  1,  1,  1,  1, -3,  0; ...
Production matrix begins
  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.
- _Paul Barry_, Apr 08 2011

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
H. Belbachir and F. Bencherif, On some properties of bivariate Fibonacci and Lucas Polynomials, JIS 11 (2008) 08.2.6.
Hacene Belbachir and Athmane Benmezai, Expansion of Fibonacci and Lucas Polynomials: An Answer to Prodinger's Question, Journal of Integer Sequences, Vol. 15 (2012), #12.7.6.
Emeric Deutsch, L. Ferrari and S. Rinaldi, Production Matrices, Advances in Applied Mathematics, 34 (2005) pp. 101-122.
Kyu-Hwan Lee and Se-jin Oh, Catalan triangle numbers and binomial coefficients, arXiv:1601.06685 [math.CO], 2016.
H. Prodinger, On the expansion of Fibonacci and Lucas Polynomials, JIS 12 (2009) 09.1.6.

Cf. A174294, A174295, A174296, A174297. - Mats Granvik, Mar 15 2010
Cf. A072547 (central terms), A112555 (reversed rows), A112465, A052953, A112556, A112739, A119258.
See A279006 for another version.

T:= function(n,k)
    if k=0 or k=n then return 1;
    else return T(n-1,k-1) - T(n-1,k);
    fi;
  end;
Flat(List([0..12], n-> List([0..n], k-> T(n,k) ))); # G. C. Greubel, Nov 13 2019

a112468 n k = a112468_tabl !! n !! k
a112468_row n = a112468_tabl !! n
a112468_tabl = iterate (\xs -> zipWith (-) ([2] ++ xs) (xs ++ [0])) [1]
-- Reinhard Zumkeller, Jan 03 2014

function T(n,k)
  if k eq 0 or k eq n then return 1;
  else return T(n-1,k-1) - T(n-1,k);
  end if;
  return T;
end function;
[T(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Nov 13 2019

T := (n,k,m) -> (1-m)^(-n+k)-m^(k+1)*pochhammer(n-k,k+1)*hypergeom( [1,n+1],[k+2],m)/(k+1)!; A112468 := (n,k) -> T(n,n-k,-1);
seq(print(seq(simplify(A112468(n,k)),k=0..n)),n=0..10); # Peter Luschny, Jul 25 2014

T[n_, 0] = 1; T[n_, n_] = 1; T[n_, k_ ]:= T[n, k] = T[n-1, k-1] - T[n-1, k]; Table[T[n, k], {n, 0, 12}, {k, 0, n}]//Flatten (* Jean-François Alcover, Mar 06 2013 *)

{T(n,k)=local(m=1,x=X+X*O(X^n),y=Y+Y*O(Y^k)); polcoeff(polcoeff((1+(m-1)*x)*(1+m*x)/(1+m*x-x*y)/(1-x),n,X),k,Y)} \\ Paul D. Hanna, Jan 20 2006

T(n,k) = if(k==0 || k==n, 1, T(n-1, k-1) - T(n-1, k)); \\ G. C. Greubel, Nov 13 2019

@CachedFunction
def T(n, k):
    if (k<0 or n<0): return 0
    elif (k==0 or k==n): return 1
    else: return T(n-1, k-1) - T(n-1, k)
[[T(n, k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Nov 13 2019

Triangle T(n,k) read by rows: T(n,0)=1, T(n,k) = T(n-1,k-1) - T(n-1,k). - Mats Granvik, Mar 15 2010
Number triangle T(n, k)= Sum_{j=0..n-k} C(n-j-1, n-k-j)*(-1)^(n-k-j).
G.f. of matrix power T^m: (1+(m-1)*x)*(1+m*x)/(1+m*x-x*y)/(1-x). G.f. of matrix log: x*(1-2*x*y+x^2*y)/(1-x*y)^2/(1-x). - Paul D. Hanna, Jan 20 2006
T(n, k) = R(n,n-k,-1) where R(n,k,m) = (1-m)^(-n+k)-m^(k+1)*Pochhammer(n-k,k+1)*hyper2F1([1,n+1],[k+2],m)/(k+1)!. - Peter Luschny, Jul 25 2014

A116466 Unsigned row sums of triangle A114700.

Original entry on oeis.org

1, 2, 2, 4, 2, 4, 4, 8, 10, 20, 32, 64, 112, 224, 408, 816, 1514, 3028, 5680, 11360, 21472, 42944, 81644, 163288, 311896, 623792, 1196132, 2392264, 4602236, 9204472, 17757184, 35514368, 68680170, 137360340, 266200112, 532400224, 1033703056

Offset: 0

Paul D. Hanna, Feb 19 2006

nonn

Both triangles A112555 and A114700 have the property that the m-th matrix power of the triangles satisfy T^m = I + m*(T - I). So it is curious that the row squared sums of A112555 is a bisection of the unsigned row sums of A114700.

Cf. A112556, A112555, A114700.

CoefficientList[Series[(1 + 2*x)*(2*(1 + x^2)/(1 - x^2) + x^2/(1 - 4*x^2)^(1/2))/(2 + x^2), {x, 0, 50}], x] (* Wesley Ivan Hurt, Feb 21 2017 *)

a(n)=local(x=X+X*O(X^n)); polcoeff((1+2*x)*(2*(1+x^2)/(1-x^2)+x^2/(1-4*x^2)^(1/2))/(2+x^2),n,X)

/* a(n) as the unsigned row sums of A114700 */ a(n)=sum(k=0,n,abs(polcoeff(polcoeff(1/(1-x*y)+ x*(1+x-2*x^2*y)/(1-x)/(1+x+x*y+x*O(x^n)+y*O(y^k))/(1-x*y),n,x),k,y)))

G.f.: (1+2*x)*( 2*(1+x^2)/(1-x^2) + x^2/(1-4*x^2)^(1/2) )/(2+x^2). Also, a(2*n+1) = 2*a(2*n), a(2*n) = A112556(n), where A112556 equals the row squared sums of triangle A112555.

cp's OEIS Frontend

A112555 Triangle T, read by rows, such that the m-th matrix power satisfies T^m = I + m*(T - I) and consequently the matrix logarithm satisfies log(T) = T - I, where I is the identity matrix.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Sage

Formula

A112468 Riordan array (1/(1-x), x/(1+x)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A116466 Unsigned row sums of triangle A114700.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

PARI

PARI

Formula