基数排序——尚未补完的坑QAQ

基数排序复杂度是(n+b)logn/logb

我们找一个基数 每次处理一部分位 从低位到高位处理

t是出现次数 s是这个桶管辖的起点 然后就可以写了

不过我这里是指针版的 有点难看

#include<cstdio>
#include<cstring>
#include<algorithm>
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
const int N=2e5+,B=;
int *s[],*mp;
int n,t[],b[N],a[N];
int main(){
    n=read(); for(int i=;i<n;i++) a[i]=read();
    memset(t,,sizeof(t)); mp=b;
    for(int i=;i<n;i++) ++t[a[i]&];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[a[i]&]++=a[i];
    memset(t,,sizeof(t)); mp=a;
    for(int i=;i<n;i++) ++t[b[i]>>];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[b[i]>>]++=b[i];
    for(int i=;i<n;i++) printf("%d ",a[i]);
    return ;
}

其实排结构体也差不多呢

#include<cstdio>
#include<cstring>
#include<algorithm>
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
const int N=2e5+,B=;
struct pos{int v;};
pos *s[],*mp;
pos b[N],a[N];
int n,t[];
int main(){
    n=read(); for(int i=;i<n;i++) a[i].v=read();
    memset(t,,sizeof(t)); mp=b;
    for(int i=;i<n;i++) ++t[a[i].v&];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[a[i].v&]++=a[i];
    memset(t,,sizeof(t)); mp=a;
    for(int i=;i<n;i++) ++t[b[i].v>>];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[b[i].v>>]++=b[i];
    for(int i=;i<n;i++) printf("%d ",a[i].v);
    return ;
}

基数排序其实也可以做多关键字排序 不过要按关键字顺序倒着来

就是先处理第n关键字倒着处理到第1关键字

代码待填

#include<cstdio>
#include<cstring>
#include<algorithm>
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
const int N=2e5+,B=;
struct pos{int l,r;};
pos *s[],*mp;
pos b[N],a[N];
int n,t[];
int main(){
    n=read(); for(int i=;i<n;i++) a[i].l=read(),a[i].r=read();
    memset(t,,sizeof(t)); mp=b;
    for(int i=;i<n;i++) ++t[a[i].r&];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[a[i].r&]++=a[i];
    memset(t,,sizeof(t)); mp=a;
    for(int i=;i<n;i++) ++t[b[i].r>>];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[b[i].r>>]++=b[i];
    memset(t,,sizeof(t)); mp=b;
    for(int i=;i<n;i++) ++t[a[i].l&];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[a[i].l&]++=a[i];
    memset(t,,sizeof(t)); mp=a;
    for(int i=;i<n;i++) ++t[b[i].l>>];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[b[i].l>>]++=b[i];
    for(int i=;i<n;i++) printf("[%d %d]\n",a[i].l,a[i].r);
    return ;
}

#include<cstdio>
#include<cstring>
#include<algorithm>
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
const int N=2e5+,B=;
struct pos{int l,r;};
pos *s[],*mp;
pos b[N],a[N];
int n,t[];
int main(){
    n=read();
    for(int i=;i<n;i++) a[i].l=read(),a[i].r=read();
    memset(t,,sizeof(t)); mp=b;
    for(int i=;i<n;i++) ++t[a[i].r&];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[a[i].r&]++=a[i];
    memset(t,,sizeof(t)); mp=a;
    for(int i=;i<n;i++) ++t[b[i].r>>];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[b[i].r>>]++=b[i];
    memset(t,,sizeof(t)); mp=b;
    for(int i=;i<n;i++) ++t[a[i].l&];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[a[i].l&]++=a[i];
    memset(t,,sizeof(t)); mp=a;
    for(int i=;i<n;i++) ++t[b[i].l>>];
    for(int i=;i<B;i++) s[i]=mp,mp+=t[i];
    for(int i=;i<n;i++) *s[b[i].l>>]++=b[i];
    for(int i=;i<n;i++) printf("[%d %d]\n",a[i].l,a[i].r);
    return ;
}

存一下ccz的代码

大小为1e6范围内的

#include<bits/stdc++.h>
const int N=1e8,B=;
int a[N],*ls[],*rs[],t[],b[N],*mp;
int main(){
    for(int i=;i<N;++i)a[i]=(rand()^rand()<<)%;//generate data
    int tt=clock();
    //radix sort int in range[0,1e6)
    memset(t,,sizeof(t));
    mp=b;
    for(int i=;i<N;++i)++t[a[i]&];
    for(int i=;i<B;++i)ls[i]=rs[i]=mp,mp+=t[i];
    for(int i=;i<N;++i)*rs[a[i]&]++=a[i];
    memset(t,,sizeof(t));
    mp=a;
    for(int i=;i<N;++i)++t[b[i]>>];
    for(int i=;i<B;++i)ls[i]=rs[i]=mp,mp+=t[i];
    for(int i=;i<N;++i)*rs[b[i]>>]++=b[i];
    //the result is in a
    printf("%g s",(clock()-tt)*./CLOCKS_PER_SEC);
    return ;
}

大小为1e9范围内的

#include<bits/stdc++.h>
const int N=1e8,B=;
int a[N],*ls[],*rs[],t[],b[N],*mp;
int main(){
    for(int i=;i<N;++i)a[i]=(rand()^rand()<<)%;//generate data
    int tt=clock();
    //radix sort int in range[0,1e9)
    memset(t,,sizeof(t));
    mp=b;
    for(int i=;i<N;++i)++t[a[i]&];
    for(int i=;i<B;++i)ls[i]=rs[i]=mp,mp+=t[i];
    for(int i=;i<N;++i)*rs[a[i]&]++=a[i];
    memset(t,,sizeof(t));
    mp=a;
    for(int i=;i<N;++i)++t[b[i]>>&];
    for(int i=;i<B;++i)ls[i]=rs[i]=mp,mp+=t[i];
    for(int i=;i<N;++i)*rs[b[i]>>&]++=b[i];
    memset(t,,sizeof(t));
    mp=b;
    for(int i=;i<N;++i)++t[a[i]>>];
    for(int i=;i<B;++i)ls[i]=rs[i]=mp,mp+=t[i];
    for(int i=;i<N;++i)*rs[a[i]>>]++=a[i];
    //the result is in b
    printf("%g s",(clock()-tt)*./CLOCKS_PER_SEC);
    return ;
}