POJ1226 - Substrings(KMP+二分)

题目大意

给定n个字符串，字符串可逆序可顺序，求它们的最长公共子串

题解

在输入的过程中记录一下最短的那个字符串，然后枚举起点，然后进行二分求出子串末位置，然后再验证是否是公共子串，记录最长的公共子串就是最终答案~~~~时间复杂度为O(N^3*logn)

代码：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
using namespace std;
#define  MAXN 105
#define INF 0x3f3f3f3f
char T[MAXN][MAXN];
int f[MAXN];
void getfail(char *p,int len)
{
    int j;
    f[0]=j=-1;
    for(int i=1;i<len;i++)
    {
        while(j>=0&&p[j+1]!=p[i]) j=f[j];
        if(p[j+1]==p[i])j++;
        f[i]=j;
    }
}
bool find(char *s,int len,int n)
{
    getfail(s,len);
    for(int i=0;i<n;i++)
    {
        int j=-1,lt=strlen(T[i]);
        bool flag1=false,flag2=false;
        for(int t=0;t<lt;t++)
        {
            while(j>=0&&s[j+1]!=T[i][t]) j=f[j];
            if(s[j+1]==T[i][t]) j++;
            if(j+1==len) 
            {
                flag1=true;
                break;
            }
        }
        char ss[MAXN];
        strcpy(ss,T[i]);
        reverse(ss,ss+lt);
        j=-1;
        for(int t=0;t<lt;t++)
        {
            while(j>=0&&s[j+1]!=ss[t]) j=f[j];
            if(s[j+1]==ss[t]) j++;
            if(j+1==len) 
            {
                flag2=true;
                break;
            }
        }
        if(!flag1&&!flag2) return false;
    }
    return true;

}
int main()
{
    int cases;
    scanf("%d",&cases);
    while(cases--)
    {
        int n,maxlen=-INF,len=INF;
        char s[MAXN];
        scanf("%d",&n);
        for(int i=0;i<n;i++)
        {
            scanf("%s",T[i]);
            if(len>strlen(T[i]))
            {
                len=strlen(T[i]);
                strcpy(s,T[i]);
            }
        }
        for(int i=0;i<len;i++)
        {
            int l=i,r=len-1,lens=-INF;
            while(l<=r)
            {
                int mid=(l+r)>>1;
                bool flag=find(s+i,mid-i+1,n);
                if(flag)
                {
                    l=mid+1;
                    lens=mid-i+1;
                }
                else
                    r=mid-1;
                if(lens!=-INF&&lens>maxlen) maxlen=lens;
            }
        }
        if(maxlen!=-INF)
            printf("%d\n",maxlen);
        else
        {
            printf("0\n");
        }
    }
    return 0;
}