LIS,LCS,LICS 学习笔记
1.最长上升子序列（LIS）
子序列： 1.可以不连续 2.相对位置不变
dp[i][j] 表示前i位置，最大值为j的LIS长度
1. dp[i-1][j] 前i-1位置，最大值为j的LIS长度 (没有考虑a[i])
2. dp[i][j]=dp[i-1][k]+1 (j==a[i] k < j)
ans=max(dp[n][i])
DP复杂度：状态数量*单个状态转移复杂度
O(n^2) 空间 O(n^2)
序列： 前i个位置，以第i个位置结尾。
f[i] 以第i个位置结尾的LIS长度
f[i] <- f[j]+1 (j< i a[j]< a[i])
ans=max(f[i])
O(n^2) 空间 O(n)
for(int i=1;i<=n;i++){
f[i]=1;
for(int j=1;j<i;j++){
if(a[j]<a[i])
f[i]=max(f[j]+1,f[i]);
}
}O(nlogn):
1. 用一个数组(栈)来维护最可能成为LIS的序列 （和DP没有关系）
2. 用树状数组来优化第二种DP（有推广意义）。
1 3 5 6 4 7 8
[1,3,5,6] 4
[1,3,4(5),6] 5
[1,3,4(5),5(6),6]
向前查找位置（二分或STL）nlogn
upper_bound: “元素值>查找值”的第一个元素的位置
lower_bound: “元素值>=查找值”的第一个元素的位置
#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
int a[40005];
int d[40005];
int main()
{
int n;
scanf("%d",&n);
for (int i=1;i<=n;i++) scanf("%d",&a[i]);
if (n==0)
{
printf("0\n");
return 0;
}
d[1]=a[1];
int len=1;
for (int i=2;i<=n;i++)
{
if (a[i]>d[len]) d[++len]=a[i];
else {
int j=lower_bound(d+1,d+len+1,a[i])-d; //找到第一个>=它的d的下标
d[j]=a[i];
}
}
printf("%d\n",len);
return 0;
}树状数组： 1. 求前缀和， 2.单点加减
int ask(int pos){
int ret=0;
while(pos>0){
ret+=c[pos];
pos-=lowbit(pos);
}
return ret;
}
void add(int pos,int w){
while(pos<=n){
c[pos]+=w;
pos+=lowbit(pos);
}
}树状数组 1. 求前缀最大值， 2.单点修改(往大里改)
int ask(int pos){
int ret=0;
while(pos>0){
ret=max(ret,c[pos]);
pos-=lowbit(pos);
}
return ret;
}
void modify(int pos,int w){
while(pos<=n){
c[pos]=max(c[pos],w);
pos+=lowbit(pos);
}
}for(int i=1;i<=n;i++){
f[i]=1+max(0,b[0],b[1],b[2],...,b[a[i]-1]);
b[a[i]]=f[i];
}for(int i=1;i<=n;i++){
f[i]=ask(a[i]-1)+1;
modify(a[i],f[i]);
}2.最长公共子序列
a 1 4 5 2 3
b 1 5 2 4 3
1 5 2 3
1) 前…个元素
f[i][j] a串前i个元素，b串前j个元素的LCS长度
a[i] != b[j] f[i][j] <- f[i-1][j] f[i][j-1]
a[i] == b[j] f[i][j] <- f[i-1][j-1]+1
O(1)*O(n^2)
f[n][m]
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
if(a[i] != b[j]) f[i][j]=max( f[i-1][j] , f[i][j-1]);
else f[i][j]=f[i-1][j-1]+1;
}
}
2) 以…结尾
f[i][j] a串以i结尾，b串以j结尾的LCS长度
a[i] != b[j] f[i][j] = 0
a[i] == b[j] f[i][j] <- f[k][l] (k< i l< j)
O(n^2)*O(n^2)
ans=max(f[i][j])
3.LICS(LCIS):
1) 前…个元素
f[i][j] a串前i个元素，b串前j个元素的LICS长度
无法转移
f[i][j][k] a串前i个元素，b串前j个元素的LICS长度最大值为k
a[i] != b[j] f[i][j][k] <- f[i-1][j][k] f[i][j-1][k]
a[i] == b[j] && a[i]==k f[i][j][k] <- f[i-1][j-1][l] l< k
O(n^3) 空间 (空间可以滚动数组优化) n^3 时间
ans=max(f[n][m][i])
2) 以…结尾
f[i][j] a串以i结尾，b串以j结尾的LICS长度
a[i] != b[j] f[i][j] = 0
a[i] == b[j] f[i][j] < - f[k][l]+1
(k < i l< j a[k]==b[l] a[k]< a[i] )
O(n^2)*O(n^2)
O(n^2)空间 O(n^4)时间
ans=max(f[i][j])
3)
f[i][j] a串前i个元素，b串以j结尾的LICS长度
a[i] != b[j] f[i][j] <- f[i-1][j]
a[i] == b[j] f[i][j] <- f[i-1][k] +1 (k< j b[k]< b[j])
for(int i=1;i< =n;i++){
for(int j=1;j< =m;j++){
if(a[i]!=b[j]) f[i][j]=f[i-1][j];
else{
f[i][j]=1;
for(int k=1;k< j;k++){
if(b[k]< b[j]) f[i][j]=max(f[i][j],f[i-1][k]+1);
}
}
}
}O(n^2)空间 O(n^3)时间
ans=max(f[n][i])
O(n^2logn)
int mx[];
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){ // LIS 树状数组优化
mx[j] = ask(b[j]-1);
modify(b[j],f[i-1][j]);
}
for(int j=1;j<=m;j++){// LCS
if(a[i]!=b[j]) f[i][j]=f[i-1][j];
else{
f[i][j]=mx[j]+1;
}
}
}O(n^2)? 思考
//By Menteur_Hxy
#include<cstdio>
#include<iostream>
using namespace std;
const int MAX=3010;
int n,m,top;
int a[MAX],b[MAX],f[MAX][MAX];
int main() {
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
// scanf("%d",&m);
for(int i=1;i<=n;i++) scanf("%d",&b[i]);
for(int i=1;i<=n;i++){
int maxn=0;
for(int j=1;j<=n;j++) {
if(a[i]!=b[j]) f[i][j]=f[i-1][j];
else f[i][j]=maxn+1;
if(a[i]>b[j]) maxn=max(maxn,f[i-1][j]);
}
}
int ans=0;
for(int i=1;i<=n;i++)
ans=max(ans,f[n][i]);
printf("%d",ans);
return 0;
}输出方案
f[i] <- max( f[j]+1) j< i
g[i] j
f[n] g[n] f[g[n]] g[f[g[n]]]
f[i][j]=max(f[i+1][j],f[i+1][j+1])+a[i][j];
int g[][]// 记录转移
if(f[i+1][j]>f[i+1][j+1]){
f[i][j]=f[i+1][j]+a[i][j];
g[i][j]=j;
}else{
f[i][j]=f[i+1][j+1]+a[i][j];
g[i][j]=j+1;
}