Download Document

Document related concepts

Probability interpretations wikipedia , lookup

Ars Conjectandi wikipedia , lookup

Transcript
The Zone Theorem
The Cutting Lemma Revisited
Tom Jurgenson
1
The Zone Theorem
2
Definitions reminders
 Is a sub-space of d-1 dimensions.
 Is a partition of
into relatively open convex sets.
 Are 0/1/(d-1)-dimension faces (respectively) in
 Are the connected components of
3
.
also called d-faces
Example d=2 (Plane)
 Vertices are points, edges are segments, and facets are also
segments since 2-1=1.
 Cells are 2-dimesional connected components in
 The blue lines are the Hyperplanes in H
 Their arrangement creates:
 Vertices – in red
 Edges\Facets – blue segments(excluding red)
 Cells – light blue areas.
4
The Zone Theorem - Motivation
 Motivation – we would like to bound the number of faces
that can see any part of a given hyperplane. For that we’ll
need to define what ‘see’ means, and define the set of objects
(faces) that ‘see’ g (the zone).
5
Definition for ‘sees’
 A face F can see hyperplane g of H (set of hyperplanes), if
there are points
such that the open segment
xy does not intersect any hyperplane of H.
 Important Note: It does not matter which point x we
choose, either all of them can see g or none can.
 In this example, the red line represents a segment between
the thick line – g, and the face – the cell. The blue line
exemplifies the fact that no segment connects that cell with
the hyperplane g.
6
Zone Definition
 Zone – The zone of hyperplane g, is the set of the faces of the
arrangement of H that can see g.
 Here is the zone (dark grey area) for hyperplane g (the thick line)
and hyperplane set H (all the thin lines). Each face (vertex, edge
and cell) in the zone of g can see g.
7
The Zone Theorem
The number of faces in the zone of any hyperplane in an
arrangement of n hyperplanes in
is
with the
constant of proportionality depending on d.
we will prove this theorem in an inductive manner
8
Base case, d=2 (The plane)
 Let H be a set of n lines in the plane in general position. We
consider the zone of line g.
 In this case the faces would be: vertices (0-faces), edges (1faces), and the cells are convex polygons that may or may not
be bounded(1-faces).
 The bound for the number of faces in the zone in this case
should be
9
Bound for Edges is enough
 The number of vertices is proportional to the number of
edges since in a convex polygon the number of vertices is
equal to the number of edges. And in open cells, the number
of edges is at most twice the number of vertices.
 It is clear that the number of polygons is bounded by the
number of vertices by using Euler’s formula: F=v-e+2.
Therefore proportional to the number of edges.
10
Bound for Edges
 Let g be a horizontal line.
 We’ll count the number of edges that see g that are also
above g.
 The number of edges that intersect g that are also above g is
bounded by n. Since each line (hyperplane in H) out of the n
lines can only intersect g in one point.
 The other edges are disjoint from g, and in the next slides we
are going to find a bound for these edges.
11
Bound for the disjoint edges
 Let uv be an edge that sees g but is disjoint from it.
 Let
be the line containing uv. And let a be the
intersection of g and h.
 Since u is a vertex - it is the intersection of two lines – so let
be the other line that passes through u. Let b be the
intersection of l and g.
12
Right Edge of Line l
 Definition: uv is right edge of l if b is on the right of a. uv is
left edge of l if b is on the left of a.
 We’ll show that for each such
there is at most one
right edge.
 Let’s assume uv is a right edge for l.
13
Only one right edge
 Assume by contradiction that there are two right edges for l.
Let uv and xy be the two right edges, and assume that uv lies
below xy.
 By definition xy should see some point of g.
 However, every point to the left of a on g is obscured by the
line l. And every point to the right of a on g is obscured by
the line h.
 This contradiction shows that only one right edge can exist
for each line.
14
Only one right edge (2)
15
Conclusion for d=2
 Each line in H has at most one right edge, and by symmetry
