BZOJ4278 [ONTAK2015]Tasowanie 后缀数组 贪心

题目链接
还是权限题，于是链接是darkbzoj的。

题意：
给你两个数字序列，一个长度为$n$，一个长度为$m$，你每次可以从两个序列的开头拿走一个数字放到新序列的队尾，要求新形成的这个长度为$n+m$的序列字典序最小。输出这个新的序列。$n,m<=2e5$，权值$<=1000$。

题解：
首先由于是要字典序最小，于是我们考虑贪心，如果当前能拿的数一个大一个小，那么我们一定拿这个小的。但是问题是，如果当前的两个数大小一样怎么办呢，特别是，如果接下来两个序列连续的很多个数都是相同的，该怎么办呢？直接暴力找下一个不一样的位置并比较大小的复杂度是爆炸的，那么我们要另找办法。其实可以发现，如果我们把每个数看作一个带权值的字符，那么我们就相当于比较两个字符串的大小。而我们每次比较的字符串其实都是一个后缀，于是我们用后缀数组做一遍后缀排序。这样我们每次在决策的时候只需要比较两个后缀的rank数组的大小，然后贪心地选那个rank靠前的，一定不会更劣。就口胡一个并不严谨的证明，如果一个串它的rank比另一个靠前，那么我们在选出前缀的若干位之后，一定会在某一位处的值比选另一个串优，于是先从这个rank小的串选是一定不会更劣的。

于是这样就做完了。复杂度$O(nlogn)$。

现在都9102年了啊，我学OI时间也不短了，像我这样现在都不会SA的都找不到了啊。感觉都快要退役了，还是什么都不会，也不想多说什么了，好好珍惜与享受自己的OI生涯吧。

代码：

#include <bits/stdc++.h>
using namespace std;

int n,m,sa[400010],rk[400010],s[400010];
int a[400010],b[400010],c[400010];
inline int read()
{
	int x=0;
	char s=getchar();
	while(s>'9'||s<'0')
	s=getchar();
	while(s>='0'&&s<='9')
	{
		x=x*10+s-'0';
		s=getchar();
	}
	return x;
}
int main()
{
	n=read();
	for(int i=1;i<=n;++i)
	a[i]=read();
	m=read();
	a[n+1]=1001;
	for(int i=n+2;i<=n+m+1;++i)
	a[i]=read();
	for(int i=1;i<=n+m+1;++i)
	rk[i]=a[i];
	for(int i=1;i<=n+m+1;++i)
	s[rk[i]]++;
	for(int i=1;i<=1001;++i)
	s[i]+=s[i-1];
	for(int i=m+n+1;i>=1;--i)
	sa[s[rk[i]]--]=i;
	int len=1,p=0,ji=1001;
	while(p<n+m+1)
	{
		int k=0;
		for(int i=n+m+1-len+1;i<=n+m+1;++i)
		b[++k]=i;
		for(int i=1;i<=n+m+1;++i)
		{
			if(sa[i]>len)
			b[++k]=sa[i]-len;
		}
		for(int i=1;i<=n+m+1;++i)
		c[i]=rk[b[i]];
		memset(s,0,sizeof(s));
		for(int i=1;i<=n+m+1;++i)
		s[c[i]]++;
		for(int i=1;i<=ji;++i)
		s[i]+=s[i-1];
		for(int i=n+m+1;i>=1;--i)
		sa[s[c[i]]--]=b[i];
		for(int i=1;i<=n+m+1;++i)
		c[i]=rk[i];
		p=1;
		rk[sa[1]]=1;
		for(int i=2;i<=n+m+1;++i)
		{
			if(!(c[sa[i]]==c[sa[i-1]]&&c[sa[i]+len]==c[sa[i-1]+len]))
			++p;
			rk[sa[i]]=p;
		}
		ji=p;
		len*=2;
	}
	int l=1,r=n+2;
	for(int i=1;i<=n+m;++i)
	{
		if((rk[l]>rk[r]||l>n)&&r<=n+m+1)
		printf("%d ",a[r++]);
		else
		printf("%d ",a[l++]);		
	}
	printf("\n");
	return 0;
}