从零学算法5

5.给你一个字符串 s，找到 s 中最长的回文子串。
如果字符串的反序与原始字符串相同，则该字符串称为回文字符串。
示例 1：
输入：s = “babad”
输出：“bab”
解释：“aba” 同样是符合题意的答案。
示例 2：
输入：s = “cbbd”
输出：“bb”

• 暴力解法。无非就是双重循环，截取出所有子串，判断是否为回文串，然后取最大值。稍微优化一下，如果此时子串长度小于当前最长的回文子串 max 就直接跳过，当前子串长度为 j - i +1

•   String s;public String longestPalindrome(String s) {this.s = s;int n = s.length();// start：最长回文子串起始下标int max = 0,start = 0;for(int i=0;i<n;i++){for(int j=i;j<n;j++){if(j-i+1 < max)continue;if(isPalindrome(i,j)){max = Math.max(max, j-i+1);start = i;}}}return s.substring(start, start+max);}public boolean isPalindrome(int i,int j){// 根据回文串对称的性质从左右两端往中间找，有不一样的就不是while(i<j){if(s.charAt(i++)!=s.charAt(j--))return false;}return true;}
  

• 中心扩散法。遍历 s 的每个下标，以其为中心点向两边扩散，看是否为相同字符，就能得到以该点为中心的回文串。但是回文串不一定为奇数长度，即中心点处长度不一定为 1，比如 abccba，中心点处其实为 cc，所以就在扩散时定义初始左右起点为这两个字符的下标即可。

•   public String longestPalindrome(String s) {int n = s.length();int max = 0,start = 0;for(int mid=0;mid<n;mid++){if(n-mid<max/2)break;int i=mid,j=mid;//有相同的点就移动初始右端点// 为什么不需要移动左端点?看下面 mid = j 部分while(j<n-1 && s.charAt(j)==s.charAt(j+1))j++;// i~j 这一段的点作为中心处已经在这一轮考虑过，所以之后直接跳过即可// 同时这也保证了下一轮 i 初始化为 mid 时,i 的左边不会和 i 处字符重复// 相当于在上一轮就移动了左端点，遍历顺序为从左到右，所以首轮左处无端点无需处理mid = j;// 向两端延伸求最大回文串长度while(i>0 && j<n-1 && s.charAt(i-1)==s.charAt(j+1)){i--;j++;}// 更新最大值和最大回文子串起始下标if(j-i+1>max){max = j-i+1;start=i;}}return s.substring(start, start+max);}
  

• 动态规划法。根据上面中心扩散法判断回文串的方式，可以大致看到动态规划的雏形。假定 boolean dp[left][right] 表示子串 s[left:right] 是否为回文串。那么从 dp[left][right] 为 true 开始，如果 s.charAt(left-1)==s.charAt(right+1)，就可以得到 dp[left-1][right+1] 为 true，即 dp[l][r] = dp[l+1][r-1] && s.charAt(l-1)==s.charAt(r+1)。得到递推公式以后还需要确定一下边界条件。以下 s.charAt(left) 我就简写成 s[left] 了，如果 s[left]!=s[right]，说明当前子串左右两端字符不相同，那么你怎么都不会是回文串；如果相同，那么当子串长度小于等于 3 时，你一定是回文串，比如 aba,bb,a，否则就根据 dp[left + 1][right - 1] 递推。
• 考虑到dp[left][right] 依赖 dp[left + 1][right - 1] 的特性，比如 dp[0][3]，要先知道 dp[1][2] ，那 left 在外层遍历是肯定不可能先得到 dp[1][x] 再得到 dp[0][x] 的，所以遍历顺序需要调整为外层为 right，内层为 left

•   public String longestPalindrome(String s) {int n = s.length();boolean[][] dp = new boolean[n][n];int max = 0,start = 0;for(int j=0;j<n;j++){for(int i=0;i<=j;i++){if(s.charAt(i)!=s.charAt(j))continue;if(j-i<3)dp[i][j]=true;else dp[i][j]=dp[i+1][j-1];if(dp[i][j] && j-i+1>max){max = j-i+1;start=i;}}}return s.substring(start, start+max);}
  

• 长度为 1 的子串肯定是回文串，所以还可以根据这点稍微优化一下

•   int max = 1,start = 0;for(int j=1;j<n;j++){for(int i=0;i<j;i++){
  

• 马拉车算法（Manacher）：回文串有奇数 aba 和偶数 abba 两种形式，使用马拉车算法，会在字符串的每两个字符间以及字符首尾加上一个特殊字符，这样最终字符串长度都会为奇数，同时，原本的每个回文串的长度也会都变成奇数，以上面两个字符串为例就可以证明
• aba->*a*b*a*
• abba->*a*b*b*a*
• 其实就是奇数加偶数必定为奇数
• 再引用一个变量回文半径，表示回文串左边或右边到中心点的长度，比如 *a*b*a* 回文半径为 *a*b 或者 b*a* 的长度也就是 4， *a*b*b*a* 回文半径为 5
• 铺垫完毕，进入正题，中心扩散法每次换一个中心点，都要重新计算此时以 i 点为中心点的回文串长度，而马拉车在中心扩散法的基础上改进，可以利用之前的结果。
• 因为遍历顺序是从左到右，所以我们取之前的回文串结果中右边界最大的一个，利用它推出当前以 i 为中心的回文串长度。我们就假设之前的那个结果中心点为 maxCenter，左右端点为 left 和 maxRight。我们根据 i 的位置划分出不同情况下怎么推出当前结果

  1. i < maxRight，那么以 maxCenter 为对称中心对称过去的点 j 肯定在 left 到 maxCenter，如果以 j 为中心点得到的回文串也在 left 到 maxCenter 之内，那根据对称的性质，i 对应的结果和 j 对应的一样，下图顺便说明了得到 j 点的计算公式，把 left 带入 0 为例子比较容易推导， p[] 为回文串半径数组
     

  2. i < maxRight，但是以 j 为中心的回文串长度超过了 left，还是根据对称，我起码可以肯定我的回文串半径最小也等于 j - left ，也就是 maxRight - i，我们以此为基础用中心扩散法继续扩散即可
     

  3. i > maxRight，那没办法了，老老实实从头开始中心扩散法吧

• i 为什么没有小于 maxCenter 的情况，因为 maxCenter 就是我们之前不断遍历得到的某个 i，新的 i 当然大于之前的 i 了

•   int charLen = s.length();//源字符串的长度int length = charLen * 2 + 1;//添加特殊字符之后的长度char[] chars = s.toCharArray();//源字符串的字符数组char[] res = new char[length];//添加特殊字符的字符数组int index = 0;//添加特殊字符for (int i = 0; i < res.length; i++) {res[i] = (i % 2) == 0 ? '#' : chars[index++];}//新建p数组 ，p[i]表示以res[i]为中心的回文串半径int[] p = new int[length];//maxRight(某个回文串延伸到的最右边下标)//maxCenter(maxRight所属回文串中心下标),//resCenter（记录遍历过的最大回文串中心下标）//resLen（记录遍历过的最大回文半径）int maxRight = 0, maxCenter = 0, resCenter = 0, resLen = 0;//遍历字符数组resfor (int i = 0; i < length; i++) {if (i < maxRight) {//情况一，i没有超出范围[left,maxRight]//2 * maxCenter - i其实就是j的位置，实际上是判断p[j]<maxRight - i// maxRight - i 其实就是 j-left，也就是说 j 的回文半径在 left ~ jif (p[2 * maxCenter - i] < maxRight - i) {//j的回文半径没有超出范围[left,maxRight]，直接让p[i]=p[j]即可p[i] = p[2 * maxCenter - i];} else {//情况二，j的回文半径已经超出了范围[left,maxRight]，我们可以确定p[i]的最小值//是maxRight - i，至于到底有多大，后面还需要在计算// i - p[i] 表示以 i 为中心， p[i] 为半径的回文串的左端点，同理 i + p[i] 为右端点// 半径不断增加就等于中心点不断向外扩散p[i] = maxRight - i;while (i - p[i] >= 0 && i + p[i] < length && res[i - p[i]] == res[i + p[i]])p[i]++;}} else {//情况三，i超出了范围[left,maxRight]，就没法利用之前的已知数据，而是要一个个判断了p[i] = 1;while (i - p[i] >= 0 && i + p[i] < length && res[i - p[i]] == res[i + p[i]])p[i]++;}// 匹配完之后，如果右边界i + p[i]超过maxRight，那么就更新maxRight和maxCenter// 因为主要看 i 相对右边界来推导结果，所以能更新右边界就更新if (i + p[i] > maxRight) {maxRight = i + p[i];maxCenter = i;}//记录最长回文串的半径和中心位置if (p[i] > resLen) {resLen = p[i];resCenter = i;}}//计算最长回文串的长度和开始的位置resLen = resLen - 1;int start = (resCenter - resLen) >> 1;//截取最长回文子串return s.substring(start, start + resLen);
  

• 如下图所示，有特殊字符的回文半径 - 1 其实就是原始回文串的长度，这也就是 resLen = resLen - 1; 的由来
  
• 上面三种情况还能合并，把他们都看成确定半径以及扩散这两步即可，情况 1 确定完半径以后不满足扩散条件所以不会扩散，情况 2 和 3 也就是初始半径不同，都会扩散

•  int charLen = s.length();//源字符串的长度int length = charLen * 2 + 1;//添加特殊字符之后的长度char[] chars = s.toCharArray();//源字符串的字符数组char[] res = new char[length];//添加特殊字符的字符数组int index = 0;//添加特殊字符for (int i = 0; i < res.length; i++) {res[i] = (i % 2) == 0 ? '#' : chars[index++];}//新建p数组 ，p[i]表示以res[i]为中心的回文串半径int[] p = new int[length];//maxRight(某个回文串延伸到的最右边下标)//maxCenter(maxRight所属回文串中心下标),//resCenter（记录遍历过的最大回文串中心下标）//resLen（记录遍历过的最大回文半径）int maxRight = 0, maxCenter = 0, resCenter = 0, resLen = 0;//遍历字符数组resfor (int i = 0; i < length; i++) {//合并后的代码p[i] = maxRight > i ? Math.min(maxRight - i, p[2 * maxCenter - i]) : 1;//上面的语句只能确定i~maxRight的回文情况，至于maxRight之后的部分是否对称，//就看是否需要一个个去匹配了，匹配的时候首先数组不能越界while (i - p[i] >= 0 && i + p[i] < length && res[i - p[i]] == res[i + p[i]])p[i]++;//匹配完之后，如果右边界i + p[i]超过maxRight，那么就更新maxRight和maxCenterif (i + p[i] > maxRight) {maxRight = i + p[i];maxCenter = i;}//记录最长回文串的半径和中心位置if (p[i] > resLen) {resLen = p[i];resCenter = i;}}//计算最长回文串的长度和开始的位置resLen = resLen - 1;int start = (resCenter - resLen) >> 1;//截取最长回文子串return s.substring(start, start + resLen);