Cod sursa(job #551789)
#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 i;
int cnt[256],ind[256];
memset(cnt,0,sizeof(cnt));
memset(ind,0,sizeof(ind));
for(i=0;i<N;i++) cnt[(a[i]&(255<<byte))>>byte]++;
for(i=1;i<256;i++) ind[i]=ind[i-1]+cnt[i-1];
for(i=0;i<N;i++) b[ind[(a[i]&(255<<byte))>>byte]++]=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(8,x,v);
radix(16,v,x);
radix(24,x,v);
for(i=0;i<N-1;i++) out<<v[i]<<' ';
out<<v[N-1]<<'\n';
return 0;
}