51nod 1732 51nod婚姻介绍所

江老板说是后缀数组模板题，然后把模板给了我。然后发现 k1==k2 时要特判。然后交上去超时，然后加了读入优化和输出优化。。。。。

板子是江老板的，压缩饼干模式。

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

const int N=2010;
int n,m,Rank[N],sa[N],h[N],tmp[N],cnt[N],ans;char s[N],t[N];
void suffixarray(int n,int m){
    int i,j,k;n++;
    for(i=0;i<2*n+5;i++)Rank[i]=sa[i]=h[i]=tmp[i]=0;
    for(i=0;i<m;i++)cnt[i]=0;
    for(i=0;i<n;i++)cnt[Rank[i]=s[i]]++;
    for(i=1;i<m;i++)cnt[i]+=cnt[i-1];
    for(i=0;i<n;i++)sa[--cnt[Rank[i]]]=i;
    for(k=1;k<=n;k<<=1){
        for(i=0;i<n;i++){
            j=sa[i]-k;
            if(j<0)j+=n;
            tmp[cnt[Rank[j]]++]=j;
        }sa[tmp[cnt[0]=0]]=j=0;
        for(i=1;i<n;i++){
            if(Rank[tmp[i]]!=Rank[tmp[i-1]]||Rank[tmp[i]+k]!=Rank[tmp[i-1]+k])cnt[++j]=i;
            sa[tmp[i]]=j;
        }memcpy(Rank,sa,n*sizeof(int));
        memcpy(sa,tmp,n*sizeof(int));
        if(j>=n-1)break;
    }for(j=Rank[h[i=k=0]=0];i<n-1;i++,k++)
    while(~k&&s[i]!=s[sa[j-1]+k])h[j]=k--,j=Rank[sa[j]+1];
}
/*
    sa[i]表示排名第i的后缀开始的位置
    h[i]表示sa[i]和sa[i-1]的最长公共前缀
    Rank[i]求从第i为开始的后缀的排名
*/
int f[N][30],lg2[N];
void rmq_init(int n){  
	lg2[0]=lg2[1]=0;
    for(int i=2;i<=n;i++)lg2[i]=lg2[i/2]+1;
    for(int i=1;i<=n;i++)f[i][0]=h[i];
    for(int j=1;(1<<j)<=n;j++)  
    for(int i=1;i+(1<<j)-1<=n;i++)
        f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
} 
int rmq_min(int l,int r){
    if(l>r)swap(l,r);l++;
    int k=lg2[r-l+1];
    return min(f[l][k],f[r-(1<<k)+1][k]);
}

void read(int&a){
    char ch;while(!((ch=getchar())>='0')&&(ch<='9'));
    a=ch-'0';while(((ch=getchar())>='0')&&(ch<='9'))a*=10,a+=ch-'0';
}

inline void prin_d(int x)
{
    if (x > 9)
    {
        prin_d(x / 10);
    }
    putchar(x % 10 + '0');
    return ;
}

int main()
{
	int q,k1,k2;
	while(scanf("%d",&n)!=EOF)
	{
		scanf(" %s",s);
		suffixarray(n,'z');
		rmq_init(n);
		read(q);
		while(q--)
		{
			read(k1);
			read(k2);
			if(k1==k2)
				prin_d(n-k1);
			else
				prin_d(rmq_min(Rank[k1],Rank[k2]));
			puts("");
		}
	}
}