manacher算法 - 回文串

manacher算法

历史总是惊人的相似，乍一看还以为没学这里，原来是又双叒叕忘记，那就在这里复习一下

学习的话可以看这里：
互联网侦察微信公众号 - 关于manacher算法的漫画讲解
CDSN博客 - 这是一个很随便的名字
CSDN博客 - dyx心心

代码：

#include <iostream>
#include <algorithm>
#include<string.h>
using namespace std;

const int N = 1e5+10;
int len[N];

char a[N], b[N];// aasfaagabagaatu

int manacher(char *a){

    int mxn, id, k = 0, res = 0;

    int n = strlen(a);

    b[k++] = '@';
    for(int i = 0; i < n; i++){
        b[k++] = '#';
        b[k++] = a[i];
    }
    b[k++] = '#';
    b[k++] = '!';

    for(int i = 0; i < k; i++){
        if(mxn > i) len[i] = min(mxn - i, len[2*id-i]);
        else len[i] = 1;

        while(b[i-len[i]] == b[i+len[i]]) len[i]++;
        if(len[i] + i > mxn) mxn = len[i] + i, id = i;

        res = max(res, len[i]);
    }

    return (res - 1); 
}

int main(){

    cin>> a;
    cout<< manacher(a) <<endl;

    return 0;
}

实现原理

利用回文串的对称性 减少重复匹配的时间

首先是将字符串改造，变成无论原字符串长度是奇数还是偶数都能一样匹配的字符串

例如：
char a[15] = “aasfaagabagaatu”

改造为：char b[15] = “@#a#a#s#f#a#a#g#a#b#a#g#a#a#t#u#!”
与原字符串比较 a = " a a s f a a g a b a g a a t u "

收尾插入不一样的字符 防止越界
相邻两字符之间插入相同的字符，不影响匹配

构造字符串b的len[ i ]数组（以 i 为 中心的回文字符串的半径）， 如图：

其中 ：

1. p₀为某一回文串的中心坐标， p 为p₀回文串匹配到的最远的右坐标（保证p为匹配到的最右坐标）
2. i 为 当前匹配点，j为i关于回文串p₀对称的坐标点 ( j 已经匹配完成)
3. len[j] 为 以j 为中心的回文串的半径长度，len[j] -1为 j 回文串的原长度

当匹配到 i 时， 若 p > i , 说明 i 包含在回文串p₀中，并且关于p₀对称，对称点为 j，那么有两种情况：

1. 在回文串p₀内，len[i] 的长度肯定等于len[j] （如上图1）
2. 若 len[j] > p - i 的长度（即 回文串半径 > i到p的坐标距离长度）, 此时大于的部分无法保证对称性，因此需要一个一个的进行匹配

若 p < i , 也需要一个一个再匹配