Bipartite Diameter and Other Measures Under Translation Boris - - PowerPoint PPT Presentation

bipartite diameter and other measures under translation
SMART_READER_LITE
LIVE PREVIEW

Bipartite Diameter and Other Measures Under Translation Boris - - PowerPoint PPT Presentation

Bipartite Diameter and Other Measures Under Translation Boris Aronov, Omrit Filtser , Matthew J. Katz, and Khadijeh Sheikhan March 14, 2019 Similarity between two sets of points Goal : Determining the similarity between two sets of points. ? B.


slide-1
SLIDE 1

Bipartite Diameter and Other Measures Under Translation

Boris Aronov, Omrit Filtser, Matthew J. Katz, and Khadijeh Sheikhan March 14, 2019

slide-2
SLIDE 2

Similarity between two sets of points

Goal: Determining the similarity between two sets of points.

?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 2 / 22

slide-3
SLIDE 3

Similarity between two sets of points

Goal: Determining the similarity between two sets of points.

◮ A well investigated problem in computational geometry.

?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 2 / 22

slide-4
SLIDE 4

Similarity between two sets of points

Goal: Determining the similarity between two sets of points.

◮ A well investigated problem in computational geometry. ◮ Problem: Sometimes, a bipartite measure is meaningless,

unless one of the sets undergoes some transformation.

?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 2 / 22

slide-5
SLIDE 5

Similarity between two sets of points

Goal: Determining the similarity between two sets of points.

◮ A well investigated problem in computational geometry. ◮ Problem: Sometimes, a bipartite measure is meaningless,

unless one of the sets undergoes some transformation. This paper: Find a translation which minimizes some bipartite measure.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 2 / 22

slide-6
SLIDE 6

Similarity between two sets of points

Goal: Determining the similarity between two sets of points.

◮ A well investigated problem in computational geometry. ◮ Problem: Sometimes, a bipartite measure is meaningless,

unless one of the sets undergoes some transformation. This paper: Find a translation which minimizes some bipartite measure.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 2 / 22

slide-7
SLIDE 7

Bipartite measures under translation

A = {a1, . . . , an} and B = {b1, . . . , bm} – two sets of points in Rd.

Problem

Find a translation t∗ that minimizes some bipartite measure of A and B + t over all translations t.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 3 / 22

slide-8
SLIDE 8

Bipartite measures under translation

A = {a1, . . . , an} and B = {b1, . . . , bm} – two sets of points in Rd.

Problem

Find a translation t∗ that minimizes some bipartite measure of A and B + t over all translations t.

Remarks

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 3 / 22

slide-9
SLIDE 9

Bipartite measures under translation

A = {a1, . . . , an} and B = {b1, . . . , bm} – two sets of points in Rd.

Problem

Find a translation t∗ that minimizes some bipartite measure of A and B + t over all translations t.

Remarks

◮ For the sake of simplicity, we assume that m = n.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 3 / 22

slide-10
SLIDE 10

Bipartite measures under translation

A = {a1, . . . , an} and B = {b1, . . . , bm} – two sets of points in Rd.

Problem

Find a translation t∗ that minimizes some bipartite measure of A and B + t over all translations t.

Remarks

◮ For the sake of simplicity, we assume that m = n. ◮ This class of problems naturally extends to other types of

transformations, such as rotations, rigid motions, etc.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 3 / 22

slide-11
SLIDE 11

Some bipartite measure?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 4 / 22

slide-12
SLIDE 12

Some bipartite measure?

When comparing two sets of points A and B of the same size:

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 4 / 22

slide-13
SLIDE 13

Some bipartite measure?

When comparing two sets of points A and B of the same size:

◮ Congruence testing: decide if there exists a transformation that

maps A exactly or approximately into B.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 4 / 22

slide-14
SLIDE 14

Some bipartite measure?

When comparing two sets of points A and B of the same size:

◮ Congruence testing. ◮ RMS distance: minimize the sum of squares of distances in a

perfect matching between A and B.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 4 / 22

slide-15
SLIDE 15

Some bipartite measure?

When comparing two sets of points A and B of the same size:

◮ Congruence testing. ◮ RMS distance.

When comparing two sets of points A and B of different sizes:

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 4 / 22

slide-16
SLIDE 16

Some bipartite measure?

