最长不下降子序列是动态规划的必做入门题
普通做法如下:
#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;
int a[1001][3]; //a[i][1]用来记录输入的值,a[i][2]用来记录从最后一个到编号为i的数最长不下降子序列的长度,a[i][3]用来记录上一个数,便于打印;
int main(){
int i,j,k,m,n;
cin>>m;
for(i=1;i<=m;i++){
cin>>a[i][1];
a[i][2]=1;//初始化长度为1,包括他本身;
a[i][3]=0; //0为没有数
}
for(i=m-1;i>=1;i--){
k=0;//k记录从后往前枚举到i最长长度;
n=0;//n记录上一个数的编号;
for(j=i+1;j<=m;j++){
if((a[i][1]<a[j][1])&&(a[j][2]>k)){//若当前值比上一个值小且比当前最长长度小;
k=a[j][2];//更新
n=j;//记录
}
}
if(k>0){
a[i][2]=k+1;//加上当前值的长度1
a[i][3]=n;
}
}
n=1;//从后往前搜索看长度最长的是多少
for(j=2;j<=m;j++){
if(a[j][2]>a[n][2])n=j;
}
cout<<a[n][2]<<endl;
return 0;
}
若要打印序列,则还要加上如下几行代码:
while(n!=0){
cout<<a[n][1]<<' ';
n=a[n][3];
}
接下来介绍O(n)的二分算法
这个算法的弊端是只能算出不下降子序列的长度而不能打印序列
这个算法思想就是把输入的数放入一个栈里,
从第一个数开始,一个一个和后面的比较;
若后面的数比栈内队尾的数大,添加此元素进栈,队尾加一
若位于中间位置,则把栈内刚好比他大的数删除,并添加至此元素被删的位置,因为在序列长度不会变长的情况下,小的数更有前途
记住比较时是二分比较
最后输出栈长即可
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int a[100010],b[100010],top;
int dp(int j){
int f=0,l=top;
while(f<l){
int mid=(f+l)/2;
if(b[mid]>a[j])l=mid;
else f=mid+1;
}
return f;
}
int main(){
int i,j,k,m,n,p;
scanf("%d",&p);
for(k=1;k<=p;k++){
top=0;
memset(b,0,sizeof(b));
memset(a,0,sizeof(a));
scanf("%d",&m);
for(i=1;i<=m;i++)
scanf("%d",&a[i]);
b[++top]=a[1];
for(j=2;j<=m;j++){
if(a[j]>b[top])b[++top]=a[j];
else{
n=dp(j);
b[n]=a[j];
}
}
printf("%d\n",top);
}
return 0;
}
因篇幅问题不能全部显示,请点此查看更多更全内容