Poj 1837 Balance【题解报告|dp】

题目大意：

有一个天平，天平左右两边各有若干个钩子，总共有C个钩子，有G个钩码，求将钩码全部挂到钩子上使天平平衡的方法的总数。

其中可以把天枰看做一个以x轴0点作为平衡点的横轴

与这道题有点类似，都是将某个元素看作体积，另一个元素看作价值进行0-1背包，而且都需要处理正负数，建议一试。

输入：

2 4 //C 钩子数 与 G钩码数

-2 3 //负数：左边的钩子距离天平中央的距离；正数：右边的钩子距离天平中央的距离c[k]

3 4 5 8 //G个重物的质量w[i]

dp思路：

每向天平中方一个重物，天平的状态就会改变，而这个状态可以由若干前一状态获得。

首先定义一个平衡度j的概念

当平衡度j=0时，说明天枰达到平衡，j>0，说明天枰倾向右边（x轴右半轴），j<0则相反

那么此时可以把平衡度j看做为衡量当前天枰状态的一个值

因此可以定义一个 状态数组$dp[i][j]$，意为在挂满前i个钩码时，平衡度为j的挂法的数量。(这里为什么不以位置为体积，即dp[i]意为使用完前i个位置时平衡度为j的挂法的数量，读者可以思考一下。提示：重复使用)

由于距离c[i]的范围是$-15-15$，钩码重量的范围是1~25，钩码数量最大是20

因此最极端的平衡度是所有物体都挂在最远端，因此平衡度最大值为$j=15\ast 20\ast 25=7500$。原则上就应该有$dp[120][-75007500]$。

因此为了不让下标出现负数，做一个处理，使使得数组开为 $dp[120][015000]$，则当$j=7500$时天枰为平衡状态

那么每次挂上一个钩码后，对平衡状态的影响因素就是每个钩码的 力臂

力臂=重量 *臂长 = $w[i]\ast c[k]$

那么若在挂上第i个砝码之前，天枰的平衡度为j

(换言之把前i-1个钩码全部挂上天枰后，天枰的平衡度为j)

则挂上第i个钩码后，即把前i个钩码全部挂上天枰后，天枰的平衡度 $j=j+w[i]\ast c[k]$

其中$c[k]$为天枰上钩子的位置，代表第i个钩码挂在不同位置会产生不同的平衡度

不难想到，假设 $dp[i-1][j]$ 的值已知，设$dp[i-1][j]=num$

（即已知把前i-1个钩码全部挂上天枰后得到状态j的方法有num次）

那么$dp[i][j+w[i]\ast c[k]]=dp[i-1][j]=num$

(即以此为前提，在第k个钩子挂上第i个钩码后，得到状态$j+w[i]\ast c[k]$的方法也为num次)

想到这里，利用递归思想，不难得出 状态方程$\mathbf{d}\mathbf{p}[\mathbf{i}][\mathbf{j}+\mathbf{w}[\mathbf{i}]\ast \mathbf{c}[\mathbf{k}]]=\sum （\mathbf{d}\mathbf{p}[\mathbf{i}-1][\mathbf{j}]）$

有些前辈推导方式稍微有点不同，得到的 状态方程为$dp[i][j]=\sum (dp[i-1][j-c[i]\ast w[i]])$

其实两条方程是等价的，这个可以简单验证出来，而且若首先推导到第二条方程，也必须转化为第一条方程，这是为了避免下标出现负数

结论：

最终转化为01背包问题

状态方程$dp[i][j+w[i]\ast c[k]]=\sum （dp[i-1][j]）$

初始化：$dp[0][7500]=1;$ //不挂任何重物时天枰平衡，此为一个方法

复杂度$O(C\ast G\ast 15000)$ 完全可以接受

#define ll long long
#define vec vector<int>
#define P pair<int,int>
#define inf 0x3f3f3f3f
#define MAX 15005

int dp[25][15005];//dp[i][j]:前i个重物能凑出来的值
int C, G, c[25], g[25];

int main() {
	while (scanf("%d %d", &C, &G) != EOF) {
		for (int i = 1; i <= C; i++)scanf("%d", &c[i]);
		for (int i = 1; i <= G; i++)scanf("%d", &g[i]);

		memset(dp, 0, sizeof(dp)); dp[0][7500] = 1;
		for (int i = 1; i <= G; i++) {//每一个重物
			for (int j = 1; j <= C; j++) {//每一个位置
				if (c[j] >= 0) {//右侧的位置
					for (int k = 0; k < MAX - c[j] * g[i]; k++)
						if (dp[i - 1][k])
							dp[i][k + c[j] * g[i]] += dp[i - 1][k];
				}
				else//挂在左侧某个位置
					for (int k = MAX - 1; k >= -c[j] * g[i]; k--)
						if (dp[i - 1][k])
							dp[i][k + c[j] * g[i]] += dp[i - 1][k];
			}
		}
		
		printf("%d\n", dp[G][7500]);
	}
}