hihoCoder#1032 最长回文子串——Manacher算法

这道题在之前的Leet上也做过，后来看到过一些比较形象的图，这里就摘过来了，原地址

首先,在字符串s中,用cnt[i]表示第i个字符的回文半径,即cnt[i]尽可能大,且满足:
s[i-cnt[i],i-1]=s[i+1,i+cnt[i]]
很明显,求出了所有的cnt,就求出了所有的长度为奇数的回文子串.
（至于偶数的怎么求,相比大家也都知道，就是在字符串中间添加特殊字符串，使其变成奇数，那么也就方便用奇数的方法求解了）

假设现在求出了cnt[1..i-1],现在要求后面的cnt值,并且通过前面的操作,得知了当前字符i的cnt值至少为j.现在通过试图扩大j来扫描,求出了cnt[i].再假设现在有个指针k,从1循环到cnt[i],试图通过某些手段来求出[i+1,i+cnt[i]]的cnt值.
根据定义,黑色的部分是一个回文子串,两段红色的区间全等.
因为之前已经求出了cnt[i-k],所以直接用它.有3种情况:

①cnt[i]-k < cnt[i-k]
如图,cnt[i-k]的范围为青色.因为黑色的部分是回文的,且青色的部分超过了黑色的部分,所以cnt[i+k]肯定至少为cnt[i]-k,即橙色的部分.那橙色以外的部分就不是了吗?这是肯定的.因为如果橙色以外的部分也是回文的,那么根据青色和红色部分的关系,可以证明黑色部分再往外延伸一点也是一个回文子串,这肯定不可能,因此cnt[i+k]=cnt[i]-k.为了方便下文,这里的cnt[i+k]=cnt[i]-k=min(cnt[i]-k,cnt[i-k]).

②cnt[i]-k > cnt[i-k]
如图,cnt[i-k]的范围为青色.因为黑色的部分是回文的,且青色的部分在黑色的部分里面,根据定义,很容易得出:cnt[i+k]=cnt[i-k].为了方便下文,这里的cnt[i+k]=cnt[i-k]=min(cnt[i]-k,cnt[i-k]).

根据上面两种情况,可以得出结论:当cnt[i]-k!=cnt[i-k]的时候,cnt[i+k]=min(cnt[i]-k,cnt[i-k]).
注意:当cnt[i]-k==cnt[i-k]的时候,就不同了,这是第三种情况:

如图,通过和第一种情况对比之后会发现,因为青色的部分没有超出黑色的部分,所以即使橙色的部分全等,也无法像第一种情况一样引出矛盾,因此橙色的部分是有可能全等的,但是,根据已知的信息,我们不知道橙色的部分是多长,因此就把i指针移到i+k的位置,j=cnti-k,等下次循环的时候再做了.

———-
整个算法就像上面的那样，下面是AC的代码：

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

const int L =2000050;
char s1[L],s2[L];
int cnt[L];

int main() {
    //freopen("F:/in.txt","r",stdin);
    int num;
    cin >> num;
    int l1,l2;
    for(int t = 0; t < num; ++t) {
    cin >> s1;
        l1 = strlen(s1);
        l2 = 0;
        s2[l2++] = '*';
        s2[l2++] = '#';
        for(int i = 0;i < l1; ++i){
            s2[l2++] = s1[i];
            s2[l2++] = '#';
        }
        s2[l2] = '\0';
        for(int i = 1,j = 0,k; i< l2; i+=k) {
            while(s2[i-j-1] == s2[i+j+1]) ++j;
            cnt[i] = j;
            for(k = 1; k <= cnt[i]&&cnt[i-k]!=cnt[i]-k;++k) {
                cnt[i+k] = min(cnt[i-k], cnt[i]-k);
            }
            j = max(j-k, 0);
        }
        int ans = 0;
        for(int i = 0; i < l2; ++i) {
            if(ans < cnt[i]) ans = cnt[i];
        }

        cout << ans << endl;
    }

    //fclose(stdin);
    return 0;
}