Michel Waldschmidt. Linear recurrence sequences, exponential polynomials and Diophantine approximation. Institut de Mathématiques de Jussieu Paris VI

Size: px
Start display at page:

Download "Michel Waldschmidt. Linear recurrence sequences, exponential polynomials and Diophantine approximation. Institut de Mathématiques de Jussieu Paris VI"

Transcription

1 8th ICSMEDC December 6, 2015 The Eighth International Conference on Science and Mathematics Education in Developing Countries Yangon University, Yangon, The Republic of the Union of Myanmar. Linear recurrence sequences, exponential polynomials and Diophantine approximation Michel Waldschmidt Institut de Mathématiques de Jussieu Paris VI

2 Abstract Linear recurrence sequences are ubiquitous. They occur in biology, economics, computer science (analysis of algorithms), digital signal processing. We give a survey of this subject, including connections with linear combinations of powers and with exponential polynomials, with an emphasis on arithmetic questions. This lecture will include new results, arising from a joint work with Claude Levesque, involving families of Diophantine equations, with explicit examples related to some units of L. Bernstein and H. Hasse.

3 Leonardo Pisano (Fibonacci) Fibonacci sequence (F n ) n 0 0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, is defined by F 0 =0,F 1 =1, F n+2 = F n+1 + F n (n 0). Leonardo Pisano (Fibonacci) ( )

4 Fibonacci rabbits Fibonacci considers the growth of a rabbit population. Anewlybornpairofrabbits, one male, one female, are put in a field. Rabbits are able to mate at the age of one month so that at the end of its second month a female can produce another pair of rabbits ; rabbits never die and amatingpairalwaysproduces one new pair (one male, one female) every month from the second month on. The puzzle that Fibonacci posed was : how many pairs will there be in one year? Answer : F 12 =144.

5 Fibonacci squares

6 Geometric construction of the Fibonacci sequence

7 Fibonacci numbers in nature Ammonite (Nautilus shape)

8 Reflections of a ray of light Consider three parallel sheets of glass and a ray of light which crosses the first sheet. Each time it touches one of the sheets, it can cross it or reflect on it. Denote by p n the number of di erent paths with the ray going out of the system after n reflections. p 0 =1, p 1 =2, p 2 =3, p 3 =5. In general, p n = F n+2.

9 Reflections of a ray of light Consider three parallel sheets of glass and a ray of light which crosses the first sheet. Each time it touches one of the sheets, it can cross it or reflect on it. Denote by p n the number of di erent paths with the ray going out of the system after n reflections. p 0 =1, p 1 =2, p 2 =3, p 3 =5. In general, p n = F n+2.

10 Levels of energy of an electron of an atom of hydrogen An atom of hydrogen can have three levels of energy, 0 at the ground level when it does not move, 1 or 2. Ateachstep, alternatively, it gains or losses 1 or 2 levels of energy, without going below 0 nor above 2. Let`n. bethenumberofdi erent possible histories of this electron are there after n steps. We have `0 =1(initial state level 0) `1 =2:state1 or 2, histories 01 or 02. `2 =3:histories010, 021 or 020. `3 =5:histories0101, 0102, 0212, 0201 or In general, `n = F n+2.

11 Rhythmic patterns The Fibonacci sequence appears in Indian mathematics, in connection with Sanskrit prosody. Several Indian scholars, Pingala (200 BC), Virahanka (c. 700 AD), Gopāla (c. 1135), and the Jain scholar Hemachandra (c. 1150) studied rhythmic patterns that are formed from one-beat notes (or short syllables, ti in Morse Alphabet) : and two-beat notes (or long syllables, ta ta in Morse) :. 1 beat, 1 pattern : 2 beats, 2 patterns : and 3 beats, 3 patterns :, and 4 beats, 5 patterns :,,,, n beats, F n+1 patterns

12 Rhythmic patterns The Fibonacci sequence appears in Indian mathematics, in connection with Sanskrit prosody. Several Indian scholars, Pingala (200 BC), Virahanka (c. 700 AD), Gopāla (c. 1135), and the Jain scholar Hemachandra (c. 1150) studied rhythmic patterns that are formed from one-beat notes (or short syllables, ti in Morse Alphabet) : and two-beat notes (or long syllables, ta ta in Morse) :. 1 beat, 1 pattern : 2 beats, 2 patterns : and 3 beats, 3 patterns :, and 4 beats, 5 patterns :,,,, n beats, F n+1 patterns

13 Rhythmic patterns The Fibonacci sequence appears in Indian mathematics, in connection with Sanskrit prosody. Several Indian scholars, Pingala (200 BC), Virahanka (c. 700 AD), Gopāla (c. 1135), and the Jain scholar Hemachandra (c. 1150) studied rhythmic patterns that are formed from one-beat notes (or short syllables, ti in Morse Alphabet) : and two-beat notes (or long syllables, ta ta in Morse) :. 1 beat, 1 pattern : 2 beats, 2 patterns : and 3 beats, 3 patterns :, and 4 beats, 5 patterns :,,,, n beats, F n+1 patterns

14 Rhythmic patterns The Fibonacci sequence appears in Indian mathematics, in connection with Sanskrit prosody. Several Indian scholars, Pingala (200 BC), Virahanka (c. 700 AD), Gopāla (c. 1135), and the Jain scholar Hemachandra (c. 1150) studied rhythmic patterns that are formed from one-beat notes (or short syllables, ti in Morse Alphabet) : and two-beat notes (or long syllables, ta ta in Morse) :. 1 beat, 1 pattern : 2 beats, 2 patterns : and 3 beats, 3 patterns :, and 4 beats, 5 patterns :,,,, n beats, F n+1 patterns

15 Rhythmic patterns The Fibonacci sequence appears in Indian mathematics, in connection with Sanskrit prosody. Several Indian scholars, Pingala (200 BC), Virahanka (c. 700 AD), Gopāla (c. 1135), and the Jain scholar Hemachandra (c. 1150) studied rhythmic patterns that are formed from one-beat notes (or short syllables, ti in Morse Alphabet) : and two-beat notes (or long syllables, ta ta in Morse) :. 1 beat, 1 pattern : 2 beats, 2 patterns : and 3 beats, 3 patterns :, and 4 beats, 5 patterns :,,,, n beats, F n+1 patterns

16 Rhythmic patterns The Fibonacci sequence appears in Indian mathematics, in connection with Sanskrit prosody. Several Indian scholars, Pingala (200 BC), Virahanka (c. 700 AD), Gopāla (c. 1135), and the Jain scholar Hemachandra (c. 1150) studied rhythmic patterns that are formed from one-beat notes (or short syllables, ti in Morse Alphabet) : and two-beat notes (or long syllables, ta ta in Morse) :. 1 beat, 1 pattern : 2 beats, 2 patterns : and 3 beats, 3 patterns :, and 4 beats, 5 patterns :,,,, n beats, F n+1 patterns

17 Fibonacci sequence and the Golden ratio For n 0, thefibonacci number F n is the nearest integer to 1 p 5 n, where is the Golden Ratio : = 1+p 5 2 which satisfies = lim n!1 F n+1 F n = =1+ 1

18 Binet s formula For n 0, F n = n ( ) n p 5 Jacques Philippe Marie Binet (1843) = (1 + p 5) n (1 2 np 5 p 5) n, = 1+p 5, 2 1 = 1 p 5, 2 X 2 X 1=(X )(X + 1 ).

19 The so called Binet Formula Formula of A. De Moivre (1718, 1730), Daniel Bernoulli (1726), L. Euler (1728, 1765), J.P.M. Binet (1843) : for n 0, F n = p 1 1+ p! n p p! n 5. Abraham de Moivre ( ) Daniel Bernoulli ( ) Leonhard Euler ( ) Jacques P.M. Binet ( )

20 Generating series AsingleseriesencodesalltheFibonacci sequence : X F n X n = X + X 2 +2X 3 +3X 4 +5X F n X n + n 0 Fact : this series is the Taylor expansion of a rational fraction : X F n X n X = 1 X X 2 n 0 Proof :theproduct (X + X 2 +2X 3 +3X 4 +5X 5 +8X 6 + )(1 X X 2 ) is a telescoping series X + X 2 +2X 3 +3X 4 +5X 5 +8X 6 + = X. X 2 X 3 2X 4 3X 5 5X 6 X 3 X 4 2X 5 3X 6

21 Generating series AsingleseriesencodesalltheFibonacci sequence : X F n X n = X + X 2 +2X 3 +3X 4 +5X F n X n + n 0 Fact : this series is the Taylor expansion of a rational fraction : X F n X n X = 1 X X 2 n 0 Proof :theproduct (X + X 2 +2X 3 +3X 4 +5X 5 +8X 6 + )(1 X X 2 ) is a telescoping series X + X 2 +2X 3 +3X 4 +5X 5 +8X 6 + = X. X 2 X 3 2X 4 3X 5 5X 6 X 3 X 4 2X 5 3X 6

22 Generating series AsingleseriesencodesalltheFibonacci sequence : X F n X n = X + X 2 +2X 3 +3X 4 +5X F n X n + n 0 Fact : this series is the Taylor expansion of a rational fraction : X F n X n X = 1 X X 2 n 0 Proof :theproduct (X + X 2 +2X 3 +3X 4 +5X 5 +8X 6 + )(1 X X 2 ) is a telescoping series X + X 2 +2X 3 +3X 4 +5X 5 +8X 6 + = X. X 2 X 3 2X 4 3X 5 5X 6 X 3 X 4 2X 5 3X 6

23 Generating series of the Fibonacci sequence Remark. The denominator 1 X X 2 in the right hand side of X + X 2 +2X 3 +3X F n X n + = X 1 X X 2 is X 2 f(x 1 ),wheref(x) =X 2 X 1 is the irreducible polynomial of the Golden ratio.

24 Fibonacci and powers of matrices The Fibonacci linear recurrence relation F n+2 = F n+1 + F n for n 0 can be written Fn Fn =. F n F n+1 By induction one deduces, for n 0, n Fn = F n+1 An equivalent formula is, for n 1, n 0 1 Fn = 1 F n. 1 1 F n F n+1

25 Fibonacci and powers of matrices The Fibonacci linear recurrence relation F n+2 = F n+1 + F n for n 0 can be written Fn Fn =. F n F n+1 By induction one deduces, for n 0, n Fn = F n+1 An equivalent formula is, for n 1, n 0 1 Fn = 1 F n. 1 1 F n F n+1

26 Fibonacci and powers of matrices The Fibonacci linear recurrence relation F n+2 = F n+1 + F n for n 0 can be written Fn Fn =. F n F n+1 By induction one deduces, for n 0, n Fn = F n+1 An equivalent formula is, for n 1, n 0 1 Fn = 1 F n. 1 1 F n F n+1

27 Characteristic polynomial The characteristic polynomial of the matrix 0 1 A = 1 1 is det(xi X 1 A) =det = X 2 X 1, 1 X 1 which is the irreducible polynomial of the Golden ratio.

