Hello people…! In this post I will talk about another sorting algorithm, the Merge Sort, which is considered to be the fastest sorting technique by many. Merge Sort was invented by a famous computer scientist John von Neumann. It is a stable sort, so the relative position of the elements are preserved. It uses O(n) auxiliary memory to sort the elements in O(n log n) time. It is the classic example of Divide-and-Conquer Algorithms. It is very easy to understand and code. So, let’s get started..!

In the Merge Sort, we divide the array into halves recursively until we cannot divide further. Then, we merge them into a single array which would be sorted. The partitioning is similar to that in Quick Sort. We can put the merge part of the algorithm in a step-by-step process –

- Create an array which is large enough to accommodate all the elements of the two arrays to be merged.
- Compare the first elements of both arrays. Pick the least one and put it into the array we just created.
- Now we compare the first element of one array and the second element of the other array (from which the element was picked). We choose the smaller element and put it into the large array.
- We do this until we run out of elements. Then, our larger array is sorted.

For a better understanding of the merge procedure, look at the sketch below –

This is the way we merge two arrays. As you can see the prerequisite for this procedure to work is that the two smaller arrays that we are given must be sorted. How do we ensure this in our algorithm. Remember, at the lowest level, the partitions stop when the array has only one element, which is considered to be sorted. Then, we would merge two arrays of size 1 into a sorted array of size 2. This helps in merging two arrays of size 2 to form an array of size 4. In this way the recursion goes. So we will have the smaller arrays sorted initially.

As we have gained the knowledge about the working of the algorithm. Let’s put it in terms of pseudo-code. It comprises of two functions –

mergeSort(Array, low, high) if low < high mid = (low + high) / 2 for i = low to mid LeftPartition[i - low] = Array[i] for i = mid + 1 to high RightPartition[i - mid - 1] = Array[i] mergeSort(LeftPartition, low, mid) mergeSort(RightPartition, mid + 1, high) mergePartitions(Array, low, high, LeftPartition, RightPartition) mergePartitions(Array, low, high, LeftPartition, RightPartition) i = j = k = 0 mid = (low + high) / 2 while (i <= low - mid && j < high - mid) if (LeftPartition[i] < RightPartition[j]) Array[k] = LeftPartition[i] ++i else Array[k] = RightPartition[j] ++j ++k // If any array is left out if i > low - mid && j < high - mid while j < high - mid Array[k] = RightPartition[j] ++j ++k else while i <= low - mid Array[k] = RightPartition[i] ++i ++k

The pseudo-code is clear and self-explanatory. If you have any doubts, feel free to comment them. Looking at the pseudo-code try to code Merge Sort algorithm. If you get any doubts, keep reading the explanation again and again. I’m sure you’ll get it if you give your best. If you didn’t, don’t fret, I posted my code below –

/* * Merge Sort Algorithm * * Authored by, * Vamsi Sangam */ #include <stdio.h> // Merges two sorted partitions into a bigger sorted array void mergePartitions(int Array[], int low, int high, int LeftPartition[], int RightPartition[]) { int i, j, k; int mid = (low + high + 1) / 2; // Initialise carefully i = j = 0; k = low; while (i < mid - low && j <= high - mid) { if (LeftPartition[i] < RightPartition[j]) { Array[k] = LeftPartition[i]; ++i; } else { Array[k] = RightPartition[j]; ++j; } ++k; } // Copying any leftover elements if (i == mid - low && j <= high - mid) { while (j <= high - mid) { Array[k] = RightPartition[j]; ++j; ++k; } } else { while (i < mid - low) { Array[k] = LeftPartition[i]; ++i; ++k; } } } // Partitions the array into two // Does nothing if array has only one element void mergeSort(int Array[], int low, int high) { if (low < high) { int mid = (low + high + 1) / 2; int LeftPartition[mid - low]; int RightPartition[high - mid + 1]; int i; // Copying elements for (i = low; i < mid; ++i) { LeftPartition[i - low] = Array[i]; } // Copying elements for (i = mid; i <= high; ++i) { RightPartition[i - mid] = Array[i]; } // Recursive call mergeSort(LeftPartition, 0, mid - low - 1); mergeSort(RightPartition, 0, high - mid); // Merge the two partitions mergePartitions(Array, low, high, LeftPartition, RightPartition); } } int main() { int n, i; printf("Enter the size of input -\n"); scanf("%d", &n); int arr[n]; for (i = 0; i < n; ++i) { scanf("%d", &arr[i]); } mergeSort(arr, 0, n - 1); printf("\nThe Sorted Array -\n"); for (i = 0; i < n; ++i) { printf("%d ", arr[i]); } printf("\n"); return 0; }

I hope you get the intuition behind the Merge Sort Algorithm. It is the best example to explain Divide-and-Conquer technique. It has a runtime complexity of O(n log n). Which can be proved easily by using a Recursion Tree –

**Optimisation** – Merge Sort can be optimised to perform better. The idea behind this is the asymptotic complexities of various sorting algorithms. When we say that Merge Sort is better than Insertion Sort because it has a better worst case of O(n log n) which is less than O(n^{2}), this is difference can be seen only for large values of ‘n’. This is because, the terms like O(n log n), O(n) and all are valid only when n → ∞ . So, for small values of ‘n’, like 100, we may not find any significant performance differences between Merge Sort and Insertion Sort. For smaller ‘n’ an optimised Bubble Sort can actually out perform the Merge Sort, because Merge Sort comes with a lot of overheads, the partition, the merging, the recursive calls and so. Whereas, the Bubble Sort has no such over heads. So, how about we apply an optimised Bubble Sort for smaller sized partitions inside the algorithm of Merge Sort…? That is how we optimise the Merge Sort.

But how do we decide how long this “smaller sized partitions” should be..? One way is you could do an analysis. You can just test your optimised Merge Sort for various lengths of cut-off’s. You need to run your code for the same input and keep varying the ‘threshold’ value below which the smaller sized partitions are sorted by optimised Bubble Sort and not by Merge Sort partitions. You could actually plot a graph between threshold and runtime, where you will find a minimum runtime for a certain value of threshold. But remember, if you want to see significant results you need to input as many as 10^{5} or 10^{6} values.

The next best thing you can do is make a guess by looking at your resources. The idea is to look at your cache. Let’s say your cache is 2MB. Then, your cache can accommodate an integer array of size 5 ✗ 10^{5}. But, the whole of cache will not be available to you, so let’s say when you operate on an array of size 100, it is most likely that it’ll end up being in the cache. So, if we do an optimised Bubble Sort for arrays of size 100, we are most likely doing that in cache and as our process has less overheads, it actually increases the performance.

Its okay if your don’t get the whole threshold concept. But the gist of the discussion is that we will use an optimised Bubble Sort for smaller sized partitions in Merge Sort. Try to code this weird practice. You only have to make minor adjustments to your code. I’m sure you’ll get it if you give a little serious thought. I will post the code very shortly, until then keep trying !

This is the Merge Sort Algorithm. If you have any doubts or anything to say, feel free to comment them. I hope my post helped you in getting an idea of the Merge Sort. If it did, let me know by commenting…! Keep practicing… Happy Coding…! 🙂