洛谷P2066 机器分配【DP+求方案】【黄】

Date：2022.02.25
题目描述
总公司拥有高效设备M台，准备分给下属的N个分公司。各分公司若获得这些设备，可以为国家提供一定的盈利。问：如何分配这M台设备才能使国家得到的盈利最大？求出最大盈利值。其中M≤15，N≤10。分配原则：每个公司有权获得任意数目的设备，但总台数不超过设备数M。
输入格式
第一行有两个数，第一个数是分公司数N，第二个数是设备台数M。
接下来是一个N*M的矩阵，表明了第 I个公司分配 J台机器的盈利。
输出格式
第1行为最大盈利值
第2到第n为第i分公司分x台
P.S.要求答案的字典序最小
输入输出样例
输入 #1复制
3 3
30 40 50
20 30 50
20 25 30
输出 #1复制
70
1 1
2 1
3 1

思路：$f[i][j]:$给前$i$个公司总共分配$j$台机器的最大盈利。
状态转移：$f[i][j]=max(f[i-1][j-k]+a[i][k],f[i][j])$;【其中所枚举的$k$为第$i$个公司选多少台。】
接着不妨总结一个求dp中方案数的通法，根据状态的结果反推。能力有限，给不出严谨证明，但我们不妨一想：每个上一步的状态都是从起点走到上一步的“局部最优解”；类比过来，走到当前状态的最优解即为：上一步的局部最优解+这一步的某一种选择。因此，根据当前最优解求上一步的局部最优解方式显然是“等式匹配”：即只有当上一步的局部最优解+这一步的某种选择==当前最优解，才能保证这一步的某种选择是包含在最优解中的，通俗一点就是递归反推。
代码如下：

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 20,mod=1e6+7;
typedef long long LL;
LL n,m,k,a[N][N],f[N][N];
void print(LL i,LL j,LL maxx)
{
    if(i==0) return;
    for(int k=0;k<=j;k++)
        if(maxx==f[i-1][k]+a[i][j-k])
        {
            print(i-1,k,f[i-1][k]);
            cout<<i<<' '<<j-k<<endl;
            break;
        }
}
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++)
            for(int k=0;k<=j;k++)
                f[i][j]=max(f[i][j],f[i-1][j-k]+a[i][k]);
    cout<<f[n][m]<<endl;
    print(n,m,f[n][m]);
    return 0;
}