28 Fibonacci sequence and the Golden ratio (continued) For n 1, n 2 Z[ ] = Z + Z is a linear combination of 1 and with integer coe cients, namely n = F n 1 + F n

29 Fibonacci sequence and the Golden ratio (continued) For n 1, n 2 Z[ ] = Z + Z is a linear combination of 1 and with integer coe cients, namely n = F n 1 + F n

30 Fibonacci sequence and Hilbert s 10thproblem Yuri Matiyasevich (1970) showed that there is a polynomial P in n, m, andanumberofothervariablesx, y, z,... having the property that n = F 2m i there exist integers x, y, z,... such that P (n, m, x, y, z,... )=0. This completed the proof of the impossibility of the tenth of Hilbert s problems (does there exist a general method for solving Diophantine equations?) thankstothe previous work of Hilary Putnam, Julia Robinson and Martin Davis.

31 The Fibonacci Quarterly The Fibonacci sequence satisfies a lot of very interesting properties. Four times a year, the Fibonacci Quarterly publishes an issue with new properties which have been discovered.

32 Lucas sequence The Lucas sequence (L n ) n 0 satisfies the same recurrence relation as the Fibonacci sequence, namely L n+2 = L n+1 + L n (n 0), only the initial values are di erent : L 0 =2,L 1 =1. The sequence of Lucas numbers starts with 2, 1, 3, 4, 7, 11, 18, 29, 47, 76, 123, 199, AclosedforminvolvingtheGoldenratio is L n = n +( ) n, from which it follows that for n to n. 2, L n is the nearest integer

33 Lucas sequence The Lucas sequence (L n ) n 0 satisfies the same recurrence relation as the Fibonacci sequence, namely L n+2 = L n+1 + L n (n 0), only the initial values are di erent : L 0 =2,L 1 =1. The sequence of Lucas numbers starts with 2, 1, 3, 4, 7, 11, 18, 29, 47, 76, 123, 199, AclosedforminvolvingtheGoldenratio is L n = n +( ) n, from which it follows that for n to n. 2, L n is the nearest integer

34 François Édouard Anatole Lucas ( ) Edouard Lucas is best known for his results in number theory. He studied the Fibonacci sequence and devised the test for Mersenne primes still used today. Mathematicians/Lucas.html

35 Generating series of the Lucas sequence The generating series of the Lucas sequence X L n X n =2+X +3X 2 +4X L n X n + n 0 is nothing else than 2 X 1 X X 2

36 The Lucas sequence and power of matrices From the linear recurrence relation L n+2 = L n+1 + L n one deduces, (as we did for the Fibonacci sequence), for n 0, Ln Ln =, 1 1 L n+2 L n+1 hence Ln L n+1 = n 2. 1 Any one of the three sequences (F n ) n 0, (L n ) n 0, ( n ) n 0 can be written as a linear combination of the two others.

37 The Lucas sequence and power of matrices From the linear recurrence relation L n+2 = L n+1 + L n one deduces, (as we did for the Fibonacci sequence), for n 0, Ln Ln =, 1 1 L n+2 L n+1 hence Ln L n+1 = n 2. 1 Any one of the three sequences (F n ) n 0, (L n ) n 0, ( n ) n 0 can be written as a linear combination of the two others.

38 Perrin sequence The Perrin sequence (also called skiponacci sequence) is the linear recurrence sequence defined by with the initial conditions P n+3 = P n+1 + P n for n 0, P 0 =3,P 1 =0,P 2 =2. It starts with 3, 0, 2, 3, 2, 5, 5, 7, 10, 12, 17, 22, 29, 39, 51, François Olivier Raoul Perrin :

39 Perrin sequence The Perrin sequence (also called skiponacci sequence) is the linear recurrence sequence defined by with the initial conditions P n+3 = P n+1 + P n for n 0, P 0 =3,P 1 =0,P 2 =2. It starts with 3, 0, 2, 3, 2, 5, 5, 7, 10, 12, 17, 22, 29, 39, 51, François Olivier Raoul Perrin :

40 Plastic (or silver) constant The ratio of successive terms in the Perrin sequence approaches the plastic number %, whichistheminimal Pisot Vijayaraghavan number, real root of x 3 x 1 which has a value of approximately This constant is equal to % = p p p p 69 6

41 Plastic (or silver) constant The ratio of successive terms in the Perrin sequence approaches the plastic number %, whichistheminimal Pisot Vijayaraghavan number, real root of x 3 x 1 which has a value of approximately This constant is equal to % = p p p p 69 6

42 Perrin sequence and the plastic constant Decompose the polynomial X 3 X 1 into irreducible factors over C and over R X 3 X 1=(X %)(X )(X ). X 3 X 1=(X %)(X 2 + %X + % 1 ). Hence and are the roots of X 2 + %X + % 1.Then,for n 0, P n = % n + n + n, It follows that, for n 0, P n is the nearest integer to % n

43 Perrin sequence and the plastic constant Decompose the polynomial X 3 X 1 into irreducible factors over C and over R X 3 X 1=(X %)(X )(X ). X 3 X 1=(X %)(X 2 + %X + % 1 ). Hence and are the roots of X 2 + %X + % 1.Then,for n 0, P n = % n + n + n, It follows that, for n 0, P n is the nearest integer to % n

44 Perrin sequence and the plastic constant Decompose the polynomial X 3 X 1 into irreducible factors over C and over R X 3 X 1=(X %)(X )(X ). X 3 X 1=(X %)(X 2 + %X + % 1 ). Hence and are the roots of X 2 + %X + % 1.Then,for n 0, P n = % n + n + n, It follows that, for n 0, P n is the nearest integer to % n

45 Perrin sequence and the plastic constant Decompose the polynomial X 3 X 1 into irreducible factors over C and over R X 3 X 1=(X %)(X )(X ). X 3 X 1=(X %)(X 2 + %X + % 1 ). Hence and are the roots of X 2 + %X + % 1.Then,for n 0, P n = % n + n + n, It follows that, for n 0, P n is the nearest integer to % n

46 Generating series of the Perrin sequence The generating series of the Perrin sequence X P n X n =3+2X 2 +3X 3 +2X P n X n + n 0 is nothing else than 3 X 2 1 X 2 X 3 The denominator 1 X 2 X 3 is X 3 f(x 1 ) where f(x) =X 3 X 1 is the irreducible polynomial of %.

47 Generating series of the Perrin sequence The generating series of the Perrin sequence X P n X n =3+2X 2 +3X 3 +2X P n X n + n 0 is nothing else than 3 X 2 1 X 2 X 3 The denominator 1 X 2 X 3 is X 3 f(x 1 ) where f(x) =X 3 X 1 is the irreducible polynomial of %.

48 Perrin sequence and power of matrices From we deduce P n+3 = P n+1 + P n P n P P n+2 A 0 0 P n+1 A P n P n+2 Hence P n 0 1 P n+1 A 0 0 1A P n A 2 n 0

49 Perrin sequence and power of matrices From we deduce P n+3 = P n+1 + P n P n P P n+2 A 0 0 P n+1 A P n P n+2 Hence P n 0 1 P n+1 A 0 0 1A P n A 2 n 0

50 Characteristic polynomial The characteristic polynomial of the matrix A 0 0 1A is det(xi 0 X A) = 0 X 1A = X 3 X 1, 1 1 X which is the irreducible polynomial of the plastic constant %.

51 Perrin pseudoprimes If p is prime, then p divides P p. The smallest composite n such that n divides P n is For n either = or = , the number n divides P n. Jon Grantham has proved that there are infinitely many Perrin pseudoprimes. The number c of decimal digits of P satisfies 10 c = % ,hencec =271441(log%)/(log 10) The website maintained by Richard Turk is devoted to Perrin numbers. See OEISA

52 Perrin pseudoprimes If p is prime, then p divides P p. The smallest composite n such that n divides P n is For n either = or = , the number n divides P n. Jon Grantham has proved that there are infinitely many Perrin pseudoprimes. The number c of decimal digits of P satisfies 10 c = % ,hencec =271441(log%)/(log 10) The website maintained by Richard Turk is devoted to Perrin numbers. See OEISA

53 Perrin pseudoprimes If p is prime, then p divides P p. The smallest composite n such that n divides P n is For n either = or = , the number n divides P n. Jon Grantham has proved that there are infinitely many Perrin pseudoprimes. The number c of decimal digits of P satisfies 10 c = % ,hencec =271441(log%)/(log 10) The website maintained by Richard Turk is devoted to Perrin numbers. See OEISA

54 Perrin pseudoprimes If p is prime, then p divides P p. The smallest composite n such that n divides P n is For n either = or = , the number n divides P n. Jon Grantham has proved that there are infinitely many Perrin pseudoprimes. The number c of decimal digits of P satisfies 10 c = % ,hencec =271441(log%)/(log 10) The website maintained by Richard Turk is devoted to Perrin numbers. See OEISA

55 Padovan sequence The Padovan sequence p n satisfies the same recurrence p n+3 = p n+1 + p n as the Perrin sequence but has di erent initial values : p 0 =1, p 1 = p 2 =0. It starts with 1, 0, 0, 1, 0, 1, 1, 1, 2, 2, 3, 4, 5, 7, 9, 12, 16, 21, 28, 37, 49, Richard Padovan

56 Generating series and power of matrices 1+X 3 + X p n X n + = 1 X 2 1 X 2 X 3 For n 0, p n 0 1 p n+1 A 0 0 1A p n A 0 n 0

57 Generating series and power of matrices 1+X 3 + X p n X n + = 1 X 2 1 X 2 X 3 For n 0, p n 0 1 p n+1 A 0 0 1A p n A 0 n 0

58 Padovan triangles p n = p n 2 + p n 3 p n 1 = p n 3 + p n 4 Hence p n 2 = p n 4 + p n 5 p n p n 1 = p n 5 p n = p n 1 + p n 5

59 Padovan triangles

60 Padovan triangles vs Fibonacci squares

61 Narayana sequence Narayana sequence is defined by the recurrence relation C n+3 = C n+2 + C n with the initial values C 0 =2, C 1 =3, C 2 =4. It starts with 2, 3, 4, 6, 9, 13, 19, 28, 41, 60, 88, 129, 189, 277, 406, 595,... Real root of x 3 x 2 1 s p s p =

62 Narayana sequence Narayana sequence is defined by the recurrence relation C n+3 = C n+2 + C n with the initial values C 0 =2, C 1 =3, C 2 =4. It starts with 2, 3, 4, 6, 9, 13, 19, 28, 41, 60, 88, 129, 189, 277, 406, 595,... Real root of x 3 x 2 1 s p s p =

63 Generating series and power of matrices 2+3X +4X 2 +6X C n X n + = 2+X + X2 1 X X 3 For n 0, C n 0 1 C n+1 A 0 0 1A C n A 4 n 0

64 Generating series and power of matrices 2+3X +4X 2 +6X C n X n + = 2+X + X2 1 X X 3 For n 0, C n 0 1 C n+1 A 0 0 1A C n A 4 n 0

