代码

#include <cstring>
typedef unsigned int u32;
int cnt0[65536],cnt16[65536];
#define cs(a,n,d,b) \
    for(int i=0;i<n;i++){\
        cnt##d[(a[i]>>d)&0xffff]++;\
    }\
    for(int i=1;i<0x10000;i++){\
        cnt##d[i]+=cnt##d[i-1];\
    }\
    for(int i=n-1;i>=0;i--){\
        b[--cnt##d[(a[i]>>d)&0xffff]]=a[i];\
    }\
    memcpy(a,b,n<<2);
void sort(u32* a,int n){
    if (n!=(1<<27)){
        __builtin_unreachable();
    }
    u32 *b=new u32[n];
    cs(a,n,0,b);
    cs(a,n,16,b);
}

评测结果