1837

/*
起初我实现的DP比较慢,原因是我忽略了一个很重要的事实,那就是所有的砝码都必须用上,这就是说,每一个砝码都会出现的,既然如果,不管
放的顺序是什么,都可以看成从第一个开始放,然后第二个,然后第三个。这个发现可以帮助我们消去一个循环

现在说下我的最初的DP: DP[i][j][k], 到第i时,差为j,选择k个的个数
DP[i+1][j+(i+1号砝码所有可能的放置)][k+1] += DP[i][j][k],很显然这里当k<i的这种情况在以后是用不到的。试想如果出现了DP[i][j][i-1],也就是说少了第i个砝码,
以后还会在添加进去吗?不会了,因为砝码就是我的阶段。所以这纯粹是浪费时间和内存。如果题目中没有说明全部都用的话,这个方法还是可以的

下面是针对题目的改进
DP[i][j],i既是个数,也是阶段,j是差值
这样的转移方程是就是,对每个i,枚举其所有可能性,转移到i+1.最后找出等于0的情况就行了
*/

// include file
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cctype>
#include <ctime>

#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <bitset>

#include <algorithm>
#include <string>
#include <vector>
#include <queue>
#include <set>
#include <list>
#include <functional>

using namespace std;

// typedef
typedef long long LL;
typedef unsigned long long ULL;

// 
#define read freopen("in.txt","r",stdin)
#define write freopen("out.txt","w",stdout)
#define FORi(a,b,c) for(int i=(a);i<(b);i+=c)
#define FORj(a,b,c) for(int j=(a);j<(b);j+=c)
#define FORk(a,b,c) for(int k=(a);k<(b);k+=c)
#define FORp(a,b,c) for(int p=(a);p<(b);p+=c)
#define FORii(a,b,c) for(int ii=(a);ii<(b);ii+=c)
#define FORjj(a,b,c) for(int jj=(a);jj<(b);jj+=c)
#define FORkk(a,b,c) for(int kk=(a);kk<(b);kk+=c)

#define FF(i,a)    for(int i=0;i<(a);i++)
#define FFD(i,a)   for(int i=(a)-1;i>=0;i--)

#define Z(a) (a<<1)
#define Y(a) (a>>1)

const double eps = 1e-6;
const double INFf = 1e100;
const int INFi = 1000000000;
const LL INFll = (LL)1<<62;
const double Pi = acos(-1.0);

template<class T> inline T sqr(T a){return a*a;}
template<class T> inline T TMAX(T x,T y)
{
	if(x>y) return x;
	return y;
}
template<class T> inline T TMIN(T x,T y)
{
	if(x<y) return x;
	return y;
}
template<class T> inline void SWAP(T &x,T &y)
{
	T t = x;
	x = y;
	y = t;
}
template<class T> inline T MMAX(T x,T y,T z)
{
	return TMAX(TMAX(x,y),z);
}
template<class T> inline T MMIN(T x,T y,T z)
{
	return TMIN(TMIN(x,y),z);
}


// code begin
int DP[22][15010];
int N,M;
int X[22];
int Y[22];
int zero;
int main()
{
	read;
	write;
	int Xmax,Ymax;
	while(scanf("%d %d",&N,&M)==2)
	{
		Xmax = 0;
		FORi(1,N+1,1) 
		{
			scanf("%d",X+i);
			if(abs(X[i])>Xmax) Xmax = abs(X[i]);
		}
		Ymax = 0;
		FORi(1,M+1,1) 
		{
			scanf("%d",Y+i);
			Ymax+=Y[i];
		}
		zero = Ymax*Xmax;

		memset(DP,0,sizeof(DP));
		DP[0][zero] = 1;
		FORi(0,M,1)
		{
			FORj(0,2*zero+1,1)
			{
				if(DP[i][j])
				{
					FORk(1,N+1,1)
					{
						DP[i+1][j+X[k]*Y[i+1]] += DP[i][j];
					}
				}
			}
		}

		printf("%d\n",DP[M][zero]);
	}
	return 0;
}


/*
int DP[2][15010][22];
int N,M;
int X[22];
int Y[22];
int zero;
int main()
{
	read;
	write;
	int Xmax,Ymax;
	while(scanf("%d %d",&N,&M)==2)
	{
		Xmax = 0;
		FORi(1,N+1,1) 
		{
			scanf("%d",X+i);
			if(abs(X[i])>Xmax) Xmax = abs(X[i]);
		}
		Ymax = 0;
		FORi(1,M+1,1) 
		{
			scanf("%d",Y+i);
			Ymax+=Y[i];
		}
		zero = Ymax*Xmax;
		memset(DP,0,sizeof(DP));
		//FORi(0,M+1,1) 
		DP[0][zero][0] = 1;

		FORi(1,M+1,1)
		{
			memset(DP[i&1],0,sizeof(DP[i&1]));

			FORk(0,2*zero+1,1)
			{
				FORp(0,M,1)
				{
					if(DP[(i+1)&1][k][p]!=-1)
					{
						//对每一个砝码
						FORj(1,N+1,1)
						{
							if(X[j]<0)
							{
								if(k+X[j]*Y[i]>=0) DP[i&1][k+X[j]*Y[i]][p+1] += DP[(i+1)&1][k][p];
							}
							else
							{
								if(k+X[j]*Y[i]<=2*zero) DP[i&1][k+X[j]*Y[i]][p+1] += DP[(i+1)&1][k][p];
							}
						}
					}
				}
			}
			

			FORk(0,2*zero+1,1)
			{
				FORp(0,M+1,1)
				{
					DP[i&1][k][p] = TMAX(DP[i&1][k][p],DP[(i+1)&1][k][p]);
					//printf("[%d %d %d %d]\n",i,k,p,DP[i&1][k][p]);
				}
			}
			//printf("\n");
		}

		printf("%d\n",DP[M&1][zero][M]);
	}
	return 0;
}
*/

posted @ 2011-04-12 15:58  AC2012  阅读(259)  评论(0编辑  收藏  举报