专题训练之dp

一、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;
}