LIS:
最长公共子序列,定义不再过多赘述,主要说一下它的状态方程,dp[i][j]代表的是x[1....i],和y[1.....j]的最长公共子序列,模板就是对x和y进行遍历,当想x[i]==y[j]时,dp[i][j]=dp[i-1][j-1]+1,这个很好理解,因为dp里面还没有包括x[i]==y[j]的情况,所以要加1,当x[i]和y[i]不相等的时候 ,dp[i][j]=max(dp[i][j-1],dp[i-1][j]),因为dp[i][j-1]和dp[i-1][j]的值必定大于dp【i-1】【j-1】的值,所以不需要考虑这个,也可以这样理解 x[1.....i]和y[1.....j-1]的最长公共子序列不会小与x[1...i-1]与y[1...j-1]的最长公共子序列,另一个也是一样
模板如下:
for(int i = 1; i <= A_len; ++i) {
for(int j = 1; j <= B_len; ++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]);
}
}
LCS(二分优化)
对于相同长度的子序列,尽量使用较小的那个,因为这样可以构造更长的序列
int cnt=1; //dp初始化为INF
for(int i = 0; i < length; ++i) {
int res = lower_bound(dp + 1, dp + cnt, a[i]) - dp;//返回第一个大于等于a[i]的地址
if(a[i] > dp[cnt]) dp[++cnt] = a[i];
else dp[res] = a[i];//代替掉较大元素
}
printf("%d\n", cnt);