HDU 5340 Three Palindromes (manacher算法)

题意：给一个20000长度的字符串，问能否拆成3个非空的回文

分析：可以用后缀数组做，但是感觉比较麻烦，这里采用新学的一个manacher算法

(不会的可以移步http://blog.csdn.net/sin_xf/article/details/47271713)

manacher有一个好处就是会记录每个位置为中心的最长回文，那么我们只需要用这个找出所有第一个字符串和第三个字符串，然后用枚举第一个字符串和第三的字符串的所有组合中，找到剩下的中间的字符串是回文就好了，那么怎么判断中间剩下的部分是不是回文呢，只要判断截取前后两端的字符串的中间点的回文半径是不是大于该字符串的一半就好了。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
#include<vector>
#pragma comment(linker,"/STACK:1024000000,1024000000")
using namespace std;
const int maxn = 2e4+5;
char s[maxn*2];                                 //2倍数组
int p[maxn*2],st[maxn*2],ed[maxn*2];
void for_max(int n)                             //manacher算法
{
    for(int j=n;j>=0;j--)                       //字符填充
    {
        s[j*2+2]=s[j];
        s[j*2+1]='#';
    }
    s[0]='$';                                   //防止数组越界填充第一个字符和最后一个字符为‘$’和‘/0’
    int id = 0;                                 //id为上述i
    int maxlen = 0;
    p[0] = 1;
    for(int j=2;j<2*n+2;j++)                    //j为上述i+k,所以k=j-id
    {
        if(p[id]+id<=j) p[j]=1;                 //第一种情况
        else p[j]=min(p[2*id-j],p[id]-(j-id));  //第二种情况的两种小情况
        while(s[j+p[j]]==s[j-p[j]]) p[j]++;     //继续往后找
        if(p[j]+j>p[id]+id) id=j;               //更新最靠近字符串后面的点
    }
}
void solve(int n)
{
    int flag=0;
    int cnt1=0;                                 //记录满足的第一个字符串的个数
    int cnt2=0;                                 //记录满足第三个字符串的个数
    for(int i=2;i<2*n-2;i++){                   //记录所有满足第一个字符串的下标
        if(p[i]==i) st[++cnt1]=i;              
    }
    for(int i=2*n;i>4;i--){                     //记录所有满足第三个字符串的下标
        if(p[i]+i==2*n+2) ed[++cnt2]=i;
    }           
    for(int i=1;i<=cnt1;i++){                   //判断除去第一个字符串和第三个字符串后中间的是否为回文
        for(int j=1;j<=cnt2;j++){
            int l2=p[st[i]]+st[i];
            int r2=ed[j]-p[ed[j]];
            if(l2>r2) continue;
            int mid=(l2+r2)/2;
            if(p[mid]>(r2-l2+1)/2){
                flag=1;
                break;
            }
        }
        if(flag) break;
    }
    if(flag) printf("Yes\n");
    else printf("No\n");
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%s",s);
        int n=strlen(s);
        memset(p,0,sizeof(p));
        for_max(n);
        solve(n);
    }
}