Pascal Eigenspaces and Invariant Sequences of the First or Second Kind

Size: px
Start display at page:

Download "Pascal Eigenspaces and Invariant Sequences of the First or Second Kind"

Transcription

1 Pascal Eigenspaces and Invariant Sequences of the First or Second Kind I-Pyo Kim a,, Michael J Tsatsomeros b a Department of Mathematics Education, Daegu University, Gyeongbu, 38453, Republic of Korea b Department of Mathematics and Statistics, Washington State University, Pullman, WA 99164, USA Abstract An infinite real sequence a n } is called an invariant sequence of the first (resp, second ind if a n = n ( n = ( 1 a (resp, a n = ( =n n ( 1 a We review and investigate invariant sequences of the first and second ind, and study their relationships using similarities of Pascal-type matrices and their eigenspaces Keywords: Invariant sequence, Pascal matrix, Eigenvalue, Eigenvector 1 MSC: 15B18, 11B39, 11B65 1 Introduction Inverse relations play an important role in combinatorics [11] The binomial inversion formula, which states that for sequences a n } and b n } (n =, 1,,, a n = n = ( n ( 1 b if and only if b n = n = ( n ( 1 a, (11 is a typical inverse relation considered in [6, 8, 1, 14, 15, 16] Specifically, (11 motivated Sun [14] to investigate the following sequences Definition 11 Let a n } (n =, 1,, be a sequence such that ( 1 s 1 a n = n = ( n ( 1 a ; s = 1 or s = (1 5 We refer to a n } as an invariant sequence (when s = 1 or an inverse invariant sequence (when s = of the first ind Research supported by Daegu University Research Grant 13 Corresponding author addresses: imipyo@daeguacr (I-Pyo Kim, tsat@mathwsuedu (Michael J Tsatsomeros Preprint submitted to Journal of LATEX Templates September 1, 17

2 Several examples of invariant sequences of the first ind can be found in [14], including 1 n }, nf n 1}, L n }, ( 1 n B n } (n, where F 1 = and F n }, L n }, and B n } are the Fibonacci sequence, Lucas sequence, and Bernoulli numbers [7], respectively In this paper, we will establish (see Lemma 1 the modified binomial inversion formula, that is, a n = =n ( ( 1 b if and only if b n = n Motivated by (13, we will introduce and consider the following sequences =n Definition 1 Let a n } (n =, 1,, be a sequence such that ( 1 s 1 a n = =n ( ( 1 a (13 n ( ( 1 a ; s = 1 or s = (14 n We refer to a n } as an invariant sequence (when s = 1 or an inverse invariant sequence (when s = of the second ind 1 15 Naturally arising are the questions of existence, identification, and construction of (inverse invariant sequences of the second ind, as well as the problem of characterizing such sequences and examining their relationship to their counterparts of the first ind Invariant sequences, which are also called self-inverse sequences in [15], have indeed been studied by several authors [6, 8, 14, 15] They are naturally connected to involutory (also nown as involution or self-invertible matrices [9] and to Riordan involutions [5] Involutory matrices find use in numerical methods for differential equations [, 9] They are also useful in cryptography, information theory, and computer security by providing convenient encryption and decryption methods [1] Motivated by Shapiro s open questions [1], Riordan involutions have been intensely investigated as a combinatorial concept [4, 5] In this paper, we investigate invariant sequences by means of the eigenspaces of P D and P T D, where P is the Pascal matrix and D an infinite diagonal matrix with alternating diagonal entries in 1, 1} (see Sections, 3 In fact, P D and P T D are involutory matrices and P D is a Riordan involution Our investigation follows the ideas and connections of invariant sequences to the eigenspaces of P D and P T D developed in Choi et al [6] This will allow us to associate (inverse invariant sequences of the first and second ind, as well as identify and construct such sequences (Section 4 Notation and preliminaries 5 The following notation and conventions are used throughout the manuscript The infinite matrices in this paper have infinite numbers of rows i and columns, with i,, 1,, }

3 E λ (A denotes the eigenspace of a (finite or infinite matrix A corresponding to its eigenvalue λ For a matrix A with columns A ( =, 1,, [ and ] with denoting the vector of zeros in R, A denotes the matrix whose th column is A, where is by convention vacuous For a matrix A, its (possibly infinite row and column index sets are J and K, respectively For J J, K K, let A(J K denote the matrix obtained from A by deleting rows in J and columns in K, and let A[J K ] denote the matrix A(J K, where J = J \J, K = K \K For brevity, write A( K and A(J in place of A( K and A(J, respectively Further, for m, n, 1,, }, we let A m,n = A[, 1,,, m}, 1,,, n}]; A m,m is abbreviated by A m The binomial coefficient ( i choose is denoted by ( i with the convention that it equals when i < or < [ (i P = ] (i, =, 1,, denotes the (infinite Pascal matrix D = diag(1, 1, 1, 1, Infinite real sequences x n } are identified with the infinite dimensional real vector space R consisting of column vectors x = [x, x 1, x, ] T Notice that as a consequence of the binomial inversion formula (11, we have that [ P 1 = DP D = ( 1 i ( i ] (i, =, 1,, Thus (P D 1 = P D and (11 can be converted [6] into a vector equation for x = [a, a 1, a, ] T and y = [b, b 1, b, ] T R, as follows: P Dx = y if and only if P Dy = x (1 Lemma 1 Let P and D be the Pascal matrix and the diagonal matrix defined above, and let x, y R Then P T Dx = y if and only if P T Dy = x, ( and the modified binomial inversion formula (13 holds 45 Proof As P 1 = DP D, we have (P T D 1 = D(P T 1 = P T D As a consequence, ( holds Letting x = [a, a 1, a, ] T and y = [b, b 1, b, ] T implies (13 Let F = [F, F 1, F, ] T and L = [L, L 1, L, ] T denote the vectors in R whose entries are the members of the Fibonacci and Lucas sequences, respectively; that is F =, F 1 = 1, F n = F n 1 + F n (n, L =, L 1 = 1, L n = L n 1 + L n (n 3

4 x The generating functions of F and L are h 1 (x = 1 x x and h (x = x 1 x x, respectively [3, 7] The following fact is nown, however, we include a proof for completeness Lemma P D F = F and P D L = L 5 Proof Let g(x = 1 x 1 x and f(x = 1 x For =, 1,,, the generating function of the th column of P D is g(xf(x [13] Thus the generating function of P DF is [g(x, g(xf(x, g(xf(x, ][F, F 1, F, ] T = g(x(f + F 1 f(x + F f(x + = g(xh 1 (f(x = h 1 (x, which implies that P D F = F The proof of P D L = L is similar That is, 1 and 1 are eigenvalues of P D and consequently of P T D In fact, these are the only eigenvalues of P D and P T D; see [6] The corresponding eigenspaces are infinite dimensional Indeed, if we consider the Pascal-type matrices P and Q constructed via the Pascal matrix P and the matrix [ ] 1 T Q = P +, P then, as shown in [6], the columns of [ T P ] and form bases for E 1 (P D and E 1 (P D, respectively The following observation follows directly from the definitions and properties mentioned above: Observation 3 The entries of x R form an invariant sequence of the first ind if and only if x E 1 (P D; Q 55 an inverse invariant sequence of the first ind if and only if x E 1 (P D; an invariant sequence of the second ind if and only if x E 1 (P T D; an inverse invariant sequence of the second ind if and only if x E 1 (P T D Based on Observation 3, our goal is to study the eigenspaces E λ (P D and E λ (P T D (λ 1, 1} and discover their relationships Our approach entails showing the existence of an infinite invertible matrix N such that N(P T DN 1 = (P T 1 D 1 (P T 1 D 1 (P T 1 D 1 (3 6 and D(N 1 T D (P D DN T D = (P 1 D 1 (P 1 D 1 (P 1 D 1, (4 which are infinite direct sums of copies of P1 T D 1 and P 1 D 1, respectively This result will be applied to characterize E λ (P D and E λ (P T D Extending the wor in [6], we will also show that the columns of P T and Q T ( form bases for E λ (P T D This will indeed allow us to investigate the relationships between invariant sequences of the first and second ind 4

5 3 The Eigenspaces of P T D and P D Let B = [b i ] (i, =, 1,, be the matrix defined by ( 1 i, if i, b i =, if i >, (31 and let J(a denote the infinite Jordan bloc of the form a 1 a 1 O a 1 a 1 O a (3 It readily follows that B 1 = J(1 65 In the next two lemmas, we will construct an infinite matrix N and its inverse M = N 1, which will give rise to similarity transformations of P T D and P D into direct sums as in (3 and (4 Lemma 31 Let m be a positive integer and let N (m = H (m H (m 1 H (1 with H (l = I l B (l = 1,,, m, where B is the matrix defined in (31 Then lim N (m = N = [n i ] m is the infinite matrix defined by n = 1, n = and n i = for i, = 1,,, and n ( 1 i( i 1 + i i =, if 1 i, i 1, if i > 1 Proof Let m be a positive integer and let N (m = H (m H (m 1 H (1 = [n m i ] with nm = 1, n m = and n m i = for i, = 1,, We will prove that n m ( 1 i( + i i =, if i,, if i > by induction on m, where = i 1, if i = 1,,, m; = 1,,, m 1, if i = m + 1, m +, ; = 1,, The claim is clear for m = 1 For m =, by the construction of H (1 and H (, we have n ( 1 i( + i i =, if i,, if i >, 5

6 where = i 1, if i = 1,, 3, 4; = 1,,, 1, if i = 5, 6, ; = 1,,, since for each pair of indices i and with i = 3, 4, and = 1,,, n i = ( l ( l ( 1 l i ( 1 l = ( 1 i l=i l=i ( 1 + i = ( 1 i 1 Let now m 3 Then by the construction of H (m, we have n m i = nm 1 i for all i = 1,,, m and all = 1,, By the induction hypothesis, n m i = ( 1 i( + i, if i,, if i >, where because = n m i = i 1, if i = 1,,, m; = 1,,, m 1, if i = m + 1, m +, ; = 1,,, ( m + l ( 1 l i ( 1 l m l=i ( m + l = ( 1 i m l=i ( m 1 + i = ( 1 i m 1 (33 for each pair of indices i and with i = m 1, m, and = 1,, Thus, by (33, we obtain lim m N (m = N = [n i ] given by n ( 1 i( i 1 + i i =, if i, i 1, if i > for each i, = 1,, Clearly, we have n = 1, n = and n i construction of H (l (l = 1,,, and the proof is complete = for i, = 1,, by the The difference sequence a = [ a, a 1, a, ] T of a sequence a = [a, a 1, a, ] T is defined by a i = a i+1 a i for each i =, 1,, Let a = [ a, a 1, a, ] T ( =, 1,, be the th difference sequence defined inductively by a = ( 1 a, where a = a The infinite matrix a a 1 a a a 1 a a a 1 a 6

7 is called the difference matrix of a It is well nown [3] that for each n =, 1,,, ( ( ( ( n n n n a n = a + a + a + + n a, (34 1 n which is used in the proof of the following lemma Lemma 3 Let M = [m i ] be the matrix with m = 1, m = and m i and ( m i = i, if i, (i, = 1,,, if i > = for i, = 1,,, 7 Then M = N 1, where N is the matrix defined in Lemma 31 Proof Let n i denote the (i, entry of N (i, =, 1,, We would lie to show that n il m l = δ i, l= the Kronecer delta For i = or =, there is nothing to show, so let i 1 and 1 If i = resp i >, then clearly n il m l = l= i ( i 1 n il m li = ( i i 1 = 1 resp l=i n il m l = So it is enough to show that if = i + r with r >, then l= n il m i 1 l = Let = [ and (+r ( consider the sequence z =, +r 1 (, +r (,, ( 1 r 1 +1 ( ] T, ( 1 r,, We can construct the difference matrix A = [a i ] having z as its first column as follows: A = ( +r ( +r 1 ( 1 r ( +r 1 ( +r 1 ( r 1 ( +r ( +r 3 ( 1 r Y ( 1 r ( +1 1 ( 1 r ( ( 1 r 1( 1 ( 1 r 1( 1 ( 1 r( ( 1 r( 1 1 ( 1 r( ( 1 r ( + ( 1 r 1( +1 l=, 7

8 where Y is the (r + 1 (r + 1 matrix given by Y = Thus a = for = + 1,, + r Since + 1 i+r i 1+r+r = + r, by using (34 we obtain where = i 1 n il m l = l= = i+r ( ( + l i n il m l = ( 1 l i t i + r l l=i t ( ( + r s t ( 1 r+s = ( 1 r a t =, s l=i s= i+r and t =, completing the proof Let M (m = U (1 U ( U (m, where U (l = I l J(1 (l = 1,,, m where J(1 is the matrix defined in (3 Then, by Lemmas 31 and 3, we have lim m M (m = M, which is the matrix defined in Lemma 3 In fact we have N =

9 and M = N 1 = We can now state and prove the similarity transformations of P D and P T D claimed in (3 and (4 [ (i Theorem 33 Let P = ] (i, =, 1,, and D = diag(1, 1, 1, 1, Then, 75 (a N(P T DM = (P T 1 D 1 (P T 1 D 1 (P T 1 D 1, (b (DM T D (P D (DN T D = (P 1 D 1 (P 1 D 1 (P 1 D 1, where P 1 resp D 1 are the leading submatrices of P resp D, and N resp M are the matrices defined in Lemmas 31 resp 3 Proof (a For each = 1,,, let H ( = I B and U( = I J(1, where B is the matrix defined in (31 Let m, n be arbitrary positive integers with n = m + 1 First, we will show by induction on m that ( H(m H (m 1 H (1 P T D ( U (1 U ( U (m = Z (P T D, where Z is an n n matrix such that Z = m }} (P T 1 D 1 (P T 1 D 1 (P T 1 D 1 When m = 1, the first row of H (1 P T DU (1 is clearly [1, 1,,, ] Since for a pair of indices i and with i, = 1,,, (( ( ( ( 1 (H (1 P T D i = ( ( 1 = ( 1, i i + 1 i 1 the second row of H (1 P T DU (1 is [, 1, 1, 1, 1, ]U (1 = [, 1,,,, ] For a pair of indices i and with i,, we have ( ( ( 1 (H (1 P T DU (1 i = ( ( 1 = ( 1, i 1 i 1 i and so H (1 P T DU (1 = (P T 1 D 1 (P T D By induction on m, it follows N (m P T DU (m = m }} (P T 1 D 1 (P T 1 D 1 (P T 1 D 1 (P T D, 9

10 where N (m = H (m H (m 1 H (1 and M (m = U (1 U ( U (m Since m was an arbitrary positive integer, NP T DM = (P T 1 D 1 (P T 1 D 1 (P T 1 D 1, where N = lim m N (m and M = lim m M (m (b It follows directly from (a that (DM T DP D(DN T D = (P 1 D 1 (P 1 D 1 (P 1 D 1, completing the proof 8 Let e i denote the ith column of the identity matrix I (i =, 1, Then E λ (P T D and E λ (P D (λ 1, 1} can be characterized via Theorem 33 as follows: Theorem 34 Let N resp M = N 1 be the matrices in Lemma 31 resp Lemma 3 Then the following hold: (a Me, Me, Me 4, } is a basis for E 1 (P T D 85 (b M(e + e 1, M(e + e 3, } is a basis for E 1 (P T D (c DN T D(e + e 1, DN T D(e + e 3 } is a basis for E 1 (P D (d DN T De 1, DN T De 3, DN T De 5, } is a basis for E 1 (P D 9 Proof (a Let B (1 = e, e, e 4, } and x = [x, x 1, x, ] T E 1 (NP T DM Then, by Theorem 33, we have ( [ ] (NP T 1 DM Ix = x =, i=1 which implies that x i = t i, x i+1 = for each i =,, 4, and t i R So B (1 spans E 1 (NP T DM and since e, e, e 4, are linearly independent, B (1 is a basis for E 1 (NP T DM Therefore Me, Me, Me 4, } is a basis for E 1 (P T D (b Let B ( 1 = e + e 1, e + e 3, } and y = [y, y 1, y, ] T E 1 (NP T DM Then, by Theorem 33, we have ( [ ] (NP T 1 DM + Iy = y =, i=1 [ ] [ ] y i 1 which implies that = s i for each i =,, 4, and s i R So B ( 1 spans y i+1 E 1 (NP T DM and since e + e 1, e + e 3, are linearly independent, B ( 1 is a basis for E 1 (NP T DM Therefore M(e + e 1, M(e + e 3, } is a basis for E 1 (P T D Clauses (c and (d can be proven similarly 1

11 Consider now P T = The following is a matrix expression of Theorem 34; the last two clauses appeared in [6] [ ] 1 T [ (i Corollary 35 Let Q = P + and D = diag(( 1, ( 1 1,, where P = P ] (i, =, 1,, Then the following hold: 95 (a The columns of P T form a basis for E 1 (P T D (b The columns of Q T ( form a basis for E 1 (P T D [ ] (c The columns of T P form a basis for E 1 (P D (d The columns of Q form a basis for E 1 (P D Proof (a Let (i, be a pair of indices i and with i, The ith component of Me equals ( }} i when i, and equals, otherwise; Me = [,,, ( (, ( 1,,,,, ] T is the th column of P T, which implies (a (b The ith component of M(e +e +1 equals ( ( ( i + i 1 = +1 ( i + i 1, when i +1, and equals otherwise; }} ( ( ( M(e + e +1 = [,,,,,,,,, ] T }} ( ( ( + [,,,,,,,,, ] T 1 [ ] T is the th column of P T 1 T ( + ( = Q T (, which implies (b Using the fact P T that for a pair of indices i and with i 1, ( 1 (DN T D i = ( 1 i ( 1 i + i 1 ( 1, clauses (c and (d can be proven similarly 11

12 1 4 Invariant sequences of two inds: Relations and examples We begin by some basic examples of (inverse invariant sequences Example 41 It follows from Corollary 35 (c and (d that the Fibonacci sequence F is an inverse invariant sequence of the first ind and the Lucas sequence L is an invariant sequence of the first ind The ith row of P T (i =, 1, is i }} ( i [,,, ( ( i i 1 i,,, 1, ( i,,, ], (41 15 from which we can get that J(F is an invariant sequence of the second ind and J(L is an inverse invariant sequence of the second ind Recall that J( is the infinite Jordan bloc with in the main diagonal, which is the matrix defined in (3 In the following theorem, we provide a general mechanism for transforming invariant sequences into inverse invariant sequences, and vice versa 11 Theorem 4 Let J(λ denote the matrix defined in (3 Then the following hold: (a If x is an invariant sequence of the second ind, then (J(1 + J( T x is an inverse invariant sequence of the second ind (b If x is an inverse invariant sequence of the second ind, then J( 1 J(x is an invariant sequence of the second ind (c If x is an invariant sequence of the first ind, then ( J(J( 1 T x is an inverse invariant 115 sequence of the first ind (d If x is an inverse invariant sequence of the first ind, then (J( 1 + J(x is an invariant sequence of the first ind Proof (a If x is an invariant sequence of the second ind, then by Corollary 35 (a, there exists b R such that P T b = x Let L = [l i ] (i, =, 1,, be the infinite lower triangular matrix defined by ( 1 i+, if i, l i =, if i < Then LP T ( = P T because for a pair of indices i and with i,, the (i, entry of LP T ( equals ( ( ( ( ( 1 i+ + ( 1 i ( 1 i++i = 1 i i 1

13 1 when i, and equals otherwise, which coincides with the (i, entry of P T Thus, [ ] LQ T ( = P T + L T P T From the fact that L 1 = J(1 T and by Corollary 35 (b, it follows that [ ] Q T ( b = J(1 T P T b + T P T (4 b = J(1 T x + [, x T ] T = (J(1 + J( T x is an inverse invariant sequence of the second ind (b Let x be an inverse invariant sequence of the second ind From (4 and J(J(1 T = J(1, it readily follows that J( 1 J(Q T ( = P T and so by Corollary 35 (b, we conclude that J( 1 J(x is an invariant sequence of the second ind 1 T T (c Since Q = P + T, we obtain J( T Q = J( T P 1 T + P T For nonnegative P integers i and with i, =, 1,,, since ( ( i i + i ( i i ( i+ = i+1 i++1, it implies that ([ ] [ ] [ ] T Q T P T P Ω =, where Ω is the infinite (, 1-matrix with 1 s everywhere on and [ ] [ ] below its main diagonal Since J( 1 T = Ω 1, we get (I J( 1 T 1 T T =, which implies that ( J(J( 1 T 15 x is an inverse invariant sequence of the first ind Clause (d easily follows from (c similarly Q P Let τ 1 = 1+ 5 and τ = 1 5 It is well nown that F n = 1 5 τ n τ n and L n = τ n 1 + τ n where F n and L n are the nth terms of F and L, respectively (n =, 1,, [3] Since J(F is an invariant sequence of the second ind by (41, it follows from Theorem 4 (a that J(1 T J(F + J( T J(F = J(1 T J(F + F is an inverse invariant sequence of the second ind In fact, the nth term (n =, 1,, of J(1 T J(F + F is 1 τ1 n+ 1 τ n+ + 1 τ1 n 1 τ n = L n+1, which implies that J(1 T J(F + F = J(L On the other hand, since for a pair of nonnegative integers i and with i, =, 1,,, (J( 1 ( 1 i ( 1 i = i+1, if i,, if i >, 13

14 and the n-th term of J( 1 J( L is ( 1 (1/ +1 (τ1 ++n + τ ++n = 1 τ1 n+1 1 τ n+1 = F n+1, 5 5 = we have J( 1 J( L = J(F, which is an invariant sequence of the second ind by Theorem 4 (b, as stated in Example 41 It directly follows from Theorem 4 (d that (J(+J( 1F, namely L, is an invariant sequence of the first ind Since for a pair of nonnegative integers i and with i, =, 1,,, ( J( 1 ( 1 i = i+1, if i,, if i >, we obtain ( J( 1 T [, L T ] T = F, which is an inverse invariant sequence of the first ind by Theorem 4 (c, because for n =, 1,,, n 1 n 1 (1/ n L = (1/ n ((τ 1 + (τ = 1 τ n 1 τ n = F n 5 5 = = 13 The sequence B = (B, B 1, B, T defined by B = 1 and n ( n+1 = B = (n 1 comprises the Bernoulli numbers and DB is an invariant sequence of the first ind [14], which also follows directly from the fact that P DDB = DB A new inverse invariant sequence of the first ind from the Bernoulli numbers B is provided next See Table 1 for explicit members of these sequences Corollary 43 Let B = (B, B 1, B, T be the Bernoulli numbers Then the sequence K = (K, K 1, K, T defined by n 1 K =, K n = (1/ n ( 1 B (n = 1,, is an inverse invariant sequence of the first ind Proof It follows from Theorem 4 (c that ( J( 1 T J( T DB = ( J( 1 T [, DB T ] T = 135 is an inverse invariant sequence of the first ind Notice now that the sequence K in the statement is indeed equal to ( J( 1 T [, DB T ] T By Theorem 4 (d, we get DB = (J(+J( 1K, since the first component of (J(+J( 1K is clearly ( 1 B, and for each positive integer i with i = 1,,, ith component of (J(+J( 1K is (1/ i ( 1 B + i 1 = i (1/ i+1 ( 1 B = ( 1 i B i = By Corollary 35 and Theorem 4, we can directly get more (inverse invariant sequences of the first and second ind as follows: 14

15 Table 1: (Inverse invariant sequences of the first ind associated with the Bernoulli numbers n B DB 1 1 K Corollary 44 For a positive integer n, let Φ n = S ln S ln 1 S l3 S l S l1, where P T, if l i = 1, S li = J(1 T + J( T, if l i = 1, and let Φ n = S ln Sln 1 S l3 Sl Sl1 where J( S 1 J(, if l i = 1, li = Q T (, if l i = 1 for i = 1,,, n Then for x R, we have the following: 14 (a If n is odd (even and l i = ( 1 i+1 for i = 1,,, n, then Φ n x is an invariant (inverse invariant sequence of the second ind (b If n is odd (even and l i = ( 1 i for i = 1,,, n, then Φ n x is an inverse invariant (invariant sequence of the second ind Example 45 It follows from (41 and Corollary 44 (a that for x = [x, x 1, x, ] T R, Φ x = y is an inverse invariant sequence of the second ind, where y = [y, y 1, y, ] T with y i = i t= i (( t + i 1 t ( t + 1 x t (i =, 1,, i t For example, let x = [,,,,,,, 1,,, ] T This results to the inverse invariant sequence of the second ind y = [y, y 1, y, ] T, where ( ( 7 8 y i = + (i =, 1,, i 1 7 i 7 That is, y = = y 6 =, y 16 = y 17 = =, and by direct calculation, one can compute the nonzero components y 7, y 8,, y 15 ; eg, y 1 = ( ( = 77 Corollary 46 For a positive integer n, let Ψ n = T ln T ln 1 T l3 T l T l1, where T li = Q, if l i = 1, ( J( 1 T J( T, if l i = 1 15

16 and let Ψ n = T ln Tln 1 T l3 Tl Tl1, where J( + J( 1, if l i = 1, [ ] T li = T, if l i = 1 P 145 for i = 1,,, n Then for x R, we have the following: (a If n is odd (even and l i = ( 1 i+1 for i = 1,,, n, then Ψ n x is an invariant (inverse invariant sequence of the first ind (b If n is odd (even and l i = ( 1 i for i = 1,,, n, then Ψ n x is an inverse invariant (invariant sequence of the first ind [ ] T Example 47 For i = 1,,, the ith row of is [ ( ( i 1, i ( i 1 1,,,,, ] It follows i 1 P that for x = [x, x 1, x, ] T R, Ψ x = y is an invariant sequence of the first ind, where by Corollary 46 (b, y = (y, y 1, y, T satisfies i (( i 1 t y i = + t 1 t= ( i t t x t For example, let x = [,,,,,,, 1,,, ] T This results to the invariant sequence of the first ind y = [y, y 1, y, ] T, where ( ( i 1 7 i 7 y i = + (i = 1,, Thus y = y 1 = = y 13 = and one can find by direct calculation, the nonzero components y 14, y 15, ; eg, y 14 = ( ( = In the next two theorems, we obtain direct relationships among (inverse invariant sequences of the first ind and (inverse invariant sequences of the second ind 155 Theorem 48 Let x and y be, respectively, either (i an invariant sequence of the first ind and an inverse invariant sequence of the second ind, or (ii an inverse invariant sequence of the first ind and an invariant sequence of the second ind Then x T Dy = 16 Proof From P Dx = λx and P T Dy = λy, we have x T DP T Dy = λx T Dy, which implies that λx T Dy = for λ 1, 1} The following theorem is a useful tool for proving the final theorem Theorem 49 For every positive integer n, 16

17 (a the columns of (P +D n are invariant sequences of the first ind, and the columns of (P D n are inverse invariant sequences of the first ind 165 (b the columns of (P T + D n are invariant sequences of the second ind, and the columns of (P T D n are inverse invariant sequences of the second ind Proof Let n be a positive integer Then and P D(P + D n = P D(P + D(P + D n 1 = (P + D n P T D(P T D n = P T D(P T D(P T D n 1 = (P T D n, which respectively imply that each column of (P +D n is an invariant sequence of the first ind, and each column of (P T D n is an inverse invariant sequence of the second ind The other assertions of the theorem follow similarly 17 The following theorem is a form of converse of Theorem 48 Theorem 41 Let x i R (i =, 1,,, y R \ }, and let X = [x, x 1, x, ] Then for each i =, 1,,, 175 (a if x i T Dy = and X = P + D (X = P D, then y is an inverse invariant (invariant sequence of the second ind, and x i is an invariant (inverse invariant sequence of the first ind (b if x i T Dy = and X = P T + D (X = P T D, then y is an inverse invariant (invariant sequence of the first ind and x i is an inverse invariant (invariant sequence of the second ind 18 Proof If x T i Dy = and X = P + D, then since X T Dy = (P T + DDy =, we have P T Dy = y and by Theorem 49, P Dx i = x i Thus y is an inverse invariant sequence of the second ind and x i is an invariant sequence of the first ind This proves the first case of part (a For the case of X = P D and part (b, the results can be shown similarly References 185 [1] B Acharya, SK Panigrahy, SK Patra, G Panda, Image encryption using advanced Hill Ciper Algorithm, Int J Recent Trends in Engineering, 1(1 ( [] AL Andrew, Eigenvectors of certain matrices, Linear Algebra Appl, 7 ( [3] RA Brualdi, Introductory combinatorics, Fourth ed, Prentice-Hall, Upper Saddle River, 4 17

18 19 [4] NT Cameron, A Nwanta, On some(pseudo involutions in the Riordan group, J Integer Seq 8, (5 Article 537 [5] G-S Cheon, H Kim, Simple proofs of open problems about the structure of involutions in the Riordan group, Linear Algebra Appl, 48 ( [6] G-S Choi, S-G Hwang, I-P Kim BL Shader, (±1-Invariant sequences and truncated Fibonacci sequences, Linear Algebra Appl, 395 ( [7] L Comtet, Advanced Combinatorics, Reidel, Dordrecht, 1974 [8] R Donaghey, Binomial self-inverse sequences and tangent coefficient, J Combin Theory Ser A, 1 (1976, [9] RA Horn, CR Johnson, Matrix Analysis Second ed, Cambridge University Press, New Yor, 13 [1] H Prodinger, Some information about the binomial transform, Fibonacci Quarterly, 3 ( [11] J Riordan, Combinatorial Identities, Wiley, New Yor, 1968 [1] LW Shapiro, Some open questions about random wals, involutions, limiting distributions and generating functions, Adv Appl Math, 7 ( [13] LW Shapiro, S Getu, W-J Woan, LC Woodson, The Riordan group, Discrete Appl Math, 34 ( [14] ZH Sun, Invariant sequences under binomial transformation, Fibonacci Quarterly, 39 ( [15] Y Wang, Self-inverse sequences related to a binomial inverse pair, Fibonacci Quarterly, 43 ( [16] K Woci, Binomial transform and Dole sequences, J Integer Seq, 18 (15 Article

Applications of Riordan matrix functions to Bernoulli and Euler polynomials

Applications of Riordan matrix functions to Bernoulli and Euler polynomials Illinois Wesleyan University From the SelectedWors of Tian-Xiao He 06 Applications of Riordan matrix functions to Bernoulli and Euler polynomials Tian-Xiao He Available at: https://worsbepresscom/tian_xiao_he/78/

More information

Discrete Applied Mathematics

Discrete Applied Mathematics Discrete Applied Mathematics 157 (2009 1696 1701 Contents lists available at ScienceDirect Discrete Applied Mathematics journal homepage: www.elsevier.com/locate/dam Riordan group involutions and the -sequence

More information

MAPPING AND PRESERVER PROPERTIES OF THE PRINCIPAL PIVOT TRANSFORM

MAPPING AND PRESERVER PROPERTIES OF THE PRINCIPAL PIVOT TRANSFORM MAPPING AND PRESERVER PROPERTIES OF THE PRINCIPAL PIVOT TRANSFORM OLGA SLYUSAREVA AND MICHAEL TSATSOMEROS Abstract. The principal pivot transform (PPT) is a transformation of a matrix A tantamount to exchanging

More information

Yi Wang Department of Applied Mathematics, Dalian University of Technology, Dalian , China (Submitted June 2002)

Yi Wang Department of Applied Mathematics, Dalian University of Technology, Dalian , China (Submitted June 2002) SELF-INVERSE SEQUENCES RELATED TO A BINOMIAL INVERSE PAIR Yi Wang Department of Applied Mathematics, Dalian University of Technology, Dalian 116024, China (Submitted June 2002) 1 INTRODUCTION Pairs of

More information

Bounds for Levinger s function of nonnegative almost skew-symmetric matrices

Bounds for Levinger s function of nonnegative almost skew-symmetric matrices Linear Algebra and its Applications 416 006 759 77 www.elsevier.com/locate/laa Bounds for Levinger s function of nonnegative almost skew-symmetric matrices Panayiotis J. Psarrakos a, Michael J. Tsatsomeros

More information

Geometric Mapping Properties of Semipositive Matrices

Geometric Mapping Properties of Semipositive Matrices Geometric Mapping Properties of Semipositive Matrices M. J. Tsatsomeros Mathematics Department Washington State University Pullman, WA 99164 (tsat@wsu.edu) July 14, 2015 Abstract Semipositive matrices

More information

The generalized order-k Fibonacci Pell sequence by matrix methods

The generalized order-k Fibonacci Pell sequence by matrix methods Journal of Computational and Applied Mathematics 09 (007) 33 45 wwwelseviercom/locate/cam The generalized order- Fibonacci Pell sequence by matrix methods Emrah Kilic Mathematics Department, TOBB University

More information

On identities with multinomial coefficients for Fibonacci-Narayana sequence

On identities with multinomial coefficients for Fibonacci-Narayana sequence Annales Mathematicae et Informaticae 49 08 pp 75 84 doi: 009/ami080900 http://amiuni-eszterhazyhu On identities with multinomial coefficients for Fibonacci-Narayana sequence Taras Goy Vasyl Stefany Precarpathian

More information

arxiv: v1 [math.co] 11 Jul 2015

arxiv: v1 [math.co] 11 Jul 2015 arxiv:50703055v [mathco] Jul 05 Duals of Bernoulli Numbers and Polynomials and Euler Number and Polynomials Tian-Xiao He and Jinze Zheng Department of Mathematics Illinois Wesleyan University Bloomington,

More information

Math 113 Homework 5. Bowei Liu, Chao Li. Fall 2013

Math 113 Homework 5. Bowei Liu, Chao Li. Fall 2013 Math 113 Homework 5 Bowei Liu, Chao Li Fall 2013 This homework is due Thursday November 7th at the start of class. Remember to write clearly, and justify your solutions. Please make sure to put your name

More information

Central Groupoids, Central Digraphs, and Zero-One Matrices A Satisfying A 2 = J

Central Groupoids, Central Digraphs, and Zero-One Matrices A Satisfying A 2 = J Central Groupoids, Central Digraphs, and Zero-One Matrices A Satisfying A 2 = J Frank Curtis, John Drew, Chi-Kwong Li, and Daniel Pragel September 25, 2003 Abstract We study central groupoids, central

More information

Intrinsic products and factorizations of matrices

Intrinsic products and factorizations of matrices Available online at www.sciencedirect.com Linear Algebra and its Applications 428 (2008) 5 3 www.elsevier.com/locate/laa Intrinsic products and factorizations of matrices Miroslav Fiedler Academy of Sciences

More information

The k-fibonacci matrix and the Pascal matrix

The k-fibonacci matrix and the Pascal matrix Cent Eur J Math 9(6 0 403-40 DOI: 0478/s533-0-0089-9 Central European Journal of Mathematics The -Fibonacci matrix and the Pascal matrix Research Article Sergio Falcon Department of Mathematics and Institute

More information

On the adjacency matrix of a block graph

On the adjacency matrix of a block graph On the adjacency matrix of a block graph R. B. Bapat Stat-Math Unit Indian Statistical Institute, Delhi 7-SJSS Marg, New Delhi 110 016, India. email: rbb@isid.ac.in Souvik Roy Economics and Planning Unit

More information

ON QUADRAPELL NUMBERS AND QUADRAPELL POLYNOMIALS

ON QUADRAPELL NUMBERS AND QUADRAPELL POLYNOMIALS Hacettepe Journal of Mathematics and Statistics Volume 8() (009), 65 75 ON QUADRAPELL NUMBERS AND QUADRAPELL POLYNOMIALS Dursun Tascı Received 09:0 :009 : Accepted 04 :05 :009 Abstract In this paper we

More information

Math 3191 Applied Linear Algebra

Math 3191 Applied Linear Algebra Math 9 Applied Linear Algebra Lecture 9: Diagonalization Stephen Billups University of Colorado at Denver Math 9Applied Linear Algebra p./9 Section. Diagonalization The goal here is to develop a useful

More information

Definition (T -invariant subspace) Example. Example

Definition (T -invariant subspace) Example. Example Eigenvalues, Eigenvectors, Similarity, and Diagonalization We now turn our attention to linear transformations of the form T : V V. To better understand the effect of T on the vector space V, we begin

More information

RESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices

RESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices Linear and Multilinear Algebra Vol. 00, No. 00, Month 200x, 1 15 RESEARCH ARTICLE An extension of the polytope of doubly stochastic matrices Richard A. Brualdi a and Geir Dahl b a Department of Mathematics,

More information

On sum of powers of the Laplacian and signless Laplacian eigenvalues of graphs

On sum of powers of the Laplacian and signless Laplacian eigenvalues of graphs On sum of powers of the Laplacian and signless Laplacian eigenvalues of graphs Saieed Akbari 1,2 Ebrahim Ghorbani 1,2 Jacobus H. Koolen 3,4 Mohammad Reza Oboudi 1,2 1 Department of Mathematical Sciences

More information

The Perron eigenspace of nonnegative almost skew-symmetric matrices and Levinger s transformation

The Perron eigenspace of nonnegative almost skew-symmetric matrices and Levinger s transformation Linear Algebra and its Applications 360 (003) 43 57 www.elsevier.com/locate/laa The Perron eigenspace of nonnegative almost skew-symmetric matrices and Levinger s transformation Panayiotis J. Psarrakos

More information

SOME RESULTS ON q-analogue OF THE BERNOULLI, EULER AND FIBONACCI MATRICES

SOME RESULTS ON q-analogue OF THE BERNOULLI, EULER AND FIBONACCI MATRICES SOME RESULTS ON -ANALOGUE OF THE BERNOULLI, EULER AND FIBONACCI MATRICES GERALDINE M. INFANTE, JOSÉ L. RAMÍREZ and ADEM ŞAHİN Communicated by Alexandru Zaharescu In this article, we study -analogues of

More information

EVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS

EVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS EVALUATION OF A FAMILY OF BINOMIAL DETERMINANTS CHARLES HELOU AND JAMES A SELLERS Abstract Motivated by a recent work about finite sequences where the n-th term is bounded by n, we evaluate some classes

More information

Remark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.

Remark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero. Sec 5 Eigenvectors and Eigenvalues In this chapter, vector means column vector Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called

More information

On the closures of orbits of fourth order matrix pencils

On the closures of orbits of fourth order matrix pencils On the closures of orbits of fourth order matrix pencils Dmitri D. Pervouchine Abstract In this work we state a simple criterion for nilpotentness of a square n n matrix pencil with respect to the action

More information

Representation of doubly infinite matrices as non-commutative Laurent series

Representation of doubly infinite matrices as non-commutative Laurent series Spec. Matrices 217; 5:25 257 Research Article Open Access María Ivonne Arenas-Herrera and Luis Verde-Star* Representation of doubly infinite matrices as non-commutative Laurent series https://doi.org/1.1515/spma-217-18

More information

ON THE MATRIX EQUATION XA AX = X P

ON THE MATRIX EQUATION XA AX = X P ON THE MATRIX EQUATION XA AX = X P DIETRICH BURDE Abstract We study the matrix equation XA AX = X p in M n (K) for 1 < p < n It is shown that every matrix solution X is nilpotent and that the generalized

More information

Some results on the existence of t-all-or-nothing transforms over arbitrary alphabets

Some results on the existence of t-all-or-nothing transforms over arbitrary alphabets Some results on the existence of t-all-or-nothing transforms over arbitrary alphabets Navid Nasr Esfahani, Ian Goldberg and Douglas R. Stinson David R. Cheriton School of Computer Science University of

More information

Linear transformations preserving the strong q-log-convexity of polynomials

Linear transformations preserving the strong q-log-convexity of polynomials Linear transformations preserving the strong q-log-convexity of polynomials Bao-Xuan Zhu School of Mathematics and Statistics Jiangsu Normal University Xuzhou, PR China bxzhu@jsnueducn Hua Sun College

More information

Matrix functions that preserve the strong Perron- Frobenius property

Matrix functions that preserve the strong Perron- Frobenius property Electronic Journal of Linear Algebra Volume 30 Volume 30 (2015) Article 18 2015 Matrix functions that preserve the strong Perron- Frobenius property Pietro Paparella University of Washington, pietrop@uw.edu

More information

Animals and 2-Motzkin Paths

Animals and 2-Motzkin Paths 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 8 (2005), Article 0556 Animals and 2-Motzkin Paths Wen-jin Woan 1 Department of Mathematics Howard University Washington, DC 20059 USA wwoan@howardedu

More information

The Jordan Canonical Form

The Jordan Canonical Form The Jordan Canonical Form The Jordan canonical form describes the structure of an arbitrary linear transformation on a finite-dimensional vector space over an algebraically closed field. Here we develop

More information

1. General Vector Spaces

1. General Vector Spaces 1.1. Vector space axioms. 1. General Vector Spaces Definition 1.1. Let V be a nonempty set of objects on which the operations of addition and scalar multiplication are defined. By addition we mean a rule

More information

Homework sheet 4: EIGENVALUES AND EIGENVECTORS. DIAGONALIZATION (with solutions) Year ? Why or why not? 6 9

Homework sheet 4: EIGENVALUES AND EIGENVECTORS. DIAGONALIZATION (with solutions) Year ? Why or why not? 6 9 Bachelor in Statistics and Business Universidad Carlos III de Madrid Mathematical Methods II María Barbero Liñán Homework sheet 4: EIGENVALUES AND EIGENVECTORS DIAGONALIZATION (with solutions) Year - Is

More information

Investigating Geometric and Exponential Polynomials with Euler-Seidel Matrices

Investigating Geometric and Exponential Polynomials with Euler-Seidel Matrices 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 14 (2011), Article 11.4.6 Investigating Geometric and Exponential Polynomials with Euler-Seidel Matrices Ayhan Dil and Veli Kurt Department of Mathematics

More information

Frame Diagonalization of Matrices

Frame Diagonalization of Matrices Frame Diagonalization of Matrices Fumiko Futamura Mathematics and Computer Science Department Southwestern University 00 E University Ave Georgetown, Texas 78626 U.S.A. Phone: + (52) 863-98 Fax: + (52)

More information

Determinants of Partition Matrices

Determinants of Partition Matrices journal of number theory 56, 283297 (1996) article no. 0018 Determinants of Partition Matrices Georg Martin Reinhart Wellesley College Communicated by A. Hildebrand Received February 14, 1994; revised

More information

Takao Komatsu School of Mathematics and Statistics, Wuhan University

Takao Komatsu School of Mathematics and Statistics, Wuhan University Degenerate Bernoulli polynomials and poly-cauchy polynomials Takao Komatsu School of Mathematics and Statistics, Wuhan University 1 Introduction Carlitz [6, 7] defined the degenerate Bernoulli polynomials

More information

RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA

RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA Discussiones Mathematicae General Algebra and Applications 23 (2003 ) 125 137 RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA Seok-Zun Song and Kyung-Tae Kang Department of Mathematics,

More information

ACI-matrices all of whose completions have the same rank

ACI-matrices all of whose completions have the same rank ACI-matrices all of whose completions have the same rank Zejun Huang, Xingzhi Zhan Department of Mathematics East China Normal University Shanghai 200241, China Abstract We characterize the ACI-matrices

More information

(a) II and III (b) I (c) I and III (d) I and II and III (e) None are true.

(a) II and III (b) I (c) I and III (d) I and II and III (e) None are true. 1 Which of the following statements is always true? I The null space of an m n matrix is a subspace of R m II If the set B = {v 1,, v n } spans a vector space V and dimv = n, then B is a basis for V III

More information

Maximizing the numerical radii of matrices by permuting their entries

Maximizing the numerical radii of matrices by permuting their entries Maximizing the numerical radii of matrices by permuting their entries Wai-Shun Cheung and Chi-Kwong Li Dedicated to Professor Pei Yuan Wu. Abstract Let A be an n n complex matrix such that every row and

More information

A q-analogue OF THE GENERALIZED FACTORIAL NUMBERS

A q-analogue OF THE GENERALIZED FACTORIAL NUMBERS J. Korean Math. Soc. 47 (2010), No. 3, pp. 645 657 DOI 10.4134/JKMS.2010.47.3.645 A q-analogue OF THE GENERALIZED FACTORIAL NUMBERS Seok-Zun Song, Gi-Sang Cheon, Young-Bae Jun, and LeRoy B. Beasley Abstract.

More information

Hessenberg Pairs of Linear Transformations

Hessenberg Pairs of Linear Transformations Hessenberg Pairs of Linear Transformations Ali Godjali November 21, 2008 arxiv:0812.0019v1 [math.ra] 28 Nov 2008 Abstract Let K denote a field and V denote a nonzero finite-dimensional vector space over

More information

Two applications of the theory of primary matrix functions

Two applications of the theory of primary matrix functions Linear Algebra and its Applications 361 (2003) 99 106 wwwelseviercom/locate/laa Two applications of the theory of primary matrix functions Roger A Horn, Gregory G Piepmeyer Department of Mathematics, University

More information

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector

More information

A Simple Proof of Fiedler's Conjecture Concerning Orthogonal Matrices

A Simple Proof of Fiedler's Conjecture Concerning Orthogonal Matrices University of Wyoming Wyoming Scholars Repository Mathematics Faculty Publications Mathematics 9-1-1997 A Simple Proof of Fiedler's Conjecture Concerning Orthogonal Matrices Bryan L. Shader University

More information

A LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS

A LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS Applicable Analysis and Discrete Mathematics, 1 (27, 371 385. Available electronically at http://pefmath.etf.bg.ac.yu A LINEAR BINOMIAL RECURRENCE AND THE BELL NUMBERS AND POLYNOMIALS H. W. Gould, Jocelyn

More information

2 b 3 b 4. c c 2 c 3 c 4

2 b 3 b 4. c c 2 c 3 c 4 OHSx XM511 Linear Algebra: Multiple Choice Questions for Chapter 4 a a 2 a 3 a 4 b b 1. What is the determinant of 2 b 3 b 4 c c 2 c 3 c 4? d d 2 d 3 d 4 (a) abcd (b) abcd(a b)(b c)(c d)(d a) (c) abcd(a

More information

Improved Upper Bounds for the Laplacian Spectral Radius of a Graph

Improved Upper Bounds for the Laplacian Spectral Radius of a Graph Improved Upper Bounds for the Laplacian Spectral Radius of a Graph Tianfei Wang 1 1 School of Mathematics and Information Science Leshan Normal University, Leshan 614004, P.R. China 1 wangtf818@sina.com

More information

Catalan triangle numbers and binomial coefficients

Catalan triangle numbers and binomial coefficients Contemporary Mathematics Volume 73, 208 https://doiorg/0090/conm/73/435 Catalan triangle numbers and binomial coefficients Kyu-Hwan Lee and Se-jin Oh Abstract The binomial coefficients and Catalan triangle

More information

CATALAN TRIANGLE NUMBERS AND BINOMIAL COEFFICIENTS arxiv:6006685v2 [mathco] 7 Oct 207 KYU-HWAN LEE AND SE-JIN OH Abstract The binomial coefficients and Catalan triangle numbers appear as weight multiplicities

More information

Kernels of Directed Graph Laplacians. J. S. Caughman and J.J.P. Veerman

Kernels of Directed Graph Laplacians. J. S. Caughman and J.J.P. Veerman Kernels of Directed Graph Laplacians J. S. Caughman and J.J.P. Veerman Department of Mathematics and Statistics Portland State University PO Box 751, Portland, OR 97207. caughman@pdx.edu, veerman@pdx.edu

More information

arxiv: v1 [math.ra] 8 Aug 2016

arxiv: v1 [math.ra] 8 Aug 2016 Each n-by-n matrix with n > 1 is a sum of 5 coninvolutory matrices arxiv:1608.02503v1 math.ra 8 Aug 2016 Ma. Nerissa M. Abara b, Dennis I. Merino 1,, Vyacheslav I. Rabanovich c, Vladimir V. Sergeichuk

More information

Spectrally arbitrary star sign patterns

Spectrally arbitrary star sign patterns Linear Algebra and its Applications 400 (2005) 99 119 wwwelseviercom/locate/laa Spectrally arbitrary star sign patterns G MacGillivray, RM Tifenbach, P van den Driessche Department of Mathematics and Statistics,

More information

Upper triangular matrices and Billiard Arrays

Upper triangular matrices and Billiard Arrays Linear Algebra and its Applications 493 (2016) 508 536 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Upper triangular matrices and Billiard Arrays

More information

AN ALGEBRAIC INTERPRETATION OF THE q-binomial COEFFICIENTS. Michael Braun

AN ALGEBRAIC INTERPRETATION OF THE q-binomial COEFFICIENTS. Michael Braun International Electronic Journal of Algebra Volume 6 (2009) 23-30 AN ALGEBRAIC INTERPRETATION OF THE -BINOMIAL COEFFICIENTS Michael Braun Received: 28 February 2008; Revised: 2 March 2009 Communicated

More information

A hyperfactorial divisibility Darij Grinberg *long version*

A hyperfactorial divisibility Darij Grinberg *long version* A hyperfactorial divisibility Darij Grinberg *long version* Let us define a function H : N N by n H n k! for every n N Our goal is to prove the following theorem: Theorem 0 MacMahon We have H b + c H c

More information

Research Article The Adjacency Matrix of One Type of Directed Graph and the Jacobsthal Numbers and Their Determinantal Representation

Research Article The Adjacency Matrix of One Type of Directed Graph and the Jacobsthal Numbers and Their Determinantal Representation Applied Mathematics Volume 20, Article ID 423163, 14 pages doi:101155/20/423163 Research Article The Adjacency Matrix of One Type of Directed Graph and the Jacobsthal Numbers and Their Determinantal Representation

More information

Eigenvectors and Reconstruction

Eigenvectors and Reconstruction Eigenvectors and Reconstruction Hongyu He Department of Mathematics Louisiana State University, Baton Rouge, USA hongyu@mathlsuedu Submitted: Jul 6, 2006; Accepted: Jun 14, 2007; Published: Jul 5, 2007

More information

Lecture 15, 16: Diagonalization

Lecture 15, 16: Diagonalization Lecture 15, 16: Diagonalization Motivation: Eigenvalues and Eigenvectors are easy to compute for diagonal matrices. Hence, we would like (if possible) to convert matrix A into a diagonal matrix. Suppose

More information

Generalized Akiyama-Tanigawa Algorithm for Hypersums of Powers of Integers

Generalized Akiyama-Tanigawa Algorithm for Hypersums of Powers of Integers 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 16 (2013, Article 1332 Generalized Aiyama-Tanigawa Algorithm for Hypersums of Powers of Integers José Luis Cereceda Distrito Telefónica, Edificio Este

More information

Counting Matrices Over a Finite Field With All Eigenvalues in the Field

Counting Matrices Over a Finite Field With All Eigenvalues in the Field Counting Matrices Over a Finite Field With All Eigenvalues in the Field Lisa Kaylor David Offner Department of Mathematics and Computer Science Westminster College, Pennsylvania, USA kaylorlm@wclive.westminster.edu

More information

Some Formulas for the Principal Matrix pth Root

Some Formulas for the Principal Matrix pth Root Int. J. Contemp. Math. Sciences Vol. 9 014 no. 3 141-15 HIKARI Ltd www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.4110 Some Formulas for the Principal Matrix pth Root R. Ben Taher Y. El Khatabi and

More information

A Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)!

A Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)! 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 8 (2005), Article 05.2.3 A Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)! Ira M. Gessel 1 and Guoce Xin Department of Mathematics Brandeis

More information

Definition: An n x n matrix, "A", is said to be diagonalizable if there exists a nonsingular matrix "X" and a diagonal matrix "D" such that X 1 A X

Definition: An n x n matrix, A, is said to be diagonalizable if there exists a nonsingular matrix X and a diagonal matrix D such that X 1 A X DIGONLIZTION Definition: n n x n matrix, "", is said to be diagonalizable if there exists a nonsingular matrix "X" and a diagonal matrix "D" such that X X D. Theorem: n n x n matrix, "", is diagonalizable

More information

0.1 Rational Canonical Forms

0.1 Rational Canonical Forms We have already seen that it is useful and simpler to study linear systems using matrices. But matrices are themselves cumbersome, as they are stuffed with many entries, and it turns out that it s best

More information

Wavelets and Linear Algebra

Wavelets and Linear Algebra Wavelets and Linear Algebra () (05) 49-54 Wavelets and Linear Algebra http://wala.vru.ac.ir Vali-e-Asr University of Rafsanjan Schur multiplier norm of product of matrices M. Khosravia,, A. Sheikhhosseinia

More information

ELA ON A SCHUR COMPLEMENT INEQUALITY FOR THE HADAMARD PRODUCT OF CERTAIN TOTALLY NONNEGATIVE MATRICES

ELA ON A SCHUR COMPLEMENT INEQUALITY FOR THE HADAMARD PRODUCT OF CERTAIN TOTALLY NONNEGATIVE MATRICES ON A SCHUR COMPLEMENT INEQUALITY FOR THE HADAMARD PRODUCT OF CERTAIN TOTALLY NONNEGATIVE MATRICES ZHONGPENG YANG AND XIAOXIA FENG Abstract. Under the entrywise dominance partial ordering, T.L. Markham

More information

The Jordan Normal Form and its Applications

The Jordan Normal Form and its Applications The and its Applications Jeremy IMPACT Brigham Young University A square matrix A is a linear operator on {R, C} n. A is diagonalizable if and only if it has n linearly independent eigenvectors. What happens

More information

m We can similarly replace any pair of complex conjugate eigenvalues with 2 2 real blocks. = R

m We can similarly replace any pair of complex conjugate eigenvalues with 2 2 real blocks. = R 1 RODICA D. COSTIN Suppose that some eigenvalues are not real. Then the matrices S and are not real either, and the diagonalization of M must be done in C n. Suppose that we want to work in R n only. Recall

More information

MATH JORDAN FORM

MATH JORDAN FORM MATH 53 JORDAN FORM Let A,, A k be square matrices of size n,, n k, respectively with entries in a field F We define the matrix A A k of size n = n + + n k as the block matrix A 0 0 0 0 A 0 0 0 0 A k It

More information

Lemma 8: Suppose the N by N matrix A has the following block upper triangular form:

Lemma 8: Suppose the N by N matrix A has the following block upper triangular form: 17 4 Determinants and the Inverse of a Square Matrix In this section, we are going to use our knowledge of determinants and their properties to derive an explicit formula for the inverse of a square matrix

More information

arxiv: v1 [math.nt] 9 May 2017

arxiv: v1 [math.nt] 9 May 2017 The spectral norm of a Horadam circulant matrix Jorma K Merikoski a, Pentti Haukkanen a, Mika Mattila b, Timo Tossavainen c, arxiv:170503494v1 [mathnt] 9 May 2017 a Faculty of Natural Sciences, FI-33014

More information

JORDAN NORMAL FORM. Contents Introduction 1 Jordan Normal Form 1 Conclusion 5 References 5

JORDAN NORMAL FORM. Contents Introduction 1 Jordan Normal Form 1 Conclusion 5 References 5 JORDAN NORMAL FORM KATAYUN KAMDIN Abstract. This paper outlines a proof of the Jordan Normal Form Theorem. First we show that a complex, finite dimensional vector space can be decomposed into a direct

More information

e j = Ad(f i ) 1 2a ij/a ii

e j = Ad(f i ) 1 2a ij/a ii A characterization of generalized Kac-Moody algebras. J. Algebra 174, 1073-1079 (1995). Richard E. Borcherds, D.P.M.M.S., 16 Mill Lane, Cambridge CB2 1SB, England. Generalized Kac-Moody algebras can be

More information

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v ) Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O

More information

An Analytic Approach to the Problem of Matroid Representibility: Summer REU 2015

An Analytic Approach to the Problem of Matroid Representibility: Summer REU 2015 An Analytic Approach to the Problem of Matroid Representibility: Summer REU 2015 D. Capodilupo 1, S. Freedman 1, M. Hua 1, and J. Sun 1 1 Department of Mathematics, University of Michigan Abstract A central

More information

Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems

Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Linear Algebra and its Applications 430 (2009) 532 543 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: wwwelseviercom/locate/laa Computing tight upper bounds

More information

QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY

QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY D.D. Olesky 1 Department of Computer Science University of Victoria Victoria, B.C. V8W 3P6 Michael Tsatsomeros Department of Mathematics

More information

arxiv: v1 [cs.sc] 29 Jul 2009

arxiv: v1 [cs.sc] 29 Jul 2009 An Explicit Construction of Gauss-Jordan Elimination Matrix arxiv:09075038v1 [cssc] 29 Jul 2009 Yi Li Laboratory of Computer Reasoning and Trustworthy Computing, University of Electronic Science and Technology

More information

Two finite forms of Watson s quintuple product identity and matrix inversion

Two finite forms of Watson s quintuple product identity and matrix inversion Two finite forms of Watson s uintuple product identity and matrix inversion X. Ma Department of Mathematics SuZhou University, SuZhou 215006, P.R.China Submitted: Jan 24, 2006; Accepted: May 27, 2006;

More information

Computationally, diagonal matrices are the easiest to work with. With this idea in mind, we introduce similarity:

Computationally, diagonal matrices are the easiest to work with. With this idea in mind, we introduce similarity: Diagonalization We have seen that diagonal and triangular matrices are much easier to work with than are most matrices For example, determinants and eigenvalues are easy to compute, and multiplication

More information

0.1 Tangent Spaces and Lagrange Multipliers

0.1 Tangent Spaces and Lagrange Multipliers 01 TANGENT SPACES AND LAGRANGE MULTIPLIERS 1 01 Tangent Spaces and Lagrange Multipliers If a differentiable function G = (G 1,, G k ) : E n+k E k then the surface S defined by S = { x G( x) = v} is called

More information

A property concerning the Hadamard powers of inverse M-matrices

A property concerning the Hadamard powers of inverse M-matrices Linear Algebra and its Applications 381 (2004 53 60 www.elsevier.com/locate/laa A property concerning the Hadamard powers of inverse M-matrices Shencan Chen Department of Mathematics, Fuzhou University,

More information

The Fibonacci Identities of Orthogonality

The Fibonacci Identities of Orthogonality The Fibonacci Identities of Orthogonality Kyle Hawins, Ursula Hebert-Johnson and Ben Mathes January 14, 015 Abstract In even dimensions, the orthogonal projection onto the two dimensional space of second

More information

Some bounds for the spectral radius of the Hadamard product of matrices

Some bounds for the spectral radius of the Hadamard product of matrices Some bounds for the spectral radius of the Hadamard product of matrices Guang-Hui Cheng, Xiao-Yu Cheng, Ting-Zhu Huang, Tin-Yau Tam. June 1, 2004 Abstract Some bounds for the spectral radius of the Hadamard

More information

On the spectra of striped sign patterns

On the spectra of striped sign patterns On the spectra of striped sign patterns J J McDonald D D Olesky 2 M J Tsatsomeros P van den Driessche 3 June 7, 22 Abstract Sign patterns consisting of some positive and some negative columns, with at

More information

A note on restricted Lie supertriple systems

A note on restricted Lie supertriple systems Advances in Theoretical and Applied Mathematics. ISSN 0973-4554 Volume 10, Number 1 (2015), pp. 1 13 c Research India Publications http://www.ripublication.com/atam.htm A note on restricted Lie supertriple

More information

Chapter 2 The Group U(1) and its Representations

Chapter 2 The Group U(1) and its Representations Chapter 2 The Group U(1) and its Representations The simplest example of a Lie group is the group of rotations of the plane, with elements parametrized by a single number, the angle of rotation θ. It is

More information

WAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME)

WAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME) WAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME ITAI BENJAMINI, GADY KOZMA, LÁSZLÓ LOVÁSZ, DAN ROMIK, AND GÁBOR TARDOS Abstract. We observe returns of a simple random wal on a finite graph to a fixed node,

More information

Recall : Eigenvalues and Eigenvectors

Recall : Eigenvalues and Eigenvectors Recall : Eigenvalues and Eigenvectors Let A be an n n matrix. If a nonzero vector x in R n satisfies Ax λx for a scalar λ, then : The scalar λ is called an eigenvalue of A. The vector x is called an eigenvector

More information

L(C G (x) 0 ) c g (x). Proof. Recall C G (x) = {g G xgx 1 = g} and c g (x) = {X g Ad xx = X}. In general, it is obvious that

L(C G (x) 0 ) c g (x). Proof. Recall C G (x) = {g G xgx 1 = g} and c g (x) = {X g Ad xx = X}. In general, it is obvious that ALGEBRAIC GROUPS 61 5. Root systems and semisimple Lie algebras 5.1. Characteristic 0 theory. Assume in this subsection that chark = 0. Let me recall a couple of definitions made earlier: G is called reductive

More information

Markov Chains, Stochastic Processes, and Matrix Decompositions

Markov Chains, Stochastic Processes, and Matrix Decompositions Markov Chains, Stochastic Processes, and Matrix Decompositions 5 May 2014 Outline 1 Markov Chains Outline 1 Markov Chains 2 Introduction Perron-Frobenius Matrix Decompositions and Markov Chains Spectral

More information

Key words. Strongly eventually nonnegative matrix, eventually nonnegative matrix, eventually r-cyclic matrix, Perron-Frobenius.

Key words. Strongly eventually nonnegative matrix, eventually nonnegative matrix, eventually r-cyclic matrix, Perron-Frobenius. May 7, DETERMINING WHETHER A MATRIX IS STRONGLY EVENTUALLY NONNEGATIVE LESLIE HOGBEN 3 5 6 7 8 9 Abstract. A matrix A can be tested to determine whether it is eventually positive by examination of its

More information

Remark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.

Remark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero. Sec 6 Eigenvalues and Eigenvectors Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called an eigenvalue of A if there is a nontrivial

More information

MULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng

MULTISECTION METHOD AND FURTHER FORMULAE FOR π. De-Yin Zheng Indian J. pure appl. Math., 39(: 37-55, April 008 c Printed in India. MULTISECTION METHOD AND FURTHER FORMULAE FOR π De-Yin Zheng Department of Mathematics, Hangzhou Normal University, Hangzhou 30036,

More information

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution

More information

Elementary Row Operations on Matrices

Elementary Row Operations on Matrices King Saud University September 17, 018 Table of contents 1 Definition A real matrix is a rectangular array whose entries are real numbers. These numbers are organized on rows and columns. An m n matrix

More information

and let s calculate the image of some vectors under the transformation T.

and let s calculate the image of some vectors under the transformation T. Chapter 5 Eigenvalues and Eigenvectors 5. Eigenvalues and Eigenvectors Let T : R n R n be a linear transformation. Then T can be represented by a matrix (the standard matrix), and we can write T ( v) =

More information