You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
141 lines
2.4 KiB
141 lines
2.4 KiB
/* Sort.c -- Sort functions |
|
2014-04-05 : Igor Pavlov : Public domain */ |
|
|
|
#include "Precomp.h" |
|
|
|
#include "Sort.h" |
|
|
|
#define HeapSortDown(p, k, size, temp) \ |
|
{ for (;;) { \ |
|
size_t s = (k << 1); \ |
|
if (s > size) break; \ |
|
if (s < size && p[s + 1] > p[s]) s++; \ |
|
if (temp >= p[s]) break; \ |
|
p[k] = p[s]; k = s; \ |
|
} p[k] = temp; } |
|
|
|
void HeapSort(UInt32 *p, size_t size) |
|
{ |
|
if (size <= 1) |
|
return; |
|
p--; |
|
{ |
|
size_t i = size / 2; |
|
do |
|
{ |
|
UInt32 temp = p[i]; |
|
size_t k = i; |
|
HeapSortDown(p, k, size, temp) |
|
} |
|
while (--i != 0); |
|
} |
|
/* |
|
do |
|
{ |
|
size_t k = 1; |
|
UInt32 temp = p[size]; |
|
p[size--] = p[1]; |
|
HeapSortDown(p, k, size, temp) |
|
} |
|
while (size > 1); |
|
*/ |
|
while (size > 3) |
|
{ |
|
UInt32 temp = p[size]; |
|
size_t k = (p[3] > p[2]) ? 3 : 2; |
|
p[size--] = p[1]; |
|
p[1] = p[k]; |
|
HeapSortDown(p, k, size, temp) |
|
} |
|
{ |
|
UInt32 temp = p[size]; |
|
p[size] = p[1]; |
|
if (size > 2 && p[2] < temp) |
|
{ |
|
p[1] = p[2]; |
|
p[2] = temp; |
|
} |
|
else |
|
p[1] = temp; |
|
} |
|
} |
|
|
|
void HeapSort64(UInt64 *p, size_t size) |
|
{ |
|
if (size <= 1) |
|
return; |
|
p--; |
|
{ |
|
size_t i = size / 2; |
|
do |
|
{ |
|
UInt64 temp = p[i]; |
|
size_t k = i; |
|
HeapSortDown(p, k, size, temp) |
|
} |
|
while (--i != 0); |
|
} |
|
/* |
|
do |
|
{ |
|
size_t k = 1; |
|
UInt64 temp = p[size]; |
|
p[size--] = p[1]; |
|
HeapSortDown(p, k, size, temp) |
|
} |
|
while (size > 1); |
|
*/ |
|
while (size > 3) |
|
{ |
|
UInt64 temp = p[size]; |
|
size_t k = (p[3] > p[2]) ? 3 : 2; |
|
p[size--] = p[1]; |
|
p[1] = p[k]; |
|
HeapSortDown(p, k, size, temp) |
|
} |
|
{ |
|
UInt64 temp = p[size]; |
|
p[size] = p[1]; |
|
if (size > 2 && p[2] < temp) |
|
{ |
|
p[1] = p[2]; |
|
p[2] = temp; |
|
} |
|
else |
|
p[1] = temp; |
|
} |
|
} |
|
|
|
/* |
|
#define HeapSortRefDown(p, vals, n, size, temp) \ |
|
{ size_t k = n; UInt32 val = vals[temp]; for (;;) { \ |
|
size_t s = (k << 1); \ |
|
if (s > size) break; \ |
|
if (s < size && vals[p[s + 1]] > vals[p[s]]) s++; \ |
|
if (val >= vals[p[s]]) break; \ |
|
p[k] = p[s]; k = s; \ |
|
} p[k] = temp; } |
|
|
|
void HeapSortRef(UInt32 *p, UInt32 *vals, size_t size) |
|
{ |
|
if (size <= 1) |
|
return; |
|
p--; |
|
{ |
|
size_t i = size / 2; |
|
do |
|
{ |
|
UInt32 temp = p[i]; |
|
HeapSortRefDown(p, vals, i, size, temp); |
|
} |
|
while (--i != 0); |
|
} |
|
do |
|
{ |
|
UInt32 temp = p[size]; |
|
p[size--] = p[1]; |
|
HeapSortRefDown(p, vals, 1, size, temp); |
|
} |
|
while (size > 1); |
|
} |
|
*/
|
|
|