BZOJ 2565 最长双回文串 Manacher

BZOJ 2565 最长双回文串 Manacher

---

题目

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

Input

一行由小写英文字母组成的字符串S。

Output

一行一个整数，表示最长双回文子串的长度。

Sample Input

baacaabbacabb

Sample Output

12

题解

这道题是我写了一整天的一道题，也就是说我几乎一天都在写这题。。。由于我非常的弱，我后来才发现可以在Manacher的过程中求出一些很有用的东西，下文将会予以论述，而我要是和Manacher一起求出来的话那么感觉就会非常的好

就是，这道题其实非常的水，我们只要先Manacher一下，然后求出一个位置的左边的右端点为这个字符的最长回文串的长度，不妨存入 Left 数组，右侧同理，存入 Right 数组，然后我们求 Ans=max(Ans,Left[i]+Right[i+1]),i∈(0,len−1) 即可

代码

让注释记录我思维无力的挣扎

#include
#include
#include
#include
#define maxn 100005
using namespace std;
int f[maxn*2],len1,pre[maxn*2],suf[maxn*2],len;
char txt[maxn],s[maxn*2];
void Manacher(){
    for(int i=0;i1]='@';
        s[(i<<1)+1]=txt[i];
    }
    s[(len1)<<1]='@';
    len=2*len1+1;
    int mx=0,p=0;
    for(int i=0;iif(mx>i) f[i]=min(mx-i,f[2*p-i]);
        else f[i]=0;
        while(i+f[i]+11>=0&&s[i+f[i]+1]==s[i-f[i]-1]) f[i]++;
        if(f[i]+i>mx){
            mx=f[i]+i;
            p=i;
        }
    }
    return;
}
void process(){
    for(int i=0;iif(i&1){
            f[i]=((f[i]>>1)<<1)+1;
        }else{
            f[i]=((f[i]+1)>>1)<<1;
        }
    }
    int now=-1;
    for(int i=1;i1;i++){
        if(i&1){
            if((i>>1)>now){
                pre[++now]=1;
            }
            while(((i>>1)+(f[i]>>1))>now){
                now++;
                pre[now]=((now-(i>>1))<<1)+1;
            }
        }
        else{
            if((i>>1)>now){
                pre[++now]=1;
            }
            while(((i>>1)+(f[i]>>1)-1)>now){
                now++;
                pre[now]=(now-(i>>1)+1)<<1;
            }
        }
    }
    now=len1;
    for(int i=len-2;i>0;i--){
        if(i&1){
            if((i>>1)1;
            }
            while(((i>>1)-(f[i]>>1))>1)-now)<<1)+1;
            }
        }
        else{
            if((i>>1)1;
            }
            while(((i>>1)-(f[i]>>1))>1))<<1;
            }
        }
    }
    /*int now=0;
    for(int i=0;if[2*i+1]+i){
            while(f[2*i]-1+i>now){
                pre[now]=((now-i)<<1)+1;
                now++;
            }
        }
        else{
            while(f[2*i+1]+i>now){
                pre[now]=(now-i)<<1;
                now++;
            }
        }
    }
    pre[0]=1;
    now=len1-1;
    for(int i=len1-1;i>=0;i--){
        if(i-f[2*i]+1
    /*int now=0;
    for(int i=0;inow){
            if(i>=now){
                pre[now]=1;
                now++;
                continue;
            }
            if(!(i&1))pre[now]=(((now-i)>>1)<<1)+1;
            else pre[now]=(((now-i+1)>>1)<<1);
            now++;
        }
    }
    now=len-1;
    for(int i=len-1;i>=0;i--){
        while(i-f[i]>1)<<1)+1;
            else suf[now]=(((i-now+1)>>1)<<1);
            now--;
        }
    }
    suf[0]=1;*/
    /*int now=0;
    for(int i=0;inow){
            now++;
            pre[now]=((now-i)<<1)+1;
        }
    }
    for(int i=0;i>=1;
        else pre[i]=(pre[i]+1)>>1;
    }
    now=len-1;
    for(int i=len-1;i>=0;i--){
        while(i-f[i]=0;i--){
        if(!(i&1))suf[i]>>=1;
        else suf[i]=(suf[i]+1)>>1;
    }*/
    int ans=0;
    for(int i=0;i1;i++){
        ans=max(ans,pre[i]+suf[i+1]);
    }
    printf("%d",ans);
    return;
}
int main(){
    /*freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);*/
    scanf("%s",txt);
    len1=strlen(txt);
    Manacher();
    process();
    return 0;
}

细节

我实在不知道怎么说，满屏的细节QAQ

总结

要灵活地运用Manacher算法
我真是太弱了