LCS
- 最长公共子序列问题,十分经典,最简单的状态转移方程(板子)为:
if(x[i - 1] == y[j - 1])
lcs[i][j] = lcs[i - 1][j - 1] + 1
else
lcs[i][j] = max(lcs[i - 1][j],lcs[i][j - 1]);
- 观察上面的式子,可以看到对于lcs[i][jx],只与lcs[i - 1][jy]有关,其中jy = jx - 1 or jx,于是可以将这n * m的二维数组改成两个一维数组,一个存储目前的最大长度lcs[maxn],一个存储上一次的最大长度last[maxn]:
if(x[i - 1] == y[j - 1])
lcs[j] = last[j - 1] + 1;
else
lcs[j] = max(last[j],lcs[j - 1]);
同时在内层循环结束后,要将lcs copy到last中:
for(int j = 1;j <= m;j++){
last[j] = lcs[j];
}
- 再优化,注意到,lcs数组与last数组实际上在每次内层循环结束时最后一个元素是一样的,而计算过程中就只用到最后一个元素,于是只用记录那最后一个元素就好了,如下:
if(x[i - 1] == y[j - 1])
lcs[j] = pre + 1;
else
lcs[j] = max(pre,lcs[j]);
同时在外层循环开始要初始化pre = 0,内层循环结束要给pre赋值lcs[j]。
LIS
记录一种上课时听到的精妙思路(但空间复杂度较高),注意到,一个单调递增的序列的任意子序列也是单调递增的,于是问题可以进行转化:设输入的序列为X,将其副本进行排序(O(nlogn))得到序列Y,对X和Y求LCS(O(n2)),即得到答案。
void LIS(){
my_copy();//将x拷贝到y
sort(y,y + strlen(y),less<int>());//给y进行题目所要求的排序
LCS();//LCS问题
}