NEGATIVELY SUBSCRIPTED FIBONACCI AND LUCAS NUMBERS AND THEIR COMPLEX FACTORIZATIONS E. KILIC 1 AND D. TASCI2
Abstract. In this paper, we …nd families of (0; 1; 1) tridiagonal matrices whose determinants and permanents equal to the negatively subscripted Fibonacci and Lucas numbers. Also we give complex factorizations of these numbers by the …rst and second kinds of Chebyshev polynomials.
1. Introduction The well-known Fibonacci sequence, fFn g ; is de…ned by the recurrence relation, for n 2 Fn+1 = Fn + Fn 1 (1.1) where F1 = F2 = 1: The Lucas Sequence, fLn g ; is de…ned by the recurrence relation, for n 2 Ln+1 = Ln + Ln
(1.2)
1
where L1 = 1; L2 = 3: Rules (1.1) and (1.2) can be used to extend the sequence backward, respectively, thus F L
= F1 = L1
1 1
F0 ; L0;
F L
= F0 F L 2 = L0
2
1 1; : : : ;
and so on. Clearly F L
n n
= F = L
n+2 n+2
F L
n+1
= ( 1) Fn ; n n+1 = ( 1) Ln :
n+1
(1.3) (1.4)
In [9] and [5], the authors give complex factorizations of the Fibonacci numbers by considering the roots of Fibonacci polynomials as follows Fn =
nQ1
1
2i cos
k=1
k n
;
n
2:
(1.5)
2000 Mathematics Subject Classi…cation. 11B37, 15A15, 15A36. Key words and phrases. Fibonacci and Lucas numbers, tridiagonal matrices, determinant, permanent, factorization. 1
E. KILIC 1 AND D. TASCI 2
2
In [10] and [11], the authors establish the following forms: Fn
= in
Ln
=
n cos 1 sin cos 1
1 sin
2in cos n cos
1
i 2 i 2
i 2
; n ; n
1
(1.6) 1:
In [3], the authors prove (1.5) by considering how to the Fibonacci numbers can be connected to Chebyshev polynomials by determinants of a sequence of matrices, and then show that a connection between the Lucas numbers and Chebyshev polynomials by using a slightly di¤erent sequence of matrices as follows ! n Q k 21 ; n 1: Ln = 1 2i cos n k=1
There are many connections between permanents or determinants of tridiagonal matrices and the Fibonacci and Lucas numbers. For example, Minc [12] de…nes a n n super diagonal matrix F (n; k) for n > k 2; and shows that the permanent of F (n; k) equals the generalized orderk Fibonacci numbers. In [14], the author gives the same result of Minc by the same matrix F (n; k) and using di¤erent a computing method of permanent, contraction. In particular, when k = 2; the matrix F (n; k) is reduced to the tridiagonal toeplitz matrix 3 2 1 1 0 7 6 7 6 1 1 ... 7 F (n; 2) = 6 6 7 .. .. 4 . 1 5 . 0 1 1
and perF (n; 2) = Fn+1 : In [15], Lehmer proves a very general result on permanents of tridiagonal matrices whose main diagonal and super-diagonal elements are ones and whose subdiagonal entries are somewhat arbitrary. Also in [16] and [17], the authors de…ne the n n tridiagonal matrix Mn and show that the determinant of M (n) is the Fibonacci number F2n+2 . In [2] and [3], the authors de…ne the n n tridiagonal matrix H (n) and show that the determinant of H (n) is the Fibonacci number Fn : In a similar family of matrices, the (1; 1) element of H (n) is replaced with a 3, thus the determinants, [18], now generate the Lucas sequence Ln : Recently, in [7], the authors …nd families of square matrices such that (i) each matrix is the adjacency matrix of a bipartite graph; and (ii) the permanent of the matrix is a sum of consecutive Fibonacci or Lucas numbers.” Also, in [8], the authors de…ne two tridiagonal matrices and then give the
NEGATIVELY SUBSCRIPTED FIBONACCI AND LUCAS NUM BERS
3
relationships between the permanents and determinants of these matrices, and the terms of second order linear recurrences. In this paper, we consider negatively subscripted Fibonacci and Lucas numbers and …nd associated families of tridiagonal matrices whose determinants or permanents equal to these numbers. Then we give the complex factorizations of these numbers by Chebyshev polynomial. The permanent of an n-square matrix A = (aij ) is de…ned by perA =
n X Y
ai
(i)
2Sn i=1
where the summation extends over all permutations of the symmetric group Sn : Also one can …nd more applications of permanents in [13]. A matrix is said to be a ( 1; 0; 1)-matrix if each of its entries are 1; 0 or 1: Let A = [aij ] be an m n real matrix row vectors 1 ; 2 ; : : : ; m : We say A is contractible on column (resp. row:) k if column (resp. row:) k contains exactly two nonzero entries. Suppose A is contractible on column k with aik 6= 0 6= ajk and i 6= j: Then the (m 1) (n 1) matrix Aij:k obtained from A by replacing row i with ajk i + aik j and deleting row j and column k is called the contraction of A on column k relative to rows i and j: If A is contractible on row k with aki 6= 0 6= akj and i 6= j; then the h iT matrix Ak:ij = ATij:k is called the contraction of A on row k relative to columns i and j: Every contraction used in this paper will be on the …rst column using the …rst and second rows. We say that A can be contracted to a matrix B if either B = A or exist matrices A0 ; A1 ; : : : At (t 1) such that A0 = A; At = B; and Ar is a contraction of Ar 1 for r = 1; 2; : : : ; t: Let we consider the following result (see [1]): Let A be a nonnegative integral matrix of order n > 1 and let B be a contraction of A. Then perA = perB:
(1.7)
2. Negatively Subscripted Fibonacci and Lucas numbers In this section, we de…ne families of tridiagonal matrices and then show that the determinants and permanents of these matrices equal to the negatively subscripted Fibonacci and Lucas numbers. We start with negatively subscripted Fibonacci numbers. Now we de…ne a n n tridiagonal toeplitz (0; 1; 1) matrix An = [aij ] with aii = 1 for
E. KILIC 1 AND D. TASCI 2
4
1
i
n; ai;i+1 = ai+1;i = 1 for 1 i 2 1 1 6 6 1 1 An = 6 6 . .. 4 0
Then we give following Theorem.
n
1 and 0 otherwise. That is, 3 0 7 .. 7 . 7: (2.1) 7 .. 5 . 1 1 1
Theorem 1. Let the matrix An have the form (2.1). Then, for n perAn = F where F
n
1
(n+1)
is the nth negatively subscripted Fibonacci number.
Proof. If n = 1; then perA1 = per [ 1] = F 2 = If n = 2; then 1 1 A2 = 1 1 and hence perA2 = F 3 = 2: Let Apn be pth contraction of An ; 1 p n An ; the matrix An can be contracted on column 2 2 1 6 1 1 1 6 6 . 1 1 1 .. An = 6 6 6 .. .. 4 . . 1 Since the matrix 2; 2 3 6 1 6 6 2 An = 6 6 6 4
1:
2: From the de…nition of 1 so that 3 7 7 7 7: 7 7 1 5 1
A1n can be contracted on column 1 and F 2 1 1
3
1 1 .. .
.. ..
.
. 1
2
F 4 7 6 1 7 6 7 6 7=6 7 6 7 6 1 5 4 1
Continuing this process, we obtain 2 F (r+2) F 6 1 6 6 r An = 6 6 6 4
F
1 1
1 ..
1 .
..
..
1 .. .
..
.
. 1
.
. 1
7 7 7 7 7 7 1 5 1
3; F
3
3
1
3
1 ..
=
3
(r+1)
1
4
7 7 7 7: 7 7 1 5 1
=
NEGATIVELY SUBSCRIPTED FIBONACCI AND LUCAS NUM BERS
for 3
r
n
4: Hence, Ann
3
2
=4
which, by contraction of Ann Ann
2
5
F
=
F
(n 2)
F
F
(n 1)
1 0 4
(n 2)
1 1
on column 1, gives F
(n 1)
1
(n 1)
=
1
3 0 1 5 1 F
n
F
1
(n 1)
1
:
By the Eq. (1.7) and the de…nition of the negatively subscripted Fibonacci numbers, we obtain perAn = perAnn
2
=F
F
(n 1)
n
=F
(n+1) :
So the proof is complete. Second, we de…ne a n n tridiagonal (0; 1; 1) matrix Bn = [bij ] with bii = 1 for 2 i n; bi;i+1 = bi+1;i = 1 for 1 i n 1; b11 = 21 and 0 otherwise. That is, 3 2 1 1 0 2 7 6 . 7 6 1 1 .. 7: 6 (2.2) Bn = 6 7 . . .. .. 4 1 5 0 1 1 Now we give following Theorem.
Theorem 2. Let the matrix Bn has the form (2.2). Then perBn = where L
n
L n 2
is the nth negatively subscripted Lucas number.
Proof. If n = 1; then perB1 = per
1 =L 2
1 =2
If n = 2; then B2 = and hence perB2 = L
2 =2
= 3=2:
1 2
1
1 1
=
1=2:
E. KILIC 1 AND D. TASCI 2
6
Let Bnp be pth contraction of Bn ; 1 p n 2: From the de…nition of Bn ; the matrix Bn can be contracted on column 1 so that 2 3 3 1 0 2 2 6 1 7 1 1 6 7 6 7 . 1 . 6 7: . 1 1 Bn = 6 7 6 7 .. .. 4 . . 1 5 0 1 1
Since the matrix Bn1 can be contracted on column 1 and L 3; 3 2 L 3 2 4 L 2 3 6 6 6 2 Bn = 6 6 6 4
2
1
2
1 1
7 6 7 6 7 6 7=6 7 6 7 6 1 5 4 1
1
1 .. .
..
.
..
. 1
Continuing this process, we obtain 2 L (r+1)
for 3
r
n
2
6 6 6 r Bn = 6 6 6 4
1
Bnn
3
L
=4
2
=
L
(n 3)
L 1
1
(n 2)
1
..
1 .. .
L
.
..
. 1
(n 3) =2
on column 1, gives (n 2) =2
1
4; L
1 .. .
=
3
7 7 7 7: 7 7 1 5 1
1 ..
2
.
..
. 1
7 7 7 7 7 7 1 5 1
1 1
=2 L
1
=
3
1
(n 2) =2
4
2
1
1 0
which, by contraction of Bnn Bnn
L r 2
1
4: Hence, 2
2
3
=
3 0 1 5 1 L
(n 1) =2
1
L
(n 2) =2
1
By the Eq. (1.7) and the de…nition of the negatively subscripted Lucas numbers, we obtain perBn = perBnn
2
= L
(n 2)
L
(n 1)
=2 = L
n =2:
So the proof is complete. A matrix A is called convertible if there is an n n (1; 1) matrix H such that perA = det (A H) ; where A H denotes the Hadamard product of A and H: Such a matrix H is called a converter of A:
:
NEGATIVELY SUBSCRIPTED FIBONACCI AND LUCAS NUM BERS
7
Let S be a (1; 1) matrix of order n, de…ned by 2
6 6 6 S=6 6 4
1 1 1 .. .
1 ::: 1 ::: 1 ::: .. .
1
1 :::
3 1 1 1 1 7 7 1 1 7 7: .. .. 7 . . 5 1 1
Then we have that F (n+1) = det (An S) and L n =2 = det (Bn S) where F n and L n are the nth negatively subscripted Fibonacci and Lucas number, respectively. Let we denote the matrices An S and Bn S by Cn and Dn ; respectively. Thus 3 2 1 1 0 7 6 .. 7 6 1 . 1 7 Cn = 6 7 6 .. .. 4 . . 1 5 0 1 1 and
2
6 6 Dn = 6 6 4
1 2
1
1
1 .. .
0 ..
7 7 7 7 1 5 1
.
..
. 1
0
3
Also it is clear that the value of following determinant is independent of x : (see p.105, [19]) a 1 x
x
0
a .. .
..
.
..
.
1 x
0
: x a
Using the above result and considering the following matrices 2
1 6 p 6 1 C^n = 6 6 4 0
p
1
1 .. .
0 ..
.
.. . p 1
3
7 7 7 p 7 1 5 1
E. KILIC 1 AND D. TASCI 2
8
and
2
p
1 2
6 p 6 1 ^ Dn = 6 6 4
1
1 .. .
0
we can write that
det C^n ^n det D
3
0 ..
.
.. . p 1
7 7 7 p 7; 1 5 1
=
det Cn = perAn = F
(n+1) ;
=
det Dn = perBn = L
n =2:
Furthermore, from [13]; we have that let A be a tridiagonal matrix, and let Ap = (aij ) be de…ned by ast = iast if s 6= t and ass = ass ; for all s and t i= 1 : Then we have per (A) = det A : Also let we de…ne the following matrices; 2 1 i 6 . 6 i 1 .. 6 Cn = 6 .. .. 4 . . 0 i
and
2
6 6 Dn = 6 6 4
1 2
i
i
1 .. .
0
0
..
. . i
(2.3)
3
(2.4)
7 7 7 7 i 5 1 0
..
3
7 7 7: 7 i 5 1
Thus we have following Corollaries without proof. Corollary 1. Let the n Then, for n 1
n tridiagonal toeplitz matrix Cn as in (2.3). det Cn = F
Corollary 2. Let the n for n 1
(n+1) :
n tridiagonal matrix Dn be as in (2.4). Then, det Dn = L
n =2:
NEGATIVELY SUBSCRIPTED FIBONACCI AND LUCAS NUM BERS
9
3. Complex Factorization of the negatively subscripted Fibonacci numbers In [3], the authors consider the relationships between the certain tridiagonal determinants, and, the usual Fibonacci and Lucas numbers. Then using the eigenvalues of these tridiagonal matrices, the authors give the complex factorizations of the usual Fibonacci and Lucas numbers. Following the method of [3], we …nd the eigenvalues of the two tridiagonal matrices whose determinants associated with the negatively subscripted Fibonacci and Lucas numbers. Therefore, we give the complex factorizations of the negatively subscripted Fibonacci and Lucas numbers. There are variety of ways of computing matrix determinants (see [4] and [6] for more details). In addition to the method of cofactor expansion, the determinant of a matrix can be found by taking the product of its eigenvalues. Therefore, we will compute the spectrum of Cn to …nd an alternative representation of det Cn : Now we de…ne another n n tridiagonal toeplitz matrix Vn = [vij ] with vii = 0 for 1 i n and vi;i 1 = vi 1;i = 1 for 2 i n and 0 otherwise. Clearly 3 2 0 1 0 7 6 7 6 1 0 ... 7: (3.1) Vn = 6 7 6 .. .. 4 . . 1 5 0 1 0 So it is clear that Cn =
Theorem 3. Let F Then, for n 1 F
n
I + iVn : Then we give following Theorem.
be the nth negatively subscripted Fibonacci number.
(n+1)
=
n Q
1
j n+1
2i cos
j=1
:
Proof. Let j ; j = 1; 2; : : : ; n; be the eigenvalues of Vn with respect to eigenvectors xi : Then, for all j Cn xj = ( I + iVn ) xj = Therefore, for n 1
j
=
1+i
j;
Ixj + iVn xj =
xj + i
j xj
=( 1+i
j ) xj :
(3.2) j = 1; 2; : : : ; n; are the eigenvalues of Cn : Hence,
det Cn =
n Q
( 1+i
j) :
(3.3)
j=1
To compute the j ’s, we recall that each j is a zero of the characteristic polynomial pn ( ) = jVn Ij : Since Vn I is a tridiagonal toeplitz
E. KILIC 1 AND D. TASCI 2
10
matrix, i.e.,
0
B B 1 I=B B @
Vn
1
1 .. ..
.
.
..
. 1
C C C; C 1 A
(3.4)
we can establish a recursive formula for the characteristic polynomials Vn : p1 ( )
=
; 2
p2 ( ) = pn ( ) =
1; pn
1
(3.5) ( )
pn
2
( ):
This family of characteristic polynomials can be transformed into another family fUn (x) ; n 1g by taking 2x : U1 (x)
=
2x;
U2 (x) Un (x)
= =
4x2 1; 2xUn 1 (x)
(3.6) Un
2
(x) :
The family fUn (x) ; n 1g is the set of Chebyshev polynomials of second kind. It is a well-known fact (see [10]) that de…ning x cos allows the Chebyshev polynomials of the second kind to be written as: Un (x) =
sin [(n + 1) ] : sin
(3.7)
From (3.7), we can see that the roots of Un (x) = 0 are given by k = k k n+1 ; k = 1; 2; : : : ; n; or xk = cos k = cos n+1 ; k = 1; 2; : : : ; n: Applying the transformation 2x; we have the eigenvalues of Vn : k
=
2 cos
k n+1
; k = 1; 2; : : : ; n:
(3.8)
Considering Corollary 1, the Eqs. (3.3) and (3.8), we obtain F
(n+1)
= det Cn =
n Q
1
j n+1
2i cos
j=1
which is desired. Theorem 4. Let F Then, for n 1 F
n
be nth negatively subscripted Fibonacci number.
(n+1)
= in
sin (n + 1) cos 1 sin cos 1 2i
i 2
:
NEGATIVELY SUBSCRIPTED FIBONACCI AND LUCAS NUM BERS
11
Proof. From (3.4), we can think of Chebyshev polynomials of the second kind as being generated by determinants of successive matrices of the form 2 3 2x 1 6 7 6 1 2x . . . 7 6 7; (3.9) Kn (x) = 6 7 . . .. .. 1 5 4 1 2x where Kn (x) is n
i 2
n: If we denote that Cn = iKn i 2
det Cn = in det Kn
; then we obtain:
i 2
= in Un
:
(3.10)
Combining the result of Corollary 1, the Eqs. (3.7) and (3.10) yields, for n 1 sin (n + 1) cos 1 2i : F (n+1) = in sin cos 1 2i So the proof is complete. Theorem 5. Let L for n 1 L
n
n
be nth negatively subscripted Lucas number. Then,
=
n Q
1
2i cos
k=1
k n
1 2
!
:
Proof. From Corollary 2, we have that 2 det Dn = L n : We will not compute the spectrum of Dn directly. Instead, we will note that the following: (det I + e1 eT1 = 2) det Dn =
1 det 2
I + e1 eT1 Dn ;
(3.11)
where ej is the jth column of the identity matrix. Thus we can write that the right-side of (3.11) as follows 1 det 2
I + e1 eT1 Dn =
1 det 2
I + i Vn + e1 eT2
(3.12)
where the matrix Vn is given by (3.1). Let j ; j = 1; 2; : : : ; n; be the eigenvalues of Vn + e1 eT2 with respect to eigenvectors yj : Then, for all j I + i Vn + e1 eT2
yj =
Iyj +i Vn + e1 eT2 yj =
yj +i j yj =
1+i
Thus 1 det 2
I + i Vn + e1 eT2
=
n 1 Q 2 k=1
1+i
j
:
(3.13)
j
yj :
E. KILIC 1 AND D. TASCI 2
12
To compute the j ’s, we recall that all is a zero of the characteristic polynomial tn ( ) = det Vn + e1 eT2 I : Since det I 21 e1 eT1 = 12 ; we can alternately write the characteristic polynomial as tn ( ) = 2 det
T 1 2 e1 e1
I
Vn + e1 eT2
I
:
(3.14)
Since tn ( ) is twice the determinant of a tridiagonal matrix, that is, 2 3 1 2 6 7 .. 6 1 7 . 6 7; tn ( ) = 2 det 6 (3.15) 7 .. .. 4 . . 1 5 1
one can derive a recursive formula for t1 ( ) 2 t2 ( ) 2 tn ( ) 2
= =
tn ( ) 2
:
2 2
1
2
=
tn
1
( )
tn
2
( ):
This family of polynomials can be transformed into another family fTn (x) ; n by taking 2x : T1 (x)
= x;
T2 (x) = 2x2 1; Tn (x) = 2xTn 1 (x)
Tn
2
(x) :
The family fTn (x) ; n 1g is the set of Chebyshev polynomials of …rst kind. In [10], Rivlin presents that de…ning x cos allows the Chebyshev polynomials of the …rst kind to be written as Tn (x) = cos n :
(3.16)
From the Eq. (3.16), one can see that the roots of Tn (x) = 0 are given by k 21 or xk = cos k = cos n Applying the transformation 2x and (3.14) are also roots of det Vn + e1 eT2 I T of Vn + e1 e2 : k
=
k 12 for k = 1; 2; : : : ; n: n considering the roots of the = 0; we have the eigenvalues
k 12 for k = 1; 2; : : : ; n: n From Corollary 2, the Eqs. (3.17) and (3.13), we obtain ! n Q k 21 L n= 1 2i cos : n k=1 k
=
2 cos
So the proof is complete.
(3.17)
1g
NEGATIVELY SUBSCRIPTED FIBONACCI AND LUCAS NUM BERS
Theorem 6. Let L for n 1
n
L
13
be nth negatively subscripted Lucas number. Then, n
= 2in cos n cos
1
i 2
:
Proof. From (3.15), we think of Chebyshev polynomials of the …rst kind as being generated by determinants of successive matrices of the form 2 3 x 1 6 7 6 1 2x . . . 7 6 7 Gn (x) = 6 : 7 . . .. .. 1 5 4 1 2x n n
We note that det Dn = iGn
i 2
; thus
det Dn = in det Gn
i 2
= in Tn
i 2
:
(3.18)
From Corollary 2, the Eqs. (3.16) and (3.18), we obtain L
n
= 2in cos n cos
1
i 2
:
So the proof is complete.
Acknowledgements The authors would like to thank to the anonymous referee for helpful suggestions. References [1] R. A. Brualdi & P. M. Gibson. "Convex polyhedra of Doubly Stochastic matrices I: Applications of the Permanent Function." J. Combin. Theory A 22 (1977): 194230. [2] N.D. Cahill, J.R. D’Ericco, D.A. Narayan and J.Y. Narayan. "Fibonacci determinants." College Math. J. 3.3 (2002): 221-225. [3] N. Cahill, J. R. D’Ericco and J. P. Spence. "Complex factorizations of the Fibonacci and Lucas numbers." Fib. Quart. 41.1 (2003): 13-19. [4] G. Golub & C. Van Loan. Matrix computations. 3rd ed., pp. 50-51, 310. Baltimore: Johns Hopkins University Press, 1996. [5] V. Hoggatt, Jr., & C. Long. "Divisibility properties of generalized Fibonacci polynomials." Fib. Quart. 12.2 (1974): 113-120. [6] R.A. Horn and C.R. Johnson. Matrix analysis. Cambridge University Press, Cambridge-New York, 1990. [7] E. Kilic and D. Tasci. "On families of bipartite graphs associated with sums of Fibonacci and Lucas numbers." Ars Combin. (to appear). [8] E. Kilic and D. Tasci. "On the second order linear recurrences by tridiagonal matrices" Ars Combin. (to appear). [9] W. A. Webb & E. A. Parberry. "Divisibility properties of Fibonacci polynomials." Fib. Quart. 7.5 (1969): 457-463.
E. KILIC 1 AND D. TASCI 2
14
[10] T. Rivlin. Chebyshev Polynomials : From Approximation Theory to Algebra and Number Theory. 2nd ed. New York: Wiley & Sons, 1990. [11] J. Morgado. "Note on the Chebyshev polynomials and applications to the Fibonacci numbers." Port. Math. 52 (1995): 363-78. [12] H. Minc. "Permanents of (0; 1) Circulants." Canad. Math. Bull., 7.2 (1964), 253263. [13] H. Minc. Permanents, Encyclopedia of Mathematics and its Applications. AddisonWesley, New York, 1978. [14] G. -Y. Lee and S. -G. Lee. "A note on generalized Fibonacci numbers." Fib. Quart. 33 (1995): 273-278. [15] D. Lehmer. "Fibonacci and related sequences in periodic tridiagonal matrices." Fib. Quart. 13 (1975): 150-158. [16] G. Strang. Introduction to Linear algebra, 3rd Edition, Wellesley MA, WellesleyCambridge, 2003. [17] G. Strang and K. Borre. Linear algebra. Geodesy and GPS. Wellesley MA, WellesleyCambridge, 1997, pp. 555-557. [18] P.F. Byrd. "Problem B-12: A Lucas determinant." Fib. Quart. 1.4 (1963): 78. [19] F. Zhang. Matrix theory. Springer-Verlag, New York-Berlin-Heidelberg, 1999. 1 TOBB
Economics and Technology University Mathematics Department 06560 Ankara Turkey E-mail address :
[email protected] 2 Gazi
University, Department of Mathematics, 06500 Ankara Turkey