BUPT Summer Journey #test9 C

 

C. 新来的小妹妹 2014新生暑假个人排位赛09

时间限制 6000 ms     内存限制 65536 KB    

题目描述

Mays王国的女王大人每天过着自由自在的生活，她最大的乐趣就是给邻国的帅气王子写信。而负责给她送信的就是皇家小妹妹快递公司。
今天负责给女王大人送信的是一个新来的小妹妹，她非常好奇女王大人的信，于是悄悄的把它拆开来看了！但是机智的女王大人早就想到了会有这种情况发生，她和邻国帅气王子的信都是加密过的～
小妹妹研究了一路，她感觉，里面重复比较多的内容应该是有用信息。为了安慰自己的智商，小妹妹希望找到信的一个最长连续的子串，这个子串出现2次或以上。为了能找到的子串尽可能长，小妹妹认为即便出现的2次有一部分重叠也是可以的。

输入格式

输入第一行为数据组数T(T<=10)，每组一行字符串str，str中只包含小写字母，且长度不超过2000。

 

输出格式

每组答案输出一行。

输入样例

2
aabbaabbaa
abcde

输出样例

6
0

 

思路：原题即为求原串的最长可重叠的重复子串。可用后缀数组来解决。

代码：

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#define maxn 100005
//#define LOCAL
using namespace std;
string s;
void radix(int *str,int *a,int *b,int n,int m)
{
    static int count[maxn];
    memset(count,0,sizeof(count));
    for(int i=0;i<n;++i) ++count[str[a[i]]];
    for(int i=1;i<=m;++i) count[i]+=count[i-1];
    for(int i=n-1;i>=0;--i) b[--count[str[a[i]]]]=a[i];
}
void suffix_array(int *str,int *sa,int n,int m)
{
    static int rank[maxn],a[maxn],b[maxn];
    for(int i=0;i<n;++i)rank[i]=i;
    radix(str,rank,sa,n,m);
    rank[sa[0]]=0;
    for(int i=1;i<n;++i)rank[sa[i]]=rank[sa[i-1]]+(str[sa[i]]!=str[sa[i-1]]);
    for(int i=0;1<<i <n;++i)
    {
        for(int j=0;j<n;++j)
        {
            a[j]=rank[j]+1;
            b[j]=j+(1<<i)>=n?0:rank[j+(1<<i)]+1;
            sa[j]=j;
        }
        radix(b,sa,rank,n,n);
        radix(a,rank,sa,n,n);
        rank[sa[0]]=0;
        for(int j=1;j<n;++j)
        {
            rank[sa[j]]=rank[sa[j-1]]+(a[sa[j-1]]!=a[sa[j]]||b[sa[j-1]]!=b[sa[j]]);
        }
    }
}
string duplicate_substr(string str)
{
    string rev;
    static int s[maxn],sa[maxn],rank[maxn],h[maxn];
    int n=str.length();
    copy(str.begin(),str.end(),s);
    suffix_array(s,sa,n,256);
    for(int i=0;i<n;i++)rank[sa[i]]=i;
    int k=0;
    int ans1=0,pos1=0;
    for(int i=0;i<n;++i)
    {
        k= k==0?0:k-1;
        while(rank[i]>0 && s[i+k]==s[sa[rank[i]-1]+k])++k;
        h[rank[i]]=k;
        if(h[rank[i]]>ans1)
        {
            ans1=h[rank[i]];
            pos1=i;
        }
    }
    return str.substr(pos1,ans1);
}
int main()
{
    #ifdef LOCAL
    freopen("input.txt","r",stdin);
    #endif // LOCAL
    cin>>s;
    printf("%d\n",duplicate_substr(s).length());
    return 0;
}