bzoj4278 [ONTAK2015]Tasowanie

http://www.elijahqi.win/2018/02/19/bzoj4278/
Description

给定两个数字串A和B，通过将A和B进行二路归并得到一个新的数字串T，请找到字典序最小的T。

Input

第一行包含一个正整数n(1<=n<=200000)，表示A串的长度。
第二行包含n个正整数，其中第i个数表示Ai。
第三行包含一个正整数m(1<=m<=200000)，表示B串的长度。
第四行包含m个正整数，其中第i个数表示Bi。

Output

输出一行，包含n+m个正整数，即字典序最小的T串。

Sample Input

6
1 2 3 1 2 4
7
1 2 2 1 3 4 3
Sample Output

1 1 2 2 1 2 3 1 2 3 4 3 4
感觉自己开始颓废了?每天就睡睡觉吃吃喝喝变胖 写题写错了也不想改 想改也该部队上来写完之后连续wa两发 后来发现n加的位置写错了..
做法：两个串拼起来 然后注意中间加一个大的数（今天才体会到之前做题为什么要加一个大的数了 然后求下后缀数组 类似双指针按照rk数组贪心的扫一遍即可得到答案
具体证明：http://blog.csdn.net/qq_20669971/article/details/53688333

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 440000
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,s[N],cnt[N],tmp[N],rk[N<<1],rk1[N<<1],sa[N],k;
int main(){
    //freopen("bzoj4278.in","r",stdin);
    n=read();for (int i=1;i<=n;++i) s[i]=read();s[n+1]=1001;
    m=read();for (int i=1;i<=m;++i) s[n+i+1]=read();int st=n+1;
    n+=m+1;for (int i=1;i<=n;++i) cnt[s[i]]=1;m=1001;
    for (int i=1;i<=1001;++i) cnt[i]+=cnt[i-1];
    for (int i=1;i<=n;++i) rk[i]=cnt[s[i]];
    for (int p=1;k!=n;m=k,p<<=1){
        for (int i=1;i<=m;++i) cnt[i]=0;
        for (int i=1;i<=n;++i) ++cnt[rk[i+p]];
        for (int i=1;i<=m;++i) cnt[i]+=cnt[i-1];
        for (int i=n;i;--i) tmp[cnt[rk[i+p]]--]=i;
        for (int i=1;i<=m;++i) cnt[i]=0;
        for (int i=1;i<=n;++i) ++cnt[rk[i]];
        for (int i=1;i<=m;++i) cnt[i]+=cnt[i-1];
        for (int i=n;i;--i) sa[cnt[rk[tmp[i]]]--]=tmp[i];
        memcpy(rk1,rk,sizeof(rk)>>1);rk[sa[1]]=k=1;
        for (int i=2;i<=n;++i){
            if (rk1[sa[i-1]]!=rk1[sa[i]]||rk1[sa[i-1]+p]!=rk1[sa[i]+p]) ++k;
            rk[sa[i]]=k;
        }
    }int l=1,r=st+1;
    while(l<st&&r<=n){
        if (rk[l]<rk[r]) printf("%d ",s[l++]);else printf("%d ",s[r++]);
    }
    while(l<st) printf("%d ",s[l++]);while(r<=n) printf("%d ",s[r++]);puts("");
    return 0;
}