一、dp是什么
我所理解的dp,就是要求多阶段的问题解(或可以转换为多阶段的问题最优解),求出每个阶段的解,最终推出答案。既将大问题转换为一个个小问题,并将这些小问题的答案存储下来,在解决后面的小问题或者最终的大问题时可以用到。
二、dp和贪心的区别
刚开始接触dp的时候,感觉和贪心区别不太大,两者的写法都基本相似(实际为贪心一般可用dp实现)。于是去问了师兄,师兄是这么解答的:
贪心很多时候每一步是不相互影响的,dp是很多时候每一步之间会相互影响。对于会相互影响的,可能在dp的时候就要用到贪心的思想,记录前面步骤最优的结果。然后贪心有时候也要用到dp来记录前面的结果。
我的理解是贪心着重于当前最优解,既每一个的最优;而dp则会看前面几步的答案。
三、dp的类型
本菜菜没有做太多题,但队里有巨佬总结了,这里引用一下
1.线性dp
2.背包问题
3.区间dp
4.树形dp
5.状压dp
6.换根dp
7.概论dp
8.期望dp
9.数位dp
10.虚树dp
11.斜率优化
12.单调队列优化
13.插头dp
14.凸优化dp
15.四边形不等式
四、 解决dp问题的步骤
1.状态的设置
我认为解决dp问题的重点是状态的设置,既如何构造状态去表示所有的情况;有师兄在讲课的时候说过,其实dp就是将所有暴力的情况,用很少的状态去表示出来。我们拿到dp题的时候,所要想的就是怎么在合理的时间复杂度之类,去构造状态。
构造状态不一定是直接构造出答案,既并不是说dp[n]/dp[n][n]就是最终的答案;例如专题练习中就有一题就是,求出dp[i],然后根据每个dp[i]去算出最后的答案。
2.状态的初始化/边界化
在dp问题中,我们也要注意状态的初始化,我就曾试过因为状态初始化没做好,导致一直wa
3.状态的转移
我们构造出状态,并将状态初始化之后,就要想如何实现状态直接的转换。一般而言,如果你对状态的设置比较清楚,状态的转移也不成问题。
五、本专题题解
1.送快弟
思路:01背包板子题
#include <iostream>
#include <algorithm>
using namespace std;
int main()
{
int t;
cin>>t;
while (t--)
{
int dp[1010]={},v[1010],w[1010];
int n,k,i,j;
cin>>n>>k;
for (i=0;i<n;i++)
cin>>v[i];
for (i=0;i<n;i++)
cin>>w[i];
for (i=0;i<n;i++)
for (j=k;j>=w[i];j--)
dp[j]=max(dp[j],dp[j-w[i]]+v[i]);
cout<<dp[k]<<endl;
}
return 0;
}
2.CD
题意:01背包,并输出选择了哪些
思路:01背包变型,另加一个二维数组,记录在第i件物品,在价值j的时候是否使用了。最后遍历一遍,输出所用的物品。
#include <iostream>
#include <cstdio>
using namespace std;
int main()
{
int n,m;
while (scanf("%d%d",&n,&m)!=EOF)
{
int i,j;
int v[10010],dp[10010]={};
for (i=0; i<m; i++)
cin>>v[i];
int a[25][10010]={};
for (i=0; i<n; i++)
for (j=n; j>=v[i]; j--)
{
if (dp[j-v[i]]+v[i]>=dp[j])
{
dp[j]=dp[j-v[i]]+v[i];
a[i][j]=1;
}
}
for (i=m-1,j=n;i>=0;i--)
{
if (a[i][j])
{
cout<<v[i]<<' ';
j-=v[i];
}
}
cout<<"sum:"<<dp[n]<<endl;
}
return 0;
}
3.猪钱罐
题意:01背包,要求最小价值
思路:先全部赋值极大值,再每次比较取最小
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
int main()
{
int t;
cin>>t;
while (t--)
{
int x,y,n,i,j;
cin>>x>>y;
int m=y-x;
cin>>n;
int v[5100]={},w[5100]={},dp[100100]={};
memset(dp,0x3f,sizeof(dp));
dp[0]=0;
for (i=0;i<n;i++)
{
cin>>v[i]>>w[i];
}
for (i=1;i<=m;i++)
for (j=0;j<n;j++)
{
if (w[j]>i) continue;
dp[i]=min(dp[i],dp[i-w[j]]+v[j]);
}
dp[j]=min(dp[j],dp[j-w[i]]+v[i]);
if (dp[m]!=0x3f3f3f3f) cout<<"The minimum amount of money in the piggy-bank is "<<dp[m]<<'.'<<endl;
else cout<<"This is impossible."<<endl;
}
return 0;
}
学习了一个小技巧,memset用0x3f,可以将整个数组赋值为极大值
4.划分
题意:现有面值为1、2、3、4、5、6的硬币若干枚,现在需要知道能不能将这些硬币分成等额的两堆。
思路:多重背包存在性问题 (比较复杂)
01背包变型,可以将同样的硬币分成2的n次方个一堆,不能分的放一堆,然后再用01背包来解。
#include <iostream>
#include <algorithm>
using namespace std;
int main()
{
int t=0;
while (1)
{
int i,p=0,sum=0,j,v[100010]={},dp[100010]={};
int a[7]={};
for (i=1;i<=6;i++)
{
cin>>a[i];
if (a[i]!=0)
p=1;
sum+=a[i]*i;
}
if (!p) break;
t++;
cout<<"Collection #"<<t<<':'<<endl;
if (sum%2)
{
cout<<"Can't be divided."<<endl<<endl;
continue;
}
sum/=2;
int k=0;
for (i=1;i<=6;i++)
{
for (j=1;j<=a[i];j*=2)
{
v[k++]=i*j;
a[i]-=j;
}
if (a[i]>0)
v[k++]=a[i]*i;
}
for (i=0;i<k;i++)
for (j=sum;j>=v[i];j--)
dp[j]=max(dp[j],dp[j-v[i]]+v[i]);
if (dp[sum]==sum)
cout<<"Can be divided."<<endl<<endl;
else cout<<"Can't be divided."<<endl<<endl;
}
return 0;
}