When comparing two sets of points A and B of the same size:

◮ Congruence testing. ◮ RMS distance.

When comparing two sets of points A and B of different sizes:

◮ Hausdorff distance: the maximum of the distances from a point in

each of the sets to the nearest point in the other set. Huttenlocher,Kedem, Sharir: ˜ O(n3) in 2D.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 4 / 22

slide-17
SLIDE 17

Some bipartite measure?

When comparing two sets of points A and B of the same size:

◮ Congruence testing. ◮ RMS distance.

When comparing two sets of points A and B of different sizes:

◮ Hausdorff distance: ˜

O(n3) in 2D.

◮ Maximum overlap between the convex hulls of the sets A and B.

de Berg et al.:O(n log n) in 2D, Ahn et al.: ˜ O(n3) in 3D.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 4 / 22

slide-18
SLIDE 18

Some bipartite measure?

When comparing two sets of points A and B of the same size:

◮ Congruence testing. ◮ RMS distance.

When comparing two sets of points A and B of different sizes:

◮ Hausdorff distance: ˜

O(n3) in 2D.

◮ Maximum overlap between the convex hulls of the sets A and B.

de Berg et al.:O(n log n) in 2D, Ahn et al.: ˜ O(n3) in 3D.

All the above measures (under various geometric transformations) were widely investigated in the literature.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 4 / 22

slide-19
SLIDE 19

Our results

The main bipartite measures that we consider are:

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 5 / 22

slide-20
SLIDE 20

Our results

The main bipartite measures that we consider are:

◮ diameter – the distance between the farthest bichromatic

pair, i.e. max{a − b | (a, b) ∈ A × B}.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 5 / 22

slide-21
SLIDE 21

Our results

The main bipartite measures that we consider are:

◮ diameter – max{a − b | (a, b) ∈ A × B}. ◮ uniformity – the difference between the bipartite diameter

and the distance between the closest bichromatic pair, i.e. diam(A, B) − min{a − b | (a, b) ∈ A × B}.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 5 / 22

slide-22
SLIDE 22

Our results

The main bipartite measures that we consider are:

◮ diameter – max{a − b | (a, b) ∈ A × B}. ◮ uniformity – diam(A, B) − min{a − b | (a, b) ∈ A × B}. ◮ union width – the width of A ∪ B, where the width of a set

  • f points in the plane is the smallest distance between a pair
  • f parallel lines, such that the closed strip between the lines

contains the entire set.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 5 / 22

slide-23
SLIDE 23

Our results

The main bipartite measures that we consider are:

◮ diameter – max{a − b | (a, b) ∈ A × B}. ◮ uniformity – diam(A, B) − min{a − b | (a, b) ∈ A × B}. ◮ union width – the width of A ∪ B. ◮ red-blue width – ...

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 5 / 22

slide-24
SLIDE 24

Our results

The main bipartite measures that we consider are:

◮ diameter – max{a − b | (a, b) ∈ A × B}. ◮ uniformity – diam(A, B) − min{a − b | (a, b) ∈ A × B}. ◮ union width – the width of A ∪ B. ◮ red-blue width – ...

Surprisingly, all of these measures (under translation) were not investigated previously in the literature.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 5 / 22

slide-25
SLIDE 25

Our results

measure dimension running time d = 2 O(n log n) d = 3 O(n log2 n) diameter d > 3 (fixed) O(n2) uniformity d = 2 O(n9/4+ε) union width d = 2 O(n log n) d = 3 O(n2)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 6 / 22

slide-26
SLIDE 26

Our results

measure dimension running time d = 2 O(n log n) d = 3 O(n log2 n) diameter d > 3 (fixed) O(n2) uniformity d = 2 O(n9/4+ε) union width d = 2 O(n log n) d = 3 O(n2)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 6 / 22

slide-27
SLIDE 27

Diameter

A = {a1, . . . , an} and B = {b1, . . . , bm} – two sets of points in Rd. diam(A, B) = max{a − b | (a, b) ∈ A × B}

Problem (Bipartite Diameter under Translation)

Find a translation t∗ such that for any translation t, diam(A, B + t∗) ≤ diam(A, B + t).

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 7 / 22

slide-28
SLIDE 28

Diameter

P = {a − b | (a, b) ∈ A × B}

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 8 / 22