65 Narayana s cows Narayana was an Indian mathematician in the 14th century who proposed the following problem : Acowproducesonecalfeveryyear.Beginninginitsfourth year each calf produces one calf at the beginning of each year. How many calves are there altogether after, for example, 17 years?

66 Music : In working this out, Tom Johnson found a way to translate this into a composition called Narayana s Cows. Music : Tom Johnson Saxophones : Daniel Kientzy

67 Narayana s cows Year Original Cow Second generation Third generation Fourth generation Fifth generation Sixth generation Total

68 Jean-Paul Allouche and Tom Johnson bibliorecente.html

69 Cows, music and morphisms Jean-Paul Allouche and Tom Johnson Narayana s Cows and Delayed Morphisms In 3èmes Journées d Informatique Musicale (JIM 96), Ile de Tatihou, Les Cahiers du GREYC (1996 no. 4), pages 2-7, May Finite automata and morphisms in assisted musical composition, Journal of New Music Research, no. 24 (1995), nl/jnmr/vol24_2.html

70 Linear recurrence sequences : definitions A linear recurrence sequence is a sequence of numbers u =(u 0,u 1,u 2,...) for which there exist a positive integer d together with numbers a 1,...,a d with a d 6=0such that, for n 0, (?) u n+d = a 1 u n+d a d u n. Here, a number means an element of a field K of zero characteristic. Given a =(a 1,...,a d ) 2 K d,thesete a of linear recurrence sequences u =(u n ) n 0 satisfying (?) isak vector subspace of dimension d of the space K N of all sequences. The characteristic (or companion) polynomial of the linear recurrence is f(x) =X d a 1 X d 1 a d.

71 Linear recurrence sequences : definitions A linear recurrence sequence is a sequence of numbers u =(u 0,u 1,u 2,...) for which there exist a positive integer d together with numbers a 1,...,a d with a d 6=0such that, for n 0, (?) u n+d = a 1 u n+d a d u n. Here, a number means an element of a field K of zero characteristic. Given a =(a 1,...,a d ) 2 K d,thesete a of linear recurrence sequences u =(u n ) n 0 satisfying (?) isak vector subspace of dimension d of the space K N of all sequences. The characteristic (or companion) polynomial of the linear recurrence is f(x) =X d a 1 X d 1 a d.

72 Linear recurrence sequences : definitions A linear recurrence sequence is a sequence of numbers u =(u 0,u 1,u 2,...) for which there exist a positive integer d together with numbers a 1,...,a d with a d 6=0such that, for n 0, (?) u n+d = a 1 u n+d a d u n. Here, a number means an element of a field K of zero characteristic. Given a =(a 1,...,a d ) 2 K d,thesete a of linear recurrence sequences u =(u n ) n 0 satisfying (?) isak vector subspace of dimension d of the space K N of all sequences. The characteristic (or companion) polynomial of the linear recurrence is f(x) =X d a 1 X d 1 a d.

73 Linear recurrence sequences : definitions A linear recurrence sequence is a sequence of numbers u =(u 0,u 1,u 2,...) for which there exist a positive integer d together with numbers a 1,...,a d with a d 6=0such that, for n 0, (?) u n+d = a 1 u n+d a d u n. Here, a number means an element of a field K of zero characteristic. Given a =(a 1,...,a d ) 2 K d,thesete a of linear recurrence sequences u =(u n ) n 0 satisfying (?) isak vector subspace of dimension d of the space K N of all sequences. The characteristic (or companion) polynomial of the linear recurrence is f(x) =X d a 1 X d 1 a d.

74 Linear recurrence sequence : examples Constant sequence : u n = u 0. Linear recurrence sequence of order 1 : u n+1 = u n. Characteristic polynomial : f(x) =X 1. Generating series : X n 0 X n = 1 1 X Geometric progression : u n = u n 0. Linear recurrence sequence of order 1 : u n = u n 1. Characteristic polynomial f(x) =X. Generating series : X n 0 u 0 n X n = u 0 1 X

75 Linear recurrence sequence : examples Constant sequence : u n = u 0. Linear recurrence sequence of order 1 : u n+1 = u n. Characteristic polynomial : f(x) =X 1. Generating series : X n 0 X n = 1 1 X Geometric progression : u n = u n 0. Linear recurrence sequence of order 1 : u n = u n 1. Characteristic polynomial f(x) =X. Generating series : X n 0 u 0 n X n = u 0 1 X

76 Linear recurrence sequence : examples u n = n. Linearrecurrencesequenceoforder2 : Characteristic polynomial n +2=2(n +1) n. f(x) =X 2 2X +1=(X 1) 2. Generating series X nx n = n X + X 2 Power of matrices : n 0 1 = 1 2 n +1 n. n n+1

77 Linear recurrence sequence : examples u n = f(n), f polynomial of degree d. Linearrecurrence sequence of order d +1. Proof. The sequences (f(n)) n 0, (f(n +1)) n 0,, (f(n + k)) n 0 are K linearly independent in K N for k = d dependent for k = d. 1 and linearly

78 Linear recurrence sequence : examples u n = f(n), f polynomial of degree d. Linearrecurrence sequence of order d +1. Proof. The sequences (f(n)) n 0, (f(n +1)) n 0,, (f(n + k)) n 0 are K linearly independent in K N for k = d dependent for k = d. 1 and linearly

79 Linear sequences which are ultimately recurrent The sequence 1, 0, 0,... is not a linear recurrence sequence. The condition is satisfied only for n 1. u n+1 = u n The relation u n+2 = u n+1 +0u n with d =2, a d =0does not fulfill the requirement a d 6=0.

80 Linear sequences which are ultimately recurrent The sequence 1, 0, 0,... is not a linear recurrence sequence. The condition is satisfied only for n 1. u n+1 = u n The relation u n+2 = u n+1 +0u n with d =2, a d =0does not fulfill the requirement a d 6=0.

81 Linear sequences which are ultimately recurrent The sequence 1, 0, 0,... is not a linear recurrence sequence. The condition is satisfied only for n 1. u n+1 = u n The relation u n+2 = u n+1 +0u n with d =2, a d =0does not fulfill the requirement a d 6=0.

82 Order of a linear recurrence sequence If u =(u n ) n 0 satisfies the linear recurrence, the characteristic polynomial of which is f, then,foranymonicpolynomial g 2 K[X], thissequenceualso satisfies the linear recurrence, the characteristic polynomial of which is fg. Example :forg(x) =X,from (?) u n+d a 1 u n+d 1 a d u n =0 we deduce u n+d+1 a 1 u n+d a d u n+1 + (u n+d a 1 u n+d 1 a d u n )=0. The order of a linear recurrence sequence is the smallest d such that (?) holdsforalln 0.

83 Order of a linear recurrence sequence If u =(u n ) n 0 satisfies the linear recurrence, the characteristic polynomial of which is f, then,foranymonicpolynomial g 2 K[X], thissequenceualso satisfies the linear recurrence, the characteristic polynomial of which is fg. Example :forg(x) =X,from (?) u n+d a 1 u n+d 1 a d u n =0 we deduce u n+d+1 a 1 u n+d a d u n+1 + (u n+d a 1 u n+d 1 a d u n )=0. The order of a linear recurrence sequence is the smallest d such that (?) holdsforalln 0.

