Reachability Problems for Continuous Linear Dynamical Systems

Size: px
Start display at page:

Download "Reachability Problems for Continuous Linear Dynamical Systems"

Transcription

1 Reachability Problems for Continuous Linear Dynamical Systems James Worrell Department of Computer Science, Oxford University (Joint work with Ventsislav Chonev and Joël Ouaknine) FICS 2015 September 12th, 2015

2 Reachability for Continuous-Time Markov Chains Stagnant market Bull market 0.02 Bear market 0.3

3 Reachability for Continuous-Time Markov Chains Stagnant market Bull market 0.02 Bear market 0.3 Distribution P(t) at time t satisfies P (t) = P(t)Q, where Q = is the rate matrix.

4 Reachability for Continuous-Time Markov Chains Stagnant market Bull market 0.02 Bear market 0.3 Distribution P(t) at time t satisfies P (t) = P(t)Q, where Q = is the rate matrix. Is it ever more likely to be a Bear market than a Bull market? t (P(t) Bear P(t) Bull )

5 Reachability - Some Basic Insights Is it ever more likely to be a Bear market than a Bull market? t (P(t) Bear P(t) Bull )

6 Reachability - Some Basic Insights Is it ever more likely to be a Bear market than a Bull market? t (P(t) Bear P(t) Bull ) Reduce to the time-bounded case by computing the stationary distribution: π = (0.885, 0.071, 0.044)

7 Reachability - Some Basic Insights Is it ever more likely to be a Bear market than a Bull market? t (P(t) Bear P(t) Bull ) Reduce to the time-bounded case by computing the stationary distribution: π = (0.885, 0.071, 0.044) Reduction does not work if π is on boundary of the target set.

8 Reachability - Some Basic Insights Is it ever more likely to be a Bear market than a Bull market? t (P(t) Bear P(t) Bull ) Reduce to the time-bounded case by computing the stationary distribution: π = (0.885, 0.071, 0.044) Reduction does not work if π is on boundary of the target set. Semantic shift consider whether the closure of the orbit {P(t) : t 0} meets the closure of the target set.

9 Cyber-Physical Systems To analyze a cyber-physical system, such as a pacemaker, we need to consider the discrete software controller interacting with the physical world, which is typically modelled by differential equations Rajeev Alur (CACM, 2013)

10 Hybrid Automata: Various Continuous Dynamics Hybrid automaton = states + variables x R k

11 Hybrid Automata: Various Continuous Dynamics Hybrid automaton = states + variables x R k ẋ = 1 timed automata

12 Hybrid Automata: Various Continuous Dynamics Hybrid automaton = states + variables x R k ẋ = 1 timed automata ẋ = c rectangular hybrid automata

13 Hybrid Automata: Various Continuous Dynamics Hybrid automaton = states + variables x R k ẋ = 1 timed automata ẋ = c rectangular hybrid automata ẋ = Ax linear hybrid automata

14 Hybrid Automata: Various Continuous Dynamics Hybrid automaton = states + variables x R k ẋ = 1 timed automata ẋ = c rectangular hybrid automata ẋ = Ax linear hybrid automata...

15 Hybrid Automata: Various Continuous Dynamics Hybrid automaton = states + variables x R k ẋ = 1 timed automata ẋ = c rectangular hybrid automata ẋ = Ax linear hybrid automata... o-minimal flows + strong resets reachability decidable

16 Hybrid Automata: Various Continuous Dynamics Hybrid automaton = states + variables x R k ẋ = 1 timed automata ẋ = c rectangular hybrid automata ẋ = Ax linear hybrid automata... o-minimal flows + strong resets reachability decidable Is this location a trap? x := 2 y := 4 ẋ = 3x y ẏ = x 5y x 10 y 2?

17 Reachability for Continuous Linear Dynamical Systems Is this location a trap? x := 2 y := 4 ẋ = 3x y ẏ = x 5y x 10 y 2? Stagnant market Is ever more likely to be a Bear Bull market Bear market market than a Bull market: t (P(t) Bear P(t) Bull )?

18 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax

19 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0)

20 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0)

21 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0)

22 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0)

23 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0)

24 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0)

25 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0) u

26 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0) f (t) = u T x(t)

27 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0) f (t) = u T x(t) f (k) (t) + a k 1 f (k 1) (t) a 1 f (t) + a 0 f (t) = 0

28 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0) f (t) = u T x(t) f (k) (t) + a k 1 f (k 1) (t) a 1 f (t) + a 0 f (t) = 0 m f (t) = P j (t)e λ j t j=1

29 Reachability for Continuous Linear Dynamical Systems x : R 0 R k ẋ = Ax x(t) = exp(at)x(0) f (t) = u T x(t) f (k) (t) + a k 1 f (k 1) (t) a 1 f (t) + a 0 f (t) = 0 m f (t) = P j (t)e λ j t j=1 Note the λ j are complex in general.

30 Reachability for Continuous Linear Dynamical Systems Let f : R 0 R be given as above, with all coefficients algebraic.

31 Reachability for Continuous Linear Dynamical Systems Let f : R 0 R be given as above, with all coefficients algebraic. BOUNDED-ZERO Problem Instance: f and bounded interval [a, b] Question: Is there t [a, b] such that f (t) = 0?

32 Reachability for Continuous Linear Dynamical Systems Let f : R 0 R be given as above, with all coefficients algebraic. BOUNDED-ZERO Problem Instance: f and bounded interval [a, b] Question: Is there t [a, b] such that f (t) = 0? ZERO Problem Instance: f Question: Is there t R 0 such that f (t) = 0?

33 Reachability for Continuous Linear Dynamical Systems Let f : R 0 R be given as above, with all coefficients algebraic. BOUNDED-ZERO Problem Instance: f and bounded interval [a, b] Question: Is there t [a, b] such that f (t) = 0? ZERO Problem Instance: f Question: Is there t R 0 such that f (t) = 0? Decidability open! [Bell, Delvenne, Jungers, Blondel 2010]

34 Related Work A lot of work since 1920s on the zeros of exponential polynomials f (z) = m P j (z)e λ j z j=1 (Polya, Ritt, Tamarkin, Kac, Voorhoeve, van der Poorten,... ) but mostly on distribution of complex zeros.

35 Related Work A lot of work since 1920s on the zeros of exponential polynomials f (z) = m P j (z)e λ j z j=1 (Polya, Ritt, Tamarkin, Kac, Voorhoeve, van der Poorten,... ) but mostly on distribution of complex zeros. CONTINUOUS-ORBIT Problem The problem of whether the trajectory x(t) = e At x(0) reaches a given target point was shown to be decidable by Hainry (2008) and in PTIME by Chen, Han and Yu (2015).

36 Reachability for Continuous Linear Dynamical Systems Theorem (Bell, Delvenne, Jungers, Blondel 2010) In dimension 2, BOUNDED-ZERO and ZERO are decidable.