slide-29
SLIDE 29

Diameter

P = {a − b | (a, b) ∈ A × B}

◮ The set of all possible translations taking a point of B to a

point of A.

b1 b2 a1 a2 a1 − b2 a2 − b2 a1 − b1 a2 − b1

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 8 / 22

slide-30
SLIDE 30

Diameter

P = {a − b | (a, b) ∈ A × B}

◮ The set of all possible translations taking a point of B to a

point of A.

◮ Clearly, |P| = O(n2).

b1 b2 a1 a2 a1 − b2 a2 − b2 a1 − b1 a2 − b1

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 8 / 22

slide-31
SLIDE 31

Diameter

P = {a − b | (a, b) ∈ A × B}

Claim

Given a point t, the radius of the minimum enclosing ball of P centered at t is equal to diam(A, B + t).

t r =diam(A, B + t) P

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 9 / 22

slide-32
SLIDE 32

Diameter

P = {a − b | (a, b) ∈ A × B}

Claim

Given a point t, the radius of the minimum enclosing ball of P centered at t is equal to diam(A, B + t).

Proof.

This radius is at most max

(a−b)∈P (a − b) − t =

max

(a,b)∈A×B a − (b + t) = diam(A, B + t).

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 9 / 22

slide-33
SLIDE 33

Diameter

P = {a − b | (a, b) ∈ A × B}

Claim

Given a point t, the radius of the minimum enclosing ball of P centered at t is equal to diam(A, B + t).

Corollary

The optimal translation t∗ minimizing the bipartite diameter coincides with the center of the minimum enclosing ball of P.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 9 / 22

slide-34
SLIDE 34

Diameter: Algorithm (naive implementation)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 10 / 22

slide-35
SLIDE 35

Diameter: Algorithm (naive implementation)

◮ Compute the set of translations P.

b1 b2 a1 a2 a1 − b2 a2 − b2 a1 − b1 a2 − b1

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 10 / 22

slide-36
SLIDE 36

Diameter: Algorithm (naive implementation)

◮ Compute the set of translations P. ◮ Find the center c of the minimum enclosing ball of P.

b1 b2 a1 a2 a1 − b2 a2 − b2 a1 − b1 a2 − b1

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 10 / 22

slide-37
SLIDE 37

Diameter: Algorithm (naive implementation)

◮ Compute the set of translations P. ◮ Find the center c of the minimum enclosing ball of P. ◮ Translating B by c minimizes the diameter.

b1 b2 a1 a2 a1 − b2 a2 − b2 a1 − b1 a2 − b1

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 10 / 22

slide-38
SLIDE 38

Diameter: Running time

The minimum enclosing ball can be computed in:

◮ linear time using Megiddo’s (’83) algorithm, or ◮ expected linear time using Welzl’s (’91) simpler randomized

algorithm.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 11 / 22

slide-39
SLIDE 39

Diameter: Running time

The minimum enclosing ball can be computed in:

◮ linear time using Megiddo’s (’83) algorithm, or ◮ expected linear time using Welzl’s (’91) simpler randomized

algorithm. ⇒ O(n2)-time solution.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 11 / 22

slide-40
SLIDE 40

Diameter: Running time

The minimum enclosing ball can be computed in:

◮ linear time using Megiddo’s (’83) algorithm, or ◮ expected linear time using Welzl’s (’91) simpler randomized

algorithm. ⇒ O(n2)-time solution. BUT, in 2D and 3D, we can do better! In fact, computing the minimum enclosing ball of P in 2D and 3D (without computing P explicitly) can be done in near-linear time...

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 11 / 22

slide-41
SLIDE 41

Diameter

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P implicitly.

B A P

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 12 / 22

slide-42
SLIDE 42

Diameter

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P implicitly.

Observations

B A P

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 12 / 22

slide-43
SLIDE 43

Diameter

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P implicitly.

Observations

  • 1. We only need to look at the convex hull (CH) of P.

B A P CH(P)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 12 / 22

slide-44
SLIDE 44

Diameter

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P implicitly.

Observations

  • 1. We only need to look at the convex hull (CH) of P.
  • 2. P is the Minkowski sum of A and −B, i.e. P = A ⊕ −B.

B A P

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 12 / 22

slide-45
SLIDE 45