also at most one left edge. There is also an edge that
intersects g. The same can be argued for edges below g, for
the lower right and lower left edges, and an edge that
intersects. Altogether there are at most 6 edges for each
line in H in the zone.
 There are n lines in H so the number of edges is at most 6n.
Linear in n.
 We already argued that the total number of faces is linear in
the number of edges therefore the number of faces in the
zone when d=2 is:
as needed.
16
Inductive Step from d-1 to d
 Assume that the total number of faces of a zone in
is
 Mark
as the max. number of (d-1) faces in the zone in
an arrangement of n hyperplanes in
.
 Let H be an arrangement of hyperplanes and g a base
hyperplane s.t
is attained for them.
17
Motivation: number of blue facets
 We color one hyperplane in H red, let it be




18
. The rest
n-1 hyperplanes in H we color blue.
A Blue Facet is a facet that lies in a blue hyperplane.
Since every facet has a chance of
becoming blue
(choose any blue hyperplane out of n blue hyperplanes), the
expected number of blue facets is
Now we bound the expected number of blue facets in a
different way and use it to estimate
Example in the next slide…
Blue and Red Facets
19
Bound for Expected Blue Facets
 Consider the arrangement of blue hyperplanes only. Since
there are n-1 such hyperplanes the number of blue facets is
bounded by
 In the next stage we add the red hyperplane h and see how it
affects the number of blue facets in the zone of g.
20
The increase in number of blue facets
caused by h
 The red hyperplane h may split an existing blue facet F into
two new facets – F1 and F2. If both new facets can see g then
we get an increase in the number of facets.
 Claim:
is visible from
(we’ll prove it later)
 If we intersect all blue hyperplanes with g and with h we get
a d-1 dimensional arrangement in which
is a facet in
a zone of the d-2 dimensional hyperplane
. Using the
induction assumption we get
such “splits”. The
total number of facets becomes:
21
Split Example
22
Proof for the above claim
 Reminder: We wanted to show the following claim:
is visible from
 Proof: Let C be a cell of the zone in the arrangement of blue
hyperplanes having F on the boundary.
 Because F1 and F2 see g we get that
sees
and
sees
 The interior of
is contained in C and we get that
the intersection of
with the hyperplane h
contains a segment witnessing the visibility of
from
23
Split Example (II)
24
Conclusion for facets
 On the one hand we calculated the expected number of blue
facets:
 On the other hand we found an upper bound for the number
of blue facets:
 We combine the two and get:
25
Conclusion for facets (II)
 And as we can see:
as needed.
 Next we are going to discuss d-k faces and not just d-1 faces (or
facets)
26
Expected number of d-k faces
 Denote




27
the max. possible number of j-faces in the
zone for an arrangement of n hyperplanes of dimension d.
Let H be an arrangement of n hyperplanes where
is attained.
Again a random hyperplane h in H is colored red and the rest
blue.
Reminder: a d-k face is the intersection of k hyperplanes.
A d-k face is blue if its relative interior is disjoint from the
red hyperplane h.
Expected number of d-k faces (2)
 The probability for a d-k face to be blue is
:
[For the selection of first blue hyperplane, then the second on
so on until we select the k-th hyperplane]
 As in the facet case we get that the expected number of blue
d-k faces in the required zone is:
28
Upper bound for d-k faces
 By adding the red hyperplane h the number of blue d-k faces
can increase by at most
by the inductive
hypothesis.
 Using the same arguments in lower dimensions as the facets
case we get to the conclusion that:
 And in the next slide we’ll show:
29
Algebraic proof
30
Bound for edges and vertices
 For the case k=d-1 (edges) by using this method we only get
the bound:
 So the number of edges and vertices must be bound
separately:
 Vertices: Each vertex is contained in some3-face of the zone.
Within such 3-face, the number of vertices is at most 3 times
the number of 2-faces, because the 3-face is a 3D convex
polyhedron. Since H is simple each 2-face is contained in a
bounded number of 3-face. It follows that the total number of
vertices is at most proportional to
31
Final details
 Edges are proportional to vertices.
 In conclusion, for every k=1…d-2 we get
 For vertices and edges we get
 In conclusion we get that the zone has
