线性动态规划:最长上升子序列
最长上升子序列
给定一个长度为 的数列,求数值严格单调递增的子序列的长度最长是多少。(子序列不一定要连续)
例如:
数列 3 1 2 1 8 5 6 中的最长上升子序列为1 2 5 6,长度为 4 。
算法 1:动态规划
数列中的第i项为a[i],对应f[i]。
-
状态表示:
f[i]表示从第一个数开始,以a[i]结尾的上升子序列最大长度。 -
状态计算:先将
f[i]初始化为 , 循环从 到 ,if (a[i] > a[j]) f[i] = max(f[i], f[j] + 1) -
分析:对于每个数,它本身就是一个以它结尾的子序列,所以如果在前面的数里面没有以它结尾的上升子序列,那
f[i]就是1。如果有,则把它加上,再判断是不是最长的,所以取max。
for (int i=0;i<n;i++)
{
f[i] = 1;
for (int j=0;j<i;j++)
{
if (a[i] > a[j]) f[i] = max(f[i], f[j]+1);
}
ans = max(ans, f[i]);
}算法2:动态规划+二分(贪心)
数列中的第i项为a[i],对应f[i]。
-
状态表示:
f[i]表示长度为i的最长上升子序列,末尾的最小数字。 -
状态计算:遍历a数组,在
f数组中查找小于a[i]的最大的数,将a[i]插到这个数的后面一个位置。
- 分析:容易看出,f数组是单调递增的,所以可以用二分查找。如果cnt表示f数组的长度(即答案),二分的范围是 0 ~ cnt。这里注意,每次二分完r+1的值,就是a[i]要插入的位置,也是cnt可能要更新的值,如果r+1大于cnt就更新cnt为r+1。
for (int i=0;i<n;i++)
{
int l = 0, r = cnt;
while (l < r)
{
int mid = l + r + 1>> 1;
if (f[mid] < a[i]) l = mid;
else r = mid - 1;
}
cnt = max(cnt, r + 1);
f[r+1] = a[i];
}例题:
最长公共子序列
给定两个长度分别为 和 的字符串 和 ,求既是 的子序列又是 的子序列的字符串长度最长是多少。
例如:
acbd 和 abedc 的最长公共子序列是 abd 长度为 。
动态规划
dp[i][j] 表示 的前 个字符和 的前 个字符中,最长公共子序列的长度。
状态转移:
- 如果
a[i] == b[j],dp[i][j] = dp[i-1][j-1] + 1,即把不包含这对字符的最长公共子序列加上这一对。 - 如果
a[i] != b[j]则dp[i][j]由dp[i-1][j]或者dp[i][j-1]转移来,可以理解为这两个序列中一定有一个的最后一个字符不在最长公共子序列中,所以对分别去掉最后一个字符的状态取两者的最大值作为当前的状态。
for (int i=1;i<=n;i++)
{
for (int j=1;j<=m;j++)
{
if (a[i] == b[j]) dp[i][j] = dp[i-1][j-1] + 1;
else dp[i][j] = max(dp[i-1][j], dp[i][j-1]);
}
}例题:
最短编辑距离
给定两个长度分别为 和 的字符串 和 ,现在要将 经过若干操作变为 ,可进行的操作有:
- 删除 – 将字符串 中的某个字符删除。
- 插入 – 在字符串 的某个位置插入某个字符。
- 替换 – 将字符串 中的某个字符替换为另一个字符。
现在请你求出,将 变为 至少需要进行多少次操作。
动态规划
dp[i][j] 表示将 的前 个字符变成 的前 个字符所用的最少的操作数。
状态转移: 对三个操作进行分析:
- 删除操作:把 删除之后 的前 个字符与 的前 个字符匹配,所以在这之前要满足 的前 个字符与 的前 个字符已经匹配,所以转移的是
dp[i-1][j] + 1。 - 插入操作:添加一个字母之后 的前 个和 的前 个变得相同,说明没有添加前 的前 个已经和 的前 个已经相同,即
dp[i][j-1] + 1。 - 替换操作:把 改成 之后 的前 个与 的前 个相同,所以之前 的前 个与 的前 个已经匹配,此时再判断 是否等于 即
dp[i-1][j-1] + (a[i] != b[j])。
所以 dp[i][j] 是从上面三个状态转移过来的,因为要求最短的操作数,所以上面三个值取最小值作为 dp[i][j] 的值。
初始化:
dp[0][i] = i如果a初始长度就是0,那么只能用插入操作让它变成b
dp[i][0] = i同样地,如果b的长度是0,那么a只能用删除操作让它变成b
for (int i=0;i<=n;i++) dp[i][0] = i;
for (int i=0;i<=m;i++) dp[0][i] = i;
for (int i=1;i<=n;i++)
{
for (int j=1;j<=m;j++)
{
dp[i][j] = min(dp[i-1][j] + 1, dp[i][j-1] + 1);
dp[i][j] = min(dp[i][j], dp[i-1][j-1] + (a[i] != b[j]));
}
}例题: