manacher算法（对付最大回文子串的利器），O(n)级别，快狠准

 

https://segmentfault.com/a/1190000008484167?utm_source=tag-newest

https://blog.csdn.net/liuwei0604/article/details/50414542

 

 

 

 

 

p[i] = min(p[2 * id - i], mx - i);

s[ ]    # a # b # a # b # a #

p[i]    1 2 1 4 1 6 1 4 1 2 1

值也是对称存在的，即：p[i] = p[2 * id - i]   同时  存在  p[i] = mx - i 的情况，所以推理为此。

 

i	0	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19
s_new[i]	$	#	a	#	b	#	b	#	a	#	h	#	o	#	p	#	x	#	p	#
p[i]		1	2	1	2	5	2	1	2	1	2	1	2	1	2	1	4	1	2	1

可以看出，p[i] - 1正好是原字符串中最长回文串的长度。

接下来的重点就是求解 p 数组，如下图：

设置两个变量，mx 和 id 。mx 代表以 id 为中心的最长回文的右边界，也就是mx = id + p[id]。

假设我们现在求p[i]，也就是以 i 为中心的最长回文半径，如果i < mx，如上图，那么：

if (i < mx)  
    p[i] = min(p[2 * id - i], mx - i);

2 * id - i为 i 关于 id 的对称点，即上图的 j 点，而p[j]表示以 j 为中心的最长回文半径，因此我们可以利用p[j]来加快查找

。

由于回文串范围内前后的p[i]的值相等    

   

我们可以根据它的 对称的值 或者 mx-i的长度值 取得其的最小值，

再在最小值得基础上进行p[i]++；向外扩充它的长度值   

之后进行mx的更新和求取最大值！

（mx值开始是0，刚刚开始时p[1]为0，mx才为2，进行求大值；

之后每次的id，为上一次的或者之前的i值，因此i总是在id之前）

 

#include <iostream>
#include <cstring>
#include <algorithm>
#include<stdio.h>

using namespace std;

char s[1000];
char s_new[2000];
int p[2000];

int Init()
{
    int len = strlen(s);
    s_new[0] = '$';
    s_new[1] = '#';
    int j = 2;

    for (int i = 0; i < len; i++)
    {
        s_new[j++] = s[i];
        s_new[j++] = '#';
    }

    s_new[j] = '\0';  // 别忘了哦

    return j;  // 返回 s_new 的长度
}

int Manacher()
{
    int len = Init();  // 取得新字符串长度并完成向 s_new 的转换
    int max_len = -1;  // 最长回文长度

    int id;
    int mx = 0;

    for (int i = 1; i < len; i++)
    {
        if (i < mx)
            p[i] = min(p[2 * id - i], mx - i);  // 需搞清楚上面那张图含义, mx 和 2*id-i 的含义
        else
            p[i] = 1;

        while (s_new[i - p[i]] == s_new[i + p[i]])  // 不需边界判断，因为左有'$',右有'\0'
            p[i]++;

        // 我们每走一步 i，都要和 mx 比较，我们希望 mx 尽可能的远，这样才能更有机会执行 if (i < mx)这句代码，从而提高效率
        if (mx < i + p[i])
        {
            id = i;
            mx = i + p[i];
        }

        max_len = max(max_len, p[i] - 1);
    }

    return max_len;
}

int main()
{
    while (printf("请输入字符串：\n"))
    {
        scanf("%s", s);
        printf("最长回文长度为 %d\n\n", Manacher());
    }
    return 0;
}