142 lines
		
	
	
		
			2.5 KiB
		
	
	
	
		
			C
		
	
	
		
			Executable File
		
	
	
	
	
			
		
		
	
	
			142 lines
		
	
	
		
			2.5 KiB
		
	
	
	
		
			C
		
	
	
		
			Executable File
		
	
	
	
	
| /* 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);
 | |
| }
 | |
| */
 |