算法-最长回文子串

问题概述

---

题目：算法-最长回文子串

问题描述：
输入一个字符串，求出其中最长的回文子串。子串的含义是：
在原串中连续出现的字符串片段。回文的含义是：正着看和倒着看相同。
如abba和yyxyy。在判断回文时，应该忽略所有标点符号和空格，且忽略大小写，但输出应保持原样（在回文串的首部和尾部不要输出多余字符）。
输入字符串长度不超过5000，且占据单独的一行。
应该输出最长的回文串，如果有多个，输出起始位置最靠左的。

输入：
一行字符串，字符串长度不超过5000。
输出：
字符串中的最长回文子串。
样例输入：
Confuciuss say:Madam,I’m Adam.
样例输出：
Madam,I’m Adam

---

1. 问题分析

---

1.1 算法核心

若是使用枚举算法，算法复杂度可达O_n³，在面对一些n较大的题中很明显是不能满足时间要求的，所以本文分享的是动态规划实现求最长回文字符串，选择是最容易理解的一种动态规划方式，实际上解决这个问题的方法有很多种，最优秀的当属复杂度为O_n的Manacher算法，这个就不在此介绍了。

---

1.2 算法分析

令dp[i][j]表示s[i]到[j]所表示的子串是否是回文子串，是则为1，不是则为0，这样根据s[i]是否等于s[j]，可以把转移情况分为两类：

1. 若s[i] == s[j]，那么只要s[i +1][j - 1]是回文子串，s[i]和s[j]就是回文子串；若s[i + 1][j - 1]不是回文子串，那么s[i]和s[j]也不是回文子串；
2. 若s[i] != s[j]，那么s[i]到s[j]一定不是回文子串；

由此可以写出状态转移方程：

边界：dp[i][j] = 1, dp[i][i + 1] = (s[i] == s[i + 1])? 1 : 0;
到这一步还有一个问题没有解决，那就是如果按照i和j从小到大的顺序来枚举子串的两个端点，然后更新dp[i][j]，会无法保证dp[i + 1][j - 1]已经被计算过，从而无法得到正确的dp[i][j];
针对这个问题，我们采取按子串的长度和子串的初始位置进行枚举的方法来解决，即第一遍将长度为3的子串的dp值全部求出，第二遍通过第一遍结果计算出长度为4的子串的dp值…这样就可以避免状态无法转移的问题。

---

2. 解决方案

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <string.h>
using namespace std;

#define maxlength 5001

int dp[maxlength][maxlength] = { 0 }, t[maxlength], n, j = 0, res;
char A[maxlength], B[maxlength];

int main(void)
{
	fgets(A, 5001, stdin);
	n = strlen(A) - 1;
	for (int i = 0; i < n; i++)
	{
		if (isalpha(A[i]) || isdigit(A[i]))
		{
			if (isalpha(A[i]))
				B[j] = toupper(A[i]);
			else
				B[j] = A[i];
			t[j++] = i;
		}
	}
	n = strlen(B);
	for (int i = 0; i < n - 1; i++)
	{
		dp[i][i] = 1;
		if (B[i] == B[i + 1])
		{
			dp[i][i + 1] = 1;
			res = 2;
		}
	}
	dp[n - 1][n - 1] = 1;

	int maxi = 0, maxj = 0;

	for (int L = 3; L <= n; L++)
	{
		for (int j = n - 1; j - L + 1 >= 0; j--)
		{
			int i = j - L + 1;
			if (B[i] == B[j] && dp[i + 1][j - 1] == 1)
			{
				dp[i][j] = 1;
				res = L;
				maxi = i;
				maxj = j;
			}
		}
	}
	int maxij = 0;
	for (int i = 0; i < n; i++)
	{
		for (int j = n - 1; j >= 0; j--)
		{
			if (dp[i][j] == 1 && (j - i) > maxij)
			{
				maxi = i;
				maxj = j;
				maxij = (j - i);
			}
		}
	}
	for (int i = t[maxi]; i <= t[maxj]; i++)
		printf("%c", A[i]);
	printf("\n");

	return 0;
}

3. 资源分享

旗木白哉のGitHub：
https://github.com/YangMengHeng
GitHub/Gitee源代码下载地址：
GitHub
https://github.com/YangMengHeng/myCode/tree/master/VSCode
https://github.com/YangMengHeng/myCode/tree/master/Java
Gitee
https://gitee.com/QMBZ/myCode/tree/master/VSCode
https://gitee.com/QMBZ/myCode/tree/master/Java
旗木白哉のblog源代码下载地址：
https://github.com/YangMengHeng/myCode/tree/master/%E6%97%97%E6%9C%A8%E7%99%BD%E5%93%89%E3%81%AEblog
https://gitee.com/QMBZ/myCode/tree/master/%E6%97%97%E6%9C%A8%E7%99%BD%E5%93%89%E3%81%AEblog

---