最长回文子串
原创
©著作权归作者所有:来自51CTO博客作者软糖酱八号机的原创作品,请联系作者获取转载授权,否则将追究法律责任
最长回文子串问题是求一个序列中最长的回文子串,这个子串可以不连续,如{1, 3, 4, 2, 1}的最长回文子串为{1, 2, 1},典型的区间dp
状态转移方程如下:
当i > j时,dp[i,j]= 0。
当i = j时,dp[i,j] = 1。
当i < j并且str[i] == str[j]时,dp[i][j] = dp[i+1][j-1]+2;
当i < j并且str[i] ≠ str[j]时,dp[i][j] = max(dp[i][j-1],dp[i+1][j]);
由于dp[i][j]依赖i+1,所以循环计算的时候,第一维必须倒过来计算,从len-1到0。
最后,s的最长回文子序列长度为dp[0][len-1]。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 1000 + 10;
char str[maxn];
int dp[maxn][maxn];
int main()
{
while (cin >> str) {
memset(dp, 0, sizeof(dp));
int len = strlen(str);
for (int i = len - 1; i >= 0; i--) {
dp[i][i] = 1;
for (int j = i + 1; j < len; j++) {
if (str[i] == str[j]) dp[i][j] = dp[i + 1][j - 1] + 2;
else dp[i][j] = max(dp[i][j - 1], dp[i + 1][j]);
}
}
printf("%d\n", dp[0][len - 1]);
}
return 0;
}
下面来将空间复杂度优化到O(n):
我们发现计算第i行时只用到了第i + 1行,因此我们只需要两行即可,即定义dp【2】【maxn】
起初先在第0行计算dp【0】【len - 1】,然后用第0行的结果计算dp【1】【len - 2】,再用第1行的结果在第0行计算dp【0】【len - 3】,以此类推
设正在计算的哪一行为cur,那么计算第cur行时,就要用到第1-cur行的结果
当计算完成时,如果len是奇数,则结果在第0行,如果是偶数,则结果在第1行
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 1000 + 10;
char str[maxn];
int dp[2][maxn];
int main()
{
while (cin >> str) {
memset(dp, 0, sizeof(dp));
int len = strlen(str), cur = 0;
for (int i = len - 1; i >= 0; i--) {
dp[cur][i] = 1;
for (int j = i + 1; j < len; j++) {
if (str[i] == str[j]) dp[cur][j] = dp[1 - cur][j - 1] + 2;
else dp[cur][j] = max(dp[cur][j - 1], dp[1 - cur][j]);
}
}
if (len % 2) printf("%d\n", dp[0][len - 1]);
else printf("%d\n", dp[1][len - 1]);
}
return 0;
}