[leetcode] 5. 最长回文子串

题目描述

给你一个字符串 s，找到 s 中最长的回文子串。

如果字符串的反序与原始字符串相同，则该字符串称为回文字符串。

示例 1：
输入：s = “babad”
输出：“bab”
解释：“aba” 同样是符合题意的答案。

示例 2：
输入：s = “cbbd”
输出：“bb”

提示：
1 <= s.length <= 1000
s 仅由数字和英文字母组成

解题方法

Manacher算法

Mannacher算法能够在$O(N)$的时间内找到最长回文子串，所以重点介绍下。
首先我们需要将字符串字符之间以‘#’字符（可以是其它字符，只要是字符串中不存在的字符皆可）填充。
举个例子，假设原字符串为“babbac”，以‘#’填充后，则变为“#b#a#b#b#a#c#”。可以看到这个字符串的最长回文子串是“abba”，但是它的中心点在两个b之间，因为中心点不存在，所以不好通过中心点遍历的方式求出最长回文子串；而填充之后的最长回文子串就变为了“#a#b#b#a#”，这样以两个b之间的‘#’作为中心点遍历时，就可以得到最长回文子串（当然最后还需要将填充的‘#’去掉）。

将字符串填充之后，我们需要用一个数组rs来保存从中心向两边最多扩展的字符个数，即 rs[i] 表示从当前字符出发分别能往左右两边扩展 rs[i]个字符（包括i位置），还是以“#b#a#b#b#a#c#”为例，则rs数组如下图所示。

不难发现rs[i] * 2 -1即为以该字符位置为中心点时的最长回文长度。那么问题的关键点就变为了怎么求rs数组。

以下部分可以结合java代码进行理解。

我们用C代表回文字符串的中心，用R表示我们能够扩展的最大右边界坐标（边界值取不到)，那么R = C + rs[C]。我们需要从起始位置i开始遍历，计算rs[i]的值，同时不断更新C和R的值。

i从起始位置开始遍历。当i < R时, r[i]的初始值可以取以C为中心点的另一边镜像位置即rs[C * 2 - i]和R-i中的较小值。然后比较cc[i - rs[i]]和cc[i + rs[i]]两边位置的字符是否相等，相等则r[i]++，扩大两边边界比较后续字符；不相等则rs[i]的值就固定了。r[i]的值固定后。判断i+r[i]是否大于右边界R，若大于，则中心点C变为i，右边界的值也跟着变为i+r[i]；否则，R与C的值都不变。然后i++，重复上面步骤，直到遍历结束。

为何r[i]的初始值可以取以C为中心点的另一边镜像位置即rs[C * 2 - i]和R-i中的较小值，这一步可能有童鞋不太明白，我简单讲一下。

以上图为例子，此时i的镜像位置是j，因为r[j]< R - i，则rs[i]的初始值为rs[j]。因为L到R之间的字符（不包括L和R）都是以C为中心的镜像，且r[j]的长度也在L到R之间，所以rs[i]的初始长度可取rs[j]。那么可能会有童鞋说，若图上r[j]的值为5，r[j]的长度超过了左边界L，那么r[i]的初始值就只能取到R-i即4了，因为以j为中心的回文字符串有一部分的字符超过了L，不在以C为中心的镜像范围之内，r[i]初始的右边界只能在R之内，超出R(包括R)之后，不再适用于以C为中心的镜像，此时需要以i为中心扩展边界比较两端的字符。这就解释了为何r[i]的初始值可以取以C为中心点的另一边镜像位置即rs[C * 2 - i]和R-i中的较小值。

java代码

public class Solution {
    public String longestPalindrome(String s) {
        if (s == null || s.length() <= 1) {
            return s;
        }
        char[] c = s.toCharArray();
        char[] cc = getManaC(c);
        int[] rs = getManacher(cc);
        int maxLen = 0;
        int maxIndex = 0;
        for (int i = 0; i < rs.length; i++) {
            if (rs[i] > maxLen) {
                maxLen = rs[i];
                maxIndex = i;
            }
        }
        int start = maxIndex + 1 - maxLen;
        int end = start + (2 * maxLen - 1);
        StringBuilder bd = new StringBuilder();
        for (int i = start; i < end; i++) {
            if (cc[i] != '#') {
                bd.append(cc[i]);
            }
        }
        return bd.toString();
    }

    public int[] getManacher(char[] cc) {
        int[] rs = new int[cc.length];
        int C = -1;
        int R = -1;
        for (int i = 0; i < cc.length; i++) {
            rs[i] = i < R ? Math.min(rs[C * 2 - i], R - i) : 0;
            while (i - rs[i] >= 0 && i + rs[i] < rs.length) {
                if (cc[i - rs[i]] != cc[i + rs[i]]) {
                    break;
                }
                rs[i]++;
            }
            if (i + rs[i] > R) {
                C = i;
                R = i + rs[i];
            }
        }
        return rs;
    }

    public char[] getManaC(char[] c) {
        int len = c.length * 2 + 1;
        char[] cc = new char[len];
        int index = 0;
        for (int i = 0; i < len; i++) {
            if ((i & 1) == 0) {
                cc[i] = '#';
            } else {
                cc[i] = c[index++];
            }
        }
        return cc;
    }
}

时间复杂度：$O(N)$
空间复杂度：$O(N)$

---

• 个人公众号
  
• 个人小游戏