hdu 4763 扩展的KMP

一开始没有看清楚串的长度，一百万啊，我就用后缀数组做，结果肯定T了，一看是一百万的数量级，就想到了扩展的KMP。

next[i]: 以第i位置开始的子串 与 母串的公共前缀。用了这个数组，我们就可以在线性时间内求解了。

n为串的长度，那么结果k的最大值不会超过n/3;

我们首先算出next[n/3]~next[n-n/3-n/3]区间内的最大值MAX，然后枚举k的长度，每次k减一的时候只需更新MAX就可以了，

而MAX只需和next[k],next[n-k-k],next[n-k-k-1]比较。从而在线性时间内可以求解。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 1001010;
int next[N],extand[N];
int n,MAX;
void getnext(char *T){// next[i]: 以第i位置开始的子串 与 T的公共前缀
     int i,length = strlen(T);
     next[0] = length;
     for(i = 0;i<length-1 && T[i]==T[i+1]; i++);
          next[1] = i;
     int a = 1;
     for(int k = 2; k < length; k++){
        int p = a+next[a]-1, L = next[k-a];
        if( (k-1)+L >= p ){
           int j = (p-k+1)>0? (p-k+1) : 0;
           while(k+j<length && T[k+j]==T[j]) j++; // 枚举(p+1，length) 与(p-k+1,length) 区间比较
           next[k] = j, a = k;
        }
        else next[k] = L;
     }
}
bool judge(int k)
{
    MAX=max(MAX,next[k]);
    MAX=max(MAX,next[n-2*k]);
    MAX=max(MAX,next[n-2*k-1]);
    if(next[n-k]<k) return 0;
    int l=k,r=n-k-k;
    if(l>r) return 0;
    if(MAX>=k) return 1;
    return 0;
}
int work(int len)
{
    int i;
    for(i=len;i<=n-2*len;i++)
    {
        MAX=max(MAX,next[i]);
        if(MAX>=len&&next[n-len]>=len) return len;
    }
    len--;
    while(len)
    {
        if(judge(len)) break;
        else len--;
    }
    return len;
}
int main()
{
    int i,j,T,ans;
    char t[N];
    scanf("%d",&T);
    while(T--)
    {
        scanf("%s",t);
        n=strlen(t);
        if(n<3) {printf("0\n");continue;}
        getnext(t);
        MAX=0;
        ans=work(n/3);
        printf("%d\n",ans);
     }
     return 0;
}