动态规划(DP)

本文通过洛谷和LeetCode上的典型问题,展示了动态规划在解决01背包问题、最短路径、最大面积、斐波那契数列等问题中的应用。代码实例涵盖从基本到高级的动态规划技巧,如一维和二维数组优化、分组背包问题等。
摘要由CSDN通过智能技术生成

洛谷 P.1359

输入说明 :

3

5 (1->2) 15 (1->3)

7 (2 -> 3)

#include<iostream>
#include<vector>
#include<climits>
using namespace std;

vector<int> dp(201);
int  main(){
	int n ;
	cin>>n;
	vector<vector<int>> ori(201,vector<int>(201));
	for(int i = 1; i < n;i++ ){
		for(int j = i + 1; j <= n;j++){
			cin>>ori[i][j];
		}
		dp[i] = INT_MAX;
	}
	for(int i = n - 1;i > 0;i--){
		for(int j = i + 1; j <= n;j++  ){
			dp[i] = min(dp[i] , dp[j] + ori[i][j]);
		}
	}
	cout<<dp[1];
	return 0;
}

状态转移关系 : 我们要用i把n上流的中转站从大到小跑一遍。我们先记录中转站2到中转站3的最小价钱,我们要用jj跑一遍中转站2下流的所有中转站,记录ori[i][j]+dp[j]的最小价钱,记录到dp[i]里面。

 

细节  : 因为dp[i}是每次min来的,但dp各项初始化为0,因此我们要让除了最后一个中转站外设为INT_MAX(用头文件#include<climits>)

洛谷 P.1060

 

 01背包问题!

无优化代码

for(int i=1;i<=n;i++)
{
    for(int c=0;c<=m;c++)
    {
        f[i][c]=f[i-1][c];
        if(c>=w[i])
        f[i][c]=max(f[i][c],f[i-1][c-w[i]]+v[i]);
    }
}

 一维数组优化 :

for(int i=1;i<=n;i++)
{
    for(int c=m;c>=0;c--)
    {
        if(c>=w[i])
        f[c]=max(f[c],f[c-w[i]]+v[i]);
    }
}
#include<iostream>
using namespace std;
int w[30],v[30],dp[30][100000];
int main(){
	int m, n;//m为总价,n为数量
	cin>>m>>n;
	 for(int i = 1;i <=n;i++){
	 	cin>>v[i]>>w[i];
	 	w[i] *= v[i];//w[i]变为总收获 
	 }
	 for(int i = 1;i <= n;i++ ){
	 	for(int j = 0;j <= m;j++){
	 		dp[i][j] = dp[i - 1][j];
	 		if(j >= v[i]){
	 			dp[i][j] = max(dp[i][j] , dp[i - 1][j - v[i]] + w[i]);
			 }
		 }
	 }
	 cout<<dp[n][m];
	 return 0;
}
//优化版
#include<iostream>
using namespace std;
int w[30],v[30],dp[100000];
int main(){
	int m, n;//m为总价,n为数量
	cin>>m>>n;
	 for(int i = 1;i <= n;i++){
	 	cin>>v[i]>>w[i];
	 	w[i] *= v[i];//w[i]变为总收获 
	 }
	 for(int i = 1;i <= n;i++ ){
	 	for(int j = m;j >= v[i];j--){
	 			dp[j] = max(dp[j] , dp[j - v[i]] + w[i]);
			 
		 }
	 }
	 cout<<dp[m];
	 return 0;
}

洛谷P.1802

与上题相似 直接上代码‘

#include <iostream>
#include <cstdio>
#include <cstdlib>
using namespace std;
int lose[1001],win[1001],use[1001];
long long f[1001][1001];
int main(){
    int n,x;
    cin>>n>>x;
    for (int i=1;i<=n;i++)
      cin>>lose[i]>>win[i]>>use[i];
    for (int i=1;i<=n;i++)
      for (int j=0;j<=x;j++)
        if (j>=use[i]) f[i][j]=max(f[i-1][j-use[i]]+win[i],f[i-1][j]+lose[i]);
 else f[i][j]=f[i-1][j]+lose[i];
    cout<<f[n][x]*5;
    return 0;
}

 

#include <cstdio>
#include <iostream>
using namespace std;
int dp[1100];
int win[1100],lose[1100],use[1100];
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
     scanf("%d%d%d",lose+i,win+i,use+i);
    for(int i=1;i<=n;i++)
     {
         for(int j=m;j>=use[i];j--)
          dp[j]=max(dp[j]+lose[i],dp[j-use[i]]+win[i]);
         for(int j=use[i]-1;j>=0;j--)
          dp[j]+=lose[i];//其余情况不是0而要加lose[i]
     }
     printf("%lld",5ll*dp[m]);
}

