Solution1:最长不下降子序列(LIS)
题目大意:
给出m种颜色作为喜欢的颜色(同时也给出顺序),然后给出一串长度为L的颜色序列,现在要去掉这个序列中的不喜欢的颜色,然后求剩下序列的一个子序列,使得这个子序列表示的颜色顺序符合自己喜欢的颜色的顺序,不一定要所有喜欢的颜色都出现。
分析:
因为喜欢的颜色是不重复的,把喜欢的颜色的序列依次存储到数组中,
pos[i] = j表示i颜色的下标为j。先在输入的时候剔除不在喜欢的序列中的元素,然后把剩余的保存在数组a中。按照最长不下降子序列的方式做,对于从前到后的每一个i,如果它前面的所有的j,一下子找到了一个j的下标pos[j]比pos[i]小,此时就更新dp[i]使它 = max(dp[i], dp[j] + 1);并且同时再每一次遍历完成一次j后更新maxl的值为长度的最大值,最后输出maxl。
dp[i]表示以第i元素为结尾的最长不下降子序列长度,那么对于每一个dp[i]而言,初始值即为1
普通做法代码:O(n2)
#include<iostream>
#include<math.h>
using namespace std;
int n,m,l;//n为总颜色数、m为喜欢的颜色数、l为布的长度
int pos[205],like[10005];//pos为喜欢的颜色序列,pos[x]=i表示喜欢的颜色x位于第i个位置
//like为布上剔除不喜欢的颜色后剩余的颜色所处的位置
int num=0;//num为布中喜欢的颜色数量
int dp[10005];
int maxl=0;
int main(){
cin>>n>>m;
int color;//颜色
for(int i=1;i<=m;i++){//读入
cin>>color;
pos[color]=i;
}
cin>>l;
for(int i=1;i<=l;i++){//剔除不喜欢的颜色
cin>>color;
if(pos[color]>=1){
like[num++]=pos[color];
}
}
for(int i=0;i<num;i++){//dp
dp[i]=1;
for(int j=0;j<i;j++){
if(like[i]>=like[j]){
dp[i]=max(dp[i],dp[j]+1);
}
}
maxl=max(dp[i],maxl);
}
cout<<maxl;
return 0;
}
优化做法代码:O(nlogn)
#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<cmath>
#include<queue>
#include<cstring>
#include<vector>
#include<map>
#define MAX 10005
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
int n,m,l;//n为总颜色数、m为喜欢的颜色数、l为布的长度
int pos[205],like[10005];//pos为喜欢的颜色序列,pos[x]=i表示喜欢的颜色x位于第i个位置
//like为布上剔除不喜欢的颜色后剩余的颜色所处的位置
int num=0;//num为布中喜欢的颜色数量
int d[10005];//d用来存储最长不下降子序列
int maxl=0;
int main(){
cin>>n>>m;
int color;//颜色
for(int i=1;i<=m;i++){//读入
cin>>color;
pos[color]=i;
}
cin>>l;
for(int i=1;i<=l;i++){//剔除不喜欢的颜色
cin>>color;
if(pos[color]>=1){
like[num++]=pos[color];
}
}
int len=1;
d[1]=like[0];
for(int i=1;i<num;i++){
if(like[i]>=d[len]){
d[++len]=like[i];
}else{
int p=upper_bound(d+1,d+1+len,like[i])-d;
d[p]=like[i];
}
}
cout<<len;
return 0;
}
Solution2:最长公共子序列(LCS)
经典LCS模型的两个序列的元素匹配必须是一一对应的,但是本题中允许公共部分产生重复元素,例如“ABBC”与“AABC”的最长公共子序列为“AABBC”,这就需要对原模型进行一些修改。
现在考虑一个问题:在计算dp[i][j]
时,当A[i]与B[i]不相等时,dp[i][j]
应该转移到哪些状态。很容易知道,由于可以产生重复元素,因此dp[i-1][j]
、dp[i][j-1]
都是可以对dp[i][j]
产生影响的(dp[i-1][j-1]
会通过dp[i-1][j]
、dp[i][j-1]
来影响dp[i][j]
,所以不需要考虑)。这样就可以知道,
(1)当A[i]与B[j]不相等时,dp[i][j]
应该取dp[i-1][j]
、dp[i][j-1]
中的较大值
(2)当A[i]与B[j]相等时,dp[i][j]
应该取dp[i-1][j]
、dp[i][j-1]
中的较大值+1(即在A[i]与B[j]不相等的结果上+1)
代码如下:
#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<cmath>
#include<queue>
#include<cstring>
#include<vector>
#include<map>
#define MAX 100005
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
const int maxc=210;//颜色的最大种类数
const int maxn=10010;//颜色序列的最大长度
int a[maxc],b[maxn],dp[maxc][maxn];
int main(){
int n,m;
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++){
scanf("%d",&a[i]);
}
int l;
scanf("%d",&l);//长度
for(int i=1;i<=l;i++){
scanf("%d",&b[i]);
}
//边界
for(int i=0;i<=m;i++){
dp[i][0]=0;
}
for(int i=0;i<=l;i++){
dp[0][i]=0;
}
//状态转移方程
for(int i=1;i<=m;i++){
for(int j=1;j<=l;j++){
dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
if(a[i]==b[j]){
dp[i][j]++;
}
}
}
printf("%d",dp[m][l]);
return 0;
}