Manacher算法--回文查找

部分转载自https://segmentfault.com/a/1190000008484167

1.思想
首先在字符串中每隔一个字符插入一个字符（前提这个字符未出现在串里），无论是奇偶字符串都转换成奇字符串，这样跳过奇偶的分类。
定义一个辅助数组int p[]，p[i]表示以s_new[i]为中心的最长回文的半径。在求p[i]时，首先为它初始化。设置两个变量，mx 和 id 。mx 代表以s_new[id]为中心的最长回文最右边界，也就是mx=id+p[id]。假设我们现在求p[i]，也就是以s_new[i]为中心的最长回文半径，如果i小于mx，那么p[i] = min(p[2 * id - i], mx - i);

2.代码

/**
 * 
 * author 刘毅（Limer）
 * date   2017-02-25
 * mode   C++ 
 */
#include<iostream>  
#include<string.h>
#include<algorithm>  
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 maxLen = -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]++;


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

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

    return maxLen;
}

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

    return 0;
}