bzoj 2565: 最长双回文串 manachar

题意

顺序和逆序读起来完全一样的串叫做回文串。比如acbca是回文串，而abc不是（abc的顺序为“abc”，逆序为“cba”，不相同）。
输入长度为n的串S，求S的最长双回文子串T,即可将T分为两部分X，Y，（|X|,|Y|≥1）且X和Y都是回文串。
2≤|S|≤10^5

分析

直接用manachar求出f[i]表示以i为结尾的最长回文串，g[i]表示以i为开头的最长回文串，然后扫一遍即可。

代码

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

const int N=200005;

int n,w[N],f[N],a[N],g[N],len[N];
char s[N];

void manachar()
{
    for (int i=1;i<=n;i++) a[i*2]=s[i]-'a',a[i*2-1]=30;
    a[n*2+1]=30;
    int mx=0,pos=0;
    for (int i=1;i<=n*2+1;i++)
    {
        if (mx>=i) len[i]=min(mx-i+1,len[pos*2-i]);
        else len[i]=1;
        while (i-len[i]>0&&i+len[i]<=n*2+1&&a[i-len[i]]==a[i+len[i]])
        {
            len[i]++;
            if (i+len[i]-1>mx) w[i+len[i]-1]=len[i]*2-1;
        }
        if (i+len[i]-1>mx) mx=i+len[i]-1,pos=i;
    }
}

int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);
    manachar();
    for (int i=1;i<=n;i++) f[i]=w[i*2+1]/2;
    for (int i=1;i<=n/2;i++) swap(s[i],s[n-i+1]);
    manachar();
    for (int i=1;i<=n;i++) g[n-i+1]=w[i*2+1]/2;
    int ans=0;
    for (int i=1;i<n;i++) ans=max(ans,f[i]+g[i+1]);
    printf("%d",ans);
    return 0;
}