hdu4763——KMP算法变型

hdu4763题目

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int NN=1000005;
char B[NN];
int prefix[NN];
int lenB,T;
void getprefix()//即nexts数组，next[i]里面存放的是0-i个里面前缀和后缀相等的最大长度
{
    memset(prefix,0,sizeof(prefix));
    prefix[0]=-1;
    prefix[1]=0;//如果读者不理解的话可以初始化i=0和len=-1带下面的循环看看不管是什么条件都是 prefix[1]=0；在意义上理解是属于长度为1的字符串它的next始终为0
    int i=1,j;
    int len=0;
    while(i<lenB)
    {
        if(len==-1||B[i]==B[len])//如果第len位和第i位（即上一个的最大前缀的最后一位和这次的长度为i的字符串增加的后缀最后一位）
        {
            len++;
            i++;
            prefix[i]=len;
        }
        else len=prefix[len];
        //这个有点难懂，可以模拟一下，实际上就是上一个的最大前缀和后缀的长度是len，既然在上面的if中B[i]!=B[len]则长度为i+1的字符串前缀和后缀长度必定<=len
        //因此，prefix[len]就代表了从字符串头部开始的长度为len的的子串的前后缀相同的最长长度
        //------------------------------------------------------------------
        //是因为prefix[i+1]的值是 B[i]这个子串中，B[0~k-1]=B[i+1-k~i]里面最大的k（即len）
        //我们可以把B[0~k-1]分为B[0~k-2]和B[k-1]; B[i+1-k~i]分为B[i+1-k~i-1]和B[i]
        //即B[0~k-2]和B[i-k+1~i-1]比较，B[k-1]和B[i]比较
        //下面感觉不太好理解……姑且看成是从结论推向结果
        //B[0~k-1] =  B[i-k~i](prefix[i+1]=k的意义)
        //B[1~k-1] vs B[i-k+1~i-1]
        //B[1~k-1] vs B[0~k-2]
        //即得到len=prefix[len]；
    }
}
void KMP()
{
    getprefix();
    int i=0,j=0;
    bool cnt=0;
    int k;
    for(k=prefix[lenB];k&&cnt==0;k=prefix[k])//k的长度从整个字符串相等前后缀子串的最大长度开始,即满足了EAEBE中的前后E
    {
        for(j=2*k; j<=lenB-k; j++)//EAEBE中最小的EAE长度为2*k，j代表EAE的长度，满足EAE模式
        {
            if(prefix[j]==k)//如果外部的"EAEBE"首尾E最大：k 和内部EAE首尾最大prefix[j]相等就找到啦
            {
                cnt=1;
                printf("%d\n",k);
                break;
            }
        }
    }
    if(!cnt)cout<<"0"<<endl;
}
int main()
{
    scanf("%d",&T);
    for(int i=0; i<T; i++)
    {
        scanf("%s",&B);
        lenB=strlen(B);
        KMP();
    }
    return 0;
}