2251: [2010Beijing Wc]外星联络
Time Limit: 30 Sec
Memory Limit: 256 MB
Submit: 791
Solved: 475
题解:
后缀数组的题,因为要求输出按字典序排序,而后缀数组中的height数组是满足这个条件的,我们只需要利用height数组来找出所有的LCP能向后扩展多少
具体的说一下:
height可以抽象成在水平方向有长度的量,
这样我们把后缀排列后逆时针旋转90度,就
可以把height数组看成一个有高度的东西,这样对于每一个高度只需从左向右求最多能扩展多少就可以了。
图片中红色线的长度就是我们所需要输出的答案
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=3005;
int n,m,p,wa[N],wb[N],SA[N],RK[N],HE[N],c[N],pos[N];
char s[N];
void GetSA()
{
int *x=wa,*y=wb;
m=255;
for(int i=1;i<=n;i++) c[x[i]=s[i]]++;
for(int i=1;i<=m;i++) c[i]+=c[i-1];
for(int i=n;i>=1;i--) SA[c[x[i]]--]=i;
for(int k=1;p<n;k<<=1)
{
p=0;
for(int i=n-k+1;i<=n;i++) y[++p]=i;
for(int i=1;i<=n;i++)
if(SA[i]>k) y[++p]=SA[i]-k;
memset(c,0,sizeof(c));
for(int i=1;i<=n;i++) c[x[y[i]]]++;
for(int i=1;i<=m;i++) c[i]+=c[i-1];
for(int i=n;i>=1;i--) SA[c[x[y[i]]]--]=y[i];
swap(x,y);
p=1;
x[SA[1]]=p;
for(int i=2;i<=n;i++)
{
if(y[SA[i]]!=y[SA[i-1]]||y[SA[i]+k]!=y[SA[i-1]+k]) p++;
x[SA[i]]=p;
}
m=p;
}
}
void GetRK()
{
for(int i=1;i<=n;i++) RK[SA[i]]=i;
}
void GetHE()
{
int j=0;
for(int i=1;i<=n;i++)
{
if(RK[i]==1) continue;
while(s[i+j]==s[SA[RK[i]-1]+j]) j++;
HE[RK[i]]=j;
if(j) j--;
}
}
int main()
{
scanf("%d%s",&n,s+1);
GetSA();
GetRK();
GetHE();
int mx=0;
for(int i=2;i<=n;i++)
{
if(HE[i]>HE[i-1])
for(int j=mx+1;j<=HE[i];j++)
{
int k=i,now=1;
while(HE[k]>=j) now++,k++;
printf("%d\n",now);
}
mx=HE[i];
}
}
Description
小 P 在看过电影《超时空接触》(Contact)之后被深深的打动,决心致力于寻
找外星人的事业。于是,他每天晚上都爬在屋顶上试图用自己的收音机收听外星
人发来的信息。虽然他收听到的仅仅是一些噪声,但是他还是按照这些噪声的高
低电平将接收到的信号改写为由 0 和 1 构成的串, 并坚信外星人的信息就隐藏在
其中。他认为,外星人发来的信息一定会在他接受到的 01 串中重复出现,所以
他希望找到他接受到的 01 串中所有重复出现次数大于 1 的子串。但是他收到的
信号串实在是太长了,于是,他希望你能编一个程序来帮助他。
Input
输入文件的第一行是一个整数N ,代表小 P 接收到的信号串的长度。
输入文件第二行包含一个长度为N 的 01 串,代表小 P 接收到的信号串。
Output
输出文件的每一行包含一个出现次数大于1 的子串所出现的次数。输出的顺
序按对应的子串的字典序排列。
Sample Input
7
1010101
1010101
Sample Output
3
3
2
2
4
3
3
2
2
3
2
2
4
3
3
2
2
HINT
对于 100%的数据,满足 0 <= N <=3000