BZOJ2251

【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

********************************************************************************************

题意：01串中，查找所有出现次数大于1的子串出现次数（按照字典序输出）

height[i]表示排名第i和第i-1的前缀的公共前缀长度。这一部分一定次数大于是。

这样依次判断后面的height[]是不是大于height[i]就可以判断对应的串的长度。

但是height[i]对应的串一定要从height[i-1]+1开始，因为那些是height[i-1]对应的串。

看代码吧，有点胡言乱语了！

*******************************************************************************************

#include<bits/stdc++.h>
using namespace std;
const int maxn=3010;
char c[maxn];
int s[maxn],sa[maxn],rank[maxn],saf[maxn],cs[maxn],height[maxn];
int n,m;
void init()
{
    scanf("%d",&n);
    scanf("%s",c);
    for(int i=0;i<n;++i)s[i+1]=c[i];
    s[0]=s[n+1]=-1;
}
void rsort()
{
    for(int i=0;i<=m;++i)cs[i]=0;
    for(int i=1;i<=n;++i)cs[rank[saf[i]]]++;
    for(int i=1;i<=m;++i)cs[i]+=cs[i-1];
    for(int i=n;i>=1;--i)sa[cs[rank[saf[i]]]--]=saf[i];
}
int cmp(int *f,int x,int y,int w)
{
    return f[x]==f[y] && f[x+w]==f[y+w];
}
void suffix()
{
    for(int i=1;i<=n;++i)rank[i]=s[i],saf[i]=i;
    m=127;rsort();
    for(int p=1,w=1,i;p<n;w<<=1,m=p)
    {
        for(i=n-w+1,p=0;i<=n;++i)saf[++p]=i;
        for(i=1;i<=n;++i)if(sa[i]>w)saf[++p]=sa[i]-w;
        rsort();swap(rank,saf);rank[sa[1]]=p=1;
        for(i=2;i<=n;++i)rank[sa[i]]=cmp(saf,sa[i],sa[i-1],w)?p:++p;
    }
    int j,k=0;
    for(int i=1;i<=n;height[rank[i++]]=k)
        for(k=k?k-1:k,j=sa[rank[i]-1];s[i+k]==s[j+k];++k);
}
void work()
{
    for(int i=0;i<=n;++i)    //i:paiming
    {   //j:changdu
        for(int j=1+height[i];;j++)
        {   int js=1;
            for(int k=i+1;height[k]>=j;++k,++js);
            if(js>1)printf("%d
",js);
            else break;
        }
    }
}
int main()
{
    init();
    suffix();
    work();
    return 0;
}