gistfile1.cpp

void rsort(uint64_t *a, uint64_t n, uint64_t *b)
{
  uint32_t cnt[0x4001]={};

  for (size_t i=0;i<n;i++) cnt[(a[i]&0x1fff)+1]++;
  for (size_t i=0;i<0x2000;i++) cnt[i+1]+=cnt[i];
  for (size_t i=0;i<n;i++) b[cnt[a[i]&0x1fff]++]=a[i];
  std::swap(a, b);

  for (size_t i=0;i<0x2000;i++) cnt[i]=0;
  for (size_t i=0;i<n;i++) cnt[((a[i]>>13)&0x1fff)+1]++;
  for (size_t i=0;i<0x2000;i++) cnt[i+1]+=cnt[i];
  for (size_t i=0;i<n;i++) b[cnt[(a[i]>>13)&0x1fff]++]=a[i];
  std::swap(a, b);

  for (size_t i=0;i<0x2000;i++) cnt[i]=0;
  for (size_t i=0;i<n;i++) cnt[((a[i]>>26))+1]++;
  for (size_t i=0;i<0x4000;i++) cnt[i+1]+=cnt[i];
  for (size_t i=0;i<n;i++) b[cnt[(a[i]>>26)]++]=a[i];
}