Diameter in 2D

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P in 2D.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 13 / 22

slide-46
SLIDE 46

Diameter in 2D

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P in 2D.

A fact from the textbook

For points in 2D, the size of CH(A ⊕ −B) is O(n), and it can be constructed in O(n) time from CH(A) and CH(B) using the well-known rotating calipers method...

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 13 / 22

slide-47
SLIDE 47

Diameter in 2D

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P in 2D.

A fact from the textbook

For points in 2D, the size of CH(A ⊕ −B) is O(n), and it can be constructed in O(n) time from CH(A) and CH(B) using the well-known rotating calipers method... ⇒ O(n log n)-time solution for points in 2D!

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 13 / 22

slide-48
SLIDE 48

Diameter in 3D

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P in 3D.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 14 / 22

slide-49
SLIDE 49

Diameter in 3D

P = {a − b | (a, b) ∈ A × B} Goal: compute the minimum enclosing ball of P in 3D. Idea: The minimum enclosing ball is an LP-type problem ⇒ adapt Clarkson’s (’95) scheme for solving LP-type problems.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 14 / 22

slide-50
SLIDE 50

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found:

P (points)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-51
SLIDE 51

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 1: Pick a random sample R of P of size 4n.

P (points) R (random sample)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-52
SLIDE 52

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 2: Compute the minimum enclosing ball S of R ∪ X.

P (points) R (random sample)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-53
SLIDE 53

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 3: Find the set of violators V . If |V | ≥ 2n, go to 1.

P (points) R (random sample) V (violators)

|V | ≥ 2n, “bad” iteration :(

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-54
SLIDE 54

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 1: Pick a random sample R of P of size 4n.

P (points) R (random sample)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-55
SLIDE 55

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 2: Compute the minimum enclosing ball S of R ∪ X.

P (points) R (random sample)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-56
SLIDE 56

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 3: Find the set of violators V . If |V | ≥ 2n, go to 1.

P (points) R (random sample) V (violators)

|V | < 2n, “good” iteration :)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-57
SLIDE 57

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 4: If V = ∅, then X ← X ∪ V and go to 1. Else, return S.

P (points) R (random sample) V (violators) X ← X ∪ V

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-58
SLIDE 58

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 1: Pick a random sample R of P of size 4n.

P (points) X ← X ∪ V

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-59
SLIDE 59

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 1: Pick a random sample R of P of size 4n.

P (points) R (random sample) X ← X ∪ V

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-60
SLIDE 60

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 2: Compute the minimum enclosing ball S of R ∪ X.

P (points) R (random sample) X ← X ∪ V

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-61
SLIDE 61

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 3: Find the set of violators V . If |V | ≥ 2n, go to 1.

P (points) R (random sample) V (violators) X ← X ∪ V

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-62
SLIDE 62

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 4: If V = ∅, then X ← X ∪ V and go to 1. Else, return S.

P (points) R (random sample) X ← X ∪ V

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-63
SLIDE 63

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 4: If V = ∅, then X ← X ∪ V and go to 1. Else, return S.

P (points) X ← X ∪ V

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-64
SLIDE 64

Diameter in 3D: Algorithm

X – an empty set of points. Repeat until the minimum enclosing ball is found: 4: If V = ∅, then X ← X ∪ V and go to 1. Else, return S.

P (points) X ← X ∪ V

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 15 / 22

slide-65
SLIDE 65

Diameter in 3D: number of iterations

How many iterations?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 16 / 22

slide-66
SLIDE 66

Diameter in 3D: number of iterations

How many iterations?

◮ The number of good iterations cannot exceed five! (why?)

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 16 / 22

slide-67
SLIDE 67

Diameter in 3D: number of iterations

How many iterations?

◮ The number of good iterations cannot exceed five! (why?)

⇒ |X| = O(n).

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 16 / 22

slide-68
SLIDE 68

Diameter in 3D: number of iterations

How many iterations?

◮ The number of good iterations cannot exceed five! (why?)

⇒ |X| = O(n).

◮ By Clarkson’s analysis:

In each iteration the expected size of V is n.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 16 / 22

slide-69
SLIDE 69

Diameter in 3D: number of iterations

How many iterations?

◮ The number of good iterations cannot exceed five! (why?)

