今天给动态规划扫个尾,还有两题。
暴力解法
两层for循环,遍历区间起始位置和终止位置,然后还需要一层遍历判断这个区间是不是回文。所以时间复杂度:O(n^3)。
动态规划
动规五部曲:
确定dp数组(dp table)以及下标的含义
本题如果定义,dpi 为 下标i结尾的字符串有 dpi个回文串的话,会发现很难找到递归关系。
dpi 和 dpi-1 ,dpi + 1 看上去都没啥关系。
在判断字符串S是否是回文,那么如果知道 s1,s2,s3 这个子串是回文的,那么只需要比较 s0和s4这两个元素是否相同,如果相同的话,这个字符串s 就是回文串。
那么此时是不是能找到一种递归关系,也就是判断一个子字符串(字符串下标范围i,j)是否回文,依赖于,子字符串(下标范围i + 1, j - 1)) 是否是回文。
所以为了明确这种递归关系,dp数组要定义成一位二维dp数组。
布尔类型的dpij:表示区间范围i,j (注意是左闭右闭)的子串是否是回文子串,如果是dpij为true,否则为false。
确定递推公式
在确定递推公式时,就要分析如下几种情况。
整体上是两种,就是si与sj相等,si与sj不相等这两种。
当si与sj不相等,那没啥好说的了,dpij一定是false。
当si与sj相等时,这就复杂一些了,有如下三种情况
- 情况一:下标i 与 j相同,同一个字符例如a,当然是回文子串
- 情况二:下标i 与 j相差为1,例如aa,也是回文子串
- 情况三:下标:i 与 j相差大于1的时候,例如cabac,此时si与sj已经相同了,看i到j区间是不是回文子串就看aba是不是回文就可以了,那么aba的区间就是 i+1 与 j-1区间,这个区间是不是回文就看dpi + 1j - 1是否为true。
以上三种情况分析完了,那么递归公式如下:
cpp
if (s[i] == s[j]) {
if (j - i <= 1) { // 情况一 和 情况二
result++;
dp[i][j] = true;
} else if (dp[i + 1][j - 1]) { // 情况三
result++;
dp[i][j] = true;
}
}
result就是统计回文子串的数量。
注意这里没有列出当si与sj不相等的时候,因为在下面dpij初始化的时候,就初始为false。
dp数组如何初始化
dpij可以初始化为true么? 当然不行,怎能刚开始就全都匹配上了。
所以dpij初始化为false。
确定遍历顺序
首先从递推公式中可以看出,情况三是根据dpi + 1j - 1是否为true,在对dpij进行赋值true的。
dpi + 1j - 1 在 dpij的左下角。
如果这矩阵是从上到下,从左到右遍历,那么会用到没有计算过的dpi + 1j - 1,也就是根据不确定是不是回文的区间i+1,j-1,来判断了i,j是不是回文,那结果一定是不对的。
所以一定要从下到上,从左到右遍历,这样保证dpi + 1j - 1都是经过计算的。
C++代码如下:
cpp
class Solution {
public:
int countSubstrings(string s) {
vector<vector<bool>> dp(s.size(), vector<bool>(s.size(), false));
int result = 0;
for (int i = s.size() - 1; i >= 0; i--) { // 注意遍历顺序
for (int j = i; j < s.size(); j++) {
if (s[i] == s[j]) {
if (j - i <= 1) { // 情况一 和 情况二
result++;
dp[i][j] = true;
} else if (dp[i + 1][j - 1]) { // 情况三
result++;
dp[i][j] = true;
}
}
}
}
return result;
}
};
然后看第二道:516. 最长回文子序列 - 力扣(LeetCode)
回文子串是要连续的,回文子序列可不是连续的!
动规五部曲分析如下:
确定dp数组(dp table)以及下标的含义
dpij:字符串s在i, j范围内最长的回文子序列的长度为dpij。
确定递推公式
在判断回文子串的题目中,关键逻辑就是看si与sj是否相同。
如果si与sj相同,那么dpij = dpi + 1j - 1 + 2;
如果si与sj不相同,说明si和sj的同时加入 并不能增加i,j区间回文子序列的长度,那么分别加入si、sj看看哪一个可以组成最长的回文子序列。
加入sj的回文子序列长度为dpi + 1j。
加入si的回文子序列长度为dpij - 1。
那么dpij一定是取最大的,即:dpij = max(dpi + 1j, dpij - 1);
代码如下:
cpp
if (s[i] == s[j]) {
dp[i][j] = dp[i + 1][j - 1] + 2;
} else {
dp[i][j] = max(dp[i + 1][j], dp[i][j - 1]);
}
dp数组如何初始化
首先要考虑当i 和j 相同的情况,从递推公式:dpij = dpi + 1j - 1 + 2; 可以看出 递推公式是计算不到 i 和j相同时候的情况。
所以需要手动初始化一下,当i与j相同,那么dpij一定是等于1的,即:一个字符的回文子序列长度就是1。
其他情况dpij初始为0就行,这样递推公式:dpij = max(dpi + 1j, dpij - 1); 中dpij才不会被初始值覆盖。
cpp
vector<vector<int>> dp(s.size(), vector<int>(s.size(), 0));
for (int i = 0; i < s.size(); i++) dp[i][i] = 1;
确定遍历顺序
从递归公式中,可以看出,dpij 依赖于 dpi + 1j - 1 ,dpi + 1j 和 dpij - 1,
所以遍历i的时候一定要从下到上遍历,这样才能保证下一行的数据是经过计算的。
j的话,可以正常从左向右遍历。
C++代码如下:
cpp
class Solution {
public:
int longestPalindromeSubseq(string s) {
vector<vector<int>> dp(s.size(), vector<int>(s.size(), 0));
for (int i = 0; i < s.size(); i++) dp[i][i] = 1;
for (int i = s.size() - 1; i >= 0; i--) {
for (int j = i + 1; j < s.size(); j++) {
if (s[i] == s[j]) {
dp[i][j] = dp[i + 1][j - 1] + 2;
} else {
dp[i][j] = max(dp[i + 1][j], dp[i][j - 1]);
}
}
}
return dp[0][s.size() - 1];
}
};