bzoj 2565——最长双回文串（manacher）

点击打开原题目链接

题意：给一个字符串，求这个字符串中最长的双回文串的长度。（双回文串：一个能被分为两个回文串的字符串。）

题解：最长回文串当然是直接上manacher，但是这道题和普通的manacher有区别，因为要判断的是两个连续的回文串的最长长度。首先还是用最原始的manacher处理出以每个位置为中心的回文串的长度len[i]。因为要求两个连续的回文串，那么一定满足第一个回文串的最后一个位置和第二个回文串的第一个位置是挨着的，并且这两个回文串的长度之和要最大。所以我们可以在每一个'#'上维护两个值：一个是以这个位置为左端点的回文串的长度 l[i] ，一个值是以这个位置为右端点的回文串的长度 r[i] ，那么我们可以对每一个位置进行遍历，当这个位置上的 l[i] 和 r[i]都不为1时，说明这个位置的左边和右边都有回文串，那就可以组成一个双回文串，在遍历的过程中不断更新最大的双回文串的长度，最后得到的长度就是要求的答案了。至于这里为什么是在每一个'#'上面维护 l[i] 和 r[i] 的值呢？联系一下manacher的构造方式就很容易明白了，因为每两个原字符串中都有一个'#'，所以每个'#'不影响当前最长回文串的长度,所以在一个'#'上记录左右是否同时有回文串是最方便的。

这里再细讲一下维护 l[i] 和 r[i] 的过程。

对于每一个位置i，manacher维护出来的len[i]表示以 i 为中心的最长回文串的半径，不算上'#'的最长回文串的长度就是 len[i]-1 (证明很简单：首先在转换得到的字符串T中，所有的回文字串的长度都为奇数，那么对于以 i 位置为中心的最长回文字串，其长度就为2* len[i]-1,经过观察可知，T中所有的回文子串，其中 '#' 的数量一定比其他字符的数量多1，也就是有 len[i]个 '#' ，剩下 len[i]-1个字符来自原字符串，所以该回文串在原字符串中的长度就为 len[i]-1)，那么以i为中心的最长回文串的左端点就是：i - len[i]+1,右端点的位置就是： i+ len[i]-1，那么以左端点为起点的最长回文串的长度就等于以右端点为终点的最长回文串，l[i - len[i]+1]=r[i - len[i]+1]=len[i]*2-1;如果是左端点前一位当起点的最长回文串呢？那就是去掉当前左端点和当前右端点的以i为中心的那个回文串啊，长度为当前左端点为起点的长度减2,即：l[i+1]=l[i]-2;右端点也是同样的规律，按照这个规律我们可以处理一个回文串中所有子位置的回文串长度，那么如果两个回文串有重合的部分我们也可以计算他们合成的字符串的长度。如果这个位置本来就是一个回文串远的那个左端点，那么就是它本身，与它前面的字符无关，所以我们取 l[i] 和 l[i-1]-2之间最大的那个作为l[i]的值。

实现代码：

 for (int i=1; i<=len; i++)
        r[i]=max(r[i],r[i-1]-2);//左端点，由前一位推过来的，如果这个位置是在一个回文串中的话，那前面的
    for (int i=len; i>=1; i--)
        l[i]=max(l[i],l[i+1]-2);

整体代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=1000010;
char s[maxn],c[maxn];
int lens,len[maxn],l[maxn],r[maxn];
int Manachar(int lens,char s[])//manacher模板
{
    c[0]='$';
    c[1]='#';
    for (int i=0; i<lens; i++)
        c[i*2+2]=s[i],c[i*2+3]='#';
    lens=lens*2+1;
    int j,mx=1;
    len[0]==0;
    for (int i=1; i<=lens; i++)
    {
        if (mx>i)
            len[i]=min(len[j*2-i],mx-i);
        else
            len[i]=1;
        while(c[i+len[i]]==c[i-len[i]])
            len[i]++;
        if (i+len[i]>mx)
        {
            mx=i+len[i];
            j=i;
        }
        l[i+len[i]-1]=max(l[i+len[i]-1],len[i]*2-1);//维护右端点
        r[i-len[i]+1]=max(r[i-len[i]+1],len[i]*2-1);//维护左端点
    }
    for (int i=1; i<=lens; i++)
        r[i]=max(r[i],r[i-1]-2);//左端点，由前一位推过来的，如果这个位置是在一个回文串中的话，那前面的
    for (int i=lens; i>=1; i--)
        l[i]=max(l[i],l[i+1]-2);
    int maxx=0;
    for (int i=1; i<=lens; i++)
        if (l[i]!=1&&r[i]!=1)
            maxx=max(maxx,l[i]+r[i]);
    return maxx;
}
int main()
{
    scanf("%s",s);
    lens=strlen(s);
    cout<<Manachar(lens,s)/2-1<<endl;
    return 0;
}

刚学会manacher，做的第一道模板题，讲的自己都晕乎乎的，要是有讲错的地方欢迎指正。