Final Module2
Final Module2
Final Module2
Contents
1. General method
2. Recurrence equation
3. Algorithm: Binary search
4. Algorithm: Finding the maximum and minimum
5. Algorithm: Merge sort
6. Algorithm: Quick sort
7. Algorithm: Strassen’s matrix multiplication
8. Advantages and Disadvantages
9. Decrease and Conquer Approach
10. Algorithm: Topological Sort
Module 2: Divide and Conquer
1. General method
Divide and Conquer is one of the best-known general algorithm design
techniques. It works according to the following general plan:
• Given a function to compute on ‘n’ inputs the divide-and-conquer strategy suggests
splitting the inputs into ‘k’ distinct subsets, 1<k<=n, yielding ‘k’ sub problems.
• These sub problems must be solved, and then a method must be found to combine sub
solutions into a solution of the whole.
• If the sub problems are still relatively large, then the divide-and-conquer strategy can
possibly be reapplied.
• Often the sub problems resulting from a divide-and-conquer design are of the same
type as the original problem. For those cases the re application of the divide-and-
conquer principle is naturally expressed by a recursive algorithm.
Problem of size n
Where,
• T(n) is the time for divide and conquer method on any input of size n and
• g(n) is the time to compute answer directly for small inputs.
• The function f(n) is the time for dividing the problem ‘p’ and combining the solutions
to sub problems.
For divide and conquer based algorithms that produce sub problems of the same type as the
original problem, it is very natural to first describe them by using recursion.
More generally, an instance of size n can be divided into b instances of size n/b, with a of
them needing to be solved. (Here, a and b are constants; a>=1 and b > 1.). Assuming that
size n is a power of b(i.e. n=bk), to simplify our analysis, we get the following recurrence for
the running time T(n):
..... (1)
where f(n) is a function that accounts for the time spent on dividing the problem into smaller
ones and on combining their solutions.
The recurrence relation can be solved by i) substitution method or by using ii) master
theorem.
1. Substitution Method - This method repeatedly makes substitution for each
occurrence of the function T in the right-hand side until all such occurrences
disappear.
2. Master Theorem - The efficiency analysis of many divide-and-conquer algorithms is
greatly simplified by the master theorem. It states that, in recurrence equation T(n) =
aT(n/b) + f (n), If f (n)∈ Θ (nd) where d ≥ 0 then
Problems on Substitution method & Master theorem to solve the recurrence relation
3. Binary Search
Problem definition: Let ai, 1 ≤ i ≤ n be a list of elements that are sorted in non-decreasing
order. The problem is to find whether a given element x is present in the list or not. If x is
present we have to determine a value j (element’s position) such that aj=x. If x is not in the
list, then j is set to zero.
Solution: Let P = (n, ai…al , x) denote an arbitrary instance of search problem where n is the
number of elements in the list, ai…alis the list of elements and x is the key element to be
searched for in the given list. Binary search on the list is done as follows:
Step1: Pick an index q in the middle range [i, l] i.e. q=⎝(n + 1)/2] and compare x with aq.
Step 2: if x = aqi.e key element is equal to mid element, the problem is immediately solved.
Step 3: if x <aqin this case x has to be searched for only in the sub-list ai, ai+1, ……, aq-
Therefore, problem reduces to (q-i, ai…aq-1, x).
Step 4: if x >aq,x has to be searched for only in the sub-list aq+1, ...,., al . Therefore problem
reduces to (l-i, aq+1…al, x).
For the above solution procedure, the Algorithm can be implemented as recursive or non-
recursive algorithm.
Analysis
In binary search the basic operation is key comparison. Binary Search can be analyzed with
the best, worst, and average case number of comparisons. The numbers of comparisons for
the recursive and iterative versions of Binary Search are the same, if comparison counting is
relaxed slightly. For Recursive Binary Search, count each pass through the if-then-else block
as one comparison. For Iterative Binary Search, count each pass through the while block as
one comparison. Let us find out how many such key comparison does the algorithm make on
an array of n elements.
Best case –Θ(1) In the best case, the key is the middle in the array. A constant number of
comparisons (actually just 1) are required.
Worst case - Θ(log2 n) In the worst case, the key does not exist in the array at all. Through
each recursion or iteration of Binary Search, the size of the admissible range is halved. This
halving can be done ceiling (log2n ) times. Thus, ] log2 n ] comparisons are required.
Sometimes, in case of the successful search, it may take maximum number of comparisons.
] log2 n ]. So worst case complexity of successful binary search is Θ (log2 n).
Average case - Θ (log2n) To find the average case, take the sum of the product of number of
comparisons required to find each element and the probability of searching for that element.
To simplify the analysis, assume that no item which is not in array will be searched for, and
that the probabilities of searching for each element are uniform.
Space Complexity - The space requirements for the recursive and iterative versions of binary
search are different. Iterative Binary Search requires only a constant amount of space, while
Recursive Binary Search requires space proportional to the number of comparisons to
maintain the recursion stack.
Advantages: Efficient on very big list, Can be implemented iteratively/recursively.
Limitations:
• Interacts poorly with the memory hierarchy
• Requires sorted list as an input
• Due to random access of list element, needs arrays instead of linked list.
Explanation:
▪ StraightMaxMin requires 2(n-1) comparisons in the best, average & worst cases.
▪ By realizing the comparison of a[i]>max is false, improvement in a algorithm can be
done. Hence we can replace the contents of the for loop by,
If(a[i]>Max) then Max = a[i]; Else if (a[i]<min) min=a[i]
▪ On the average a[i] is > max half the time.So, the avg. no. of comparison is 3n/2-1.
Example:
Compared with the straight forward method (2n-2) this method saves 25% in comparisons.
Space Complexity
Compared to the straight forward method, the MaxMin method requires extra stack space for
i, j, max, min, max1 and min1. Given n elements there will be ⎝log2n] + 1 levels ofrecursion
and we need to save seven values for each recursive call. (6 + 1 for return address).
5. Merge Sort
Merge sort is a perfect example of a successful application of the divide-and conquer
technique. It sorts a given array A [O ... n - 1] by dividing it into two halves A [0 .. ⎝n/2]-1]
and A [ ⎝n/2] .. n-1], sorting each of them recursively, and then merging the two smaller
sorted arrays into a single sorted one.
▪ After that, the index of the smaller element is incremented to point to its immediate
successor in the array it was copied from. This operation is repeated until one of the
two given arrays is exhausted, and then the remaining elements of the other array are
copied to the end of the new array.
Example:
The operation of the algorithm on the
list 8, 3, 2, 9, 7, 1, 5, 4 is illustrated in
the figure
Analysis
Here the basic operation is key comparison. As merge sort execution does not depend on
the order of the data, best case and average case runtime are the same as worst case
runtime.
Worst case: During key comparison, neither of the two arrays becomes empty before the
other one contains just one element leads to the worst case of merge sort. Assuming for
simplicity that total number of elements n is a power of 2, the recurrence relation for the
number of key comparisons C(n) is
where, Cmerge(n) is the number of key comparisons made during the merging stage.
Let us analyze Cmerge(n), the number of key comparisons performed during the merging stage.
At each step, exactly one comparison is made, after which the total number of elements in the
two arrays still needing to be processed is reduced by 1. In the worst case, neither of the two
arrays becomes empty before the other one contains just one element (e.g., smaller elements
may come from the alternating arrays). Therefore, for the worst case, Cmerge(n) = n –1.
Now,
Advantages:
• Number of comparisons performed is nearly optimal.
• For large n, the num er of comparisons made by this algorithm in the average case
turns out to be about 0.25n less and hence is also in Θ(n log n).
• Merge sort will never degrade to O (n2)
• Another advantage of merge sort over quicksort and heapsort is its stability. (A
sorting algorithm is said to be stable if two objects with equal keys appear in the same
order in sorted output as they appear in the input array to be sorted. )
Limitations:
• The principal shortcoming of merge sort is the linear amount [O(n) ] of extra storage
the algorithm requires. Though merging can be done in-place, the resulting algorithm
is quite complicated and of theoretical interest only.
6. Quick sort
Quicksort is the other important sorting algorithm that is based on the divide-and-conquer
approach. Unlike merge sort, which divides its input elements according to their position in
the array, quicksort divides ( or partitions) them according to their value.
A partition is an arrangement of the array’s elements so that all the elements to the left of
some element A[s] are less than or equal to A[s], and all the elements to the right of A[s] are
greater than or equal to it:
Obviously, after a partition is achieved, A[s] will be in its final position in the sorted array,
and we can continue sorting the two subarrays to the left and to the right of A[s]
independently (e.g., by the same method).
In quick sort, the entire work happens in the division stage, with no work required to combine
the solutions to the subproblems.
Partitioning
We start by selecting a pivot—an element with respect to whose value we are going to divide
the subarray. There are several different strategies for selecting a pivot. We use
the sophisticated method suggested by C.A.R. Hoare, the prominent British computer
scientist who invented quicksort.
Select the subarray’s first element: p = A[l].Now scan the subarray from both ends,
comparing the subarray’s elements to the pivot.
▪ The left-to-right scan, denoted below by index pointer i, starts with the second
element. Since we want elements smaller than the pivot to be in the left part of the
subarray, this scan skips over elements that are smaller than the pivot and stops upon
encountering the first element greater than or equal to the pivot.
▪ The right-to-left scan, denoted below by index pointer j, starts with the last element of
the subarray. Since we want elements larger than the pivot to be in the right part of the
subarray, this scan skips over elements that are larger than the pivot and stops on
encountering the first element smaller than or equal to the pivot.
After both scans stop, three situations may arise, depending on whether or not the scanning
indices have crossed.
1. If scanning indices i and j have not crossed, i.e., i< j, we simply exchange A[i] and
A[j ] and resume the scans by incrementing I and decrementing j, respectively:
2. If the scanning indices have crossed over, i.e., i> j, we will have partitioned the
subarray after exchanging the pivot with A[j]:
3. If the scanning indices stop while pointing to the same element, i.e., i = j, the value
they are pointing to must be equal to p. Thus, we have the subarray partitioned, with
the split position s = i = j :
We can combine this with the case-2 by exchanging the pivot with A[j] whenever i≥j
ALGORITHM HoarePartition(A[l..r])
//Partitions a subarray by Hoare’s algorithm, using the first element as a pivot
//Input: Subarray of array A[0..n − 1], defined by its left and right indices l and r (l<r)
//Output: Partition of A[l..r], with the split position returned as this function’s value
Note that index i can go out of the subarray’s bounds in this pseudocode.
Analysis
Best Case -Here the basic operation is key comparison. Number of key comparisons made
before a partition is achieved is n + 1 if the scanning indices cross over and n if they coincide.
If all the splits happen in the middle of corresponding subarrays, we will have the best case.
The number of key comparisons in the best case satisfies the recurrence,
According to the Master Theorem, Cbest(n) ∈Θ(n log2 n); solving it exactly form = 2k yields
Cbest(n) = n log2 n.
Worst Case – In the worst case, all the splits will be skewed to the extreme: one of the two
subarrays will be empty, and the size of the other will be just 1 less than the size of the
subarray being partitioned. This unfortunate situation will happen, in particular, for increasing
arrays. Indeed, if A[0..n − 1] is a strictly increasing array and we use A[0] as the pivot, the
left-to-right scan will stop on A[1] while the right-to-left scan will go all the way to reach
A[0], indicating the split at position 0:So, after making n + 1 comparisons to get to this
partition and exchanging the pivot A[0] with itself, the algorithm will be left with the strictly
increasing array A[1..n − 1] to sort. This sorting of strictly increasing arrays of diminishing
sizes will continue until the last one A[n−2.. n−1] has been processed. The total number of
key comparisons made will be equal to
Average Case - Let Cavg(n) be the average number of key comparisons made by quicksort on
a randomly ordered array of size n. A partition can happen in any position s (0 ≤ s ≤n−1) after
n+1comparisons are made to achieve the partition. After the partition, the left and right
subarrays will have s and n − 1−s elements, respectively. Assuming that the partition split can
happen in each position s with the same probability 1/n, we get the following recurrence
relation:
Its solution, which is much trickier than the worst- and best-case analyses, turns out to be
Thus, on the average, quicksort makes only 39% more comparisons than in the best case.
Moreover, its innermost loop is so efficient that it usually runs faster than mergesort on
randomly ordered arrays of nontrivial sizes. This certainly justifies the name given to the
algorithm by its inventor.
Variations: Because of quicksort’s importance, there have been persistent efforts over the
years to refine the basic algorithm. Among several improvements discovered by researchers
are:
▪ Better pivot selection methods such as randomized quicksort that uses a random
element or the median-of-three method that uses the median of the leftmost,
rightmost, and the middle element of the array
▪ Switching to insertion sort on very small subarrays (between 5 and 15 elements for
most computer systems) or not sorting small subarrays at all and finishing the
algorithm with insertion sort applied to the entire nearly sorted array
▪ Modifications of the partitioning algorithm such as the three-way partition into
segments smaller than, equal to, and larger than the pivot
Limitations: 1. It is not stable. 2. It requires a stack to store parameters of subarrays that are
yet to be sorted. 3. While Performance on randomly ordered arrays is known to be sensitive
not only to the implementation details of the algorithm but also to both computer architecture
and data type.
Divide and Conquer method: Using the divide and conquer matrix multiplication is
achieved through multiplying the submatrices recursively.
In the above method, we do 8 multiplications for matrices of size n/2 x n/2 and 4 additions.
Addition of two matrices takes O(n2) time. So the time complexity can be written asT(n) =
8T(n/2) + O(n2) which happens
n to be O(n3); same as the direct method
Divide and Conquer through Strassen’s Method: By using divide-and-conquer approach
proposed by Strassen in 1969, we can reduce the number of multiplications.
Multiplication of 2×2 matrices :The principal in sight of the algorithm lies in the discovery
that we can find the product C of two 2 × 2 matrices A and B with just 7 multiplications as
opposed to the eight required by the brute-force algorithm. This is accomplished by using the
following formulas:
where
Thus, to multiply two 2×2 matrices, Strassen’s algorithm makes seven multiplications and 18
additions/subtractions, whereas the brute-force algorithm requires eight multiplications and
four additions.
It is not difficult to verify that one can treat these submatrices as numbers toget the correct
product. For example, C00 can be computed as M1 + M4– M5 + M7 where M1, M4, M5, and M7
are found by Strassen’s formulas, with the numbers replaced by sub the corresponding
matrices. If the seven products of n/2 × n/2 matrices are computed same recursively by the
method, we have Strassen’s algorithm for matrix multiplication.
Analysis( fromtext book T1: Levtin et al )
Here the basic operation is multiplication. If M(n) is the number of multiplications made
by Strassen’s algorithm in multiplying two n × n matrices (where n is a power of 2), we get
the following recurrence relation for it:
This implies M(n) = Θ(n2.807)which is smaller than n3 required by the brute-force algorithm.
Example:
Note that a back edge in a DFS forest of a directed graph can connect a vertex to its parent.
Whether or not it is the case, the presence of a back edge indicates that the digraph has a
directed cycle. A directed cycle in a digraph is a sequence of three or more of its vertices that
starts and ends with the same vertex and in which every vertex is connected to its immediate
predecessor by an edge directed from the predecessor to the successor. For example, a, b, a
is a directed cycle in the digraph in Figure given above. Conversely, if a DFS forest of
a digraph has no back edges, the digraph is a dag, an acronym for directed acyclic graph.
courses? The situation can be modeled by a digraph in which vertices represent courses and
directed edges indicate prerequisite requirements.
In terms of this digraph, the uestion is whether we can list its vertices in such an order that
for every edge in the graph, the vertex where the edge starts is listed before the vertex where
the edge ends. In other words, can you find such an ordering of this digraph’s vertices? This
problem is called topological sorting.
Topological Sort: For topological sorting to be possible, a digraph in question must be a
DAG. i.e., if a digraph has no directed cycles, the topological sorting problem for it has a
solution.
There are two efficient algorithms that both verify whether a digraph is a dag and, if it is,
produce an ordering of vertices that solves the topological sorting problem. The first one is
based on depth-first search; the second is based on a direct application of the decrease-by-one
technique.
Topological Sorting based on DFS
Method
1. Perform a DFS traversal and note the order in which vertices become dead-ends
2. Reversing this order yields a solution to the topological sorting problem, provided, of
course, no back edge has been encountered during the traversal. If a back edge has
been encountered, the digraph is not a DAG, and topological sorting of its vertices is
impossible.
Illustration
a) Digraph for which the topological sorting problem needs to be solved.
b) DFS traversal stack w th the subscript numbers indicating the popping off order.
c) Solution to the problem. Here we have drawn the edges of the digraph, andthey all
point from left to right as the problem’s statement requires. It is a convenient way to
check visually the correctness of a solution to an instance of the topological sorting
problem.
Note: The solution obtained by the source-removal algorithm is different from the one
obtained by the DFS-based algorithm. Both of them are correct, of course; the topological
sorting problem may have several alternative solutions.
Applications of Topological Sorting
• Instruction scheduling in program compilation
• Cell evaluation ordering in spreadsheet formulas,
• Resolving symbol dependencies in linkers.
***