后缀数组poj1226

题意：求最长的公共子串，使她后者是他的翻转串在每个串中都出现。

思路：把每个串和他的翻转串拼接中间用没有出现过的字符隔开，然后再把这些串用没有出现过的字符拼接，然后二分判断。（对字符编号的时候要注意）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<algorithm>
using namespace std;
const int maxn=30000;
int sa[maxn],height[maxn],Rank[maxn],c[maxn],idx[maxn],t[maxn],t2[maxn];
int str[maxn];
bool vis[maxn];
int n,len;

void add(int x,int pos)
{
    idx[len]=pos;
    str[len++]=x;
}

void build_sa(int m,int n)
{
    int *x=t,*y=t2;
    for(int i=0;i<m;i++)c[i]=0;
    for(int i=0;i<n;i++)c[x[i]=str[i]]++;
    for(int i=1;i<m;i++)c[i]+=c[i-1];
    for(int i=n-1;i>=0;i--)sa[--c[x[i]]]=i;
    for(int k=1;k<=n;k<<=1)
    {
        int p=0;
        for(int i=n-k;i<n;i++)y[p++]=i;
        for(int i=0;i<n;i++)if(sa[i]>=k)y[p++]=sa[i]-k;
        for(int i=0;i<m;i++)c[i]=0;
        for(int i=0;i<n;i++)c[x[y[i]]]++;
        for(int i=1;i<m;i++)c[i]+=c[i-1];
        for(int i=n-1;i>=0;i--)sa[--c[x[y[i]]]]=y[i];
        swap(x,y);
        x[sa[0]]=0;p=1;
        for(int i=1;i<n;i++)
            x[sa[i]]=((y[sa[i-1]]==y[sa[i]])&&(y[sa[i-1]+k]==y[sa[i]+k])?p-1:p++);
        if(p>=n)break;
        m=p;
    }
}
void getheight(int n)
{
    int k=0;
    for(int i=1;i<=n;i++)Rank[sa[i]]=i;
    for(int i=0;i<n;i++)
    {
        if(k)k--;
        int j=sa[Rank[i]-1];
        while(str[i+k]==str[j+k])k++;
        height[Rank[i]]=k;
    }
}

bool can(int l,int r)
{
    if(r-l<n)return false;
    int cnt=0;
    memset(vis,0,sizeof(vis));
    for(int i=l;i<r;i++)
    {
        int pos=idx[sa[i]];
        if(!vis[pos]&&pos!=n){cnt++,vis[pos]=1;}
    }
    return cnt>=n;
}
bool judge(int mid)
{
    int l=0;
    memset(vis,0,sizeof(vis));
    for(int i=2;i<=len;i++)
    {
        if(height[i]<mid||i==len)
        {
            if(can(l,i))return true;
            l=i;
        }
    }
    return false;
}

void solve()
{
    int l=0,r=101,mid;
    while(l<r)
    {
        mid=l+(r-l+1)/2;
        if(judge(mid))l=mid;
        else r=mid-1;
    }
    printf("%d\n",l);
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        len=0;
        char word[210];
        int sp=140;
        for(int i=0;i<n;i++)
        {
            scanf("%s",word);
            int tmp=strlen(word);
            for(int j=0;j<tmp;j++)add(word[j]-'A'+1,i);
            add(sp++,n);
            for(int j=0;j<tmp;j++)add(word[tmp-j-1]-'A'+1,i);
            add(sp++,n);
        }
        add(0,n);
        if(n==1)
        {
            printf("%d\n",strlen(word));
            continue;
        }
        memset(Rank,0,sizeof(Rank));
        build_sa(300,len);
        len--;
        getheight(len);
        solve();
    }
    return 0;
}