UVA 1314-Hidden Password（后缀数组）

题目大意：有一串长度为L（5<=L<=100000）的字符串，每次以第i个点起始，循环遍历数组可以得到L串字符串，求字典序最小的，若有多个求起始点最小的。输出起始点。

首先将字符串在原字符串后拼接一次得到长度为2*L的字符串。

用后缀数组计算各个后缀的长度为L的前缀就可以得出答案，计算完毕之后利用x数组可以求得起始点最小的。

首先预处理出得到L的操作，起始数为1，操作包含两种：加1和乘2，对后缀数组的倍增模板修改使得它能完成加1的操作。

#include<stdio.h>
#include<stdlib.h>
char s[200010];
int a[110];
int b[110];
int c[200010];
int e[200010];
int t1[200010];
int t2[200010];
int sa[200010];
int sa_build(int m,int n,int q);
int main(void)
{
	int i,u,p,n,pi,qi,ans;
	scanf("%d",&pi);
	for(qi=0;qi<pi;qi++)
	{
		scanf("%d%s",&n,s);
		for(i=2*n;i>=n;i--)
		{
			s[i]=s[i-n];
		}
		p=n;
		u=0;
		while(p!=1)
		{
			if(p%2==1)
			{
				a[++u]=0;
				p--;
			}
			else
			{
				a[++u]=1;
				p/=2;
			}
		}
		for(i=1;i<=u;i++)
		{
			b[i]=a[u+1-i];
		}
		ans=sa_build(128,2*n+1,u);
		printf("%d\n",ans);
	}
	return 0;
}
int sa_build(int m,int n,int q)
{
	int i,j,k,p,min,sit,*x=t1,*y=t2,*temp;
	for(i=0;i<m;i++) c[i]=0;
	for(i=0;i<n;i++) c[x[i]=s[i]]++;
	for(i=1;i<m;i++) c[i]+=c[i-1];
	for(i=n-1;i>=0;i--)
	{
		c[x[i]]--;
		e[c[x[i]]]=sa[c[x[i]]]=i;
	}
	j=1;
	for(k=1;k<=q;k++)
	{
		j=(b[k]==0)?j+1:j*2;
		p=0;
		if(b[k]==0)
		{
			for(i=n-j+1;i<n;i++) y[p++]=i;
			for(i=0;i<n;i++) if(e[i]>=j-1) y[p++]=e[i]-j+1;
		}
		else
		{
			for(i=n-j/2;i<n;i++) y[p++]=i;
			for(i=0;i<n;i++) if(sa[i]>=j/2) y[p++]=sa[i]-j/2;
		}
		for(i=0;i<m;i++) c[i]=0;
		for(i=0;i<n;i++) c[x[y[i]]]++;
		for(i=1;i<m;i++) c[i]+=c[i-1];
		for(i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
		temp=x; x=y; y=temp;
		p=1; x[sa[0]]=0;
		if(b[k]==0)
		{
			for(i=1;i<n;i++)
			{
				x[sa[i]]=((y[sa[i]]==y[sa[i-1]])&&(y[sa[i]+1]==y[sa[i-1]+1]))?p-1:p++;
			}
		}
		else
		{
			for(i=1;i<n;i++)
			{
				x[sa[i]]=((y[sa[i]]==y[sa[i-1]])&&(y[sa[i]+j/2]==y[sa[i-1]+j/2]))?p-1:p++;
			}
		}
		if(p>=n) break;
		m=p;
	}
	min=(1<<30);
	sit=0;
	for(i=0;i<n/2;i++)
	{
		if(x[i]<min)
		{
			min=x[i];
			sit=i;
		}
	}
	return sit;
}