Problem Q: 公共回文长度

http://acm.ccniit.com/problem.php?cid=1205&pid=16

Description

给定两个长度为n的字符串s1和s2，求最长公共回文串的长度（即选取一个区间，该区间在s1和s2中都是回文区间且该区间是最大的）

 

Input

第一行一个n（n<=1000000）
第二行一个长度为n的字符串s1
第三行一个长度为n的字符串s2

Output

最长公共回文串的长度(即选取的最大区间的长度)

Sample Input Copy

5
abcba
defea

Sample Output Copy

3

HINT

选取区间2-4
s1中为bcb，回文长度为3
s2中为efe，回文长度为3

题解：它这个不仅要相同最大回文串长度，还要区间也要一样，所以我们可以用马拉车来解决，不知道马拉车的可以在网上搜一下，用两个数组分别记录两个字符串每个位置的最大回文串长度，然后比较相同位置上的最大回文串长度，取小的，然后在每次选出的最小的值中找最大值，就可以找到满足条件的结果。

代码实现：

#include<iostream>
#include<cstring>
using namespace std;
const int N=1000009;
int n,ks1=2,ks2=2;
char s1[2*N],ss1[2*N],s2[2*N],ss2[2*N];
int p[2*N],p1[2*N];
void mlc1()
{
    int maxx=0,id=0;
    for(int i=2;i<ks1;i++)
    {
        if(p[id]+id>i)
        {
            p[i]=min(p[id*2-i],p[id]+id-i);
        }
        else
        {
            p[i]=1;
        }
        while(ss1[i-p[i]]==ss1[i+p[i]])
        {
                p[i]++;
        }
        if(p[i]+i>p[id]+id)
        {
            id=i;
        }
        if(maxx<p[i]-1)
        {
            maxx=p[i];
        }
    }
}
void mlc2()
{
   int maxx=0,id=0;
    for(int i=1;i<ks2;i++)
    {
        if(p1[id]+id>i)
        {
            p1[i]=min(p1[id*2-i],p1[id]+id-i);
        }
        else
        {
            p1[i]=1;
        }
        while(ss2[i-p1[i]]==ss2[i+p1[i]])
        {
                p1[i]++;
        }
        if(p1[i]+i>p1[id]+id)
        {
            id=i;
        }
        if(maxx<p1[i]-1)
        {
            maxx=p1[i];
        }
    }
}
int main()
{
    cin>>n;
    cin>>s1;
    cin>>s2;
    ss1[0]=ss1[1]='#';
    for(int i=0;i<n;i++)
    {
        ss1[ks1++]=s1[i];
        ss1[ks1++]='#';
    }
    ss2[0]=ss2[1]='#';
    for(int i=0;i<n;i++)
    {
        ss2[ks2++]=s2[i];
        ss2[ks2++]='#';
    }
    mlc1();
    mlc2();
    int ans=0;
    for(int i=2;i<n*2+2;i++)
    {
        if(p[i]<=p1[i])
        {
            ans=max(ans,p[i]);
        }
        else
        {
            ans=max(ans,p1[i]);
        }
    }
    printf("%d",ans-1);
}