leetcode-005：Longest Palindromic Substring

// Transform S into T.
// For example, S = "abba", T = "^#a#b#b#a#$".
// ^ and $ signs are sentinels appended to each end to avoid bounds checking
string preProcess(string s) {
    int n = s.length();
    if (n == 0) return "^$";
    string ret = "^";
    for (int i = 0; i < n; i++)
        ret += "#" + s.substr(i, 1);

    ret += "#$";
    return ret;
}

string longestPalindrome(string s) {
    string T = preProcess(s);
    int n = T.length();
    int *P = new int[n];
    int C = 0, R = 0;
    for (int i = 1; i < n-1; i++) {
        int i_mirror = 2*C-i; // equals to i' = C - (i-C)

        P[i] = (R > i) ? min(R-i, P[i_mirror]) : 0;

        // Attempt to expand palindrome centered at i
        while (T[i + 1 + P[i]] == T[i - 1 - P[i]])
            P[i]++;

        // If palindrome centered at i expand past R,
        // adjust center based on expanded palindrome.
        if (i + P[i] > R) {
            C = i;
            R = i + P[i];
        }
    }

    // Find the maximum element in P.
    int maxLen = 0;
    int centerIndex = 0;
    for (int i = 1; i < n-1; i++) {
        if (P[i] > maxLen) {
            maxLen = P[i];
            centerIndex = i;
        }
    }
    delete[] P;

    return s.substr((centerIndex - 1 - maxLen)/2, maxLen);//(centerIndex-(maxlen-1))/2-1;
}

 

问题描述：

    给定一个字符串S=A1A2...An，要求找出其最长回文子串（Longest Palindromic Substring）。所谓回文子串就是S的某个子串Ai...Aj为回文。例如，对字符串S=abcdcbeba，它的回文子串有：bcdcb，cdc，beb，满足题目要求的最长回文子串为bcdcb。

推理思路：

1.由于回文可能由奇数个字符组成，也可能由偶数个字符组成。对奇数回文的处理比较直观，只需要以某个字符为中心，依次向两边扩展即可。因此，我们可以通过如下方式把对偶数回文的处理转换成对奇数回文的处理：在字符边界添加特殊符号。例如，对字符串aba，预处理后变成#a#b#a#；对字符串abba，预处理后变成#a#b#b#a。可以看出，不管是奇数回文，还是偶数回文，在与处理后都变成奇数回文。在找出与预处理后字符串的最长回文后，只需要去除所有的#即为源字符串的最长回文。

2.对寻找字符串某类子串的问题，最简单直观的想法就是穷举出所有子串一一进行判别。这里也不例外，当然时间复杂度也很高，为O（n^3）。

3.对该问题，我们可以进行一定程度的简化处理。既然回文是一种特殊的字符串，我们可以以源字符串的每个字符为中心，依次寻找出最长回文子串P0, P1,...,Pn。这些最长回文子串中的最长串Pi = max(P1, P2,...,Pn)即为所求。请看源码：

view source

print ?

01	string find_lps_native(conststring &str)

02

{

03	intcenter = 0, max_len = 0;

04	for(inti = 1; i < str.length()-1; ++i)

05

{

06

intj = 1;

07

08	//以str[i]为中心，依次向两边扩展，寻找最长回文Pi

09	while(i+j < str.length() && i-j >= 0 && str[i+j] == str[i-j])

10

++j;

11

--j;

12	if(j > 1 && j > max_len)

13

{

14	center = i;

15	max_len = j;

16

}

17

}

18

19	returnstr.substr(center-max_len, (max_len << 1) + 1);

20

}

4.可以看出，上面做法的复杂度为O（n^2）。相比穷举字符串的做法，已经降低了一个量级的复杂度。但是仔细想想，上面的算法还有改进空间吗？当然有！而且改进后能够把复杂度降低到O（n）！这就是大名鼎鼎的Manacher’s Algorithm。请看下文分析：

    举例说明：对字符串S=abcdcba而言，最长回文子串是以d为中心，半径为3的子串。当我们采用上面的做法分别求出以S[1]=a, S[2]=b, S[3]=c, S[4]=d为中心的最长回文子串后，对S[5]=c,S[6]=b...还需要一一进行扩展求吗？答案是NO。因为我们已经找到以d为中心，半径为3的回文了，S[5]与S[3]，S[6]与S[2]...，以S[4]为对称中心。因此，在以S[5],S[6]为中心扩展找回文串时，可以利用已经找到的S[3],S[2]的相关信息直接进行一定步长的偏移，这样就减少了比较的次数（回想一下KMP中next数组的思想）。优化的思想找到了。