Cod sursa(job #551541)
#include <fstream>
#include <string.h>
using namespace std;
int i,N;
int v[500001],x[500001];
inline void radix(int byte,int a[],int b[])
{
int count[256];
int index[256];
int i;
memset(count,0,sizeof(count));
memset(index,0,sizeof(index));
for(i=0;i<N;i++) count[(a[i]&(255<<(byte*8)))>>(byte*8)]++;
index[0]=0;
for(i=1;i<256;i++) index[i]=index[i-1]+count[i-1];
for(i=0;i<N;i++) b[index[(a[i]&(255<<(byte*8)))>>(byte*8)]++]=a[i];
}
int main()
{
ifstream in("algsort.in");
ofstream out("algsort.out");
in>>N;
for(i=0;i<N;i++) in>>v[i];
radix(0,v,x);
radix(1,x,v);
radix(2,v,x);
radix(3,x,v);
for(i=0;i<N-1;i++) out<<v[i]<<' ';
out<<v[N-1]<<'\n';
return 0;
}