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;
}
posted @ 2017-09-12 08:31  tbhkoymiads  阅读(140)  评论(0编辑  收藏  举报