B3637 最长上升子序列
思路:
1.使用动态规划,b[i]表示以a[i]结尾的最长上升子序列
2.遍历a[i](b[i]初始定义为1)
3.对于每个i再从头遍历到i,若有a[j]比a[i]小的,代表可以在b[j]最长长度后面再加1即b[j]+1,每次与b[i]本身比较达到取最大的目的
4.在外面一层的循环,对每一个b[i]比较取最大赋值给m
代码:
#include<bits/stdc++.h>
using namespace std;
int main(){
int n,m=0;
cin>>n;
int a[5005];
int b[5005];
for(int i=0;i<n;i++){
cin>>a[i];
}
for(int i=0;i<n;i++){
b[i]=1;
for(int j=0;j<=i-1;j++){
if(a[j]<a[i]){
b[i]=max(b[i],b[j]+1);
}
}
if(b[i]>m){
m=b[i];
}
}
cout<<m;
}
P1115 最大子段和
思路:
1.使用动态规划,dp[i]表示以a[i]结尾的最大子段和
2.dp[0]初始化为a[0]
3.从i=1开始遍历,判断是保留本身大还是加上前面最大子段和大
4.比较每个dp[i]大小,最大的赋值m
ps:因为数据有负数,m需要赋值-1e9保证被换掉
代码:
#include<bits/stdc++.h>
using namespace std;
int main(){
int n,a[200005],dp[200005],m=-1e9;
cin>>n;
for(int i=0;i<n;i++){
cin>>a[i];
}
dp[0]=a[0];
for(int i=1;i<n;i++){
dp[i]=max(dp[i-1]+a[i],a[i]);
if(dp[i]>m){
m=dp[i];
}
}
cout<<m;
}
P8707 [蓝桥杯 2020 省 AB1] 走方格
思路:
1.使用动态规划,dp[i][j]表示在走到这个方格的可能方案个数
2.初始化,dp[1][1]为0,因为只能向右或向下,dp[1][j]、dp[i][1]为1,因为不能走行号和列数都是偶数的,所以这些方格赋值为0
3.遍历所有方格,dp[i][j]即其上方和左边方格的和
代码:
#include<bits/stdc++.h>
using namespace std;
int main(){
int n,m;
cin>>n>>m;
int dp[35][35];
dp[1][1]=0;
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
if(i%2==0&&j%2==0){
dp[i][j]=0;
}
else if(i>1&&j>1){
dp[i][j]=dp[i-1][j]+dp[i][j-1];
}
else{
dp[i][j]=1;
}
}
}
cout<<dp[n][m];
}
P1216 [USACO1.5] [IOI1994]数字三角形 Number Triangles
思路:
1.用二维数组存储数字三角形
2.为了消除不确定性,从倒数第二层向上遍历
3.逐层逐个判断是与左下相加大还是与右下相加大,并存储
4.最后最上面a[1][1]即为最大值
代码:
#include<bits/stdc++.h>
using namespace std;
int main(){
int r;
cin>>r;
int a[1001][1001];
for(int i=1;i<=r;i++){
for(int j=1;j<=i;j++){
cin>>a[i][j];
}
}
for(int i=r-1;i>=1;i--){
for(int j=1;j<=i;j++){
a[i][j]=max(a[i+1][j]+a[i][j],a[i+1][j+1]+a[i][j]);
}
}
cout<<a[1][1];
return 0;
}
P1020 [NOIP1999 普及组] 导弹拦截
思路:
1.第一问即是求导弹高度数列的最长不上升序列
2.第二问根据Dilworth 定理即可转换为求最长上升子序列
3.两个可以用直dp解,但要求时间复杂度为O(nlogn),即要用二分优化
4.最长不上升序列个数等价为:若遍历的a[i]小于等于已有序列最后一位数则序列+1,反之找到与a[i]最接近的比其小的序列里的数进行交换,最后得到的序列长度。因为序列是单调的所以可以用二分查找优化(利用upper_bound()这个函数查找
5.最长上升序列同理(二分因为是找比其大的,所以用lower_bound()这个函数
代码:
#include<bits/stdc++.h>
using namespace std;
int a[100005],dp1[100005],dp2[100005],m,n,l,p;
int main(){
n=m=1;
l=0;
while(scanf("%d",&a[++l])!=EOF);
l--;
dp1[m]=dp2[n]=a[1];
for(int i=2;i<=l;i++){
if(dp1[m]>=a[i]){
dp1[++m]=a[i];
}
else{
p=upper_bound(dp1+1,dp1+m+1,a[i],greater<int>())-dp1;
dp1[p]=a[i];
}
if(dp2[n]<a[i]){
dp2[++n]=a[i];
}
else{
p=lower_bound(dp2+1,dp2+n+1,a[i])-dp2;
dp2[p]=a[i];
}
}
cout<<m<<endl;
cout<<n;
return 0;
}