POJ

传送门：POJ - 1226

这个题跟POJ - 3294  和POJ - 3450 都是一样的思路，一种题型。

POJ - 3294的题解可以见：https://www.cnblogs.com/lilibuxiangtle/p/12649853.html

题意：给你n个字符串，求最长子串的长度，要求子串或子串的翻转串在n个字符串中都出现。

题解：把每个字符串和字符串的翻转串连接起来，中间用不同且没出现过的字符隔开，然后再把n个字符串和翻转的串连到一起。然后后缀数组求出sa数组和height数组，用be数组标记每个位置的字符在哪个串中。二分长度，vis数组标记符合条件的串，cnt记录个数，如果cnt>=n则说明该长度可以，反之则不行。（重点！！vis及时标记 我憨憨的wa了好几发，找了快三个小时）

#include<cstdio>
#include<algorithm>
#include<queue>
#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;

//sa:字典序中排第i位的起始位置在str中第sa[i]  sa[1~n]为有效值

//rnk:就是str第i个位置的后缀是在字典序排第几 rnk[0~n-1]为有效值

//height：字典序排i和i-1的后缀的最长公共前缀  height[2~n]为有效值，第二个到最后一个

const int INF=0x3f3f3f3f;
const int maxn = 1e5+10;
int wa[maxn],wb[maxn],wsf[maxn],wv[maxn],sa[maxn];
int rnk[maxn],height[maxn];
char s[maxn];
int r[maxn],n,L;
int be[maxn],vis[110];

int cmp(int *r,int a,int b,int k)
{
    return r[a]==r[b]&&r[a+k]==r[b+k];
}

void getsa(int *r,int *sa,int n,int m)//n为添加0后的总长
{
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0; i<m; i++)  wsf[i]=0;
    for(i=0; i<=n; i++)  wsf[x[i]=r[i]]++;
    for(i=1; i<m; i++)  wsf[i]+=wsf[i-1];
    for(i=n; i>=0; i--)  sa[--wsf[x[i]]]=i;
    p=1;
    j=1;
    for(; p<=n; j*=2,m=p){
        for(p=0,i=n+1-j; i<=n; i++)  y[p++]=i;
        for(i=0; i<=n; i++)  if(sa[i]>=j)  y[p++]=sa[i]-j;
        for(i=0; i<=n; i++)  wv[i]=x[y[i]];
        for(i=0; i<m; i++)  wsf[i]=0;
        for(i=0; i<=n; i++)  wsf[wv[i]]++;
        for(i=1; i<m; i++)  wsf[i]+=wsf[i-1];
        for(i=n; i>=0; i--)  sa[--wsf[wv[i]]]=y[i];
        swap(x,y);
        x[sa[0]]=0;
        for(p=1,i=1; i<=n; i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)? p-1:p++;
    }
}

void getheight(int *r,int n)//n为添加0后的总长
{
    int i,j,k=0;
    for(i=1; i<=n; i++)  rnk[sa[i]]=i;
    for(i=0; i<n; i++){
        if(k)
            k--;
        else
            k=0;
        j=sa[rnk[i]-1];
        while(r[i+k]==r[j+k])
            k++;
        height[rnk[i]]=k;
    }
}

bool check(int k)
{
    memset(vis,0,sizeof(vis));
    int cnt=0;
    for(int i=2;i<=L;i++){
        if(height[i]>=k){
            if(!vis[be[sa[i]]]) cnt++;
            vis[be[sa[i]]]=1;               //就是这个地方卡了我快3个小时！！！
            if(!vis[be[sa[i-1]]]) cnt++;
            vis[be[sa[i-1]]]=1;
        }
        else {
            memset(vis,0,sizeof(vis));
            cnt=0;
        }
        if(cnt>=n) return 1;
    }
    return 0;
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t;
    cin>>t;
    while(t--){
        cin>>n;
        L=0;
        for(int i=0;i<n;i++){
            cin>>s;
            int len=strlen(s);
            for(int j=0;j<len;j++){
                r[L]=s[j];
                be[L++]=i;
            }
            r[L]=i+200;
            be[L++]=i;
            for(int j=len-1;j>=0;j--){
                r[L]=s[j];
                be[L++]=i;
            }
            if(i!=n-1) r[L]=i+n+200,be[L++]=i;
        }
        if(n==1) {
            cout<<strlen(s)<<endl;
            continue;
        }
        r[L]=0;
        getsa(r,sa,L,500);
        getheight(r,L);
        int l=0,r=105,ans=0;
        while(l<=r){
            int mid=l+r>>1;
            if(check(mid)) ans=mid,l=mid+1;
            else r=mid-1;
        }
        cout<<ans<<endl;
    }
    return 0;
}