poj3181【Dollar Dayz】

**

做完这道题，心里五味陈杂，明明是最水的一道题，我却做了最长的时间。

题意是求用1～k的和表示n的方案数。

显然是个计数dp，但我不会。思考半小时未果。

然后找尹鹏哲，他给我讲了个错的dp方程，结果调试半小时未果，看了别人的代码，发现别人虽然滚了一维，但和尹鹏哲讲的十分相似。

后来我又想了想，觉得尹鹏哲的dp方程有点问题，自己改了一下，但当时并没有验证。实际上，由于对题解方程错误理解，当时我已经形成了思维定势，认为方程是对的，却没有发现方程中一个十分微小但十分致命的错误——我把完全背包和01背包的模型搞混了。正是这个错误，让我调+改又浪费了两个小时。

好了，废话说了这么多，只是想让自己记住这次教训，希望下次不要再犯。(但背包问题的模型又多，变化又复杂，所以只能希望了)，下面写题解。

先讲二维的：

设f[i][j]表示用前i个物品，刚好凑出体积j的方案数。这个状态很好想。

然后是转移，就是在这里我的思维出现了偏差。

开始，尹鹏哲告诉我的是f[i][j]+=f[i-1][j-i],这个与一维方程一模一样，显然是不对的，但我太信任尹鹏哲，所以并没有细想。

后来，受尹鹏哲的启发，我自己写出来的方程是f[i][j]=f[i-1][j-i]+f[i-1][j]，与组合数的递推关系一样，表示第i件物品选与不选。

这个方程咋看上去没问题，实际上，这道题与组合数有着本质的区别，组合数每个物品只能选一次，而这道题一个数字显然可以用多次。但傻逼的我忽视了这个差别，结果狂WA不止。

换成背包问题的计数的话，组合数实际上是01背包的方案数，而这道题实际上是完全背包的方案数！

所以，dp方程应该为：f[i][j]=f[i][j-i]+f[i-1][j]

滚成一维就很简单了：f[j]+=f[j-i]

实际上，完全背包和01背包的代码十分相似，若碰到一定要想清楚！

另：要打高精度，详见一维代码。**

二维：

#include <bits/stdc++.h>
using namespace std;

#define db double
#define ll long long
#define RG register

inline int gi()
{
	RG int ret; RG bool flag; RG char ch;
	ret=0, flag=true, ch=getchar();
	while (ch < '0' || ch > '9')
		ch == '-' ? flag=false : 0, ch=getchar();
	while (ch >= '0' && ch <= '9')
		ret=(ret<<3)+(ret<<1)+ch-'0', ch=getchar();
	return flag ? ret : -ret;
}

const db pi = acos(-1.0);
const int N = 142857, inf = 1<<30;

int f[205][2005];

int main()
{
	freopen("Dollar_Dayz.in","r",stdin);
	freopen("Dollar_Dayz.out","w",stdout);
	int n,m,i,j;
	n=gi(), m=gi();
	f[1][0]=1;
	// 两种都可
	// for (i=1; i<=m; ++i)
	// 	for (j=0; j<=n; ++j)
	// 		{
	// 			if (j >= i)
	// 				f[i][j]=f[i][j-i];
	// 			f[i][j]+=f[i-1][j];
	// 		}
	for (i=1; i<=m; ++i)
		for (j=0; j<=n; ++j)
			f[i][j+i]+=f[i][j], f[i+1][j]=f[i][j];
	printf("%d\n",f[m][n]);
	return 0;
}

一维：

#include <bits/stdc++.h>
using namespace std;

#define db double
#define ll long long
#define RG register

inline int gi()
{
	RG int ret; RG bool flag; RG char ch;
	ret=0, flag=true, ch=getchar();
	while (ch < '0' || ch > '9')
		ch == '-' ? flag=false : 0, ch=getchar();
	while (ch >= '0' && ch <= '9')
		ret=(ret<<3)+(ret<<1)+ch-'0', ch=getchar();
	return flag ? ret : -ret;
}

const db pi = acos(-1.0);
const int N = 1e3+5, inf = 1<<30;
const ll lim = 1e18;  //要开 ll

ll f[N],g[N];

int main()
{
	freopen("Dollar_Dayz.in","r",stdin);
	freopen("std.out","w",stdout);
	int n,m,i,j;
	n=gi(), m=gi();
	g[0]=1;
	for (j=1; j<=m; ++j)
		for (i=j; i<=n; ++i)
			{
				f[i]+=f[i-j], g[i]+=g[i-j];
				if (g[i] >= lim)
					f[i]+=g[i]/lim, g[i]%=lim;
			}
	if (f[n])
		printf("%lld%lld\n",f[n],g[n]);
	else
		printf("%lld\n",g[n]);
	return 0;
}