【后缀数组】【bzoj 1031】: [JSOI2007]字符加密Cipher

http://www.lydsy.com/JudgeOnline/problem.php?id=1031

我太落后了，必须赶进度了

找了好几个模板，最终决定背这个

虽然我没懂，但看起来实现得非常清晰，模块化

交换数组的时候慢了点，hzwer的代码开了两个，交换的时候只交换了下标，但也牺牲了点空间

（我会告诉你我只把最后一部分搞懂了么）

//#define _TEST _TEST
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
/************************************************
Code By willinglive    Blog:http://willinglive.cf
************************************************/
#define rep(i,l,r) for(int i=l,___t=(r);i<=___t;i++)
#define per(i,r,l) for(int i=r,___t=(l);i>=___t;i--)
#define MS(arr,x) memset(arr,x,sizeof(arr))
#define LL long long
#define INE(i,u,e) for(int i=head[u];~i;i=e[i].next)
inline const int read()
{int r=0,k=1;char c=getchar();for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;
for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';return k*r;}
/
const int N=200010;
const int K=256;
char s[N];
int n;

int c[N],sa[N],t1[N],t2[N];
/
void getsa(int *x,int *y,int n,int m)
{
	rep(i,0,m-1) c[i]=0;
	rep(i,0,n-1) c[x[y[i]]]++;
	rep(i,1,m-1) c[i]+=c[i-1];
	per(i,n-1,0) sa[--c[x[y[i]]]]=y[i];
}
void da(char *s,int n,int m)
{
	int i,j,p,*x=t1,*y=t2,*t;
	rep(i,0,n-1) x[i]=s[i],y[i]=i;
	getsa(x,y,n,m);
	for(j=p=1;p<n;j<<=1,m=p)
	{
		p=0;
		rep(i,n-j,n-1) y[p++]=i;
		rep(i,0,n-1) if(sa[i]>=j) y[p++]=sa[i]-j;
		getsa(x,y,n,m);
		t=x,x=y,y=t,p=1,x[sa[0]]=0;
		rep(i,1,n-1) x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+j]==y[sa[i-1]+j]?p-1:p++;
	}
}
/
void input()
{
    scanf("%s",s+1);
    n=strlen(s+1);
    rep(i,1,n) s[n+i]=s[i];
    n<<=1;
}
void solve()
{
	da(s,n+1,K);
	
    rep(i,1,n)
    {
    	if(sa[i]<=n/2) putchar(s[sa[i]+n/2-1]);
    }
}
/
int main()
{
    #ifndef _TEST
    freopen("std.in","r",stdin); freopen("std.out","w",stdout);
    #endif
    input(),solve();
    return 0;
}