37 Reachability for Continuous Linear Dynamical Systems Theorem (Bell, Delvenne, Jungers, Blondel 2010) In dimension 2, BOUNDED-ZERO and ZERO are decidable. Theorem (Chonev, Ouaknine, W. 2015) In dimension 3, BOUNDED-ZERO and ZERO are decidable.

38 Reachability for Continuous Linear Dynamical Systems Theorem (Bell, Delvenne, Jungers, Blondel 2010) In dimension 2, BOUNDED-ZERO and ZERO are decidable. Theorem (Chonev, Ouaknine, W. 2015) In dimension 3, BOUNDED-ZERO and ZERO are decidable. Theorem (Chonev, Ouaknine, W. 2015) Assuming Schanuel s Conjecture, BOUNDED-ZERO is decidable in all dimensions.

39 Reachability for Continuous Linear Dynamical Systems Theorem (Bell, Delvenne, Jungers, Blondel 2010) In dimension 2, BOUNDED-ZERO and ZERO are decidable. Theorem (Chonev, Ouaknine, W. 2015) In dimension 3, BOUNDED-ZERO and ZERO are decidable. Theorem (Chonev, Ouaknine, W. 2015) Assuming Schanuel s Conjecture, BOUNDED-ZERO is decidable in all dimensions. It turns out that this result (in fact, a powerful generalisation of it) had already been discovered (but never published) in the early 1990s by Macintyre and Wilkie! [Angus Macintyre, personal communication, July 2015]

40 Reachability for Continuous Linear Dynamical Systems Theorem (Chonev, Ouaknine, W. 2015) In dimension 8 or less, ZERO reduces to BOUNDED-ZERO.

41 Reachability for Continuous Linear Dynamical Systems Theorem (Chonev, Ouaknine, W. 2015) In dimension 8 or less, ZERO reduces to BOUNDED-ZERO. Theorem (Chonev, Ouaknine, W. 2015) In dimension 9 (and above), decidability of ZERO would entail major breakthroughs in Diophantine approximation the Diophantine approximation type of α would be computable to within arbitrary precision.

42 Schanuel s Conjecture Theorem (Lindemann-Weierstrass) If a 1,..., a n are algebraic numbers linearly independent over Q, then e a 1,..., e an are algebraically independent.

43 Schanuel s Conjecture Theorem (Lindemann-Weierstrass) If a 1,..., a n are algebraic numbers linearly independent over Q, then e a 1,..., e an are algebraically independent. Schanuel s Conjecture If z 1,..., z n are complex numbers linearly independent over Q then some n-element subset of {z 1,..., z n, e z 1,..., e zn } is algebraically independent.

44 Schanuel s Conjecture Theorem (Lindemann-Weierstrass) If a 1,..., a n are algebraic numbers linearly independent over Q, then e a 1,..., e an are algebraically independent. Schanuel s Conjecture If z 1,..., z n are complex numbers linearly independent over Q then some n-element subset of {z 1,..., z n, e z 1,..., e zn } is algebraically independent. Example By Schanuel s conjecture some two-element subset of {1, πi, e 1, e πi } is algebraically independent.

45 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1

46 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) t

47 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a b t

48 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a b t

49 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a t* b t non-trivial zero t transcendental

50 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a t* b t non-trivial zero t transcendental

51 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a t* b t non-trivial zero t transcendental

52 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a b t

53 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a b t

54 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a t* b t

55 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a t* b t Can this situation arise?

56 The BOUNDED-ZERO Problem m Real-valued exponential polynomial f (t) = P j (t)e λ j t j=1 f(t) a t* b t Easily! For example, f (t) = 2 + e it + e it.

57 Laurent Polynomials and Factorisation Example Write f (t) = 2 + e it + e it in the form f (t) = P(e it ) for the Laurent polynomial P(z) = 2 + z + z 1.

