代码

unsigned C[1024];
unsigned aa[100000000];
void sort(unsigned *a, int n) {
	int i;
	memset(C,0,sizeof C);
	for(i=0;i!=n;++i) {
		++C[a[i]&1023];
	}
	for(i=1;i<1024;++i) {
		C[i]+=C[i-1];
	}
	for(int i=1023;i!=0;--i) {
		C[i]=C[i-1];
	}
	C[0]=0;
	for(i=0;i!=n;++i) {
		aa[C[a[i]&1023]]=a[i];
		C[a[i]&1023]++;
	}
	memcpy(a, aa, n*sizeof(unsigned));
	
	memset(C,0,sizeof C);
	for(i=0;i!=n;++i) {
		++C[(a[i]>>10)&1023];
	}
	for(i=1;i<1024;++i) {
		C[i]+=C[i-1];
	}
	for(int i=1023;i!=0;--i) {
		C[i]=C[i-1];
	}
	C[0]=0;
	for(i=0;i!=n;++i) {
		aa[C[(a[i]>>10)&1023]]=a[i];
		C[(a[i]>>10)&1023]++;
	}
	memcpy(a, aa, n*sizeof(unsigned));
	
	memset(C,0,sizeof C);
	for(i=0;i!=n;++i) {
		++C[(a[i]>>20)&1023];
	}
	for(i=1;i<1024;++i) {
		C[i]+=C[i-1];
	}
	for(int i=1023;i!=0;--i) {
		C[i]=C[i-1];
	}
	C[0]=0;
	for(i=0;i!=n;++i) {
		aa[C[(a[i]>>20)&1023]]=a[i];
		C[(a[i]>>20)&1023]++;
	}
	memcpy(a, aa, n*sizeof(unsigned));
}

评测结果