洛谷P.1049

 

看似与01背包不同,但我们可以看作它价值与容量相等(v[i] = w[i]),即与1060一致

状态转移关系 : dp(x)为在x的最大容积下(对应1060的最大价格)最多能装下多少容积(对应1060最多价值),最后将m - dp(m)即为最小剩余容量 。。

#include<iostream>
using namespace std;

int a[31],dp[20001];

int main(){
	int m , n;
	cin>>m>>n;
	for(int i = 1;i <= n;i++){
		cin>>a[i];
	}
	for(int i = 1;i <=n;i++){
		for(int j = m; j >= a[i];j--){
			dp[j] = max(dp[j] , dp[j -a[i]] + a[i]);  
		}
	}
	cout<<m - dp[m];
	return 0;
} 

 来题几何

LeetCode.221

 

 设dp(i,j)为以(i,j)为右下角的正方形边长的最大值

边界情况 :i = 0或 j = 0并且matrix[i][j]==1时dp[i][j]为1(无法向左上角延申)

动态转移关系 :dp[i][j] = min(min(dp[i - 1][j], dp[i][j - 1]), dp[i - 1][j - 1]) + 1;(取决于左边,坐上,上面)最小的dp(画图理解,这种题主要背)


class Solution {
public:
    int maximalSquare(vector<vector<char>>& matrix) {
        int maxSide = 0;
        int rows = matrix.size(), columns = matrix[0].size();
        vector<vector<int>> dp(rows, vector<int>(columns));
        for (int i = 0; i < rows; i++) {
            for (int j = 0; j < columns; j++) {
                if (matrix[i][j] == '1') {
                    if (i == 0 || j == 0) {
                        dp[i][j] = 1;
                    } else {
                        dp[i][j] = min(min(dp[i - 1][j], dp[i][j - 1]), dp[i - 1][j - 1]) + 1;
                    }
                    maxSide = max(maxSide, dp[i][j]);
                }
            }
        }
        int maxSquare = maxSide * maxSide;
        return maxSquare;
    }
};

LeetCode : 118

弄成二维数组看更直观 :

1   1

1    2   1

1    3   3   1

1    4   6    4   1

很容易看出边界情况 :dp[m][0] = dp[m][m] = 1,

状态转移关系 : dp[m][n] = dp[m - 1][n - 1] + dp[m - 1][n],&&n > 0&&n < m

class Solution {
public:
    vector<vector<int>> generate(int numRows) {
        vector<vector<int>> dp(numRows);
        for(int i = 0; i < numRows;i++ ){
            dp[i].resize(i + 1);//每行变成i + 1列
            dp[i][0] = 1;
            dp[i][i] = 1;
            for(int j = 1;j < i;j++){
                dp[i][j] = dp[i - 1][j - 1] + dp[i - 1][j];
            }
        }
        return dp;
    }
};

LeetCode : 119

 水题放松下,和118一致

class Solution {
public:
    vector<int> getRow(int rowIndex) {
        vector<vector<int>> dp(rowIndex + 1);
        for(int i = 0; i < rowIndex + 1;i++ ){
            dp[i].resize(i + 1);//每行变成i + 1列
            dp[i][0] = 1;
            dp[i][i] = 1;
            for(int j = 1;j < i;j++){
                dp[i][j] = dp[i - 1][j - 1] + dp[i - 1][j];
            }
        }
        return dp[rowIndex];
    }
};

洛谷P.1757

 重点 :他的物品大致可分为 k 组,每组中的物品相互冲突(分组背包问题)

思路 :不同于普通01背包,第一层for循环是以小组为元素,第二层对于一个dp[j]第三层for循环要遍历一组的所有元素,以便每组中物品不冲突

#include<iostream>
using namespace std;
int v[1001],w[1001],dp[1001],level[1001],g[1001][1001];
int main(){
	int m , n , x , t;
	cin>>m>>n;
	for(int i = 1; i <= n;i++){
		cin>>v[i]>>w[i]>>x;
		t = max(t , x);
		level[x]++;
		g[x][level[x]] = i;//存储每组第几个元素的值 
	}
	for(int i = 1;i <= t; i++){//对于每个小组 
		for(int j = m ;j >=0;j--){
			for(int k = 1;k <=level[i];k++){//一组的所有元素 
				if(j >= v[g[i][k]]){
					dp[j] = max(dp[j] , dp[j - v[g[i][k]]] + w[g[i][k]]); //不让组内元素互动
				}
			}
		}
	}
	cout<<dp[m];
	return 0;
}

LeetCode.509

class Solution {
public:
    int fib(int n) {
        if(n == 0) return 0;
        if(n == 1) return 1;//特殊情况
        int begin = 0, end = 1 , res ;
        for(int i =  2;i <= n;i++){
            res = begin + end;
            begin = end;
            end = res;//依次后移
        }
        return res;
    }
};