58 Laurent Polynomials and Factorisation Example Write f (t) = 2 + e it + e it in the form f (t) = P(e it ) for the Laurent polynomial P(z) = 2 + z + z 1. Factorisation P(z) = (1 + z)(1 + z 1 ) induces a factorisation f (t) = (1 + e it ) (1 e it ) }{{}}{{} f 1 (t) f 2 (t)

59 Laurent Polynomials and Factorisation Example Write f (t) = 2 + e it + e it in the form f (t) = P(e it ) for the Laurent polynomial P(z) = 2 + z + z 1. Factorisation P(z) = (1 + z)(1 + z 1 ) induces a factorisation f (t) = (1 + e it ) (1 e it ) }{{}}{{} f 1 (t) f 2 (t) Common zeros of f 1 and f 2 are tangential zeros of f

60 Laurent Polynomials and Factorisation Example Write f (t) = 2 + e it + e it in the form f (t) = P(e it ) for the Laurent polynomial P(z) = 2 + z + z 1. Factorisation P(z) = (1 + z)(1 + z 1 ) induces a factorisation f (t) = (1 + e it ) (1 e it ) }{{}}{{} f 1 (t) f 2 (t) Common zeros of f 1 and f 2 are tangential zeros of f Idea: factorise f.

61 Laurent Polynomials and Factorisation Example Write f (t) = 2 + e it + e it in the form f (t) = P(e it ) for the Laurent polynomial P(z) = 2 + z + z 1. Factorisation P(z) = (1 + z)(1 + z 1 ) induces a factorisation f (t) = (1 + e it ) (1 e it ) }{{}}{{} f 1 (t) f 2 (t) Common zeros of f 1 and f 2 are tangential zeros of f Idea: factorise f. Noting that factors may be complex-valued!

62 The Real Case Any exponential polynomial f (t) can be written with f (t) = P(t, e a 1t,..., e amt ) P C[x, x ±1 1,..., x ±1 m ] and {a 1,..., a m } a set of real and imaginary algebraic numbers that is linearly independent over Q.

63 The Real Case Any exponential polynomial f (t) can be written with f (t) = P(t, e a 1t,..., e amt ) P C[x, x ±1 1,..., x ±1 m ] and {a 1,..., a m } a set of real and imaginary algebraic numbers that is linearly independent over Q. Lemma Assuming Schanuel s conjecture, if f is real valued and P is irreducible then f has no tangential zeros.

64 The Real Case Any exponential polynomial f (t) can be written with f (t) = P(t, e a 1t,..., e amt ) P C[x, x ±1 1,..., x ±1 m ] and {a 1,..., a m } a set of real and imaginary algebraic numbers that is linearly independent over Q. Lemma Assuming Schanuel s conjecture, if f is real valued and P is irreducible then f has no tangential zeros. Complex case requires some new ideas...

65 The Unbounded Case there are known unknowns; that is to say we know there are some things we do not know.

66 Continued Fractions Finite continued fractions: [3, 7, 15, 1, 292] =

67 Continued Fractions Finite continued fractions: [3, 7, 15, 1, 292] = =

68 Continued Fractions Finite continued fractions: [3, 7, 15, 1, 292] = 3 + Infinite continued fractions: = [a 0, a 1, a 2, a 3,...] = a a a 2 + a 3 +

69 Real Algebraic Numbers Theorem The continued fraction expansion of a real quadratic irrational number is periodic.

70 Real Algebraic Numbers Theorem The continued fraction expansion of a real quadratic irrational number is periodic. 389 = [19, 1, 2, 1, 1, 1, 1, 2, 1, 38, 1, 2, 1, 1, 1, 1, 2, 1, 38,...]

71 Real Algebraic Numbers Theorem The continued fraction expansion of a real quadratic irrational number is periodic. 389 = [19, 1, 2, 1, 1, 1, 1, 2, 1, 38, 1, 2, 1, 1, 1, 1, 2, 1, 38,...] What about numbers of degree 3?

72 Real Algebraic Numbers Theorem The continued fraction expansion of a real quadratic irrational number is periodic. 389 = [19, 1, 2, 1, 1, 1, 1, 2, 1, 38, 1, 2, 1, 1, 1, 1, 2, 1, 38,...] What about numbers of degree 3? 3 2 = [1, 3, 1, 5, 1, 1, 4, 1, 1, 8, 1, 14, 1, 10, 2, 1, 4, 12, 2, 3, 2, 1 3, 4, 1, 1, 2, 14, 3, 12, 1, 15, 3, 1, 4, 534, 1, 1, 5, 1, 1,...]

73 Real Algebraic Numbers Theorem The continued fraction expansion of a real quadratic irrational number is periodic. 389 = [19, 1, 2, 1, 1, 1, 1, 2, 1, 38, 1, 2, 1, 1, 1, 1, 2, 1, 38,...] What about numbers of degree 3? 3 2 = [1, 3, 1, 5, 1, 1, 4, 1, 1, 8, 1, 14, 1, 10, 2, 1, 4, 12, 2, 3, 2, 1 3, 4, 1, 1, 2, 14, 3, 12, 1, 15, 3, 1, 4, 534, 1, 1, 5, 1, 1,...] Lang and Trotter: no significant departure from random behaviour

74 An Open Problem [... ] no continued fraction development of an algebraic number of higher degree than the second is known. It is not even known if such a development has bounded elements. A. Khinchin

75 An Open Problem [... ] no continued fraction development of an algebraic number of higher degree than the second is known. It is not even known if such a development has bounded elements. A. Khinchin Is there an algebraic number of degree higher than two whose simple continued fraction has unbounded partial quotients? Does every such number have unbounded partial quotients? R. K. Guy, 2004

76 A Mathematical Obstacle at Dimension 9 Given x = [a 0, a 1, a 2,...], define S(x) = sup n N a n.

77 A Mathematical Obstacle at Dimension 9 Given x = [a 0, a 1, a 2,...], define S(x) = sup n N a n. Theorem (Chonev, Ouaknine, W., 2015) If the ZERO PROBLEM is decidable at dimension 9 then {x R A : S(x) < } is recursively enumerable.

78 A Mathematical Obstacle at Dimension 9 Given x = [a 0, a 1, a 2,...], define S(x) = sup n N a n. Theorem (Chonev, Ouaknine, W., 2015) If the ZERO PROBLEM is decidable at dimension 9 then {x R A : S(x) < } is recursively enumerable. Remark Perhaps this set is recursive it may even be or R A. However proving recursive enumerability would be a significant achievement.

79 Diophantine Approximation How well can one approximate a real number x with rationals? x m n

80 Diophantine Approximation How well can one approximate a real number x with rationals? x m n Theorem (Dirichlet 1842) There are infinitely many integers m, n such that x m < 1 n n 2.

81 Diophantine Approximation How well can one approximate a real number x with rationals? x m n Theorem (Dirichlet 1842) There are infinitely many integers m, n such that x m < 1 n n 2. S(x) < if and only if there exists ε > 0 such that x m < ε n n 2 has no solutions.

82 Diophantine Approximation How well can one approximate a real number x with rationals? x m n Theorem (Dirichlet 1842) There are infinitely many integers m, n such that x m < 1 n n 2. S(x) < if and only if there exists ε > 0 such that x m < ε n n 2 has no solutions. Relate this to the existence of zeros of order-9 exponential polynomial f (t) with terms e ixt and e it.

83 The ZERO Problem ZERO Problem Instance: f Question: Is there t R 0 such that f (t) = 0?

84 The ZERO Problem ZERO Problem Instance: f Question: Is there t R 0 such that f (t) = 0? Theorem In dimension 8 or less, ZERO reduces to BOUNDED-ZERO.

85 The ZERO Problem ZERO Problem Instance: f Question: Is there t R 0 such that f (t) = 0? Theorem In dimension 8 or less, ZERO reduces to BOUNDED-ZERO. In the limit f is either never zero or infinitely often zero, and we decide which is the case.

86 The ZERO Problem ZERO Problem Instance: f Question: Is there t R 0 such that f (t) = 0? Theorem In dimension 8 or less, ZERO reduces to BOUNDED-ZERO. In the limit f is either never zero or infinitely often zero, and we decide which is the case. Diophantine-approximation bounds play a key role in the proof.

87 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t

88 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

89 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

90 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

91 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

92 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

93 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

94 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

95 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

96 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

97 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

98 Illustrative Example Consider the exponential polynomial f (t) = cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

99 Illustrative Example Consider the exponential polynomial f (t) = 2 + cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2

100 Illustrative Example Consider the exponential polynomial f (t) = 2 + cos(t + ϕ 1 ) + cos( 2t + ϕ 2 ) e t Kronecker s Theorem: (t, 2t) mod 2π is dense in [0, 2π] 2 Baker s Theorem: ( t, ) 2t (π ϕ 1, π ϕ 2 ) 1 poly(t)

101 Conclusion and Perspectives

102 The Discrete Case A linear recurrence sequence is a sequence u 0, u 1, u 2,... of integers such that there exist constants a 1,..., a k, such that for all n 0. u n+k = a 1 u n+k 1 + a 2 u n+k a k u n

103 The Discrete Case A linear recurrence sequence is a sequence u 0, u 1, u 2,... of integers such that there exist constants a 1,..., a k, such that for all n 0. u n+k = a 1 u n+k 1 + a 2 u n+k a k u n Theorem (Skolem 1934; Mahler 1935, 1956; Lech 1953) The set of zeros of a linear recurrence sequence is semi-linear: {n : u n = 0} = F A 1... A l where F is finite and each A i is a full arithmetic progression.

104 The Discrete Case A linear recurrence sequence is a sequence u 0, u 1, u 2,... of integers such that there exist constants a 1,..., a k, such that for all n 0. u n+k = a 1 u n+k 1 + a 2 u n+k a k u n Theorem (Skolem 1934; Mahler 1935, 1956; Lech 1953) The set of zeros of a linear recurrence sequence is semi-linear: {n : u n = 0} = F A 1... A l where F is finite and each A i is a full arithmetic progression. Theorem (Berstel and Mignotte 1976) In Skolem-Mahler-Lech, the infinite part (arithmetic progressions A 1,..., A l ) is fully constructive.

105 The Skolem Problem Skolem Problem Does n such that u n = 0?

106 The Skolem Problem Skolem Problem Does n such that u n = 0? It is faintly outrageous that this problem is still open; it is saying that we do not know how to decide the Halting Problem even for linear automata! Terence Tao

107 The Skolem Problem Skolem Problem Does n such that u n = 0? It is faintly outrageous that this problem is still open; it is saying that we do not know how to decide the Halting Problem even for linear automata! Terence Tao... a mathematical embarrassment... Richard Lipton

108 Wrapping Things Up Continuous Skolem Problem Does t such that f (t) = 0?

109 Wrapping Things Up Continuous Skolem Problem Does t such that f (t) = 0? Even the bounded problem is hard (apparently).

110 Wrapping Things Up Continuous Skolem Problem Does t such that f (t) = 0? Even the bounded problem is hard (apparently). Formidable mathematical obstacle at dimension 9 in the unbounded case.

111 Wrapping Things Up Continuous Skolem Problem Does t such that f (t) = 0? Even the bounded problem is hard (apparently). Formidable mathematical obstacle at dimension 9 in the unbounded case. The infinite-zeros problem is also hard.

112 Wrapping Things Up Continuous Skolem Problem Does t such that f (t) = 0? Even the bounded problem is hard (apparently). Formidable mathematical obstacle at dimension 9 in the unbounded case. The infinite-zeros problem is also hard. Diophantine-approximation techniques unavoidable.

Reachability Problems for Continuous Linear Dynamical Systems

Reachability Problems for Continuous Linear Dynamical Systems Reachability Problems for Continuous Linear Dynamical Systems James Worrell Department of Computer Science, Oxford University (Joint work with Ventsislav Chonev and Joël Ouaknine) CONCUR September 2nd

More information

Decision Problems for Linear Recurrence Sequences

Decision Problems for Linear Recurrence Sequences Decision Problems for Linear Recurrence Sequences Joël Ouaknine Department of Computer Science, Oxford University (Joint work with James Worrell and Matt Daws) RP 2012 Bordeaux, September 2012 Termination

More information

Reachability problems for Markov chains

Reachability problems for Markov chains 1 Reachability problems for Markov chains S Akshay Dept of CSE, IIT Bombay 1 st CMI Alumni Conference, Chennai 10 Jan 2015 2 Markov chains: a basic model for probabilistic systems 1 3 2 3 3 1 2 1 1 2 2

More information

P -adic root separation for quadratic and cubic polynomials

P -adic root separation for quadratic and cubic polynomials P -adic root separation for quadratic and cubic polynomials Tomislav Pejković Abstract We study p-adic root separation for quadratic and cubic polynomials with integer coefficients. The quadratic and reducible

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

On the Polytope Escape Problem for Continuous Linear Dynamical Systems

On the Polytope Escape Problem for Continuous Linear Dynamical Systems On the Polytope Escape Problem for Continuous Linear Dynamical Systems Joël Ouaknine MPI-SWS and Oxford U joel@mpi-swsorg João Sousa-Pinto Oxford U jspinto@csoxacuk James Worrell Oxford U jbw@csoxacuk

More information

LINEAR FORMS IN LOGARITHMS

LINEAR FORMS IN LOGARITHMS LINEAR FORMS IN LOGARITHMS JAN-HENDRIK EVERTSE April 2011 Literature: T.N. Shorey, R. Tijdeman, Exponential Diophantine equations, Cambridge University Press, 1986; reprinted 2008. 1. Linear forms in logarithms

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

Modeling & Control of Hybrid Systems. Chapter 7 Model Checking and Timed Automata

Modeling & Control of Hybrid Systems. Chapter 7 Model Checking and Timed Automata Modeling & Control of Hybrid Systems Chapter 7 Model Checking and Timed Automata Overview 1. Introduction 2. Transition systems 3. Bisimulation 4. Timed automata hs check.1 1. Introduction Model checking

More information

Effective Positivity Problems for Simple Linear Recurrence Sequences

Effective Positivity Problems for Simple Linear Recurrence Sequences Effective Positivity Problems for Simple Linear Recurrence Sequences Joël Ouaknine and James Worrell Department of Computer Science, Oxford University, UK {joel,jbw}@cs.ox.ac.uk Abstract. We consider two

More information

NORMAL NUMBERS AND UNIFORM DISTRIBUTION (WEEKS 1-3) OPEN PROBLEMS IN NUMBER THEORY SPRING 2018, TEL AVIV UNIVERSITY

NORMAL NUMBERS AND UNIFORM DISTRIBUTION (WEEKS 1-3) OPEN PROBLEMS IN NUMBER THEORY SPRING 2018, TEL AVIV UNIVERSITY ORMAL UMBERS AD UIFORM DISTRIBUTIO WEEKS -3 OPE PROBLEMS I UMBER THEORY SPRIG 28, TEL AVIV UIVERSITY Contents.. ormal numbers.2. Un-natural examples 2.3. ormality and uniform distribution 2.4. Weyl s criterion

More information

The Orbit Problem in Zero and One Dimensions

The Orbit Problem in Zero and One Dimensions The Orbit Problem in Zero and One Dimensions Master's dissertation of Ventsislav K. Chonev Supervised by Joël Ouaknine and James Worrell 1 Contents Contents 2 1 Introduction 3 1.1 Verication of linear

More information

arxiv: v3 [cs.cc] 22 Jun 2016

arxiv: v3 [cs.cc] 22 Jun 2016 0 On the Complexity of the Orbit Problem VENTSISLAV CHONEV, Institute of Science and Technology Austria JOËL OUAKNINE, University of Oxford JAMES WORRELL, University of Oxford arxiv:1303.2981v3 [cs.cc]

More information

Effective Divergence Analysis for Linear Recurrence Sequences

Effective Divergence Analysis for Linear Recurrence Sequences Effective Divergence Analysis for Linear Recurrence Sequences Shaull Almagor Department of Computer Science, Oxford University, UK shaull.almagor@cs.ox.ac.uk Brynmor Chapman MIT CSAIL brynmor@mit.edu Mehran

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

P -ADIC ROOT SEPARATION FOR QUADRATIC AND CUBIC POLYNOMIALS. Tomislav Pejković

P -ADIC ROOT SEPARATION FOR QUADRATIC AND CUBIC POLYNOMIALS. Tomislav Pejković RAD HAZU. MATEMATIČKE ZNANOSTI Vol. 20 = 528 2016): 9-18 P -ADIC ROOT SEPARATION FOR QUADRATIC AND CUBIC POLYNOMIALS Tomislav Pejković Abstract. We study p-adic root separation for quadratic and cubic

