#### Transcript Foundations of Algorithms, Fourth Edition

Foundations of Algorithms, Fourth Edition Richard Neapolitan, Kumarss Naimipour Chapter 2 Divide-and-Conquer Divide and Conquer • In this approach a problem is divided into sub-problems and the same algorithm is applied to every subproblem ( often this is done recursively) • Examples – Binary Search (review algorithm in book) – Mergesort (review algorithm in book) – Quicksort Figure 2.1 : The steps down by a human when searching with Binary Search. (Note: x = 18) Complexity of Binary Search Since this and many other divide and conquer algorithms are recursive you will recall that we can determine their complexity using recurrence relations. For Binary Search we have T(n) = T(n/2) + 1 =[T(n/4)+1]+1 = T(n/22) + 2 =[T(n/8+1]+ 2 = T(n/23)+ 3 … =T(n/2k)+k What is k T(n/2 )+k We if we let k get larger until n=2k then we see that k = log2n. Why? Consequently the relation becomes T(n) = T(1) + log2n T(n) = log2n Since n/2k is 1 if they are equal and T(1) =1 MergeSort Recall in this algorithm we divide the array into two equal parts and sort each half prior to merging. The recurrence relation is clearly T(n) = 2T(n/2) + n Recall that Merging is O(n) right? T(n/2) T(n/4) O(n) Figure 2.2: The steps done by a human when sorting with Mergesort. T(n) = 2T(n/2) + n T(n) = 2T(n/2) + n = 2[ 2T(n/22) + n/2] + n = 22T(n/22) + 2n = 2[2T(n/23) + n/22] +2n =23T(n/23) + 3n … =2kT(n/2k) + kn If n=2k then we have T(n) = nT(1) + (log2n)n = n+ nlog2n = O(nlog2n) QuickSort Works in situ! Void quicksort(int low, int high) { int pivot; if (high > low){ partition(low, high, pivot); quicksort(low, pivot-1); quicksort(pivot+1,high); } Figure 2.3: The steps done by a human when sorting with Quicksort. The subarrays are enclosed in rectangles whereas the pivot points are free. Partition Study this carefully void partition (int low, int high, int&pivot) { int I,j, pivotitem; pivotitem = S[low]; // select left item (hmmm) j=low; for (i=low+1; i<=high; i++) if (S[i] < pivotitem){ j++; swap S[i] and S[j]; } pivot= j; swap S[low] and S[pivot]; } There are many ways to write this function! All have a complexity of O(n). Complexity of Quicksort The complexity of this algorithm depends on how good the pivot value selection is . If the value is always in the middle of array then the best case complexity is T(n) = n + 2T(n/2) Which we already have determined is T(n) = n log2 n Worst case for Quicksort This clearly will occur if each pivot value is less than (or greater) all the elements of the array. IE the array is split into 1 and n-1 size pieces. This gives a recurrence relation of T(n) = T(1) + T(n-1) + n-1 Time to sort left array right array partition Worst Case analysis T(n) = T(1) + T(n-1) + n-1 = T(n-1) + n Assume the answer is n(n-1)/2 check it out ! n(n-1)/2 = 0 + (n-1)(n-2)/2 + n-1 = (n-1)(n-2)/2 + 2(n-1)/2 =((n-1)(n-2)+ 2(n-1))/2 = (n-1)(n-2+2)/2 = n(n-1)/2 ☺ Quick Sort Analysis Quicksort’s worst case is θ(n2) Does this mean that quick sort is just as bad as say selection sort, insertion sort and/or bubble sort. No! Its all about average case performance. The average case performance for these three is θ(n2) as well. What is the average case complexity for QS? Average Case Analysis assume prob. pivotpoint is p 𝑛 𝐴 𝑛 = 1 𝐴 𝑝 − 1 + 𝐴(𝑛 − 𝑝) + 𝑛 − 1 𝑛 𝑝=1 𝑛 1 𝐴 𝑛 = 𝑛 𝐴 𝑝 − 1 + 𝐴(𝑛 − 𝑝) + 𝑛 − 1 𝑝=1 2 𝐴 𝑛 = 𝑛 𝑛 𝐴 𝑝−1 +𝑛−1 𝑝=1 See HW 22 p 86 for above conversion Average case continued Multiplying by n 𝑛 𝑛𝐴 𝑛 = 2 (𝑛 − 1)𝐴 𝑛 − 1 = 2 𝐴 𝑝−1 +𝑛 𝑛−1 𝑝=1 𝑛−1 𝑝=1 𝐴 𝑝 − 1 + (𝑛 − 1)(𝑛 − 2) Subtracting these equations we have 𝑛𝐴 𝑛 − 𝑛 − 1 𝐴 𝑛 − 1 = 2𝐴 𝑛 − 1 + 2(𝑛 − 1) 𝐴 𝑛 𝐴 𝑛−1 2 𝑛−1 = + 𝑛+1 𝑛 𝑛 𝑛+1 Average case QS continued 𝐴 𝑛 Assume 𝑎𝑛 = 𝑛+1 2 𝑛−1 𝑎𝑛 = 𝑎𝑛−1 + , 𝑛 𝑛+1 we get 𝑎0 = 0 Applying some simple math we have 𝑎𝑛 ≈ 2 ln 𝑛 Which give 𝐴(𝑛) ≈ 𝑛 + 1 lg 𝑛 = 𝑛 + 1 2((ln 2)(lg 𝑛) ≈ 1.38 𝑛 + 1 lg 𝑛 ∈ 𝜃(𝑛 lg 𝑛) Matrix Multiplication (Strassen) Lets look at the product of two 2 by 2’s 𝑐11 𝑐12 𝑎11 𝑎12 𝑏11 𝑏12 = × 𝑐21 𝑐22 𝑎21 𝑎22 𝑏21 𝑏22 Clearly after you do the homework #26 m1=(a11+a22)(b11+b22) m2=(a11+a22)b11 m3=a11(b12+b22) m4=a22(b21+b11) m5=(a11+a12)b22 m6=(a21+a11)(b11+b12) m7=(a12+a22)(b21+b22) will give the following! And the answer is 𝑚1 + 𝑚4 − 𝑚5 + 𝑚7 𝐶= 𝑚2 + 𝑚4 𝑚3 + 𝑚5 𝑚1 + 𝑚3 − 𝑚2 + 𝑚6 Original method 8 mult, four add/sub Strassen’s method 7 mult, and 18 add/sub Hmmmm! So what’s the big deal? Big Matrices 𝐶11 𝐶21 n 2 by n 2 𝑐12 𝐴11 𝐴12 𝐵11 𝐵12 = × 𝐵21 𝐵22 𝑐22 𝐴21 𝐴22 Where C11 is the upper left hand corner of the matrix of size n/2 by n/2. The others are similarly defined. Now m1=(A11+A22)(B11+B22) Is the sum and product of matrices Our function is then void Strassen(int n, A,B, C)// these are nxn mats { if (n<= threshold) computer C=AxB normally Partition A and B into eight submatrices strassen(n/2, A11+A22, B11+B22, M1); strassen(n/2, A21+A22, B11, M2) etc // making 7 recursive calls } NOT EIGHT! Complexity T(n) = 7T(n/2) + cn2 Which is T(n) = θ(nlg7) = O(n2.81) Using the general theorem. The best know is Coppersmith and Winograd with a time complexity of O(n2.376) Why am I using big O here? Recalling General Theorem See page 588 Assume 𝑇 𝑛 = 𝑎𝑇 𝑛 𝑏 + 𝑐𝑛𝑘 for n>1 and n a power of b, T(1)=d 𝜃 𝑛𝑘 𝑖𝑓 𝑎 < 𝑏 𝑘 𝑇 𝑛 ∈ 𝜃(𝑛𝑘 lg 𝑛) 𝑖𝑓 𝑎 = 𝑏 𝑘 𝜃 𝑛𝑙𝑜𝑔𝑏 𝑎 𝑖𝑓 𝑎 > 𝑏 𝑘 Just a side note Suppose we has 8 recursive calls instead of 7 in the above case. Then the recurrence relation would be T(n) = 8T(n/2) + cn2 This has a complexity of what? When not to use divide and conquer • An instance of size n is divided into two or more instances each almost of size n. • An instance of size n is divided into almost n instances of size n/c, where c is a constant