Pagini recente » Cod sursa (job #72632) | Cod sursa (job #1768247) | Cod sursa (job #1302940) | Cod sursa (job #1731540) | Cod sursa (job #455557)
Cod sursa(job #455557)
/*
* File: main.cpp
* Author: virtualdemon
*
* Created on May 13, 2010, 9:32 PM
*/
#include <cstdlib>
#include <fstream>
#include <iterator>
#define Nmax 500011
#define t16 1<<16
#define t8 1<<8
#define t4 1<<4
#define t2 1<<2
/*
*
*/
using namespace std;
int H[Nmax];
inline void DownHeap( int k, int N )
{
for( int son=2*k; son <= N; k=son, son*=2 )
{
if( son < N && H[son] < H[son+1] )
++son;
if( H[k] >= H[son] )
return;
swap( H[k], H[son] );
}
}
inline void HeapSort( int left, int right )
{
for( int i=left+(right-left)/2; i >= left; --i )
DownHeap( i, right );
while( right >= left )
{
swap( H[left], H[right] );
--right;
DownHeap( left, right );
}
}inline int medianof3( int left, int middle, int right )
{
if( H[left] > H[middle] )
swap( H[left], H[middle] );
if( H[left] > H[right] )
swap( H[left], H[right] );
if( H[middle] > H[right] )
swap( H[middle], H[right] );
swap( H[middle], H[right-1] );
return H[right-1];
}
inline void IntroSort( int left, int right, int dlimit )
{
if( left < right )
{
if( !dlimit )
{
HeapSort( left, right );
return;
}
int lleft=left-1, rright=right+1, pValue=medianof3( left, left+(right-left)/2, right );
for( ; ; )
{
while( lleft <= right && pValue > H[++lleft] );
while( rright >= 0 && pValue < H[--rright] );
if( lleft > rright )
break;
swap( H[lleft], H[rright] );
}
if( left < rright )
IntroSort( left, rright, dlimit-1 );
if( lleft < right )
IntroSort( lleft, right, dlimit-1 );
}
}
inline int Log2( int N )
{
int i, j;
for( i=1, j=0; i < N; i<<=1, ++j );
return j-1;
}
int main(int argc, char** argv)
{
int N;
ifstream in( "algsort.in" );
in>>N;
copy( istream_iterator<int>(in), istream_iterator<int>(), H+1 );
IntroSort( 1, N, 2*Log2(N) );
ofstream out( "algsort.out" );
copy( H+1, H+N+1, ostream_iterator<int>(out, " " ) );
out<<'\n';
return (EXIT_SUCCESS);
}