More information

Automata-theoretic analysis of hybrid systems

Automata-theoretic analysis of hybrid systems Automata-theoretic analysis of hybrid systems Madhavan Mukund SPIC Mathematical Institute 92, G N Chetty Road Chennai 600 017, India Email: madhavan@smi.ernet.in URL: http://www.smi.ernet.in/~madhavan

More information

DIOPHANTINE EQUATIONS AND DIOPHANTINE APPROXIMATION

DIOPHANTINE EQUATIONS AND DIOPHANTINE APPROXIMATION DIOPHANTINE EQUATIONS AND DIOPHANTINE APPROXIMATION JAN-HENDRIK EVERTSE 1. Introduction Originally, Diophantine approximation is the branch of number theory dealing with problems such as whether a given

More information

On Computably Enumerable Sets over Function Fields

On Computably Enumerable Sets over Function Fields On Computably Enumerable Sets over Function Fields Alexandra Shlapentokh East Carolina University Joint Meetings 2017 Atlanta January 2017 Some History Outline 1 Some History A Question and the Answer

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

Automorphisms of varieties and critical density. Jason Bell (U. Michigan)

Automorphisms of varieties and critical density. Jason Bell (U. Michigan) Automorphisms of varieties and critical density Jason Bell (U. Michigan) belljp@umich.edu 1 Let X be a topological space and let S be an infinite subset of X. We say that S is critically dense if every

