|
| 1 | +module mergesort; |
| 2 | + |
| 3 | +import std.stdio; |
| 4 | + |
| 5 | +pure bool isSorted(T)(T[] data) { |
| 6 | + for (int i = 0; i < data.length - 1; ++i) { |
| 7 | + if (data[i] > data[i + 1]) { |
| 8 | + return false; |
| 9 | + } |
| 10 | + } |
| 11 | + |
| 12 | + return true; |
| 13 | +} |
| 14 | + |
| 15 | +unittest { |
| 16 | + int[] data = [5, 4, 2, 6, 2, 9, 5, 2, 0, 2, 9, 7]; |
| 17 | + mergeSort(data); |
| 18 | + assert(isSorted(data)); |
| 19 | +} |
| 20 | + |
| 21 | +public pure void mergeSort(T)(T[] data) { |
| 22 | + T[] aux = new T[data.length]; |
| 23 | + mergeSort(data, aux); |
| 24 | +} |
| 25 | + |
| 26 | +private pure void mergeSort(T)(T[] data, T[] aux) { |
| 27 | + if (data.length == 1) { |
| 28 | + return; |
| 29 | + } |
| 30 | + |
| 31 | + mergeSort(data[0..$/2], aux); |
| 32 | + mergeSort(data[$/2..$], aux); |
| 33 | + merge(data[0..$/2], data[$/2..$], aux); |
| 34 | +} |
| 35 | + |
| 36 | +private pure void merge(T)(T[] left, T[] right, T[] aux) |
| 37 | +in { |
| 38 | + debug { writeln(left, right); } |
| 39 | + assert(isSorted(left)); |
| 40 | + assert(isSorted(right)); |
| 41 | +} out { |
| 42 | + debug { writeln(aux[0..(left.length + right.length)]); } |
| 43 | + assert(isSorted(aux[0..(left.length + right.length)])); |
| 44 | +} body { |
| 45 | + if (left[$ - 1] < right[0]) { |
| 46 | + return; |
| 47 | + } |
| 48 | + |
| 49 | + immutable ulong size = left.length + right.length; |
| 50 | + |
| 51 | + int i = 0; |
| 52 | + int l = 0; |
| 53 | + int r = 0; |
| 54 | + |
| 55 | + while (l < left.length && r < right.length) { |
| 56 | + if (left[l] < right[r]) { |
| 57 | + aux[i++] = left[l++]; |
| 58 | + } else { |
| 59 | + aux[i++] = right[r++]; |
| 60 | + } |
| 61 | + } |
| 62 | + |
| 63 | + if (l == left.length) { |
| 64 | + aux[i..size] = right[r..$]; |
| 65 | + } else { |
| 66 | + aux[i..size] = left[l..$]; |
| 67 | + } |
| 68 | + |
| 69 | + left[] = aux[0..left.length]; |
| 70 | + right[] = aux[left.length..size]; |
| 71 | +} |
0 commit comments