BUPT Summer Journey #test11 B

 

B. 田田背课文 2014新生暑假个人排位赛11

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

题目描述

田田被老师要求背诵一个字符串，田田刚刚背了一个子串，说他记性差他还不承认，背下这个子串后却忘了这个子串是从哪里开始的。
现在田田把这个字符串给你，并且让你告诉他至少一次背多长的子串才不可能产生歧义，即不存在另一个子串与它相等

输入格式

输入仅一行，即要求田田背诵的字符串
长度<=10^6

输出格式

输出一行，即不产生歧义串的最小串长

输入样例

abcdabcd

输出样例

5

 

思路：即最长可重叠重复子串个数+1.直接套用后缀数组来解决。

代码：

#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()+1);
    return 0;
}