More information

1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0

1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0 4 INUTES. If, ω, ω, -----, ω 9 are the th roots of unity, then ( + ω) ( + ω ) ----- ( + ω 9 ) is B) D) 5. i If - i = a + ib, then a =, b = B) a =, b = a =, b = D) a =, b= 3. Find the integral values for

More information

Day 6: 6.4 Solving Polynomial Equations Warm Up: Factor. 1. x 2-2x x 2-9x x 2 + 6x + 5

Day 6: 6.4 Solving Polynomial Equations Warm Up: Factor. 1. x 2-2x x 2-9x x 2 + 6x + 5 Day 6: 6.4 Solving Polynomial Equations Warm Up: Factor. 1. x 2-2x - 15 2. x 2-9x + 14 3. x 2 + 6x + 5 Solving Equations by Factoring Recall the factoring pattern: Difference of Squares:...... Note: There

More information

On a Question of Wintner Concerning the Sequence of Integers Composed of Primes from a Given Set

On a Question of Wintner Concerning the Sequence of Integers Composed of Primes from a Given Set On a Question of Wintner Concerning the Sequence of Integers Composed of Primes from a Given Set by Jeongsoo Kim A thesis presented to the University of Waterloo in fulfillment of the thesis requirement

More information

Recent results on Timed Systems

Recent results on Timed Systems Recent results on Timed Systems Time Petri Nets and Timed Automata Béatrice Bérard LAMSADE Université Paris-Dauphine & CNRS berard@lamsade.dauphine.fr Based on joint work with F. Cassez, S. Haddad, D.

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

CONSEQUENCES OF POWER SERIES REPRESENTATION

CONSEQUENCES OF POWER SERIES REPRESENTATION CONSEQUENCES OF POWER SERIES REPRESENTATION 1. The Uniqueness Theorem Theorem 1.1 (Uniqueness). Let Ω C be a region, and consider two analytic functions f, g : Ω C. Suppose that S is a subset of Ω that

More information

NOTES ON IRRATIONALITY AND TRANSCENDENCE

NOTES ON IRRATIONALITY AND TRANSCENDENCE NOTES ON IRRATIONALITY AND TRANSCENDENCE Frits Beukers September, 27 Introduction. Irrationality Definition.. Let α C. We call α irrational when α Q. Proving irrationality and transcendence of numbers

More information

Effective results for unit equations over finitely generated domains. Jan-Hendrik Evertse

Effective results for unit equations over finitely generated domains. Jan-Hendrik Evertse Effective results for unit equations over finitely generated domains Jan-Hendrik Evertse Universiteit Leiden Joint work with Kálmán Győry (Debrecen) Paul Turán Memorial Conference, Budapest August 22,

More information

Robust computation with dynamical systems

Robust computation with dynamical systems Robust computation with dynamical systems Olivier Bournez 1 Daniel S. Graça 2 Emmanuel Hainry 3 1 École Polytechnique, LIX, Palaiseau, France 2 DM/FCT, Universidade do Algarve, Faro & SQIG/Instituto de

More information

MA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES

MA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES MA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES 2018 57 5. p-adic Numbers 5.1. Motivating examples. We all know that 2 is irrational, so that 2 is not a square in the rational field Q, but that we can

More information

Roth s Theorem on Arithmetic Progressions

