51nod-1732 婚姻介绍所(后缀数组)

题目大意:回答任意两个子串的最长公共前缀。

题目分析:后缀数组的模板题。但是输入输出都要外挂。

代码如下:

# include<iostream>
# include<cstdio>
# include<cstring>
# include<algorithm>
using namespace std; const int N=1000; int n,m;
char s[N+5];
int cnt[N+5],SA[N+5];
int rk[N+5],tSA[N+5];
int height[N+5];
int ans[N+5][N+5]; void in(int &x)
{
x=0;
char c=getchar();
while(c<'0'||c>'9'){
c=getchar();
}
x=c-'0';
while(c=getchar()){
if(c<'0'||c>'9') break;
x=x*10+c-'0';
}
} void out(int x)
{
if(x>=10)
out(x/10);
putchar(x%10+'0');
} bool same(int *y,int i,int j,int k)
{
if(y[SA[i]]!=y[SA[j]]) return false;
if(SA[i]+k>=n&&SA[j]+k>=n) return true;
if(SA[i]+k<n&&SA[j]+k>=n) return false;
if(SA[i]+k>=n&&SA[j]+k<n) return false;
return y[SA[i]+k]==y[SA[j]+k];
} void buildSA()
{
m=130;
int *x=rk,*y=tSA;
for(int i=0;i<m;++i) cnt[i]=0;
for(int i=0;i<n;++i) ++cnt[x[i]=s[i]];
for(int i=1;i<m;++i) cnt[i]+=cnt[i-1];
for(int i=n-1;i>=0;--i) SA[--cnt[x[i]]]=i;
for(int k=1;k<=n;k<<=1){
int p=0;
for(int i=n-k;i<n;++i) y[p++]=i;
for(int i=0;i<n;++i) if(SA[i]>=k) y[p++]=SA[i]-k;
for(int i=0;i<m;++i) cnt[i]=0;
for(int i=0;i<n;++i) ++cnt[x[y[i]]];
for(int i=1;i<m;++i) cnt[i]+=cnt[i-1];
for(int i=n-1;i>=0;--i) SA[--cnt[x[y[i]]]]=y[i];
p=1;
swap(x,y);
x[SA[0]]=0;
for(int i=1;i<n;++i){
if(same(y,i,i-1,k)) x[SA[i]]=p-1;
else x[SA[i]]=p++;
}
if(p>=n) break;
m=p;
}
} void getHeight()
{
for(int i=0;i<n;++i)
rk[SA[i]]=i;
height[0]=0;
for(int i=1;i<n;++i){
int k=0;
while(s[SA[i]+k]==s[SA[i-1]+k]) ++k;
height[i]=k;
}
} void getAns()
{
for(int i=0;i<n;++i) ans[i][i]=n-i;
for(int i=0;i<n;++i){
int minn=n;
for(int j=i+1;j<n;++j){
minn=min(minn,height[j]);
ans[SA[i]][SA[j]]=ans[SA[j]][SA[i]]=minn;
}
}
} int main()
{
int q,a,b;
while(~scanf("%d",&n))
{
scanf("%s",s);
buildSA();
getHeight();
getAns();
in(q);
while(q--)
{
in(a);
in(b);
out(ans[a][b]);
puts("");
}
}
return 0;
}

  

上一篇:linux用户及组管理


下一篇:Oracle分组函数cube VS rollup