重要图片
源代码
int Manacher()
{
int len = Init(); // 取得新字符串长度并完成向 s_new 的转换
int max_len = -1; // 最长回文长度
int id; // 以id中心的最长回文
int mx = 0; // b边界
for (int i = 1; i < len; i++)
{
// 如果i在列表右边界内,那么就比较i的对称点j 和mx-i 的大小
if (i < mx)
p[i] = min(p[2 * id - i], mx - i);
else
p[i] = 1;
while (s_new[i - p[i]] == s_new[i + p[i]]) // 不需边界判断,因为左有'$',右有'\0'
p[i]++;//无法使用if (i < mx)来判断,只能比较i两边的对应值是否相等
// 我们每走一步 i,都要和 mx 比较,我们希望 mx 尽可能的远,
// 这样才能更有机会执行 if (i < mx)这句代码,从而提高效率
if (mx < i + p[i]) // 超过mx边界的i+p[i],边界mx 需要重新赋值
{
id = i;
mx = i + p[i];
}
max_len = max(max_len, p[i] - 1); // 通过比较每个元素的p[i]和max_len的大小,更新max_len的值
}
return max_len;
}
代码分析
p[i] = min(p[2 * id - i], mx - i);
mx-i 和 2*id-i是为了防止i和j 的向左右扩展 ???会不会同时向两边扩展