Roth s Theorem on Arithmetic Progressions September 25, 2014 The Theorema of Szemerédi and Roth For Λ N the (upper asymptotic) density of Λ is the number σ(λ) := lim sup N Λ [1, N] N [0, 1] The Theorema of Szemerédi and Roth For Λ N the (upper

More information

Pade Approximations and the Transcendence

Pade Approximations and the Transcendence Pade Approximations and the Transcendence of π Ernie Croot March 9, 27 1 Introduction Lindemann proved the following theorem, which implies that π is transcendental: Theorem 1 Suppose that α 1,..., α k

More information

) = nlog b ( m) ( m) log b ( ) ( ) = log a b ( ) Algebra 2 (1) Semester 2. Exponents and Logarithmic Functions

) = nlog b ( m) ( m) log b ( ) ( ) = log a b ( ) Algebra 2 (1) Semester 2. Exponents and Logarithmic Functions Exponents and Logarithmic Functions Algebra 2 (1) Semester 2! a. Graph exponential growth functions!!!!!! [7.1]!! - y = ab x for b > 0!! - y = ab x h + k for b > 0!! - exponential growth models:! y = a(

More information

HTP over Algebraic Extensions of Q: Normforms vs. Elliptic Curves

HTP over Algebraic Extensions of Q: Normforms vs. Elliptic Curves HTP over Algebraic Extensions of Q: Normforms vs. Elliptic Curves Alexandra Shlapentokh East Carolina University Number Theory and Computability, ICMS, June 2007 A Talk in Two Rounds Round I: Normforms

More information

Simultaneous Diophantine Approximation with Excluded Primes. László Babai Daniel Štefankovič

Simultaneous Diophantine Approximation with Excluded Primes. László Babai Daniel Štefankovič Simultaneous Diophantine Approximation with Excluded Primes László Babai Daniel Štefankovič Dirichlet (1842) Simultaneous Diophantine Approximation Given reals integers α,,...,, 1 α 2 α n Q r1,..., r n

More information

Modeling and Analysis of Hybrid Systems

Modeling and Analysis of Hybrid Systems Modeling and Analysis of Hybrid Systems 5. Linear hybrid automata I Prof. Dr. Erika Ábrahám Informatik 2 - LuFG Theory of Hybrid Systems RWTH Aachen University Szeged, Hungary, 27 September - 06 October

More information

Modeling and Analysis of Hybrid Systems Linear hybrid automata I Prof. Dr. Erika Ábrahám Informatik 2 - LuFG Theory of Hybrid Systems RWTH Aachen University Szeged, Hungary, 27 September - 06 October 2017

More information

Badly approximable numbers over imaginary quadratic fields

Badly approximable numbers over imaginary quadratic fields Badly approximable numbers over imaginary quadratic fields Robert Hines University of Colorado, Boulder March 11, 2018 Simple continued fractions: Definitions 1 The Euclidean algorthim a = ba 0 + r 0,

More information

The Schanuel paradigm

The Schanuel paradigm The Schanuel paradigm Jonathan Pila University of Oxford 2014 Clay Research Conference, Oxford Themes Schanuel s conjecture, analogues, in particular Functional analogues and their connections with certain

More information

On values of Modular Forms at Algebraic Points

On values of Modular Forms at Algebraic Points On values of Modular Forms at Algebraic Points Jing Yu National Taiwan University, Taipei, Taiwan August 14, 2010, 18th ICFIDCAA, Macau Hermite-Lindemann-Weierstrass In value distribution theory the exponential

More information

Notes on uniform convergence

Notes on uniform convergence Notes on uniform convergence Erik Wahlén erik.wahlen@math.lu.se January 17, 2012 1 Numerical sequences We begin by recalling some properties of numerical sequences. By a numerical sequence we simply mean

More information

FINITE GROUPS AND EQUATIONS OVER FINITE FIELDS A PROBLEM SET FOR ARIZONA WINTER SCHOOL 2016

FINITE GROUPS AND EQUATIONS OVER FINITE FIELDS A PROBLEM SET FOR ARIZONA WINTER SCHOOL 2016 FINITE GROUPS AND EQUATIONS OVER FINITE FIELDS A PROBLEM SET FOR ARIZONA WINTER SCHOOL 2016 PREPARED BY SHABNAM AKHTARI Introduction and Notations The problems in Part I are related to Andrew Sutherland

More information

CLOSED FORM CONTINUED FRACTION EXPANSIONS OF SPECIAL QUADRATIC IRRATIONALS

CLOSED FORM CONTINUED FRACTION EXPANSIONS OF SPECIAL QUADRATIC IRRATIONALS CLOSED FORM CONTINUED FRACTION EXPANSIONS OF SPECIAL QUADRATIC IRRATIONALS DANIEL FISHMAN AND STEVEN J. MILLER ABSTRACT. We derive closed form expressions for the continued fractions of powers of certain

More information

MCPS Algebra 2 and Precalculus Standards, Categories, and Indicators*

MCPS Algebra 2 and Precalculus Standards, Categories, and Indicators* Content Standard 1.0 (HS) Patterns, Algebra and Functions Students will algebraically represent, model, analyze, and solve mathematical and real-world problems involving functional patterns and relationships.

More information

Hybrid systems and computer science a short tutorial

Hybrid systems and computer science a short tutorial Hybrid systems and computer science a short tutorial Eugene Asarin Université Paris 7 - LIAFA SFM 04 - RT, Bertinoro p. 1/4 Introductory equations Hybrid Systems = Discrete+Continuous SFM 04 - RT, Bertinoro

More information

Existence of a Limit on a Dense Set, and. Construction of Continuous Functions on Special Sets

Existence of a Limit on a Dense Set, and. Construction of Continuous Functions on Special Sets Existence of a Limit on a Dense Set, and Construction of Continuous Functions on Special Sets REU 2012 Recap: Definitions Definition Given a real-valued function f, the limit of f exists at a point c R

More information

O-minimality and Diophantine Geometry

O-minimality and Diophantine Geometry O-minimality and Diophantine Geometry Jonathan Pila University of Oxford ICM 2014, Seoul Themes Rational points on definable sets in o-minimal structures and applications to some Diophantine problems comprehended

More information

Automatic Sequences and Transcendence of Real Numbers

Automatic Sequences and Transcendence of Real Numbers Automatic Sequences and Transcendence of Real Numbers Wu Guohua School of Physical and Mathematical Sciences Nanyang Technological University Sendai Logic School, Tohoku University 28 Jan, 2016 Numbers

More information

Mahler measure as special values of L-functions

Mahler measure as special values of L-functions Mahler measure as special values of L-functions Matilde N. Laĺın Université de Montréal mlalin@dms.umontreal.ca http://www.dms.umontreal.ca/~mlalin CRM-ISM Colloquium February 4, 2011 Diophantine equations

More information

Math 430/530: Advanced Calculus I Chapter 1: Sequences

Math 430/530: Advanced Calculus I Chapter 1: Sequences Math 430/530: Advanced Calculus I Chapter 1: Sequences Fall 2018 Sequences A sequence is just a list of real numbers a 1, a 2,... with a definite order. We write a = (a n ) n=1 for a sequence a with n

More information

On the Total Variation Distance of Labelled Markov Chains

On the Total Variation Distance of Labelled Markov Chains On the Total Variation Distance of Labelled Markov Chains Taolue Chen Stefan Kiefer Middlesex University London, UK University of Oxford, UK CSL-LICS, Vienna 4 July 04 Labelled Markov Chains (LMCs) a c

More information

8. Dirichlet s Theorem and Farey Fractions

8. Dirichlet s Theorem and Farey Fractions 8 Dirichlet s Theorem and Farey Fractions We are concerned here with the approximation of real numbers by rational numbers, generalizations of this concept and various applications to problems in number

More information

Root separation for irreducible integer polynomials

Root separation for irreducible integer polynomials Root separation for irreducible integer polynomials Yann Bugeaud and Andrej Dujella 1 Introduction The height H(P ) of an integer polynomial P (x) is the maximum of the absolute values of its coefficients.

More information

Ramanujan s first letter to Hardy: 5 + = 1 + e 2π 1 + e 4π 1 +

Ramanujan s first letter to Hardy: 5 + = 1 + e 2π 1 + e 4π 1 + Ramanujan s first letter to Hardy: e 2π/ + + 1 = 1 + e 2π 2 2 1 + e 4π 1 + e π/ 1 e π = 2 2 1 + 1 + e 2π 1 + Hardy: [These formulas ] defeated me completely. I had never seen anything in the least like

More information

A Note on Cyclotomic Integers

A Note on Cyclotomic Integers To the memory of Alan Thorndike, former professor of physics at the University of Puget Sound and a dear friend, teacher and mentor. A Note on Cyclotomic Integers Nicholas Phat Nguyen 1 Abstract. In this

More information

L-Polynomials of Curves over Finite Fields

L-Polynomials of Curves over Finite Fields School of Mathematical Sciences University College Dublin Ireland July 2015 12th Finite Fields and their Applications Conference Introduction This talk is about when the L-polynomial of one curve divides

More information

A remark about the positivity problem of fourth order linear recurrence sequences

A remark about the positivity problem of fourth order linear recurrence sequences ACTA ET COMMENTATIONES UNIVERSITATIS TARTUENSIS DE MATHEMATICA Volume 18, Number 1, June 2014 Available online at http://acutm.math.ut.ee A remark about the positivity problem of fourth order linear recurrence

More information

TRANSCENDENTAL NUMBERS AND PERIODS. Contents

TRANSCENDENTAL NUMBERS AND PERIODS. Contents TRANSCENDENTAL NUMBERS AND PERIODS JAMES CARLSON Contents. Introduction.. Diophantine approximation I: upper bounds 2.2. Diophantine approximation II: lower bounds 4.3. Proof of the lower bound 5 2. Periods

More information

Introduction to Techniques for Counting

Introduction to Techniques for Counting Introduction to Techniques for Counting A generating function is a device somewhat similar to a bag. Instead of carrying many little objects detachedly, which could be embarrassing, we put them all in

More information

The hyperbolic Ax-Lindemann-Weierstraß conjecture

The hyperbolic Ax-Lindemann-Weierstraß conjecture The hyperbolic Ax-Lindemann-Weierstraß conjecture B. Klingler (joint work with E.Ullmo and A.Yafaev) Université Paris 7 ICM Satellite Conference, Daejeon Plan of the talk: (1) Motivation: Hodge locus and

More information

Part IA Numbers and Sets

Part IA Numbers and Sets Part IA Numbers and Sets Definitions Based on lectures by A. G. Thomason Notes taken by Dexter Chua Michaelmas 2014 These notes are not endorsed by the lecturers, and I have modified them (often significantly)

More information

Discrete abstractions of hybrid systems for verification

Discrete abstractions of hybrid systems for verification Discrete abstractions of hybrid systems for verification George J. Pappas Departments of ESE and CIS University of Pennsylvania pappasg@ee.upenn.edu http://www.seas.upenn.edu/~pappasg DISC Summer School

More information

Ramsey theory. Andrés Eduardo Caicedo. Graduate Student Seminar, October 19, Department of Mathematics Boise State University

Ramsey theory. Andrés Eduardo Caicedo. Graduate Student Seminar, October 19, Department of Mathematics Boise State University Andrés Eduardo Department of Mathematics Boise State University Graduate Student Seminar, October 19, 2011 Thanks to the NSF for partial support through grant DMS-0801189. My work is mostly in set theory,

More information

Waiting times, recurrence times, ergodicity and quasiperiodic dynamics

Waiting times, recurrence times, ergodicity and quasiperiodic dynamics Waiting times, recurrence times, ergodicity and quasiperiodic dynamics Dong Han Kim Department of Mathematics, The University of Suwon, Korea Scuola Normale Superiore, 22 Jan. 2009 Outline Dynamical Systems

More information

DIFFERENTIABILITY OF A PATHOLOGICAL FUNCTION, DIOPHANTINE APPROXIMATION, AND A REFORMULATION OF THE THUE-SIEGEL-ROTH THEOREM

DIFFERENTIABILITY OF A PATHOLOGICAL FUNCTION, DIOPHANTINE APPROXIMATION, AND A REFORMULATION OF THE THUE-SIEGEL-ROTH THEOREM DIFFERENTIABILITY OF A PATHOLOGICAL FUNCTION, DIOPHANTINE APPROXIMATION, AND A REFORMULATION OF THE THUE-SIEGEL-ROTH THEOREM JUAN LUIS VARONA Abstract. We study the differentiability of the real function

More information

A NON-ARCHIMEDEAN AX-LINDEMANN THEOREM

A NON-ARCHIMEDEAN AX-LINDEMANN THEOREM A NON-ARCHIMEDEAN AX-LINDEMANN THEOREM François Loeser Sorbonne University, formerly Pierre and Marie Curie University, Paris Model Theory of Valued fields Institut Henri Poincaré, March 8, 2018. P. 1

More information

Critical density, automorphisms, and Harm Derksen. Jason Bell

Critical density, automorphisms, and Harm Derksen. Jason Bell Critical density, automorphisms, and Harm Derksen Jason Bell 1 Let X be a topological space and let S be an infinite subset of X. We say that S is critically dense if every infinite subset of S is dense

More information

INTRODUCTION TO O-MINIMAL STRUCTURES AND AN APPLICATION TO NEURAL NETWORK LEARNING. Contents

INTRODUCTION TO O-MINIMAL STRUCTURES AND AN APPLICATION TO NEURAL NETWORK LEARNING. Contents INTRODUCTION TO O-MINIMAL STRUCTURES AND AN APPLICATION TO NEURAL NETWORK LEARNING MARCUS TRESSL Contents 1. Definition of o-minimality and first examples 1 2. Basic real analysis and the monotonicity

More information

Notes on State Minimization

Notes on State Minimization U.C. Berkeley CS172: Automata, Computability and Complexity Handout 1 Professor Luca Trevisan 2/3/2015 Notes on State Minimization These notes present a technique to prove a lower bound on the number of

More information

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.

More information

Computing coefficients of modular forms

Computing coefficients of modular forms Computing coefficients of modular forms (Work in progress; extension of results of Couveignes, Edixhoven et al.) Peter Bruin Mathematisch Instituut, Universiteit Leiden Théorie des nombres et applications

More information

Notes on Equidistribution

Notes on Equidistribution otes on Equidistribution Jacques Verstraëte Department of Mathematics University of California, San Diego La Jolla, CA, 92093. E-mail: jacques@ucsd.edu. Introduction For a real number a we write {a} for

More information

Exponential and Weierstrass equations

Exponential and Weierstrass equations Exponential and Weierstrass equations Jonathan Kirby 8th April 2005 1 Complex functions Exponentiation is a group homomorphism C exp C. For each elliptic curve E there is a -function which gives a group

More information

COMPLEX NUMBERS WITH BOUNDED PARTIAL QUOTIENTS

COMPLEX NUMBERS WITH BOUNDED PARTIAL QUOTIENTS COMPLEX NUMBERS WITH BOUNDED PARTIAL QUOTIENTS WIEB BOSMA AND DAVID GRUENEWALD Abstract. Conjecturally, the only real algebraic numbers with bounded partial quotients in their regular continued fraction

More information

Econ Slides from Lecture 1

Econ Slides from Lecture 1 Econ 205 Sobel Econ 205 - Slides from Lecture 1 Joel Sobel August 23, 2010 Warning I can t start without assuming that something is common knowledge. You can find basic definitions of Sets and Set Operations

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

Approximation exponents for algebraic functions in positive characteristic

Approximation exponents for algebraic functions in positive characteristic ACTA ARITHMETICA LX.4 (1992) Approximation exponents for algebraic functions in positive characteristic by Bernard de Mathan (Talence) In this paper, we study rational approximations for algebraic functions

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

Factoring univariate polynomials over the rationals

Factoring univariate polynomials over the rationals Factoring univariate polynomials over the rationals Tommy Hofmann TU Kaiserslautern November 21, 2017 Tommy Hofmann Factoring polynomials over the rationals November 21, 2017 1 / 31 Factoring univariate

More information

ZEROES OF INTEGER LINEAR RECURRENCES. 1. Introduction. 4 ( )( 2 1) n

ZEROES OF INTEGER LINEAR RECURRENCES. 1. Introduction. 4 ( )( 2 1) n ZEROES OF INTEGER LINEAR RECURRENCES DANIEL LITT Consider the integer linear recurrence 1. Introduction x n = x n 1 + 2x n 2 + 3x n 3 with x 0 = x 1 = x 2 = 1. For which n is x n = 0? Answer: x n is never

More information

Math 61CM - Solutions to homework 2

Math 61CM - Solutions to homework 2 Math 61CM - Solutions to homework 2 Cédric De Groote October 5 th, 2018 Problem 1: Let V be the vector space of polynomials of degree at most 5, with coefficients in a field F Let U be the subspace of

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

1 Prologue (X 111 X 14 )F = X 11 (X 100 X 3 )F = X 11 (X 200 X 6 ) = X 211 X 17.

1 Prologue (X 111 X 14 )F = X 11 (X 100 X 3 )F = X 11 (X 200 X 6 ) = X 211 X 17. Prologue In this chapter we give a couple of examples where the method of auxiliary polynomials is used for problems that have no diophantine character. Thus we are not following Sam Goldwyn s advice to

More information

MA651 Topology. Lecture 9. Compactness 2.

MA651 Topology. Lecture 9. Compactness 2. MA651 Topology. Lecture 9. Compactness 2. This text is based on the following books: Topology by James Dugundgji Fundamental concepts of topology by Peter O Neil Elements of Mathematics: General Topology

More information

Bjorn Poonen. MSRI Introductory Workshop on Rational and Integral Points on Higher-dimensional Varieties. January 18, 2006

Bjorn Poonen. MSRI Introductory Workshop on Rational and Integral Points on Higher-dimensional Varieties. January 18, 2006 University of California at Berkeley MSRI Introductory Workshop on Rational and Integral Points on Higher-dimensional Varieties January 18, 2006 The original problem H10: Find an algorithm that solves

More information

On Reachability for Hybrid Automata over Bounded Time

On Reachability for Hybrid Automata over Bounded Time On Reachability for Hybrid Automata over Bounded Time Thomas Brihaye, Laurent Doyen 2, Gilles Geeraerts 3, Joël Ouaknine 4, Jean-François Raskin 3, and James Worrell 4 Université de Mons, Belgium 2 LSV,

More information

Notes on Complex Analysis

Notes on Complex Analysis Michael Papadimitrakis Notes on Complex Analysis Department of Mathematics University of Crete Contents The complex plane.. The complex plane...................................2 Argument and polar representation.........................

More information

Classical transcendental curves

Classical transcendental curves Classical transcendental curves Reinhard Schultz May, 2008 In his writings on coordinate geometry, Descartes emphasized that he was only willing to work with curves that could be defined by algebraic equations.

More information

INDEX UNIT 3 TSFX REFERENCE MATERIALS 2014 ALGEBRA AND ARITHMETIC

INDEX UNIT 3 TSFX REFERENCE MATERIALS 2014 ALGEBRA AND ARITHMETIC INDEX UNIT 3 TSFX REFERENCE MATERIALS 2014 ALGEBRA AND ARITHMETIC Surds Page 1 Algebra of Polynomial Functions Page 2 Polynomial Expressions Page 2 Expanding Expressions Page 3 Factorising Expressions

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

Exercises Chapter II.

Exercises Chapter II. Page 64 Exercises Chapter II. 5. Let A = (1, 2) and B = ( 2, 6). Sketch vectors of the form X = c 1 A + c 2 B for various values of c 1 and c 2. Which vectors in R 2 can be written in this manner? B y

More information

Complexity of Reachability, Mortality and Freeness Problems for Matrix Semigroups

Complexity of Reachability, Mortality and Freeness Problems for Matrix Semigroups Complexity of Reachability, Mortality and Freeness Problems for Matrix Semigroups Paul C. Bell Department of Computer Science Loughborough University P.Bell@lboro.ac.uk Co-authors for todays topics: V.

More information

Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation. Specialized values

Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation. Specialized values KleinInvariant Notations Traditional name Klein invariant modular function Traditional notation z Mathematica StandardForm notation KleinInvariant z Primary definition 09.50.0.0001.01 ϑ 0, Π z 8 ϑ 3 0,

More information

Results and open problems related to Schmidt s Subspace Theorem. Jan-Hendrik Evertse

Results and open problems related to Schmidt s Subspace Theorem. Jan-Hendrik Evertse Results and open problems related to Schmidt s Subspace Theorem Jan-Hendrik Evertse Universiteit Leiden 29 ièmes Journées Arithmétiques July 6, 2015, Debrecen Slides can be downloaded from http://pub.math.leidenuniv.nl/

More information

College Algebra. Basics to Theory of Equations. Chapter Goals and Assessment. John J. Schiller and Marie A. Wurster. Slide 1

College Algebra. Basics to Theory of Equations. Chapter Goals and Assessment. John J. Schiller and Marie A. Wurster. Slide 1 College Algebra Basics to Theory of Equations Chapter Goals and Assessment John J. Schiller and Marie A. Wurster Slide 1 Chapter R Review of Basic Algebra The goal of this chapter is to make the transition

More information

Lecture 4: Completion of a Metric Space

Lecture 4: Completion of a Metric Space 15 Lecture 4: Completion of a Metric Space Closure vs. Completeness. Recall the statement of Lemma??(b): A subspace M of a metric space X is closed if and only if every convergent sequence {x n } X satisfying

More information

1 Rings 1 RINGS 1. Theorem 1.1 (Substitution Principle). Let ϕ : R R be a ring homomorphism

1 Rings 1 RINGS 1. Theorem 1.1 (Substitution Principle). Let ϕ : R R be a ring homomorphism 1 RINGS 1 1 Rings Theorem 1.1 (Substitution Principle). Let ϕ : R R be a ring homomorphism (a) Given an element α R there is a unique homomorphism Φ : R[x] R which agrees with the map ϕ on constant polynomials

More information

Introduction to Diophantine methods Michel Waldschmidt miw/courshcmuns2007.html. Sixth course: september 21, 2007.

Introduction to Diophantine methods Michel Waldschmidt   miw/courshcmuns2007.html. Sixth course: september 21, 2007. Introduction to Diophantine methods Michel Waldschmidt http://www.math.jussieu.fr/ miw/courshcmuns2007.html Sixth course: september 2, 2007. 8 2.2.2 Diophantine approximation and applications Diophantine

More information

Zero estimates for polynomials inspired by Hilbert s seventh problem

Zero estimates for polynomials inspired by Hilbert s seventh problem Radboud University Faculty of Science Institute for Mathematics, Astrophysics and Particle Physics Zero estimates for polynomials inspired by Hilbert s seventh problem Author: Janet Flikkema s4457242 Supervisor:

More information