Adaptive heap sort

In computer science, adaptive heap sort is a comparison-based sorting algorithm of the adaptive sort family. It is a variant of heap sort that performs better when the data contains existing order. Published by Christos Levcopoulos and Ola Petersson in 1992, the algorithm utilizes a new measure of presortedness, "Osc," as the number of oscillations. Instead of putting all the data into the heap as the traditional heap sort did, adaptive heap sort only take part of the data into the heap so that the run time will reduce significantly when the presortedness of the data is high.
Heapsort.
Heap sort is a sorting algorithm that utilizes binary heap data structure. The method treats an array as a complete binary tree and builds up a Max-Heap/Min-Heap to achieve sorting. It usually involves the following four steps.
Below is a C++ implementation that builds up a Max-Heap and sorts the array after the heap is built.
 #include <iostream>
 "A C++ sample heap sort code that sort an array to an increasing order"
 using namespace std;
 void heapify(int array[], int start, int end);// A "function that build up a max-heap binary tree"
 void heapify(int array[], int start, int end)
 int parent = start;
 int child = parent * 2 + 1;
 while (child <= end)
 { if (child + 1 <= end ) "//when there are two child nodes"
 if (array[child + 1] > array[child])
 child ++; "//take the bigger child node"
 
 if (array[parent] > array[child])
 return; "//if the parent node is greater, then it's already heapified"
 if (array[parent] < array[child]) "// when child node is greater than parent node"
 swap (array[parent], array[child]); "// switch parent and child node"
 parent = child;
 child = child * 2 + 1; "//continue the loop, compare the child node and its child nodes"
 void heap_sort (int array[], int len); "// heap_sort function"
 void heap_sort (int array[], int len)
 for (int i = len/2 - 1; i >= 0; i--) //Step 1: build up the max-heap
 heapify(array, i, len);
 for (int i = len - 1; i >= 0; i--) //Step 4: repeat step 2 and 3 till finished
 swap(array[0], array[i]); // Step 2: put the max at the end of the array
 heapify (array, 0, i-1); // Step 3: remove the max from the tree and heapify again
 int main()
 int array[] = {42, 1283, 123, 654, 239847, 45, 97, 85, 763, 90, 770, 616, 328, 1444, 911, 315, 38, 5040, 1}; "//the array that will be sorted"
 int array_len = sizeof(array)/sizeof(*array); "//length of the array"
 heap_sort (array, array_len);// heap sort
 return 0;
Measures of presortedness.
Measures of presortedness measures the existing order in a given sequence. These measures of presortedness decides the number of data that will be put in to the heap during the sorting process as well as the lower bound of running time.
Oscillations ("Osc").
For sequence formula_1, "Cross"(xi) is defined as the number edges of the line plot of X that are intersected by a horizontal line through the point (i, xi). Mathematically, it is defined asformula_2, for formula_3. The oscillation("Osc") of X is just the total number of intersections, defined as formula_4.
Other measures.
Besides the original Osc measurement, other known measures include the number of inversions "Inv", the number of runs "Runs", the number of blocks "Block", and the measures "Max", "Exc" and "Rem". Most of these different measurements are related for adaptive heap sort. Some measures dominate the others: every Osc-optimal algorithm is Inv optimal and Runs optimal; every Inv-optimal algorithm is Max optimal; and every Block-optimal algorithm is Exc optimal and Rem optimal.
Algorithm.
Adaptive heap sort is a variant of heap sort that seeks optimality (asymptotically optimal) with respect to the lower bound derived with the measure of presortedness by taking advantage of the existing order in the data. In heap sort, for a data formula_1 , we put all n elements into the heap and then keep extracting the maximum (or minimum) for n times. Since the time of each max-extraction action is the logarithmic in the size of the heap, the total running time of standard heap sort is O(n log n). For adaptive heap sort, instead of putting all the elements into the heap, only the possible maximums of the data (max-candidates) will be put into the heap so that fewer runs are required when each time we try to locate the maximum(or minimum). First, a Cartesian tree is built from the input in formula_6 time by putting the data into a binary tree and making each node in the tree is greater(or smaller) than all its children nodes, and the root of the Cartesian tree is inserted into an empty binary heap. Then repeatedly extract the maximum from the binary heap, retrieve the maximum in the Cartesian tree, and add its left and right children (if any) which are themselves Cartesian trees, to the binary heap. If the input is already nearly sorted, the Cartesian trees will be very unbalanced, with few nodes having left and right children, resulting in the binary heap remaining small, and allowing the algorithm to sort more quickly than formula_7 for inputs that are already nearly sorted.
 Input: an array of n elements that need to be sorted
 Construct the Cartesian tree "l"(x)
 Insert the root of "l"(x) into a heap
 for i = from 1 to n
 Perform ExtractMax on the heap 
 if the max element extracted has any children in "l"(x)
 retrieve the children in "l"(x)
 insert the children element into the heap
Drawbacks.
Despite decades of research, there's still a gap between the theory of adaptive heap sort and its practical use. Because the algorithm makes use of Cartesian trees and pointer manipulation, it has low cache-efficiency and high memory requirements, both of which deteriorate the performance of implementations.