HYSBZ 2565 Manacher算法

题意

中文题，不解释

题解

这道题正解应该是回文树，不过用Manacher算法8s水过了。首先用Manacher算法算出来P数组，然后将所有的P数组-1，这样的话便能得到以某一个位置为中心的回文串最长长度。然后我们可以根据这个长度为两个端点赋值。dp1[i]代表从i-1开始向前的最长回文子串，dp2[i]代表从i+1开始最长回文子串。我们可以在i这个位置，得到(dp2[i-p[j]],dp1[i+p[j]])。最后我们可以从1开始遍历，寻找在哪个位置前后回文子串组成的子串长度最长。（当然，有人可能会疑惑这样搞会不会产生一些问题，因为我们可以发现这样组成的串是不包含I的。如果仔细想一下就会发现，这样是没有任何问题的。因为题目要求最长子串，所以最长子串一定是在’#’处得到，不会产生任何问题）

代码

#include <iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<cmath>
#include<queue>
#include<string>
#include<set>
#include<map>
#include<bitset>
#include<stack>
#include<string>
#define UP(i,l,h) for(int i=l;i<h;i++)
#define DOWN(i,h,l) for(int i=h-1;i>=l;i--)
#define W(a) while(a)
#define MEM(a,b) memset(a,b,sizeof(a))
#define LL long long
#define INF 0x3f3f3f3f3f3f3f3f
#define MAXN 120050
#define MOD 1000000007
#define EPS 1e-3
using namespace std;
char s[120010],ts[220010];
int p[220010];
int dp1[220010],dp2[220010];

void solve() {
    MEM(p,0);
    MEM(dp1,0);
    MEM(dp2,0);
    int n=strlen(s);
    int pos=0;
    ts[pos++]='$';
    ts[pos++]='#';
    UP(i,0,n) {
        ts[pos++]=s[i];
        ts[pos++]='#';
    }
    ts[pos]='\0';
    int mx=0,id;
    UP(i,1,pos) {
        if(mx>i) p[i]=min(p[id*2-i],mx-i);
        else p[i]=1;
        W(ts[i+p[i]]==ts[i-p[i]]) p[i]++;
        if(i+p[i]>mx) {
            mx=i+p[i];
            id=i;
        }
    }
    UP(i,1,pos) p[i]--;
    UP(i,1,pos){
        DOWN(j,p[i]+1,1){
            if(i+j>=pos) continue;
            if(dp1[i+j]>=j) break;
            dp1[i+j]=j;
        }
        DOWN(j,p[i]+1,1){
            if(i-j<=0) continue;
            if(dp2[i-j]>=j) break;
            dp2[i-j]=j;
        }
    }
    int ans=-INF;
    UP(i,1,pos){
        if(dp1[i]&&dp2[i]) ans=max(ans,dp1[i]+dp2[i]);
//        cout<<i<<" "<<dp1[i]<<" "<<dp2[i]<<endl;
    }
    printf("%d\n",ans);
}

int main() {
    W(~scanf("%s",s)) {
        solve();
    }
}