洛谷 P.1507

 太过简单,不做阐释

#include<iostream>
using namespace std;
int a[401],b[401],c[501];
int dp[401][401];

int main(){
	int v,m,n;
	cin>>v>>m>>n;
	for(int i = 1;i <= n;i++){
		cin>>a[i]>>b[i]>>c[i];
	} 
	for(int i = 1;i <= n;i++){
		for(int j = v;j >= a[i];j--){
			for(int k = m;k >= b[i];k--){
				dp[j][k] = max(dp[j][k] , dp[j - a[i]][k - b[i]] + c[i]);
			}
		}
	}
	cout<<dp[v][m];
	return 0;
}

洛谷 P.2858

 题目特色 : 入口有两个并且物品价值随着每次选择改变

用二维数组dp[left][right]表示从left 到right获得的最多钱

边界情况 : dp[i][i] = s[i] * m(只剩一份零食肯定是在最后一天)

#include<iostream>
using namespace std;
int s[2001],dp[2001][2001];

int main(){
	int m ;
	cin>>m;
	for(int i = 1; i <= m;i++ ){
		cin>>s[i];
	}
	for(int i = 1;i <= m;i++){
		dp[i][i] = m * s[i];
	}
	for(int i = 1;i <= m;i++){
		for(int left = 1;left <= m;left++){
			int right = left + i;
			if(right > m) break;
			dp[left][right] = max(dp[left + 1][right] + s[left] * (m - i),dp[left][right - 1] + s[right] * (m - i));
		}
	}
	cout<<dp[1][m];
	return 0;
}

洛谷 P .1091

 

 思路 : 找到一个数据点让它左边的递增元素数量与它右边的递减元素数量最多

#include<iostream>
using namespace std;
int a[105],b[105],c[105];

int main(){
	int n;
	cin>>n;
	for(int i = 1;i <= n;i++){
		cin>>a[i];
	}
	for(int i = 1;i <= n;i++){//细节是从第一个数开始而第二层超出第一个数以便让i = =1时候b[i] = 1 
		for(int j = 0;j < i;j++){
			if(a[i] > a[j])
			b[i] = max(b[i],b[j] + 1);//b[i]记录i左边有几个从低到高的数 (包括自己) 
		}
	}
	for(int i = n; i>= 0;i--){
		for(int j = n + 1;j > i;j--){
			if(a[i] > a[j] ){
				c[i] = max(c[i],c[j] + 1);//c[i]记录i右边有几个从高到低的数 (包括自己) 
			}
		}
	}
	int res = 0;
	for(int i = 1;i <= n;i++){
		res = max(res,b[i] + c[i] - 1);	
	} 
	cout<< n - res;
	return 0; 
}

洛谷 P.1077

 边界情况 :dp[0] = 1(摆0盆花有1种方案)

#include<iostream>
using namespace std;
const int mod = 1000007;
int a[105],dp[105];

int main(){
	int m,n;
	cin>>m>>n;
	for(int i = 1 ; i <= m;i++){
		cin>>a[i];
	}
	dp[0] = 1;
	for(int i = 1; i <= m;i++){//枚举每个品种 
		for(int j = n;j >= 1;j--){
			for(int k = 1; k <= a[i] && k <= j;k++){//相同品种的所有花 
				dp[j] = (dp[j] + dp[j - k]) % mod;
			}
		}
	}
	cout<<dp[n];
	return 0;
} 

PTA . 1020 - 月饼

2022.3.6

 特点 : 销售时允许取出一部分库存

意味着可以选择一点数据而非全选,解决方法 :结构体

//N 表示月饼的种类数
//D 表示市场最大需求量
//每种月饼的库存量 每种月饼的总售价
//贪心问题
#include<bits/stdc++.h>
using namespace std;
struct item{
	double cap;//容量 
	double price;//单价 
}moon[1001];

bool cmp(item a,item b){
	return a.price > b.price; 
}

int main(){
	int n,d;
	double res = 0.0;
	cin>>n>>d;
	for(int i = 1; i <= n;i++){
		cin>>moon[i].cap;
	}
	for(int i = 1;i <= n;i++){
		cin>>moon[i].price;
		moon[i].price /= moon[i].cap;//转化为单价
	}
	sort(moon + 1,moon + n + 1,cmp);
	for(int i = 1; i <= n;i++){
		if(d <= moon[i].cap){//供大于求
			res += moon[i].price * d;
			break;
		}
		else{//求大于供
			res += moon[i].price * moon[i].cap;
			d -= moon[i].cap;
		}
	}
	printf("%.2f",res);
	return 0;
} 

 

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值