84 Polynomial combinations of powers The sum of any two linear recurrence sequences is a linear recurrence sequence. The set [ a E a of all linear recurrence sequences with coe cients in K is a sub K algebra of K N. Given polynomials p 1,...,p` in K[X] and elements in K,thesequence 1,..., ` is a linear recurrence sequence. p 1 (n) n p`(n) ǹ n 0 Fact : any linear recurrence sequence is of this form.

85 Polynomial combinations of powers The sum of any two linear recurrence sequences is a linear recurrence sequence. The set [ a E a of all linear recurrence sequences with coe cients in K is a sub K algebra of K N. Given polynomials p 1,...,p` in K[X] and elements in K,thesequence 1,..., ` is a linear recurrence sequence. p 1 (n) n p`(n) ǹ n 0 Fact : any linear recurrence sequence is of this form.

86 Polynomial combinations of powers The sum of any two linear recurrence sequences is a linear recurrence sequence. The set [ a E a of all linear recurrence sequences with coe cients in K is a sub K algebra of K N. Given polynomials p 1,...,p` in K[X] and elements in K,thesequence 1,..., ` is a linear recurrence sequence. p 1 (n) n p`(n) ǹ n 0 Fact : any linear recurrence sequence is of this form.

87 Polynomial combinations of powers The sum of any two linear recurrence sequences is a linear recurrence sequence. The set [ a E a of all linear recurrence sequences with coe cients in K is a sub K algebra of K N. Given polynomials p 1,...,p` in K[X] and elements in K,thesequence 1,..., ` is a linear recurrence sequence. p 1 (n) n p`(n) ǹ n 0 Fact : any linear recurrence sequence is of this form.

88 Linear recurrence sequence and Brahmagupta Pell Fermat Equation Let d be a positive integer, not a square. The solutions (x, y) 2 Z Z of the Brahmagupta Pell Fermat Equation x 2 dy 2 = ±1 form a sequence (x n,y n ) n2z defined by From x n + p dy n =(x 1 + p dy 1 ) n. 2x n =(x 1 + p dy 1 ) n +(x 1 p dy1 ) n we deduce that (x n ) n 0 is a linear recurrence sequence. Same for y n,andalsoforn 0.

89 Linear recurrence sequence and Brahmagupta Pell Fermat Equation Let d be a positive integer, not a square. The solutions (x, y) 2 Z Z of the Brahmagupta Pell Fermat Equation x 2 dy 2 = ±1 form a sequence (x n,y n ) n2z defined by From x n + p dy n =(x 1 + p dy 1 ) n. 2x n =(x 1 + p dy 1 ) n +(x 1 p dy1 ) n we deduce that (x n ) n 0 is a linear recurrence sequence. Same for y n,andalsoforn 0.

90 Doubly infinite linear recurrence sequence Asequence(u n ) n2z indexed by Z is a linear recurrence sequence if it satisfies (?) u n+d = a 1 u n+d a d u n. for all n 2 Z. Recall a d 6=0. Such a sequence is determined by d consecutive values.

91 Doubly infinite linear recurrence sequence Asequence(u n ) n2z indexed by Z is a linear recurrence sequence if it satisfies (?) u n+d = a 1 u n+d a d u n. for all n 2 Z. Recall a d 6=0. Such a sequence is determined by d consecutive values.

92 Doubly infinite linear recurrence sequence Asequence(u n ) n2z indexed by Z is a linear recurrence sequence if it satisfies (?) u n+d = a 1 u n+d a d u n. for all n 2 Z. Recall a d 6=0. Such a sequence is determined by d consecutive values.

93 Linear recurrence sequence : simple roots AbasisofE a over K is obtained by attributing to the initial values u 0,...,u d 1 the values given by the canonical basis of K d. Given in K,anecessaryandsu cientconditionfora sequence ( n ) n 0 to satisfy (?) isthat is a root of the characteristic polynomial f(x) =X d a 1 X d 1 a d. If this polynomial has d distinct roots 1,..., d in K, f(x) =(X 1) (X d), i 6= j, then a basis of E a over K is given by the d sequences ( n i ) n 0, i =1,...,d.

94 Linear recurrence sequence : simple roots AbasisofE a over K is obtained by attributing to the initial values u 0,...,u d 1 the values given by the canonical basis of K d. Given in K,anecessaryandsu cientconditionfora sequence ( n ) n 0 to satisfy (?) isthat is a root of the characteristic polynomial f(x) =X d a 1 X d 1 a d. If this polynomial has d distinct roots 1,..., d in K, f(x) =(X 1) (X d), i 6= j, then a basis of E a over K is given by the d sequences ( n i ) n 0, i =1,...,d.

95 Linear recurrence sequence : simple roots AbasisofE a over K is obtained by attributing to the initial values u 0,...,u d 1 the values given by the canonical basis of K d. Given in K,anecessaryandsu cientconditionfora sequence ( n ) n 0 to satisfy (?) isthat is a root of the characteristic polynomial f(x) =X d a 1 X d 1 a d. If this polynomial has d distinct roots 1,..., d in K, f(x) =(X 1) (X d), i 6= j, then a basis of E a over K is given by the d sequences ( n i ) n 0, i =1,...,d.

96 Linear recurrence sequence : double roots The characteristic polynomial of the linear recurrence u n =2 u n 1 2 u n 2 is X 2 2 X + 2 =(X ) 2 with a double root. The sequence (n n ) n 0 satisfies n n =2 (n 1)n n 1 2 (n 2) n 2. AbasisofE a for a 1 =2, a 2 = sequences ( n ) n 0, (n n ) n 0. 2 is given by the two Given 2 K,anecessaryandsu cientconditionforthe sequence n n to satisfy the linear recurrence relation (?) is that is a root of multiplicity 2 of f(x).

97 Linear recurrence sequence : double roots The characteristic polynomial of the linear recurrence u n =2 u n 1 2 u n 2 is X 2 2 X + 2 =(X ) 2 with a double root. The sequence (n n ) n 0 satisfies n n =2 (n 1)n n 1 2 (n 2) n 2. AbasisofE a for a 1 =2, a 2 = sequences ( n ) n 0, (n n ) n 0. 2 is given by the two Given 2 K,anecessaryandsu cientconditionforthe sequence n n to satisfy the linear recurrence relation (?) is that is a root of multiplicity 2 of f(x).

98 Linear recurrence sequence : double roots The characteristic polynomial of the linear recurrence u n =2 u n 1 2 u n 2 is X 2 2 X + 2 =(X ) 2 with a double root. The sequence (n n ) n 0 satisfies n n =2 (n 1)n n 1 2 (n 2) n 2. AbasisofE a for a 1 =2, a 2 = sequences ( n ) n 0, (n n ) n 0. 2 is given by the two Given 2 K,anecessaryandsu cientconditionforthe sequence n n to satisfy the linear recurrence relation (?) is that is a root of multiplicity 2 of f(x).

99 Linear recurrence sequence : double roots The characteristic polynomial of the linear recurrence u n =2 u n 1 2 u n 2 is X 2 2 X + 2 =(X ) 2 with a double root. The sequence (n n ) n 0 satisfies n n =2 (n 1)n n 1 2 (n 2) n 2. AbasisofE a for a 1 =2, a 2 = sequences ( n ) n 0, (n n ) n 0. 2 is given by the two Given 2 K,anecessaryandsu cientconditionforthe sequence n n to satisfy the linear recurrence relation (?) is that is a root of multiplicity 2 of f(x).

100 Linear recurrence sequence : multiple roots In general, when the characteristic polynomial splits as Ỳ X d a 1 X d 1 a d = (X i) t i, i=1 abasisofe a is given by the d sequences (n k i n ) n 0, 0 apple k apple t i 1, 1 apple i apple `.

101 Generating series of a linear recurrence sequence Let u =(u n ) n 0 be a linear recurrence sequence (?) u n+d = a 1 u n+d a d u n for n 0 with characteristic polynomial f(x) =X d a 1 X d 1 a d. Denote by f the reciprocal polynomial of f : Then f (X) =X d f(x 1 )=1 a 1 X a 0 X d. 1X n=0 u n X n = r(x) f (X), where r is a polynomial of degree less than d determined by the initial values of u.

102 Generating series of a linear recurrence sequence Assume u n+d = a 1 u n+d a d u n for n 0. Then 1X n=0 u n X n = r(x) f (X) Proof. Comparingthecoe cientsofx n for n d shows that 1X f (X) u n X n n=0 is a polynomial of degree less than d

103 Generating series of a linear recurrence sequence Assume u n+d = a 1 u n+d a d u n for n 0. Then 1X n=0 u n X n = r(x) f (X) Proof. Comparingthecoe cientsofx n for n d shows that 1X f (X) u n X n n=0 is a polynomial of degree less than d

104 Taylor coe cients of rational functions Conversely, the coe cients the Taylor expansion of any rational fraction a(x)/b(x) with deg a<deg b satisfies the recurrence relation with characteristic polynomial f 2 K[X] given by f(x) =b (X). Therefore a sequence u =(u n ) n 0 satisfies the recurrence relation (?) withcharacteristicpolynomialf 2 K[X] if and only if 1X u n X n = r(x) f (X), n=0 where r is a polynomial of degree less than d determined by the initial values of u.

105 Taylor coe cients of rational functions Conversely, the coe cients the Taylor expansion of any rational fraction a(x)/b(x) with deg a<deg b satisfies the recurrence relation with characteristic polynomial f 2 K[X] given by f(x) =b (X). Therefore a sequence u =(u n ) n 0 satisfies the recurrence relation (?) withcharacteristicpolynomialf 2 K[X] if and only if 1X u n X n = r(x) f (X), n=0 where r is a polynomial of degree less than d determined by the initial values of u.

106 Matrix notation for a linear recurrence sequence The linear recurrence sequence (?) u n+d = a 1 u n+d a d u n for n 0 can be written u n+1 u n B A =. B A u n+d a d a d 1 a d 2 a 1 u n u n+1. u n+d 1 1 C A.

107 Matrix notation for a linear recurrence sequence with 0 U n = u n u n+1. u n+d 1 The determinant of I d X A) isnothingbut U n+1 = AU n C A, A =. B A a d a d 1 a d 2 a 1 A (the characteristic polynomial of f(x) =X d a 1 X d 1 a d, the characteristic polynomial of the linear recurrence sequence. By induction U n = A n U 0.

108 Matrix notation for a linear recurrence sequence with 0 U n = u n u n+1. u n+d 1 The determinant of I d X A) isnothingbut U n+1 = AU n C A, A =. B A a d a d 1 a d 2 a 1 A (the characteristic polynomial of f(x) =X d a 1 X d 1 a d, the characteristic polynomial of the linear recurrence sequence. By induction U n = A n U 0.

109 Matrix notation for a linear recurrence sequence with 0 U n = u n u n+1. u n+d 1 The determinant of I d X A) isnothingbut U n+1 = AU n C A, A =. B A a d a d 1 a d 2 a 1 A (the characteristic polynomial of f(x) =X d a 1 X d 1 a d, the characteristic polynomial of the linear recurrence sequence. By induction U n = A n U 0.

110 Powers of matrices Let A =(a ij ) 1applei,jappled 2 GL d d (K) be a d d matrix with coe cients in K and nonzero determinant. For n 0, define A n = a ij (n) 1applei,jappled. Then each of the d 2 sequences a ij (n) n 0, (1 apple i, j apple d) is a linear recurrence sequence. The roots of the characteristic polynomial of these linear recurrences are the eigenvalues of A. In particular the sequence Tr(A n ) n 0 satisfies the linear recurrence, the characteristic polynomial of which is the characteristic polynomial of the matrix A.

111 Powers of matrices Let A =(a ij ) 1applei,jappled 2 GL d d (K) be a d d matrix with coe cients in K and nonzero determinant. For n 0, define A n = a ij (n) 1applei,jappled. Then each of the d 2 sequences a ij (n) n 0, (1 apple i, j apple d) is a linear recurrence sequence. The roots of the characteristic polynomial of these linear recurrences are the eigenvalues of A. In particular the sequence Tr(A n ) n 0 satisfies the linear recurrence, the characteristic polynomial of which is the characteristic polynomial of the matrix A.

112 Powers of matrices Let A =(a ij ) 1applei,jappled 2 GL d d (K) be a d d matrix with coe cients in K and nonzero determinant. For n 0, define A n = a ij (n) 1applei,jappled. Then each of the d 2 sequences a ij (n) n 0, (1 apple i, j apple d) is a linear recurrence sequence. The roots of the characteristic polynomial of these linear recurrences are the eigenvalues of A. In particular the sequence Tr(A n ) n 0 satisfies the linear recurrence, the characteristic polynomial of which is the characteristic polynomial of the matrix A.

113 Conversely : Given a linear recurrence sequence u 2 K N,thereexistan integer d 1 and a matrix A 2 GL d (K) such that, for each n 0, u n = a 11 (n). The characteristic polynomial of A is the characteristic polynomial of the linear recurrence sequence. Everest G., van der Poorten A., Shparlinski I., Ward T. Recurrence sequences, Mathematical Surveys and Monographs (AMS, 2003), volume 104.

114 Conversely : Given a linear recurrence sequence u 2 K N,thereexistan integer d 1 and a matrix A 2 GL d (K) such that, for each n 0, u n = a 11 (n). The characteristic polynomial of A is the characteristic polynomial of the linear recurrence sequence. Everest G., van der Poorten A., Shparlinski I., Ward T. Recurrence sequences, Mathematical Surveys and Monographs (AMS, 2003), volume 104.

115 Discrete version of linear di erential equations Asequenceu 2 K N can be viewed as a linear map N! K. Define the discrete derivative D by Du : N! K n 7! u n+1 u n. Asequenceu 2 K N is a linear recurrence sequence if and only if there exists Q 2 K[T ] with Q(1) 6= 1such that Q(D)u =0. Linear recurrence sequences are a discrete version of linear di erential equations with constant coe cients. The condition Q(1) 6= 0reflects a d 6=0 otherwise one gets ultimately recurrent sequences.

116 Joint work with Claude Levesque Linear recurrence sequences and twisted binary forms. Proceedings of the International Conference on Pure and Applied Mathematics ICPAM-GOROKA South Pacific Journal of Pure and Applied Mathematics. pdf/procconfpng2014.pdf

117 Families of binary forms Consider a binary form F 0 (X, Y ) 2 C[X, Y ] which satisfies F 0 (1, 0) = 1. Wewriteitas F 0 (X, Y )=X d + a 1 X d 1 Y + + a d Y d = dy (X i Y ). Let 1,..., d be d nonzero complex numbers not necessarily distinct. Twisting F 0 by the powers n 1,..., n d (n 2 Z), we obtain the family of binary forms dy F n (X, Y )= (X i n i Y ), which we write as Therefore i=1 i=1 X d U 1 (n)x d 1 Y + +( 1) d U d (n)y d. U h (0) = ( 1) h a h (1 apple h apple d).

118 Families of binary forms Consider a binary form F 0 (X, Y ) 2 C[X, Y ] which satisfies F 0 (1, 0) = 1. Wewriteitas F 0 (X, Y )=X d + a 1 X d 1 Y + + a d Y d = dy (X i Y ). Let 1,..., d be d nonzero complex numbers not necessarily distinct. Twisting F 0 by the powers n 1,..., n d (n 2 Z), we obtain the family of binary forms dy F n (X, Y )= (X i n i Y ), which we write as Therefore i=1 i=1 X d U 1 (n)x d 1 Y + +( 1) d U d (n)y d. U h (0) = ( 1) h a h (1 apple h apple d).

119 Families of binary forms Consider a binary form F 0 (X, Y ) 2 C[X, Y ] which satisfies F 0 (1, 0) = 1. Wewriteitas F 0 (X, Y )=X d + a 1 X d 1 Y + + a d Y d = dy (X i Y ). Let 1,..., d be d nonzero complex numbers not necessarily distinct. Twisting F 0 by the powers n 1,..., n d (n 2 Z), we obtain the family of binary forms dy F n (X, Y )= (X i n i Y ), which we write as Therefore i=1 i=1 X d U 1 (n)x d 1 Y + +( 1) d U d (n)y d. U h (0) = ( 1) h a h (1 apple h apple d).

120 Families of binary forms Consider a binary form F 0 (X, Y ) 2 C[X, Y ] which satisfies F 0 (1, 0) = 1. Wewriteitas F 0 (X, Y )=X d + a 1 X d 1 Y + + a d Y d = dy (X i Y ). Let 1,..., d be d nonzero complex numbers not necessarily distinct. Twisting F 0 by the powers n 1,..., n d (n 2 Z), we obtain the family of binary forms dy F n (X, Y )= (X i n i Y ), which we write as Therefore i=1 i=1 X d U 1 (n)x d 1 Y + +( 1) d U d (n)y d. U h (0) = ( 1) h a h (1 apple h apple d).

121 Families of Diophantine equations With Claude Levesque, we consider some families of diophantine equations F n (x, y) =m obtained in the same way from a given irreducible form F (X, Y ) with coe cients in Z, when 1,..., d are algebraic units and when the algebraic numbers 1 1,..., d d are Galois conjugates with d 3. Theorem. Let K be a number field of degree d 3, S afinite set of places of K containing the places at infinity. Denote by O S the ring of S integers of K and by O S the group of S units of K. Assume 1,..., d, 1,..., d belong to K Then there are only finitely many (x, y, n) in O S O S Z satisfying F n (x, y) 2O S, xy 6= 0 and Card{ 1 n 1,..., d n d} 3.

122 Families of Diophantine equations With Claude Levesque, we consider some families of diophantine equations F n (x, y) =m obtained in the same way from a given irreducible form F (X, Y ) with coe cients in Z, when 1,..., d are algebraic units and when the algebraic numbers 1 1,..., d d are Galois conjugates with d 3. Theorem. Let K be a number field of degree d 3, S afinite set of places of K containing the places at infinity. Denote by O S the ring of S integers of K and by O S the group of S units of K. Assume 1,..., d, 1,..., d belong to K Then there are only finitely many (x, y, n) in O S O S Z satisfying F n (x, y) 2O S, xy 6= 0 and Card{ 1 n 1,..., d n d} 3.

123 Families of Diophantine equations Each of the sequences U h (n) n2z coming from the coe cients of the relation F n (X, Y )=X d U 1 (n)x d 1 Y + +( 1) d U d (n)y d is a linear recurrence sequence. For example, for n 2 Z, U 1 (n) = dx i n i, U d (n) = i=1 dy i n i. i=1 For 1 apple h apple d, thesequence combination of the sequences U h (n) n2z is a linear ( i1 ih ) n n2z, (1 apple i 1 < <i h apple d).

124 Families of Diophantine equations Each of the sequences U h (n) n2z coming from the coe cients of the relation F n (X, Y )=X d U 1 (n)x d 1 Y + +( 1) d U d (n)y d is a linear recurrence sequence. For example, for n 2 Z, U 1 (n) = dx i n i, U d (n) = i=1 dy i n i. i=1 For 1 apple h apple d, thesequence combination of the sequences U h (n) n2z is a linear ( i1 ih ) n n2z, (1 apple i 1 < <i h apple d).

125 Families of Diophantine equations Each of the sequences U h (n) n2z coming from the coe cients of the relation F n (X, Y )=X d U 1 (n)x d 1 Y + +( 1) d U d (n)y d is a linear recurrence sequence. For example, for n 2 Z, U 1 (n) = dx i n i, U d (n) = i=1 dy i n i. i=1 For 1 apple h apple d, thesequence combination of the sequences U h (n) n2z is a linear ( i1 ih ) n n2z, (1 apple i 1 < <i h apple d).

126 Units of Bernstein and Hasse Let t and s be two positive integers, D an integer c 2 { 1, +1}. Let!>1 satisfy 1, and! st = D st + c, where it is assumed that Q(!) is of degree st. Consider = D!, = D t! t. L. Bernstein and H. Hasse noticed that and are units of degree st and s respectively, and showed that these units can be obtained from the Jacobi Perron algorithm. H.-J. Stender proved that for s = t =2, {, } is a fundamental system of units of the quartic field Q(!).

127 Units of Bernstein and Hasse Let t and s be two positive integers, D an integer c 2 { 1, +1}. Let!>1 satisfy 1, and! st = D st + c, where it is assumed that Q(!) is of degree st. Consider = D!, = D t! t. L. Bernstein and H. Hasse noticed that and are units of degree st and s respectively, and showed that these units can be obtained from the Jacobi Perron algorithm. H.-J. Stender proved that for s = t =2, {, } is a fundamental system of units of the quartic field Q(!).

128 Helmut Hasse ( ) D>0, s 1, t 1, c 2 { 1, +1},!>0,! st = D st + c, = D!, = D t! t. ( D) st =( 1) st (D st + c).

129 Diophantine equations associated with some units of Bernstein and Hasse The irreducible polynomial of is F 0 (X, 1), with F 0 (X, Y )=(X DY ) st ( 1) st (D st + c)y st. For n 2 Z, thebinaryformf n (X, Y ), obtainedbytwisting F 0 (X, Y ) with the powers n of, isthehomogeneousversion of the irreducible polynomial F n (X, 1) of n.sof n depends of the parameters n, D, s, t and c. Theorem (with Claude Levesque). Suppose st 3. There exists an e ectively computable constant apple, dependingonlyon D, s and t, withthefollowingproperty.letm, a, x, y be rational integers satisfying m 2, xy 6= 0, [Q( a ):Q] =st and F n (x, y) applem. Then max{log x, log y, n } apple apple log m.

130 Diophantine equations associated with some units of Bernstein and Hasse The irreducible polynomial of is F 0 (X, 1), with F 0 (X, Y )=(X DY ) st ( 1) st (D st + c)y st. For n 2 Z, thebinaryformf n (X, Y ), obtainedbytwisting F 0 (X, Y ) with the powers n of, isthehomogeneousversion of the irreducible polynomial F n (X, 1) of n.sof n depends of the parameters n, D, s, t and c. Theorem (with Claude Levesque). Suppose st 3. There exists an e ectively computable constant apple, dependingonlyon D, s and t, withthefollowingproperty.letm, a, x, y be rational integers satisfying m 2, xy 6= 0, [Q( a ):Q] =st and F n (x, y) applem. Then max{log x, log y, n } apple apple log m.

131 Diophantine equations associated with some units of Bernstein and Hasse The irreducible polynomial of is F 0 (X, 1), with F 0 (X, Y )=(X DY ) st ( 1) st (D st + c)y st. For n 2 Z, thebinaryformf n (X, Y ), obtainedbytwisting F 0 (X, Y ) with the powers n of, isthehomogeneousversion of the irreducible polynomial F n (X, 1) of n.sof n depends of the parameters n, D, s, t and c. Theorem (with Claude Levesque). Suppose st 3. There exists an e ectively computable constant apple, dependingonlyon D, s and t, withthefollowingproperty.letm, a, x, y be rational integers satisfying m 2, xy 6= 0, [Q( a ):Q] =st and F n (x, y) applem. Then max{log x, log y, n } apple apple log m.

132 8th ICSMEDC December 6, 2015 The Eighth International Conference on Science and Mathematics Education in Developing Countries Yangon University, Yangon, The Republic of the Union of Myanmar. Linear recurrence sequences, exponential polynomials and Diophantine approximation Michel Waldschmidt Institut de Mathématiques de Jussieu Paris VI

Abstract. Michel Waldschmidt. Linear recurrence sequences : examples

Abstract. Michel Waldschmidt. Linear recurrence sequences : examples Updated: 09//207 Abstract Linear recurrence sequences, exponential polynomials and Diophantine approximation Michel Waldschmidt Institut de Mathématiques de Jussieu Paris VI http://webusersim-prgfr/~michelwaldschmidt/

More information

Introduction to Lucas Sequences

Introduction to Lucas Sequences A talk given at Liaoning Normal Univ. (Dec. 14, 017) Introduction to Lucas Sequences Zhi-Wei Sun Nanjing University Nanjing 10093, P. R. China zwsun@nju.edu.cn http://math.nju.edu.cn/ zwsun Dec. 14, 017

More information

Abstract. Roma Tre April 18 20, th Mini Symposium of the Roman Number Theory Association (RNTA) Michel Waldschmidt

Abstract. Roma Tre April 18 20, th Mini Symposium of the Roman Number Theory Association (RNTA) Michel Waldschmidt Roma Tre April 8 0, 08 4th Mini Symposium of the Roman umber Theory Association (RTA) Representation of integers by cyclotomic binary forms Michel Waldschmidt Institut de Mathématiques de Jussieu Paris

More information

Continued Fractions: Introduction and Applications

Continued Fractions: Introduction and Applications PROCEEDINGS OF THE ROMAN NUMBER THEORY ASSOCIATION Volume 2, Number, March 207, pages 6-8 Michel Waldschmidt Continued Fractions: Introduction and Applications written by Carlo Sanna The continued fraction

More information

Part II. Number Theory. Year

Part II. Number Theory. Year Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 Paper 3, Section I 1G 70 Explain what is meant by an Euler pseudoprime and a strong pseudoprime. Show that 65 is an Euler

More information

Notes on Continued Fractions for Math 4400

Notes on Continued Fractions for Math 4400 . Continued fractions. Notes on Continued Fractions for Math 4400 The continued fraction expansion converts a positive real number α into a sequence of natural numbers. Conversely, a sequence of natural

More information

A NICE DIOPHANTINE EQUATION. 1. Introduction

A NICE DIOPHANTINE EQUATION. 1. Introduction A NICE DIOPHANTINE EQUATION MARIA CHIARA BRAMBILLA Introduction One of the most interesting and fascinating subjects in Number Theory is the study of diophantine equations A diophantine equation is a polynomial

More information

#A11 INTEGERS 12 (2012) FIBONACCI VARIATIONS OF A CONJECTURE OF POLIGNAC

#A11 INTEGERS 12 (2012) FIBONACCI VARIATIONS OF A CONJECTURE OF POLIGNAC #A11 INTEGERS 12 (2012) FIBONACCI VARIATIONS OF A CONJECTURE OF POLIGNAC Lenny Jones Department of Mathematics, Shippensburg University, Shippensburg, Pennsylvania lkjone@ship.edu Received: 9/17/10, Revised:

More information

Power series and Taylor series

Power series and Taylor series Power series and Taylor series D. DeTurck University of Pennsylvania March 29, 2018 D. DeTurck Math 104 002 2018A: Series 1 / 42 Series First... a review of what we have done so far: 1 We examined series

More information

A tour through Elliptic Divisibility Sequences

A tour through Elliptic Divisibility Sequences A tour through Elliptic Divisibility Sequences Victor S. Miller CCR Princeton Princeton, NJ 08540 15 April 2010 Points and their denominators Let E : y 2 + a 1 xy + a 3 y = x 3 + a 2 x 2 + a 4 x + a 6

More information

POLYNOMIAL SOLUTIONS OF PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS

POLYNOMIAL SOLUTIONS OF PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS J. London Math. Soc. 67 (2003) 16 28 C 2003 London Mathematical Society DOI: 10.1112/S002461070200371X POLYNOMIAL SOLUTIONS OF PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS J. MCLAUGHLIN

More information

A note on the k-narayana sequence

A note on the k-narayana sequence Annales Mathematicae et Informaticae 45 (2015) pp 91 105 http://amietfhu A note on the -Narayana sequence José L Ramírez a, Víctor F Sirvent b a Departamento de Matemáticas, Universidad Sergio Arboleda,

More information

18-29 mai 2015: Oujda (Maroc) École de recherche CIMPA-Oujda Théorie des Nombres et ses Applications. Continued fractions. Michel Waldschmidt

18-29 mai 2015: Oujda (Maroc) École de recherche CIMPA-Oujda Théorie des Nombres et ses Applications. Continued fractions. Michel Waldschmidt 18-29 mai 2015: Oujda (Maroc) École de recherche CIMPA-Oujda Théorie des Nombres et ses Applications. Continued fractions Michel Waldschmidt We first consider generalized continued fractions of the form

More information

Pell s equation. Michel Waldschmidt

Pell s equation. Michel Waldschmidt Faculté des Sciences et Techniques (FAST), Bamako, Mali École de recherche CIMPA Théorie des Nombres et Algorithmique Updated: December 7, 200 Pell s equation Michel Waldschmidt This text is available

More information

Test 2. Monday, November 12, 2018

Test 2. Monday, November 12, 2018 Test 2 Monday, November 12, 2018 Instructions. The only aids allowed are a hand-held calculator and one cheat sheet, i.e. an 8.5 11 sheet with information written on one side in your own handwriting. No

More information

Some applications of Number Theory. Number Theory in Science and Communication

Some applications of Number Theory. Number Theory in Science and Communication December 13, 2005 Nehru Science Center, Mumbai Mathematics in the real life: The Fibonacci Sequence and the Golden Number Michel Waldschmidt Université P. et M. Curie (Paris VI) Alliance Française http://www.math.jussieu.fr/~miw/

More information

Chapter 1 Why Quadratic Diophantine Equations?

Chapter 1 Why Quadratic Diophantine Equations? Chapter 1 Why Quadratic Diophantine Equations? In order to motivate the study of quadratic type equations, in this chapter we present several problems from various mathematical disciplines leading to such

More information

Fundamentals of Pure Mathematics - Problem Sheet

Fundamentals of Pure Mathematics - Problem Sheet Fundamentals of Pure Mathematics - Problem Sheet ( ) = Straightforward but illustrates a basic idea (*) = Harder Note: R, Z denote the real numbers, integers, etc. assumed to be real numbers. In questions

More information

Advanced Counting Techniques. Chapter 8

Advanced Counting Techniques. Chapter 8 Advanced Counting Techniques Chapter 8 Chapter Summary Applications of Recurrence Relations Solving Linear Recurrence Relations Homogeneous Recurrence Relations Nonhomogeneous Recurrence Relations Divide-and-Conquer

More information

CSCE 222 Discrete Structures for Computing. Dr. Hyunyoung Lee

CSCE 222 Discrete Structures for Computing. Dr. Hyunyoung Lee CSCE 222 Discrete Structures for Computing Sequences and Summations Dr. Hyunyoung Lee Based on slides by Andreas Klappenecker 1 Sequences 2 Sequences A sequence is a function from a subset of the set of

More information

#A91 INTEGERS 18 (2018) A GENERALIZED BINET FORMULA THAT COUNTS THE TILINGS OF A (2 N)-BOARD

#A91 INTEGERS 18 (2018) A GENERALIZED BINET FORMULA THAT COUNTS THE TILINGS OF A (2 N)-BOARD #A91 INTEGERS 18 (2018) A GENERALIZED BINET FORMULA THAT COUNTS THE TILINGS OF A (2 N)-BOARD Reza Kahkeshani 1 Department of Pure Mathematics, Faculty of Mathematical Sciences, University of Kashan, Kashan,

More information

MY PUTNAM PROBLEMS. log(1 + x) dx = π2

MY PUTNAM PROBLEMS. log(1 + x) dx = π2 MY PUTNAM PROBLEMS These are the problems I proposed when I was on the Putnam Problem Committee for the 984 86 Putnam Exams. Problems intended to be A or B (and therefore relatively easy) are marked accordingly.

More information

QUARTIC POWER SERIES IN F 3 ((T 1 )) WITH BOUNDED PARTIAL QUOTIENTS. Alain Lasjaunias

QUARTIC POWER SERIES IN F 3 ((T 1 )) WITH BOUNDED PARTIAL QUOTIENTS. Alain Lasjaunias QUARTIC POWER SERIES IN F 3 ((T 1 )) WITH BOUNDED PARTIAL QUOTIENTS Alain Lasjaunias 1991 Mathematics Subject Classification: 11J61, 11J70. 1. Introduction. We are concerned with diophantine approximation

More information

Chapter 8: Recursion. March 10, 2008

Chapter 8: Recursion. March 10, 2008 Chapter 8: Recursion March 10, 2008 Outline 1 8.1 Recursively Defined Sequences 2 8.2 Solving Recurrence Relations by Iteration 3 8.4 General Recursive Definitions Recursively Defined Sequences As mentioned

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

PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes

PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES Notes. x n+ = ax n has the general solution x n = x a n. 2. x n+ = x n + b has the general solution x n = x + (n )b. 3. x n+ = ax n + b (with a ) can be

More information

A talk given at the City Univ. of Hong Kong on April 14, ON HILBERT S TENTH PROBLEM AND RELATED TOPICS

A talk given at the City Univ. of Hong Kong on April 14, ON HILBERT S TENTH PROBLEM AND RELATED TOPICS A talk given at the City Univ. of Hong Kong on April 14, 000. ON HILBERT S TENTH PROBLEM AND RELATED TOPICS Zhi-Wei Sun Department of Mathematics Nanjing University Nanjing 10093 People s Republic of China

More information

CHARACTERIZING INTEGERS AMONG RATIONAL NUMBERS WITH A UNIVERSAL-EXISTENTIAL FORMULA

CHARACTERIZING INTEGERS AMONG RATIONAL NUMBERS WITH A UNIVERSAL-EXISTENTIAL FORMULA CHARACTERIZING INTEGERS AMONG RATIONAL NUMBERS WITH A UNIVERSAL-EXISTENTIAL FORMULA BJORN POONEN Abstract. We prove that Z in definable in Q by a formula with 2 universal quantifiers followed by 7 existential

More information

Advanced Counting Techniques

Advanced Counting Techniques . All rights reserved. Authorized only for instructor use in the classroom. No reproduction or further distribution permitted without the prior written consent of McGraw-Hill Education. Advanced Counting

More information

Defining Valuation Rings

Defining Valuation Rings East Carolina University, Greenville, North Carolina, USA June 8, 2018 Outline 1 What? Valuations and Valuation Rings Definability Questions in Number Theory 2 Why? Some Questions and Answers Becoming

More information

RECENT RESULTS ON LINEAR RECURRENCE SEQUENCES

RECENT RESULTS ON LINEAR RECURRENCE SEQUENCES RECENT RESULTS ON LINEAR RECURRENCE SEQUENCES Jan-Hendrik Evertse (Leiden) General Mathematics Colloquium, Delft May 19, 2005 1 INTRODUCTION A linear recurrence sequence U = {u n } n=0 (in C) is a sequence

More information

Christmas Calculated Colouring - C1

Christmas Calculated Colouring - C1 Christmas Calculated Colouring - C Tom Bennison December 20, 205 Introduction Each question identifies a region or regions on the picture Work out the answer and use the key to work out which colour to

More information

Summation of certain infinite Fibonacci related series

Summation of certain infinite Fibonacci related series arxiv:52.09033v (30 Dec 205) Summation of certain infinite Fibonacci related series Bakir Farhi Laboratoire de Mathématiques appliquées Faculté des Sciences Exactes Université de Bejaia 06000 Bejaia Algeria

More information

Notes on generating functions in automata theory

Notes on generating functions in automata theory Notes on generating functions in automata theory Benjamin Steinberg December 5, 2009 Contents Introduction: Calculus can count 2 Formal power series 5 3 Rational power series 9 3. Rational power series

More information

Generalization of Fibonacci sequence

Generalization of Fibonacci sequence Generalization of Fibonacci sequence Etienne Durand Julien Chartrand Maxime Bolduc February 18th 2013 Abstract After studying the fibonacci sequence, we found three interesting theorems. The first theorem

More information

ON THE POSSIBLE QUANTITIES OF FIBONACCI NUMBERS THAT OCCUR IN SOME TYPES OF INTERVALS

ON THE POSSIBLE QUANTITIES OF FIBONACCI NUMBERS THAT OCCUR IN SOME TYPES OF INTERVALS Acta Math. Univ. Comenianae Vol. LXXXVII, 2 (2018), pp. 291 299 291 ON THE POSSIBLE QUANTITIES OF FIBONACCI NUMBERS THAT OCCUR IN SOME TYPES OF INTERVALS B. FARHI Abstract. In this paper, we show that

More information

The plastic number and its generalized polynomial

The plastic number and its generalized polynomial PURE MATHEMATICS RESEARCH ARTICLE The plastic number and its generalized polynomial Vasileios Iliopoulos 1 * Received: 18 December 2014 Accepted: 19 February 201 Published: 20 March 201 *Corresponding

More information

RANK AND PERIOD OF PRIMES IN THE FIBONACCI SEQUENCE. A TRICHOTOMY

RANK AND PERIOD OF PRIMES IN THE FIBONACCI SEQUENCE. A TRICHOTOMY RANK AND PERIOD OF PRIMES IN THE FIBONACCI SEQUENCE. A TRICHOTOMY Christian Ballot Université de Caen, Caen 14032, France e-mail: ballot@math.unicaen.edu Michele Elia Politecnico di Torino, Torino 10129,

More information

SQUARE ROOTS OF 2x2 MATRICES 1. Sam Northshield SUNY-Plattsburgh

SQUARE ROOTS OF 2x2 MATRICES 1. Sam Northshield SUNY-Plattsburgh SQUARE ROOTS OF x MATRICES Sam Northshield SUNY-Plattsburgh INTRODUCTION A B What is the square root of a matrix such as? It is not, in general, A B C D C D This is easy to see since the upper left entry

More information

Public-key Cryptography: Theory and Practice

Public-key Cryptography: Theory and Practice Public-key Cryptography Theory and Practice Department of Computer Science and Engineering Indian Institute of Technology Kharagpur Chapter 2: Mathematical Concepts Divisibility Congruence Quadratic Residues

More information

Some Results Concerning Uniqueness of Triangle Sequences

Some Results Concerning Uniqueness of Triangle Sequences Some Results Concerning Uniqueness of Triangle Sequences T. Cheslack-Postava A. Diesl M. Lepinski A. Schuyler August 12 1999 Abstract In this paper we will begin by reviewing the triangle iteration. We

More information

Discrete Math, Second Problem Set (June 24)

Discrete Math, Second Problem Set (June 24) Discrete Math, Second Problem Set (June 24) REU 2003 Instructor: Laszlo Babai Scribe: D Jeremy Copeland 1 Number Theory Remark 11 For an arithmetic progression, a 0, a 1 = a 0 +d, a 2 = a 0 +2d, to have

More information

Algebra Homework, Edition 2 9 September 2010

Algebra Homework, Edition 2 9 September 2010 Algebra Homework, Edition 2 9 September 2010 Problem 6. (1) Let I and J be ideals of a commutative ring R with I + J = R. Prove that IJ = I J. (2) Let I, J, and K be ideals of a principal ideal domain.

More information

Pell s Equation Claire Larkin

Pell s Equation Claire Larkin Pell s Equation is a Diophantine equation in the form: Pell s Equation Claire Larkin The Equation x 2 dy 2 = where x and y are both integer solutions and n is a positive nonsquare integer. A diophantine

More information

Lecture Notes. Advanced Discrete Structures COT S

Lecture Notes. Advanced Discrete Structures COT S Lecture Notes Advanced Discrete Structures COT 4115.001 S15 2015-01-13 Recap Divisibility Prime Number Theorem Euclid s Lemma Fundamental Theorem of Arithmetic Euclidean Algorithm Basic Notions - Section

More information

Table of Contents. 2013, Pearson Education, Inc.

Table of Contents. 2013, Pearson Education, Inc. Table of Contents Chapter 1 What is Number Theory? 1 Chapter Pythagorean Triples 5 Chapter 3 Pythagorean Triples and the Unit Circle 11 Chapter 4 Sums of Higher Powers and Fermat s Last Theorem 16 Chapter

More information

MULTI-VARIABLE POLYNOMIAL SOLUTIONS TO PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS

MULTI-VARIABLE POLYNOMIAL SOLUTIONS TO PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS MULTI-VARIABLE POLYNOMIAL SOLUTIONS TO PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS J. MC LAUGHLIN Abstract. Solving Pell s equation is of relevance in finding fundamental units in real

More information

PELL S EQUATION NATIONAL INSTITUTE OF TECHNOLOGY ROURKELA, ODISHA

PELL S EQUATION NATIONAL INSTITUTE OF TECHNOLOGY ROURKELA, ODISHA PELL S EQUATION A Project Report Submitted by PANKAJ KUMAR SHARMA In partial fulfillment of the requirements For award of the degree Of MASTER OF SCIENCE IN MATHEMATICS UNDER GUIDANCE OF Prof GKPANDA DEPARTMENT

More information

POLYNOMIAL SOLUTIONS TO PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS arxiv: v1 [math.nt] 27 Dec 2018

POLYNOMIAL SOLUTIONS TO PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS arxiv: v1 [math.nt] 27 Dec 2018 POLYNOMIAL SOLUTIONS TO PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS arxiv:1812.10828v1 [math.nt] 27 Dec 2018 J. MC LAUGHLIN Abstract. Finding polynomial solutions to Pell s equation

More information

Cryptography. Number Theory with AN INTRODUCTION TO. James S. Kraft. Lawrence C. Washington. CRC Press

Cryptography. Number Theory with AN INTRODUCTION TO. James S. Kraft. Lawrence C. Washington. CRC Press AN INTRODUCTION TO Number Theory with Cryptography James S Kraft Gilman School Baltimore, Maryland, USA Lawrence C Washington University of Maryland College Park, Maryland, USA CRC Press Taylor & Francis

More information

Algebra Qualifying Exam August 2001 Do all 5 problems. 1. Let G be afinite group of order 504 = 23 32 7. a. Show that G cannot be isomorphic to a subgroup of the alternating group Alt 7. (5 points) b.

More information

Introduction to Arithmetic Geometry

Introduction to Arithmetic Geometry Introduction to Arithmetic Geometry 18.782 Andrew V. Sutherland September 5, 2013 What is arithmetic geometry? Arithmetic geometry applies the techniques of algebraic geometry to problems in number theory

More information

The real line: the equational compatibility problem and one-variable variants Hilbert s 10 th Problem

The real line: the equational compatibility problem and one-variable variants Hilbert s 10 th Problem The real line: the equational compatibility problem and one-variable variants Hilbert s 10 th Problem George F. McNulty Abstract. Walter Taylor proved recently that there is no algorithm for determining

More information

Just like the ring of Gaussian integers, the ring of Eisenstein integers is a Unique Factorization Domain.

Just like the ring of Gaussian integers, the ring of Eisenstein integers is a Unique Factorization Domain. Fermat s Infinite Descent PMATH 340 Assignment 6 (Due Monday April 3rd at noon). (0 marks) Use Femtat s method of infinite descent to prove that the Diophantine equation x 3 + y 3 = 4z 3 has no solutions

More information

586 Index. vertex, 369 disjoint, 236 pairwise, 272, 395 disjoint sets, 236 disjunction, 33, 36 distributive laws

586 Index. vertex, 369 disjoint, 236 pairwise, 272, 395 disjoint sets, 236 disjunction, 33, 36 distributive laws Index absolute value, 135 141 additive identity, 254 additive inverse, 254 aleph, 465 algebra of sets, 245, 278 antisymmetric relation, 387 arcsine function, 349 arithmetic sequence, 208 arrow diagram,

More information

are the q-versions of n, n! and . The falling factorial is (x) k = x(x 1)(x 2)... (x k + 1).

are the q-versions of n, n! and . The falling factorial is (x) k = x(x 1)(x 2)... (x k + 1). Lecture A jacques@ucsd.edu Notation: N, R, Z, F, C naturals, reals, integers, a field, complex numbers. p(n), S n,, b(n), s n, partition numbers, Stirling of the second ind, Bell numbers, Stirling of the

More information

CONGRUENCES FOR BERNOULLI - LUCAS SUMS

CONGRUENCES FOR BERNOULLI - LUCAS SUMS CONGRUENCES FOR BERNOULLI - LUCAS SUMS PAUL THOMAS YOUNG Abstract. We give strong congruences for sums of the form N BnVn+1 where Bn denotes the Bernoulli number and V n denotes a Lucas sequence of the

More information

On repdigits as product of consecutive Lucas numbers

On repdigits as product of consecutive Lucas numbers Notes on Number Theory and Discrete Mathematics Print ISSN 1310 5132, Online ISSN 2367 8275 Vol. 24, 2018, No. 3, 5 102 DOI: 10.7546/nntdm.2018.24.3.5-102 On repdigits as product of consecutive Lucas numbers

More information

MATH 402 : 2017 page 1 of 6

MATH 402 : 2017 page 1 of 6 ADMINISTRATION What? Math 40: Enumerative Combinatorics Who? Me: Professor Greg Smith You: students interested in combinatorics When and Where? lectures: slot 00 office hours: Tuesdays at :30 6:30 and

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

SOLVING THE PELL EQUATION VIA RÉDEI RATIONAL FUNCTIONS

SOLVING THE PELL EQUATION VIA RÉDEI RATIONAL FUNCTIONS STEFANO BARBERO, UMBERTO CERRUTI, AND NADIR MURRU Abstract. In this paper, we define a new product over R, which allows us to obtain a group isomorphic to R with the usual product. This operation unexpectedly

More information

Sequences and Summations. CS/APMA 202 Rosen section 3.2 Aaron Bloomfield

Sequences and Summations. CS/APMA 202 Rosen section 3.2 Aaron Bloomfield Sequences and Summations CS/APMA 202 Rosen section 3.2 Aaron Bloomfield 1 Definitions Sequence: an ordered list of elements Like a set, but: Elements can be duplicated Elements are ordered 2 Sequences

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

#A5 INTEGERS 18A (2018) EXPLICIT EXAMPLES OF p-adic NUMBERS WITH PRESCRIBED IRRATIONALITY EXPONENT

#A5 INTEGERS 18A (2018) EXPLICIT EXAMPLES OF p-adic NUMBERS WITH PRESCRIBED IRRATIONALITY EXPONENT #A5 INTEGERS 8A (208) EXPLICIT EXAMPLES OF p-adic NUMBERS WITH PRESCRIBED IRRATIONALITY EXPONENT Yann Bugeaud IRMA, UMR 750, Université de Strasbourg et CNRS, Strasbourg, France bugeaud@math.unistra.fr

More information

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2. APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product

More information

198 VOLUME 46/47, NUMBER 3

198 VOLUME 46/47, NUMBER 3 LAWRENCE SOMER Abstract. Rotkiewicz has shown that there exist Fibonacci pseudoprimes having the forms p(p + 2), p(2p 1), and p(2p + 3), where all the terms in the products are odd primes. Assuming Dickson

More information

Introduction to Number Theory

Introduction to Number Theory INTRODUCTION Definition: Natural Numbers, Integers Natural numbers: N={0,1,, }. Integers: Z={0,±1,±, }. Definition: Divisor If a Z can be writeen as a=bc where b, c Z, then we say a is divisible by b or,

More information

Egyptian Fraction. Massoud Malek

Egyptian Fraction. Massoud Malek Egyptian Fraction Massoud Malek Throughout history, different civilizations have had different ways of representing numbers. Some of these systems seem strange or complicated from our perspective. The

More information

QUALIFYING EXAM IN ALGEBRA August 2011

QUALIFYING EXAM IN ALGEBRA August 2011 QUALIFYING EXAM IN ALGEBRA August 2011 1. There are 18 problems on the exam. Work and turn in 10 problems, in the following categories. I. Linear Algebra 1 problem II. Group Theory 3 problems III. Ring

More information

EIGENVALUES AND EIGENVECTORS 3

EIGENVALUES AND EIGENVECTORS 3 EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices

More information

Lecture 8: Phyllotaxis, the golden ratio and the Fibonacci sequence

Lecture 8: Phyllotaxis, the golden ratio and the Fibonacci sequence The 77 th Compton lecture series Frustrating geometry: Geometry and incompatibility shaping the physical world Lecture 8: Phyllotaxis, the golden ratio and the Fibonacci sequence Efi Efrati Simons postdoctoral

More information

Coding Theory ( Mathematical Background I)

Coding Theory ( Mathematical Background I) N.L.Manev, Lectures on Coding Theory (Maths I) p. 1/18 Coding Theory ( Mathematical Background I) Lector: Nikolai L. Manev Institute of Mathematics and Informatics, Sofia, Bulgaria N.L.Manev, Lectures

More information

TEST CODE: MMA (Objective type) 2015 SYLLABUS

TEST CODE: MMA (Objective type) 2015 SYLLABUS TEST CODE: MMA (Objective type) 2015 SYLLABUS Analytical Reasoning Algebra Arithmetic, geometric and harmonic progression. Continued fractions. Elementary combinatorics: Permutations and combinations,

More information

Europe Starts to Wake up: Leonardo of Pisa

Europe Starts to Wake up: Leonardo of Pisa Europe Starts to Wake up: Leonardo of Pisa Leonardo of Pisa, also known as Fibbonaci (from filius Bonaccia, son of Bonnaccio ) was the greatest mathematician of the middle ages. He lived from 75 to 50,

More information

THE DENOMINATORS OF POWER SUMS OF ARITHMETIC PROGRESSIONS. Bernd C. Kellner Göppert Weg 5, Göttingen, Germany

THE DENOMINATORS OF POWER SUMS OF ARITHMETIC PROGRESSIONS. Bernd C. Kellner Göppert Weg 5, Göttingen, Germany #A95 INTEGERS 18 (2018) THE DENOMINATORS OF POWER SUMS OF ARITHMETIC PROGRESSIONS Bernd C. Kellner Göppert Weg 5, 37077 Göttingen, Germany b@bernoulli.org Jonathan Sondow 209 West 97th Street, New Yor,

More information

ON THE DECIMAL EXPANSION OF ALGEBRAIC NUMBERS

ON THE DECIMAL EXPANSION OF ALGEBRAIC NUMBERS Fizikos ir matematikos fakulteto Seminaro darbai, Šiaulių universitetas, 8, 2005, 5 13 ON THE DECIMAL EXPANSION OF ALGEBRAIC NUMBERS Boris ADAMCZEWSKI 1, Yann BUGEAUD 2 1 CNRS, Institut Camille Jordan,

More information

Some remarks on diophantine equations and diophantine approximation

Some remarks on diophantine equations and diophantine approximation Some remarks on diophantine equations and diophantine approximation Claude LEVESQUE and Michel WALDSCHMIDT Dedicated to professor Hà Huy Khoái. Abstract. We first recall the connection, going back to A.

More information

Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra

Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra D. R. Wilkins Contents 3 Topics in Commutative Algebra 2 3.1 Rings and Fields......................... 2 3.2 Ideals...............................

More information

Chakravala - a modern Indian method. B.Sury

Chakravala - a modern Indian method. B.Sury Chakravala - a modern Indian method BSury Indian Statistical Institute Bangalore, India sury@isibangacin IISER Pune, India Lecture on October 18, 2010 1 Weil Unveiled What would have been Fermat s astonishment

More information

= 1 2x. x 2 a ) 0 (mod p n ), (x 2 + 2a + a2. x a ) 2

= 1 2x. x 2 a ) 0 (mod p n ), (x 2 + 2a + a2. x a ) 2 8. p-adic numbers 8.1. Motivation: Solving x 2 a (mod p n ). Take an odd prime p, and ( an) integer a coprime to p. Then, as we know, x 2 a (mod p) has a solution x Z iff = 1. In this case we can suppose

More information

NOTES ON DIOPHANTINE APPROXIMATION

NOTES ON DIOPHANTINE APPROXIMATION NOTES ON DIOPHANTINE APPROXIMATION Jan-Hendrik Evertse January 29, 200 9 p-adic Numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics

More information

GAPS IN BINARY EXPANSIONS OF SOME ARITHMETIC FUNCTIONS, AND THE IRRATIONALITY OF THE EULER CONSTANT

GAPS IN BINARY EXPANSIONS OF SOME ARITHMETIC FUNCTIONS, AND THE IRRATIONALITY OF THE EULER CONSTANT Journal of Prime Research in Mathematics Vol. 8 202, 28-35 GAPS IN BINARY EXPANSIONS OF SOME ARITHMETIC FUNCTIONS, AND THE IRRATIONALITY OF THE EULER CONSTANT JORGE JIMÉNEZ URROZ, FLORIAN LUCA 2, MICHEL

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

University of Toronto Faculty of Arts and Science Solutions to Final Examination, April 2017 MAT246H1S - Concepts in Abstract Mathematics

University of Toronto Faculty of Arts and Science Solutions to Final Examination, April 2017 MAT246H1S - Concepts in Abstract Mathematics University of Toronto Faculty of Arts and Science Solutions to Final Examination, April 2017 MAT246H1S - Concepts in Abstract Mathematics Examiners: D. Burbulla, P. Glynn-Adey, S. Homayouni Time: 7-10

More information

BALANCING NUMBERS : SOME IDENTITIES KABERI PARIDA. Master of Science in Mathematics. Dr. GOPAL KRISHNA PANDA

BALANCING NUMBERS : SOME IDENTITIES KABERI PARIDA. Master of Science in Mathematics. Dr. GOPAL KRISHNA PANDA BALANCING NUMBERS : SOME IDENTITIES A report submitted by KABERI PARIDA Roll No: 1MA073 for the partial fulfilment for the award of the degree of Master of Science in Mathematics under the supervision

More information

Chapter 8. P-adic numbers. 8.1 Absolute values

Chapter 8. P-adic numbers. 8.1 Absolute values Chapter 8 P-adic numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics 58, Springer Verlag 1984, corrected 2nd printing 1996, Chap.

More information

MATH 324 Summer 2011 Elementary Number Theory. Notes on Mathematical Induction. Recall the following axiom for the set of integers.

MATH 324 Summer 2011 Elementary Number Theory. Notes on Mathematical Induction. Recall the following axiom for the set of integers. MATH 4 Summer 011 Elementary Number Theory Notes on Mathematical Induction Principle of Mathematical Induction Recall the following axiom for the set of integers. Well-Ordering Axiom for the Integers If

More information

August 23, 2017 Let us measure everything that is measurable, and make measurable everything that is not yet so. Galileo Galilei. 1.

August 23, 2017 Let us measure everything that is measurable, and make measurable everything that is not yet so. Galileo Galilei. 1. August 23, 2017 Let us measure everything that is measurable, and make measurable everything that is not yet so. Galileo Galilei 1. Vector spaces 1.1. Notations. x S denotes the fact that the element x

More information

18.S34 linear algebra problems (2007)

18.S34 linear algebra problems (2007) 18.S34 linear algebra problems (2007) Useful ideas for evaluating determinants 1. Row reduction, expanding by minors, or combinations thereof; sometimes these are useful in combination with an induction

More information

A Short Proof of the Transcendence of Thue-Morse Continued Fractions

A Short Proof of the Transcendence of Thue-Morse Continued Fractions A Short Proof of the Transcendence of Thue-Morse Continued Fractions Boris ADAMCZEWSKI and Yann BUGEAUD The Thue-Morse sequence t = (t n ) n 0 on the alphabet {a, b} is defined as follows: t n = a (respectively,

More information

Analytic Number Theory Solutions

Analytic Number Theory Solutions Analytic Number Theory Solutions Sean Li Cornell University sxl6@cornell.edu Jan. 03 Introduction This document is a work-in-progress solution manual for Tom Apostol s Introduction to Analytic Number Theory.

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 1998 Comments to the author at krm@maths.uq.edu.au Contents 1 LINEAR EQUATIONS

More information

Narayana Sequences for Cryptographic Applications

Narayana Sequences for Cryptographic Applications Narayana Sequences for Cryptographic Applications Krishnamurthy Kirthi Abstract This paper investigates the randomness and cryptographic properties of the Narayana series modulo p, where p is a prime number.

More information

Rings and modular arithmetic

Rings and modular arithmetic Chapter 8 Rings and modular arithmetic So far, we have been working with just one operation at a time. But standard number systems, such as Z, have two operations + and which interact. It is useful to

More information

Continued Fractions New and Old Results

Continued Fractions New and Old Results Continued Fractions New and Old Results Jeffrey Shallit School of Computer Science University of Waterloo Waterloo, Ontario N2L 3G1 Canada shallit@cs.uwaterloo.ca https://www.cs.uwaterloo.ca/~shallit Joint

More information

2.5 The Fundamental Theorem of Algebra.

2.5 The Fundamental Theorem of Algebra. 2.5. THE FUNDAMENTAL THEOREM OF ALGEBRA. 79 2.5 The Fundamental Theorem of Algebra. We ve seen formulas for the (complex) roots of quadratic, cubic and quartic polynomials. It is then reasonable to ask:

More information

Mathematical Olympiad Training Polynomials

Mathematical Olympiad Training Polynomials Mathematical Olympiad Training Polynomials Definition A polynomial over a ring R(Z, Q, R, C) in x is an expression of the form p(x) = a n x n + a n 1 x n 1 + + a 1 x + a 0, a i R, for 0 i n. If a n 0,

More information

THE UNIT GROUP OF A REAL QUADRATIC FIELD

THE UNIT GROUP OF A REAL QUADRATIC FIELD THE UNIT GROUP OF A REAL QUADRATIC FIELD While the unit group of an imaginary quadratic field is very simple the unit group of a real quadratic field has nontrivial structure Its study involves some geometry

More information

Chapter 3. Rings. The basic commutative rings in mathematics are the integers Z, the. Examples

Chapter 3. Rings. The basic commutative rings in mathematics are the integers Z, the. Examples Chapter 3 Rings Rings are additive abelian groups with a second operation called multiplication. The connection between the two operations is provided by the distributive law. Assuming the results of Chapter

More information

Course 2316 Sample Paper 1

Course 2316 Sample Paper 1 Course 2316 Sample Paper 1 Timothy Murphy April 19, 2015 Attempt 5 questions. All carry the same mark. 1. State and prove the Fundamental Theorem of Arithmetic (for N). Prove that there are an infinity

More information