代码

void sort(unsigned *a,int n)
{
	unsigned *b=new unsigned[n];
	unsigned cnt0[256],cnt8[256],cnt16[256],cnt24[256];
	for(int i=0;i<n;++i){
		++cnt0[a[i]&255];
		++cnt8[a[i]>>8&255];
		++cnt16[a[i]>>16&255];
		++cnt24[a[i]>>24&255];
	}
	for(int i=1;i<=255;++i){
		cnt0[i]+=cnt0[i-1];
		cnt8[i]+=cnt8[i-1];
		cnt16[i]+=cnt16[i-1];
		cnt24[i]+=cnt24[i-1];
	}
#define countingSort(a,b,k) {\
	int tim=n>>2;\
	unsigned *now=a+n-1-4;\
    while(tim--){\
        b[--cnt##k[now[4]>>k&255]]=now[4];\
        b[--cnt##k[now[3]>>k&255]]=now[3];\
        b[--cnt##k[now[2]>>k&255]]=now[2];\
        b[--cnt##k[now[1]>>k&255]]=now[1];\
		now-=4;\
    }\
    switch(n&3){\
        case 3:b[--cnt##k[a[2]>>k&255]]=a[2];\
        case 2:b[--cnt##k[a[1]>>k&255]]=a[1];\
        case 1:b[--cnt##k[a[0]>>k&255]]=a[0];\
    }\
}
	countingSort(a,b,0);
	countingSort(b,a,8);
	countingSort(a,b,16);
	countingSort(b,a,24);
}

评测结果