【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;

}

外星