manacher 算法

今天有一题要O(n)求回文串然后再匹配一下， 所以就学了一下这个manacher 算法， 以前好像也听过的， Apio有题好像也要用这个

manacher算法:

定义数组p[i]表示以i为中心的(包含i这个字符)回文串半径长

将字符串s从前扫到后for(int i=0;i<strlen(s);++i)来计算p[i],则最大的p[i]就是最长回文串长度,则问题是如何去求p[i]?

由于s是从前扫到后的,所以需要计算p[i]时一定已经计算好了p[1]....p[i-1]

假设现在扫描到了i+k这个位置,现在需要计算p[i+k]

定义maxlen是i+k位置前所有回文串中能延伸到的最右端的位置,即maxlen=p[i]+i;//p[i]+i表示最大的

分两种情况:

1.i+k这个位置不在前面的任何回文串中,即i+k>maxlen,则初始化p[i+k]=1;//本身是回文串

然后p[i+k]左右延伸,即while(s[i+k+p[i+k]] == s[i+k-p[i+k]])++p[i+k]

2.i+k这个位置被前面以位置i为中心的回文串包含,即maxlen>i+k

这样的话p[i+k]就不是从1开始

由于回文串的性质,可知i+k这个位置关于i与i-k对称,

所以p[i+k]分为以下3种情况得出

//黑色是i的回文串范围,蓝色是i-k的回文串范围,

/上面是转载的

这篇已经讲的够清楚的了， 分类的非常清晰，

直接写出来就是这样的：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#define MAXN 110005
using namespace std;
char s1[MAXN], s[MAXN * 2];
int p[MAXN * 2], len;
int main()
{
    while(scanf("%s", s1) != EOF){
        s[0] = '*'; s[1] = '#';
        len = strlen(s1);
        for(int i = 0; i < len; i ++){
            s[i + i + 2] = s1[i];
            s[i + i + 3] = '#';    
        }
        int k = 0, ans = 0;
        for(int i = 2; i <= len * 2; i ++){
            if(k + p[k] - 1 < i)p[i] = 1;
            else{
                if(p[2 * k - i] + i - 1 > k + p[k] - 1)
                    p[i] = p[k] - i + k;
                else{
                    if(p[2 * k - i] + i - 1 < k + p[k] - 1)
                        p[i] = p[2 * k - i];
                    else p[i] = p[2 * k - i];    
                }
                //p[i] = min(p[2 * k - i], p[k] + k - i);
            }
            while(s[i + p[i]] == s[i - p[i]])p[i] ++;
            if(k + p[k] < i + p[i]) k = i;
            ans = max(ans, p[i]);
        }  cout<<ans - 1<<endl;
    }
    //system("pause");
    return 0;
}

但是还可以化简， 因为我们仔细考虑这个过程， 其实它在每次更新答案时的限定因素只有两个： 一个是p[k] + k - 1， 也就是当前访问过的最大节点， 另一个是 p[2 * k - i]， 也就是在大的回文串中与i 对应的那个点， 仔细想一想其实p[i] 的当前最大值一定是它们两个中的最小值， 因为如果大于了p[k] + k - 1, 就是没有访问过的节点， 一定不能确定， 而如果大于了 p[2 * k - i]， 一种情况会到达没有访问过的节点， 如果没有到达没有访问过的节点（也就是i + p[i] (p[i] > p[2 * k - i])   < p[k] + k - 1）， p[2 * k - i]一定会比现在的值大， 矛盾。

化简完就是我注释的那一行， 

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#define MAXN 110005
using namespace std;
char s1[MAXN], s[MAXN * 2];
int p[MAXN * 2], len;
int main()
{
    while(scanf("%s", s1) != EOF){
        s[0] = '*'; s[1] = '#';
        len = strlen(s1);
        for(int i = 0; i < len; i ++)
            s[i + i + 2] = s1[i], s[i + i + 3] = '#';    
        int k = 0, ans = 0;
        for(int i = 2; i <= len * 2; i ++){
            if(k + p[k] - 1 < i)p[i] = 1;
            else p[i] = min(p[2 * k - i], p[k] + k - i);
            while(s[i + p[i]] == s[i - p[i]])p[i] ++;
            if(k + p[k] < i + p[i]) k = i;
            ans = max(ans, p[i]);
        }  cout<<ans - 1<<endl;
    }
    return 0;
}

hdu3068  和 hdu3294,  直接模板。

poj 2974

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#define MAXN 2000005
using namespace std;
char s[MAXN];
int p[MAXN];
int main()
{
    int cnt = 0, ans;
    while(scanf("%s", s)){
        if(s[0] == 'E')break;
        int len = strlen(s);
        for(int i = len - 1; i >= 0; i --){
            s[i + i + 2] = s[i];
            s[i + i + 3] = '#';
        }s[0] = '*';s[1] = '#';
        int k = 0; ans = 0;
        for(int i = 2; i <= len + len; i ++){
            if(k + p[k] <= i)p[i] = 1;
            else p[i] = min(k + p[k] - i, p[2 * k - i]);
            while(i - p[i] > 0 && s[i + p[i]] == s[i - p[i]])p[i] ++;
            if(i + p[i] > k + p[k])k = i;
            ans = max(ans, p[i]);    
        }
        printf("Case %d: %d\n", ++ cnt, ans - 1);
    }
    return 0;
}