【题解】机器分配（DP）

题面

【题目描述】
总公司拥有高设备$M$台，准备分给下属的$N$个分公司。各分公司获得这些设备，可以为国家提供一定的盈利。问：如何分配这$M$台设备才能使国家得到的盈利最大？求出最大盈利值。其中$M\le 15,N\le 10$。分配原则：每个公司有权获得任意数目的设备，但总台数不超过设备数$M$。
【输入】
第$1$行有两个数，第$1$个数是分公司数$N$，第$2$个数是总设备台数$M$。接下来是$N\ast M$的矩阵，表明了第$I$个公司分配$J$台机器的盈利。
【输出】
第$1$行输出最大盈利值，接下来$N$行，每行$2$个数，即分公司编号和该公司获得设备台数。
【样例输入】

3 3
30 40 50
20 30 50
20 25 30

【样例输出】

70
1 1
2 1
3 1

算法分析

$DP$。
状态：
$dp[i][j]$——前$i$个公司分$j$台设备的最大盈利。

状态转移方程：
如果前$i-1$个公司已经分配完了，考虑第$i$个公司分配的机器数量，设分得$k$台，$k$可能为$0$~$j$，那么前$i-1$个公司只能分配$j-k$台。
设$a[x][y]$表示第$x$个公司分配$y$台设备的盈利。
因此，只需要考虑第$i$个公司分配了多少台机器：

$dp[i][j]=max${ $dp[i][j],dp[i-1][j-k]+a[i][k]$ }，$0\le k\le j。$

边界：
$dp[i][0]=0$，$i$个公司分配$0$台设备盈利为$0$.
$dp[0][j]=0$，$0$个公司分配$j$台设备盈利为$0$。

输出方案
最后的答案为$dp[n][m]$，即$n$个公司分$m$台设备的最优解。
通过状态转移方程可以知道，$dp[n][m]$是由$dp[n-1][j-k_n]+a[n][k_n]$得到的，即前$n-1$个公司分$m-k_n$台设备，第$n$个公司分配$k_n$台设备。数组$a$是已知的，$dp$数组也是已知的，而$k_n$是未知的，因此，可以枚举$k_n$，检验$k_n$在何时满足$dp[n][m]$，则找到了第$n$个公司分配到的设备数。
同理，继续找$dp[n-2][m-k_n-k_{n-1}]+a[n-1][k_{n-1}]$，检验第$n-1$个公司分配的设备数…
采用递归来实现。

参考程序

#include<bits/stdc++.h>
using namespace std;
int n,m;
int a[200][200],dp[200][200],num[200],mx,s[200],ans=0;
void print(int i,int j)	//递归寻找答案，前i个公司分j台 
{
	if(ans) return;		//找到答案不再找了 
    if(i==0) 
    {
    	int tmp=0;
    	for(int k=1;k<=n;k++)
    		tmp+=s[k];
    	if(tmp==m)			//m台设备全部分完 
    	{
    		for(int k=1;k<=n;k++)
    			cout<<k<<' '<<s[k]<<endl;	
			ans=1; 
		}
		return;
	}
    for(int k=0;k<=j;k++)	//前i-1个公司分j-k台，第i个公司分k台
    {
        if(dp[i][j]==dp[i-1][j-k]+a[i][k])
        {
			s[i]=k;				//记录 
            print(i-1,j-k);
        }
    }
}
int main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            cin>>a[i][j];
    for(int i=1;i<=n;i++)
    {
         
        for(int j=1;j<=m;j++)
        {
            mx=0;
            for(int k=0;k<=j;k++)	//考虑第i个公司分配了k台设备 
            {
                if(dp[i-1][j-k]+a[i][k]>mx)			//取最优解 
                    mx=(dp[i-1][j-k]+a[i][k]);
            }
            dp[i][j]=mx;	
        }
    }
    cout<<dp[n][m]<<endl;
    print(n,m);
    return 0;
}