动态规划(一)
今天第一次接触动态规划,写一点自己的心得。准确来说不是第一次,之前看到过但是因为看不懂直接跳了😂,现在有点感觉主要是对搜索有了更加深刻的理解。
从背包问题引入
有N件物品和一个容量为V的背包。第i件物品的重量是w[i],价值是v[i]。
求解将哪些物品装入背包可使这些物品的重量总和不超过背包容量,且价值总和最大。
#include <stdio.h>
#include <cmath>
using namespace std;
int n=4;
int w[4]={2,1,3,2},v[4]={3,2,4,2};
int W=5;
int rec (int i,int j)//前一个表示物体的编号,后一个表示重量余额
{
int res;
if (i==n)
{
res=0;
}else if (w[i]>j)
{
res=rec(i+1,j);//如果下一个物体的质量溢出那么跳过。
}else
{
res=max(rec(i+1,j),rec(i+1,j-w[i])+v[i]);//状态转移方程,个人感觉这只进一步理解递归的重要一步
}
return res;
}
int main()
{
printf("%d",rec(0,W));//求出挑选方案中总和的最大值
getchar();
}
观察算法我们发现,状态转移方程分叉为两条路,搜索层数是n,复杂度(在2^n),接下来用记忆化搜索提高效率。
#include <stdio.h>
#include <cmath>
#include <string.h>
using namespace std;
int n=4;
int w[4]={2,1,3,2},v[4]={3,2,4,2};
int W=5;
**int dp[100][100];//足够大就好,用于记录已经搜索的内容**
int rec (int i,int j)
{
**if (dp[i][j]>0)
{
return dp[i][j];
}**
int res;
if (i==n)
{
res=0;
}else if (w[i]>j)
{
res=rec(i+1,j);
}else
{
res=max(rec(i+1,j),rec(i+1,j-w[i])+v[i]);
}
**dp[i][j]=res;**//不能忘了返回的时候记录下来
return res;
}
int main()
{
memset(dp,-1,sizeof(dp));
printf("%d",rec(0,W));
getchar();
}
接下来我们再研究那个记忆化数组,就正式步入动态规划了。
#include <stdio.h>
#include <cmath>
#include <string.h>
using namespace std;
int n=4;
int w[4]={2,1,3,2},v[4]={3,2,4,2};
int W=5;
int dp[100][100];//足够大就好
int main()
{
for (int i = n-1; i >=0; i--)
{
for (int j = 0; j <= W; j++)
{
if (j<w[i])
{
dp[i][j]=dp[i+1][j];
}else
{
dp[i][j]=max(dp[i+1][j],dp[i+1][j-w[i]]+v[i]);
}
}
}
printf("%d",dp[0][W]);
getchar();
}
可以看到,这里数组变成了主要对象而不是辅助了,关于这种方法我还没有研究透彻,待我再研究一下届时我还会再发博客分享。
这里几个问题:
1.这样的方法看上去怎么感觉怎加了复杂度,但其实和前一种方法是一样的。
2.这样写优势在哪里