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.
142 lines
2.4 KiB
142 lines
2.4 KiB
5 years ago
|
/* 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);
|
||
|
}
|
||
|
*/
|