Evading Triangles without a Map by Braxton Carrigan A thesis submitted to the Graduate Faculty of Auburn University in partial ful llment of the requirements for the Degree of Master of Science Auburn, Alabama May 14, 2010 Keywords: shortest path, permeability, evading Copyright 2010 by Braxton Carrigan Approved by Andras Bezdek, Chair, Professor of Mathematics Wlodzimierz Kuperberg, Professor of Mathematics Chris Rodger, Professor of Mathematics George Flowers, Acting Dean of Graduate School Abstract In this thesis we will solve the following shortest path problem. Let P be an arrange- ment of equilateral non-overlapping translated triangles in the plane and two points S and T so that the segment ST is parallel to a side of each of the triangles. Assume one needs to navigate from point S to point T by evading the triangular obstacles without any previous knowledge of the location of the obstacles. The navigator becomes aware of a triangle once it is contacted along the path. We will give an algorithm which enables the navigator to reach the target point T by a path of length at mostp3(d+ 2d), where d is the length of ST. Section 4 contains the proof, which is preceded by three sections reviewing some of the main results and methods of previously considered shortest path problems. In particular we will outline three papers concerning shortest path problems. First we will address the idea of permeability of a layer mentioned by J. Pach [2] in \On the Permeability Problem" using an integration technique. Then, we will show an improvement of permeability by G. Fejes T oth [4] in his paper entitled "Evading Convex Discs" via existence of a path using the sweeping of a direction technique. Finally we will outline Chapter 3 of "Shortest Paths without a Map" by Papadimitriou and Yannakakis [3], where they show three simple heuristics of evading rectangles. ii Acknowledgments I would like to thank my advisor, Dr. Bezdek, for his academic guidance. Without his leadership, I would not have been able to come about this conclusion or culminate it into the written thesis. I also would like to extend my gratitude to the members of my committee: Dr. Wlodz- imierz Kuperberg and Dr. Chris Rodger for their time and suggestions on improving my work. Working with them in and out of the classroom has truly been a pleasure. I also would like to extend an appreciation to my parents Wanda and Ralph and my wife Catie. Throughout my life all three have been a source of strength, support, and en- couragement which has enabled me to pursue my dreams. iii Table of Contents Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Acknowledgments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii List of Illustrations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . v 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 2 Background Research . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 3 Problem and Heuristic . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 4 Proof of Theorem 3.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 5 Lower Bound for any Heuristic and Open Problems . . . . . . . . . . . . . . . . 25 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 iv List of Illustrations 2.1 Path Px in rectangle R . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 2.2 Square Regions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 2.3 Region 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 2.4 Region 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.5 Labeling Square . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 2.6 Multiple paths . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 2.7 Rectangle con guration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 2.8 graph . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 3.1 Path from S to T . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 3.2 Path from S to T . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 4.1 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 4.2 Case 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 4.3 Case 2-A . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 4.4 Maximizing Case 2-A . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 4.5 Case 2-B . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 v 4.6 Maximizing Case 2-B . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 4.7 Case 3 - A . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 4.8 Maximizing Case 3 - A . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 4.9 Case 3 - B . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 4.10 Maximizing Case 3 - B . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 5.1 Path through densest packing . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25 5.2 Layers between S and T . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 5.3 Navigating a Layer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 vi Chapter 1 Introduction The idea of evading convex discs was rst addressed in the late 1960?s by L. Fejes T oth. By the 1980?s, the topic of evading convex discs was being studied by mathematicians in a new mathematical eld called computational geometry. Numerous publications centered around the idea of evading convex shapes in the plane. The title \evading triangles without a map" refers to navigating through a plane to- ward a target point while evading unknown triangular obstacles. To visualize this problem, imagine walking through a thick forest with a compass and knowing in what direction and how far to travel. To restrict this forest we will assume there are mountains that are each translates of other, have bases of equilateral triangles, and are impassable. Finally we will also assume that the given direction is parallel to one side of the mountains? base. While moving the traveler cannot see through the thick forest, therefore he must rely on making decision once encountering a mountain. 1 Chapter 2 Background Research In the 1970?s several results were published about the notion of permeability of layers. A layer is de ned as a parallel strip of a plane with width w containing non-overlapping open domains called obstacles. The permeability of a layer was de ned by L. Fejes T oth to be winf?, where ? is the length of a path from one edge of the layer to the other which evades all open domains of the layer. J. Pach [2] in "On the Permeability Problem" proved the following: Theorem 2.1. (J. Pach) The permeability of any layer of squares is at most 23, and this bound can be achieved. Let R = x1x2x3x4 be a rectangle and R = Siji2I be a set of non-overlapping squares inside R (Figure 2.1). For any point x 2x1x2 we will de ne a path Px to be a path from x to y, where y is the only point on x3x4 such that the segment xy is perpendicular to x1x2, which evades R. In order to evade R for every segment uivi2xy which intersects a square S 2 R replace uivi by the portion of the boundary between ui and vi on Si which is the shortest; we will denote this as li(x). Let L(x) be the length of the path Px Lemma 2.1. Rx2x1 L(x)dx 23A(R), where A(R) is the area of the rectangle R. Proof of Lemma 2.1 will rst show Rx2 x1 (li(x) d(ui(x)vi(x)))dx A(Si) 2 . If i(x) denotes the length of the shortest path of the boundary of Si connecting ui(x) and 2 x1 x2 x3x4 X Y Px Figure 2.1: Path Px in rectangle R vi(x), we have Z x2 x1 li(x) d(ui(x)vi(x))dx = Z x2 x1 i(x)dx A(Si) R i(x)dx is evaluated in the following way If Si has two of its sides parallel to x1x2 the act of computing is trivially 23A(R) , thus we will assume there exist a positive angle i between one of the sides of Si and x1x2 (Figure 2.1). We can also pick a side of Si so that 0 < i < 4 . Now we will label the vertices s1, s2, s3, and s4 moving left to right with respect to x1x2. We will de ne the side length of Si to be di. We will now split Si into three regions R1, R2, and R3 by constructing lines perpendicular to x1x2 through each of the vertices of Si (Figure 2.2). For each region we will consider the leftmost vertical line to be it?s y-axis and de ne ri to be the interval spanned horizontally by region Ri. 3 ? s1 s2 s3 s4 R1 R2 R3 di r1 r2 r3 Figure 2.2: Square Regions We will compute R i(x)dx as Rr1 i(x)dx + Rr2 i(x)dx + Rr3 i(x)dx. Now by symmetry R1 = R3, thus we need only compute for R1 and R2. i = xdi sin( i)(di +di tan( i)) inside R1.(Figure 2.3) Since the horizontal distance of the region isdi sin( i) we calculateRr1 i(x)dxasRdi sin( i)0 xdi sin( i)(di+ di tan( i))dx = x2(di+di tan( i))2(di sin( i)) jdi sin( i)0 = (di sin( i))(di+di tan( i))2 . ? S1 S2 S3 S4di disin(?i) ditan(?i) R1 Figure 2.3: Region 1 Now for R2 we calculate the integral over half of the region since it is symmetric. Therefore 4 for the left half of the interval r2, (Figure 2.4) i = 2di 2( xcos( i)) with the width of r2 = dip2 cos( 4 + i). So Rr2 i(x)dx = 2R dip2cos( 4 + i) 2 0 2di 2( x cos( i))dx = 2(2dix x2 cos( i)j dip2cos( 4 + i) 20 ) = dip2 cos( 4 + i)(di di p2 cos( 4 + i) 4 cos( i) ) ? S1 S2 S3 S4 di?2cos(pi4 +?i) x Figure 2.4: Region 2 Now by adding the two previous results we get the left half of Si and by symmetry we mul- tiply this sum by 2. Therefore R i(x)dx = di sin( i)(di +di tan( i) +dip2 cos( 4 + i)(2di dip2 cos( 4 + i) 2 cos( i) = 3 2d 2 i 1+2 cos2( i) 3 cos( i) 3 2A(Si). Therefore Rx2x1 li(x) d(ui(x)vi(x))dx = R i(x)dx A(Si) 32A(Si) A(Si) = A(Si)2 Thus Rx2x1 li(x)dx = Rx2x1 li(x) d(ui(x)vi(x))+d(ui(x)vi(x))dx = Rx2x1 li(x) d(ui(x)vi(x))dx+ Rx2 x1 d(ui(x)vi(x))dx A(Si) 2 +A(Si) = 3A(Si) 2 . The bound is achieved when R lls the layer (a tiling), thus L(x) is just the sum of all li(x) which produces the bound of 32. G. Fejes T oth [4] improved this bound in his paper entitled "Evading Convex Discs". 5 Theorem 2.2. (G. Fejes T oth) Given a set of disjoint open squares with side-lengths not exceeding 1, any two points of the plane lying outside the squares at distance q from one another can be connected by a path evading the squares and having length at most 3q+12 . We let be the set of open squares of side-length at most 1, as in the previous theorem. For labeling purposes we assume X is on the side X1X2. We de ne (XY) to be the shortest path on the boundary of S from X to a point Y on S. Also de ne d(XY) = (XY) to be the distance gained by the vector !XY in direction d. We de ne the path (P;d) to be an in nite path emanating from P in direction d. (P;d) = is constructed by rst traveling on a ray !R in direction of d, then when contacting a square S 2 , evade S by traveling along the boundary of S from X to Xi such that (XXi) (XXi) is maximum, and then continue in the direction of d from point Xi. Notice that (P;d) is not uniquely determined whenever for some S2 there exist more than one point Xi at which (XXi) (XXi) is maximum. Now for some point X on we de ne l(x) to be the length of the arc of between P and X. To prove Theorem 2.2, rst show l(x) 32( (PX) + 1) To do so we need the following: Lemma 2.2. Let X be a boundary point on the square S = X1X2X3X4 with side length 1, such that the ray emanating from X in the direction d intersects S, then we have: max 1 i 4 (XXi) (XXi) 2 3 Proof of Lemma 2.2 For orientation purposes we assume that d is vertical and oriented downward. Also we can assume that contacts S between X1 and X2. Since the case where the side X1X2 is perpendicular to d is trivial, we can also by symmetry assume that X1 is the vertex which 6 is vertically the highest. Now de ne a = (X1X2) and b = (X1X4) (Figure 2.5). Obviously a2 +b2 = 1. X1 X2 X3 X4 d X a b Figure 2.5: Labeling Square Now we consider two cases Case 1: a b In this case a p2 2 , therefore (XX2) (XX2) = a p2 2 > 2 3 Case 2: 0 23 and (X2X4) (X2X4) = b a2 < 12 < 23. Also notice that as X moves from X1 to X2 the ratio (XX4) (XX4) decreases, therefore there exists some point X0 between X1 and X2 where (X0X4) (X0X4) = 23. Obviously if X 2 X0X1, then (XX4) (XX4) 2 3. Now we assume X2X0X2 and show that (XX3) (XX3) > 23. First notice that (XX3) (XX3) increases as X approaches X2, therefore it su ces to show that (X0X3) (X0X3) 2 3. We denote the distance from X0 to X1 by x therefore producing: (X0X3) = a+b ax with (X0X3) = 2 x, 7 (X0X4) = b ax with (X0X4) = 1 +x. Since (X0X4) (X0X4) = 23 = b ax1+x , we have 2(1 + x) = 3(b - ax) We want to show that 2(2 x) < 3(a+b ax). Since b> 0 this is equivalent to 3a2 4a+1 < 0, which is indeed true for p2 2