32
Zones in other arrangements
 The max. complexity of a zone can be investigated for
objects other than hyperplanes, and that leads to many
related problems that we won’t discuss here.
33
Cutting Lemma Revisited
34
General Plan
 First, we’ll prove the cutting lemma for the planar case with
a tight bound. We already saw a proof for a bound that was
not tight that uses random sampling, and a different proof
that is tight but does not use random sampling.
 The following proof is both tight and uses random
sampling in a way that makes it easy to be generalized
to higher dimensions.
 But first we’ll do a review for some cutting-lemma related
definitions.
35
Review of Definitions (I)
 A generalized triangle is the intersection of 3 half planes.
 A cutting of the plane into generalized triangles, is the
subdivision of the plane to a disjoint set generalized triangles
that their union covers the entire plane.
36
Review of Definitions (II)
 A
cutting of the plane with a set of n of lines H – is a
cutting of the plane into generalized triangles
s.t the interior of each such
is intersected
by at most
lines from H.
 Note – the triangles’ edges may or may not include lines
from H.
37
½ cutting example
 In the next example we have a ½ cutting. The red lines are
lines in H, the black lines are the lines that create the
triangulation.
38
The Cutting Lemma
 For every set H of n lines in the plane and every r>1 there
exists a 1/r cutting for H the size
 In other words, there is a subdivision of the plane into
generalized triangles s.t the interior of each triangle is
intersected by at most n/r lines of H.
39
The sampler
 Goal – sample line. Once the lines are sampled their
intersection (possibly with some modifications) would result
in a triangulation.
 In the previous lecture when we wanted to select a subgroup
S of H, we picked |S| lines out of H with repetitions. This
procedure may result with a subgroup of size smaller than
|S|.
 In order to simplify the calculation, we’ll choose S by
independent Bernoulli trials. We fix the probability of p=s/n,
and include a line from H in S using that probability.
 For this point on, let s=|S|
40
Why sampling and triangulation does
not work?
 In a previous lecture we already saw a sampling algorithm
that outputs a triangulation that with probability close to 1
non of the triangles in the outputted triangulation is
intersected by more than
where is some
constant (Weak Cutting Lemma)
 We will demonstrate that a similar statement with a
is not generally true. So our method would have to include
more than just sampling in order to achieve the desired
triangulation.
 To see that we’ll examine a simpler case than the plane – 1dimensional situation.
41
The 1-D case
 H is a set of n points.
 A generalized triangle in this case is the segment between two
points.
 We would like to count the longest consecutive count of
unselected points of H. Let k be our target number for such a
count.
 In example below, n=30 and s=15 so p=0.5. The black dots
are selected points (In this example there are exactly 15
selected dots, but it doesn’t have to be this way)
42
k consecutive unselected points
 We’ll show that for a fixed k it is very likely that k
consecutive unselected points show up in a sequence of n
points, where n is sufficiently large
 For simplicity assume n is divisible by k and divide our n
points in to blocks the length of k. Notice that this is an even
more restrictive case than a general consecutive count.
 In each such a block there is a probability of
for the entire block to be unselected. Therefore the
probability of not obtaining any block as entirely unselected
is
and this probability is
exponentially small for k proportional to
43
Why is it exponentially small?
 p is a constant, let p=0.5. The probability becomes:
 Now set k=0.5logn:
 And that exponentially small as n goes to infinity.
44
Conclusion of 1D
 A sequence of k that is proportional to log of n is very likely
to appear using this sampling method. This is why this
method is not good enough to get the 1/r cutting we require.
 Of course, in the 1D case we can define a sampler that selects
every n/s point. However it is not clear how to transform
this selection to the plane (Where the lines are not ordered).
 
45
2nd attempt – Two level decomposition
 General plan: Instead of trying to improve the sampling
