【BZOJ2251】[2010Beijing Wc]外星联络

Description

小 P 在看过电影《超时空接触》(Contact)之后被深深的打动，决心致力于寻
找外星人的事业。于是，他每天晚上都爬在屋顶上试图用自己的收音机收听外星
人发来的信息。虽然他收听到的仅仅是一些噪声，但是他还是按照这些噪声的高
低电平将接收到的信号改写为由 0 和 1 构成的串，并坚信外星人的信息就隐藏在
其中。他认为，外星人发来的信息一定会在他接受到的 01 串中重复出现，所以
他希望找到他接受到的 01 串中所有重复出现次数大于 1 的子串。但是他收到的
信号串实在是太长了，于是，他希望你能编一个程序来帮助他。

Input

输入文件的第一行是一个整数N ，代表小 P 接收到的信号串的长度。
输入文件第二行包含一个长度为N 的 01 串，代表小 P 接收到的信号串。

Output

输出文件的每一行包含一个出现次数大于1 的子串所出现的次数。输出的顺
序按对应的子串的字典序排列。

Sample Input

7
1010101

Sample Output

3
3
2
2
4
3
3
2
2

HINT

对于 100%的数据，满足 0 <= N <=3000

题解：先求出sa和height，然后只要沿着sa枚举每一个子串，在height里找出height不小于当前子串长度的范围。为了防止重复，每次枚举从height+1开始。

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int maxn=3010;
int ra[maxn],rb[maxn],st[maxn],sa[maxn],v[maxn],rank[maxn],h[maxn];
char str[maxn];
int n,m;
void work()
{int i,j,k,p,*x=ra,*y=rb;for(i=0;i<n;i++)    st[x[i]=v[i]]++;for(i=1;i<m;i++)    st[i]+=st[i-1];for(i=n-1;i>=0;i--)    sa[--st[x[i]]]=i;for(j=1,p=1;p<n;j<<=1,m=p){for(p=0,i=n-j;i<n;i++)    y[p++]=i;for(i=0;i<n;i++)    if(sa[i]>=j)    y[p++]=sa[i]-j;for(i=0;i<m;i++)    st[i]=0;for(i=0;i<n;i++)    st[x[y[i]]]++;for(i=1;i<m;i++)    st[i]+=st[i-1];for(i=n-1;i>=0;i--)    sa[--st[x[y[i]]]]=y[i];for(swap(x,y),p=1,x[sa[0]]=0,i=1;i<n;i++)x[sa[i]]=(y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+j]==y[sa[i]+j])?p-1:p++;}for(i=0;i<n;i++)    rank[sa[i]]=i;for(i=k=0;i<n-1;h[rank[i++]]=k)for(k?k--:0,j=sa[rank[i]-1];v[i+k]==v[j+k];k++);
}
int main()
{scanf("%d%s",&n,str);int i,j,r;for(i=0;i<n;i++)    v[i]=str[i]-'0'+1;m=3;n++;work();for(i=1;i<n;i++){for(j=h[i]+1;sa[i]+j<n;j++){//网上的题解这里还向下求了一下区间的左端点l，感觉不明觉厉 for(r=i+1;r<=n&&h[r]>=j;r++);if(r-i>=2)    printf("%d\n",r-i);    }}return 0;
}

转载于:https://www.cnblogs.com/CQzhangyu/p/6268708.html