200 likes | 451 Vues
Randomized Quicksort Randomized Global Min Cut. CSL758 Instructors: Naveen Garg Kavitha Telikepalli Scribe: Manish Singh Vaibhav Rastogi February 7 & 11, 2008. The Quick Sort Problem. To sort a given set of numbers
E N D
Randomized QuicksortRandomized Global Min Cut CSL758 Instructors: Naveen Garg Kavitha Telikepalli Scribe: Manish Singh Vaibhav Rastogi February 7 & 11, 2008.
The Quick Sort Problem • To sort a given set of numbers • In traditional quick sort algorithm, we pick a particular index element as pivot for splitting. • Worst Case: O(n2) • Average Case: O(n log n) • A good pivot can be selected using median finding algorithm but the total complexity will again be O(n2). • So, what if we pick a random element uniformly as pivot and do the partition. We will show that this takes expected O(n log n) time.
Randomized Quick Sort algorithm • Assuming all elements are distinct • We pick a random element x as the pivot and partition the input set S into two sets L and R such: • L = numbers less than x • R = numbers greater than x • Recursively sort L and R. • Return LxR
Analysis of Randomized Quick Sort • The running time of this algorithm is variable. • Running time = # of comparisons in this algorithm. • Expected Running Time, • Let S be the sorted sequence of the n input numbers. S = • Let Xij = 1 if Si and Sj are compared in the algo = 0 otherwise • Running time = # of comparisons = S1 S2 Si Sj Sn-1 Sn
Analysis cont… • The expected running time = = • Now, = 1 * Pr[Si and Sj are compared in our alg.] + 0 * Pr[Si and Sj are not compared] • Suppose we have a set of numbers: 2, 7, 15, 18, 19, 23, 35 • In this 18 and 19 will always be compared. • 2 and 35 will be compared only if compared at root.
Analysis cont… • Pr[Si and Sj are compared in our algo] = Pr[the first element chosen as pivot in set {si,si+1,….,sj} is either si or sj]. • To elements get compared only if they have ancestor relationship in the tree. • Pr[Picking Si or Sj] = S1 S2 Si Sj Sn-1 Sn 1st pivot 3rd pivot 2nd pivot 4th pivot 5th pivot
Analysis cont… • Thus the expected runtime: • This algorithm will always give the right answer though the running time may be different. This is an example of Las Vegas algorithms.
The global min cut problem • Input: A connected undirected multigraph G(V, E). • Output: A minimum cardinality subset of edges whose removal makes G disconnected. b a c d x
First Idea • Recall the max flow problem • An s-t min cut is the min cut which makes s and t disconnected. • Run a max flow algorithm on G for all (s, t) pairs and return the smallest of the output cuts. • Time: • Since G is undirected we could fix s and iterate t over all other vertices. • Time:
Another idea • Reduce G to G1with one vertex less such that • Contract a random edge. Each edge is contracted with probability c d a b Contract a random edge G G1 c d (a,b)
Observation 1 • Any cut of G1is also a cut of G, with a and b on the same side. • If mincut of G1< mincut of G • This mincut of G1 should be the mincut of G, by the above observation. • Thus, mincut of G1 ≥ mincut of G.
Observation 2 • Fix our favorite global mincut C, with • Let e be the edge contracted. • If e is not in C then • C is also a cut in G1. • Using Observation 1, mincut of G1 = mincut of G. e
Observation 2 cont… which is quite a high probability.
The algorithm Repeat Choose an edge of G uniformly at random and contract it. Let the resulting graph be called G. Until there are only two vertices in G. Output the set of edges between the two vertices as our candidate mincut.
Bounding the probability of correctness • Let C be our favorite mincut and F be the cut output by the algorithm. • Pr[F is C]=Pr[no edge of C is contracted in any iteration]. • Let Eibe the event that no edge of C got contracted in the ithiteration.
Bounding the probability of correctness • We know that, • Because C remains the mincut we can also bound the conditional following conditional probabilities as we bounded Thus,
Bounding the probability of correctness • is the product of these conditional probabilities. So we have, • This probability bound is too low. It can be improved by repeating the algorithm N times and returning the least cut.
Running time • Choosing a random edge. • Choose a vertex u with probability • Now choose one of u’s with probability • Time: O(n).
Running time • Contracting edges: Use adjacency matrix. • Add the rows and columns of u and v. • O(n)time. u+v u v u u+v v
Running time • The basic algorithm has iterations. • Each iteration in the basic algorithm involves selecting and contracting a random edge. O(n) time. • The basic algorithm is repeated times. • Therefore running time =