algorithm to generate a 1/r cutting from the start:
1. we’ll take a random sample (obtained as before with
probability p=s/n independent Bernoulli trials)
2. We just saw that this is not good enough, so we’ll take
every triangle that intersects more lines than required and
divide it so each new triangle will be good for our 1/r
cutting.
 Finally, we’ll show that the total number of triangles is
46
But first, more definitions
 Let T be the collection of triangles after triangulating the
random sample.
 Let
denote the set of lines of H that intersect the
triangle
and
number of lines that
intersect.
 Let the excess of
be
 Note: if
so the triangle
can be included in the 1/r cutting. If
we subdivide it
to a collection of finer triangles.
47
Subdivision of triangles with excess
greater than 1
 We consider the arrangement of
(instead of H), and
we construct a
cutting for it. (How? We’ll see later...)
 We intersect this triangulation with . This process may
produce triangles, but also quadrilaterals, pentagons, and
hexagons. Each of these convex polygons is subdivided into
triangles.
 Each triangle in the
cutting is intersected by at most
triangles, and therefore are valid for the 1/r cutting.
48
Example:
Up – the process of subdividing a triangle with excess
greater than 1.
Left – An example for a hexagon created by combining
the two triangulations.
49
But How? The suboptimal cutting
lemma.
 As of now, we can’t subdivide
because the cutting lemma is
not valid yet. So let’s use a weaker claim:
 The Suboptimal Cutting Lemma
for every finite collection of lines and any u>1, there exists a
1/u cutting consisting of at most:
triangles, where K is a suitable constant.
 We’ll prove the suboptimal cutting lemma later.
50
Using the suboptimal cutting lemma
 If we use the suboptimal cutting lemma for producing the
cuttings, we can estimate that the number of triangles
in the 1/r cutting is bounded by:
 Note that we get one if the excess is less or equal to 1. If the
excess is greater than 1 we get the bound of the yet unproven
suboptimal cutting lemma times 4.
 The times 4 is for the triangulation of a
possible hexagon as shown here:
51
Using the suboptimal cutting lemma (II)
We have two goals now:
1. Prove the suboptimal cutting lemma.
2. Show that although we typically have triangles
as large as about log(r) the expected number of triangles in
T with excess t or larger decreases exponentially as a
function of t.

52
We’ll prove the first claim by demonstrating a triangulation
method that achieves that bound – the vertical
decomposition (next slide). The second claim would be
included in the proof of the Cutting Lemma.
Vertical decomposition - reminder
 We erect vertical segments upwards and downwards from
each vertex in the arrangement of S and extends them until
they meet another line (or all the way to infinity)
 We get trapezoids, but those can be split into 2 triangles.
 Let Trpz(S) denote the set of generalized trapezoids in the
vertical decomposition of S.
53
Proposition:
Trapezoids with large excess are rare
Proposition: Let H be a fixed set of n lines in general
position, let p=r/n, where
, let S be a
random sample drawn from H by independent Bernoulli
trials with success probability of p, and let
be a real
parameter. Let
denote the set of trapezoids in
with excess at least t
Then the expected number of trapezoids in
is bounded as follows:
For a suitable constant C.
54
Notes for this proposition
 If r is not smaller than n/2 than we can easily find a cutting,
as we saw in the cutting lemma lecture it’s only gets
interesting when r is proportional to logn.
 Reminder: How to construct such a cutting? If r is
proportional to n one can select all lines in the H, as shown
before this creates the desired subdivision.
 We use this limitation to find a finer cutting only because this
is needed in the proof of this proposition (but we’ll get there
much later)
55
Proof of the suboptimal cutting lemma
 Set
for a sufficiently large constant
A.
 Assume that the trapezoid proposition holds – and choose a
sample S accordingly.
 Activate the proposition twice – once for
and once
for
56
First case
 For
by using the proposition:
for some constant B.
 Also note that
 Finally:
57
Second case
 For
by using the proposition:
for some constant C. We’ll group the constants together:
 Since:
there exists A large enough s.t
58
Proof of the Suboptimal Cutting
Lemma – Conclusion
 From these two cases and the linearity of expectation we get:
 There exist a sample were both parts of the sum are smaller