⇒ |X| = O(n).

◮ By Clarkson’s analysis:

In each iteration the expected size of V is n. ⇒ By Markov’s inequality, Pr(|V | ≥ 2n) ≤ 1

2

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 16 / 22

slide-70
SLIDE 70

Diameter in 3D: number of iterations

How many iterations?

◮ The number of good iterations cannot exceed five! (why?)

⇒ |X| = O(n).

◮ By Clarkson’s analysis:

In each iteration the expected size of V is n. ⇒ By Markov’s inequality, Pr(|V | ≥ 2n) ≤ 1

2

⇒ Expected number of bad iterations before a good one is O(1)!

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 16 / 22

slide-71
SLIDE 71

Diameter in 3D: number of iterations

How many iterations?

◮ The number of good iterations cannot exceed five! (why?)

⇒ |X| = O(n).

◮ By Clarkson’s analysis:

In each iteration the expected size of V is n. ⇒ By Markov’s inequality, Pr(|V | ≥ 2n) ≤ 1

2

⇒ Expected number of bad iterations before a good one is O(1)! ⇒ Expected number of iterations is constant!

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 16 / 22

slide-72
SLIDE 72

Diameter in 3D: implementing an iteration

What is the running time for one iteration?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-73
SLIDE 73

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 1: Pick a random sample R of P of size 4n.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-74
SLIDE 74

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 1: Pick a random sample R of P of size 4n.

◮ Repeatedly pick random points a ∈ A and b ∈ B

and return a − b.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-75
SLIDE 75

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 2: Compute the minimum enclosing ball S of R ∪ X.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-76
SLIDE 76

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 2: Compute the minimum enclosing ball S of R ∪ X.

◮ Invoke a standard minimum-ball algorithm on O(n) points,

requiring O(n) expected time.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-77
SLIDE 77

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 3: Find the set of violators V .

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-78
SLIDE 78

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 3: Find the set of violators V .

◮ ?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-79
SLIDE 79

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 3: Find the set of violators V . Idea: We solve the following problem.

Problem

Given two sets A and B, each of n points in R3, a distance r and a parameter k, report all the pairs of points a ∈ A, b ∈ B with a − b > r, if there are at most k such pairs. Otherwise, return ”TOO MANY”.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-80
SLIDE 80

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 3: Find the set of violators V . Idea: We solve the following problem.

Problem

Given two sets A and B, each of n points in R3, a distance r and a parameter k, report all the pairs of points a ∈ A, b ∈ B with a − b > r, if there are at most k such pairs. Otherwise, return ”TOO MANY”.

◮ Expected running time O((n + k) log2 n).

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-81
SLIDE 81

Diameter in 3D: implementing an iteration

What is the running time for one iteration? 3: Find the set of violators V . Idea: We solve the following problem.

Problem

Given two sets A and B, each of n points in R3, a distance r and a parameter k, report all the pairs of points a ∈ A, b ∈ B with a − b > r, if there are at most k such pairs. Otherwise, return ”TOO MANY”.

◮ Expected running time O((n + k) log2 n).

⇒ O(n log2 n)-time solution for points in 3D!

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 17 / 22

slide-82
SLIDE 82

Uniformity

A = {a1, . . . , an} and B = {b1, . . . , bm} – two sets of points in Rd. uni(A, B) = diam(A, B) − min{a − b | (a, b) ∈ A × B}

Problem (Uniformity under Translation)

Find a translation t∗ such that for any translation t, uni(A, B + t∗) ≤ uni(A, B + t).

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 18 / 22

slide-83
SLIDE 83

Uniformity

P = {a − b | (a, b) ∈ A × B}

P

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 19 / 22

slide-84
SLIDE 84

Uniformity

P = {a − b | (a, b) ∈ A × B}

Claim

The optimal translation t∗ minimizing the uniformity coincides with the center of the minimum-width annulus containing P.

P

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 19 / 22

slide-85
SLIDE 85

Uniformity in 2D

◮ Agarwal and Sharir (’96): The minimum enclosing annulus of

n points in 2D can be computed in O(n3/2+ε) expected time...

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 20 / 22

slide-86
SLIDE 86

Uniformity in 2D

◮ Agarwal and Sharir (’96): The minimum enclosing annulus of

