作用：
在线性时间内找到一个字符串的最大回文子串
原理：
奇偶变换：为处理字符串方便，现将给定的任意字符串进行处理，使所有可能的奇数/偶数长度的回文子串都转换成了奇数长度。
具体就是在每个字符的两边都插入一个特殊的符号。比如hhjj变成 #h#h#j#j#， aba变成 #a#b#a#；
为防止数组越界，可以在字符串的开始加入另一个特殊字符，比如“?#a#b#a#?” 。
定义一个辅助数组int p[]，p[i]表示以s_new[i]为中心的最长回文的半径，例如：
可以看出，p[i]-1正好是原字符串中最长回文串的长度。
Manacher 算法之所以快，就快在对 p 数组的求法上有个捷径。
if(i < mx)//mx为以id为中心的最远扩展回文串的右边界,i为扫描到的点
{p[i] = min(p[2 * id - i], mx - i);
}
注解在代码里,例题链接：http://acm.hdu.edu.cn/showproblem.php?pid=3068
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define mem(a, b) memset(a, b, sizeof(a))
using namespace std;char str[110010];
char new_s[110010 * 2]; //mannacher数组因为奇偶变换加了符号进去,数组要开原字符串2倍大小
int p[110010 * 2];void get_news(int len)//构造新字符串
{new_s[0] = '^';new_s[1] = '#';int j = 1;for(int i = 0; i < len; i ++){new_s[++ j] = str[i];new_s[++ j] = '#';}j ++;new_s[j] = '\0';
}int mannacher(int len)
{mem(p, 0);int mx = 0;int id = 0;int maxlen = -1;//答案,最长字符串长度初始化 for(int i = 1; i < len; i ++){if(i < mx)//mx为以id为中心的最远扩展回文串的右边界,i为扫描到的点 {p[i] = min(p[2 * id - i], mx - i);//需搞清楚上面那张图含义, mx和2*id-i的含义}else //i >= mx 不能利用对称性来直接求解p[i]了, 所以先p[i] = 1,然后暴力更新 p[i] = 1;while(new_s[i - p[i]] == new_s[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(scanf("%s", str)!=EOF){int len = strlen(str);get_news(len);len = strlen(new_s);int ans = mannacher(len);printf("%d\n", ans);}return 0;
}