than 2/3. Notice that the second part is at least 1 and
therefore 0.
 So there exists a sample S with both
and
59
Proof of the Suboptimal Cutting
Lemma – Conclusion (II)
 That means that there exists a 1/u cutting into
triangles.
 Suboptimal Cutting Lemma proof is complete.
60
Proof of the Cutting Lemma (I)
 Reminder: To produce a 1/r cutting we pick a sample S with
probability p=r/n as described before. We let T be the
collection of triangles of the vertical decomposition of S
(actually – trapezoids, but that has no effect). For every
triangle t in T with excess greater than 1 we refine the
triangulation with a
cutting. We’ll estimate the
expected number of triangles:
61
Proof of the Cutting Lemma (II)
 As
and
 Since each part of the sum is at least 1. Also group by
order
62
‘s
Proof of the Cutting Lemma (III)
 Notice that: if
then
so
 Also, the summing over triangles with
is less then summing over
63
Proof of the Cutting Lemma (IV)
 Using the suboptimal cutting lemma and the trapezoid proposition
for some M large enough
 The ratio test shows that the series converges to some
number. So we get that the total number of triangles is
as needed.
64
Just one more
Now we just need to prove the trapezoid proposition, for
that we’ll need some more definitions:

is the set of all trapezoids that can ever appear in the vertical
decomposition for some subset S of H (including the empty
set).

is the set of lines of H incident to at least one vertex
of . Notice that
(all the possible cases, up to
symmetry are shown in the next slide)
65
D cases
 The following shows that
66
Some Properties of Vertical
Decomposition
 C0: As we saw:
. Moreover, any subset S of H is
the defining set for at most a constant number of trapezoids
in Reg.
 C1: For any
we have
(the
defining set must be present) and
(no
intersecting lines may be present).
 C2: For any
and any
s.t
and
we have
 C3: For every
we have
67
Notes about the properties
 C3: Think of the iterative process of adding one vertical line
at a time. Each splits an existing region in two. Since the lines
are in general position the number of intersections is about
 C2 is the most interesting - it means that the vertical
decomposition is defined “locally”:
is present in
whenever it is not
excluded for simple local reasons.
68
Proof of the Trapezoid Proposition
 First case
:
From C3 we get
is the sum of independent random variables.
and so we get:
69
where |S|
First case continued:
For now we have:
for
some constant M independent of t.
On the other hand
by
combining these two and setting C=2M we get:
as required.
This concludes the case of
70
Second case
 For now we assume
. Let
be a random sample
as defined before with probability p. From C1 and C2 we get:
(a trapezoid appears iff all
his defining lines are selected into S, and no line that
intersects it are selected into S).
 Let the set of all possible trapezoids with excess at least t:
 The expected number of trapezoids can be written as:
71
Second case continued:
 We now define p’=p/t and since t>1 we get: p’<p
 Let S’ be an a sample drawn of H with probability p’
 From what we saw earlier:
 Also:
 And finally:
The last transition holds from limiting the sum.
72
Second case continued (II):
Now be expanding with p:
Where R is defined:
As we’ll see in the next slide R is bounded by:
Finally we derive:
for a sufficiently large constant C (and than choosing the max
between it and the first case) the proposition is proved!
73
Lower bound for R
 And since
 For every real p:
 Since p’<p<0.5:
74
and
we get:
Lower bound for R (II)
 By the last slide we see that:
 Since R is defined:
it will always be
75
Cutting Lemma conclusion
 Since the trapezoid decomposition lemma is proved the
entire proof is done.
 As mentioned before this Lemma can be generalized to
higher dimensions, and the proof we have just presented can
also be generalized to prove it (but it won’t be covered now).
 Cutting Lemma for Higher Dimensions: Let d>0 be a fixed
integer, let H be a set of n hyperplanes in
and let r be a
parameter
.Then there exists a 1/r cutting for H the
size of
generalized simplexes s.t the interior of each simplex
is intersected by at most n/r hyperplanes of H.
76