n points in 2D can be computed in O(n3/2+ε) expected time... ⇒ O(n3+ε)-time solution.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 20 / 22

slide-87
SLIDE 87

Uniformity in 2D

◮ Agarwal and Sharir (’96): The minimum enclosing annulus of

n points in 2D can be computed in O(n3/2+ε) expected time... ⇒ O(n3+ε)-time solution.

Claim

The minimum enclosing annulus of n points in 2D — with only O(√n) extreme points — can be computed in O(n9/8+ε) expected time.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 20 / 22

slide-88
SLIDE 88

Uniformity in 2D

◮ Agarwal and Sharir (’96): The minimum enclosing annulus of

n points in 2D can be computed in O(n3/2+ε) expected time... ⇒ O(n3+ε)-time solution.

Claim

The minimum enclosing annulus of n points in 2D — with only O(√n) extreme points — can be computed in O(n9/8+ε) expected time. ⇒ O(n9/4+ε)-time solution!

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 20 / 22

slide-89
SLIDE 89

Thank You!

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 21 / 22

slide-90
SLIDE 90

Open questions

◮ Consider other types of transformations?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-91
SLIDE 91

Open questions

◮ Consider other types of transformations? ◮ Diameter:

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-92
SLIDE 92

Open questions

◮ Consider other types of transformations? ◮ Diameter:

◮ We showed near linear time algorithms in 2D and 3D.

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-93
SLIDE 93

Open questions

◮ Consider other types of transformations? ◮ Diameter:

◮ We showed near linear time algorithms in 2D and 3D. ◮ Can we also obtain a near linear algorithm in higher

dimensions (O(n logO(d) n))?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-94
SLIDE 94

Open questions

◮ Consider other types of transformations? ◮ Diameter:

◮ We showed near linear time algorithms in 2D and 3D. ◮ Can we also obtain a near linear algorithm in higher

dimensions (O(n logO(d) n))?

◮ Uniformity:

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-95
SLIDE 95

Open questions

◮ Consider other types of transformations? ◮ Diameter:

◮ We showed near linear time algorithms in 2D and 3D. ◮ Can we also obtain a near linear algorithm in higher

dimensions (O(n logO(d) n))?

◮ Uniformity:

◮ Consider minimum ratio instead of minimum difference?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-96
SLIDE 96

Open questions

◮ Consider other types of transformations? ◮ Diameter:

◮ We showed near linear time algorithms in 2D and 3D. ◮ Can we also obtain a near linear algorithm in higher

dimensions (O(n logO(d) n))?

◮ Uniformity:

◮ Consider minimum ratio instead of minimum difference? ◮ Higher dimensions?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-97
SLIDE 97

Open questions

◮ Consider other types of transformations? ◮ Diameter:

◮ We showed near linear time algorithms in 2D and 3D. ◮ Can we also obtain a near linear algorithm in higher

dimensions (O(n logO(d) n))?

◮ Uniformity:

◮ Consider minimum ratio instead of minimum difference? ◮ Higher dimensions?

◮ Width:

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-98
SLIDE 98

Open questions

◮ Consider other types of transformations? ◮ Diameter:

◮ We showed near linear time algorithms in 2D and 3D. ◮ Can we also obtain a near linear algorithm in higher

dimensions (O(n logO(d) n))?

◮ Uniformity:

◮ Consider minimum ratio instead of minimum difference? ◮ Higher dimensions?

◮ Width:

◮ For width in 3D (without translation) there is an

O(n3/2+ǫ)-time algorithm (Agarwal and Sharir).

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22

slide-99
SLIDE 99

Open questions

◮ Consider other types of transformations? ◮ Diameter:

◮ We showed near linear time algorithms in 2D and 3D. ◮ Can we also obtain a near linear algorithm in higher

dimensions (O(n logO(d) n))?

◮ Uniformity:

◮ Consider minimum ratio instead of minimum difference? ◮ Higher dimensions?

◮ Width:

◮ For width in 3D (without translation) there is an

O(n3/2+ǫ)-time algorithm (Agarwal and Sharir).

◮ Our algorithm (for width in 3D under translation) runs in

O(n2) time. Can we do better?

  • B. Aronov, O. Filtser, M. J. Katz, K. Sheikhan

Bipartite Diameter and Other Measures Under Translation 22 / 22