LA 4794 sharing chocolate 分享巧克力

题目大意
给定一块长为x宽为y的矩形巧克力，每次操作可以沿着一条直线把一块巧克力切割成两块长宽均为整数的巧克力. 给定n个面积的巧克力试求能不能切成.n<15

分析

注意到n的规模很小。可以把n有关的子集作为动态规划的一部分. 设f(r,c,s)表示r行c列的巧克力是否可以切割成面积为S集合。
对于一块面积为（r0c0)的巧克力,状态转移为:
存在1<=r0<r和s中的集合s0 ，有f(r,c,s)=f(r-r0,c,s-s0);
存在1<=c0<c和s中的集合s0，有f(r,c,s)=f(r,c-c0,s-s0);
前者对应横着切，后者对应竖着切。
我们进一步的优化;
如果rc!=s ,我们永远且不成,显然f(r,c,s)=0；
不妨设r<c,用一个量来表示f(r,s)表示状态，这样就会降低复杂度
代码如下

#include<iostream>
#include<algorithm>
#include<string.h>
#define maxn 100+10
#define maxs (1<<16)
using namespace std;

int q[maxn],sum[maxs],vis[maxs][maxn],d[maxs][maxn];
int n;

int bitcount(int x) { return x == 0 ? x : bitcount(x / 2) + (x & 1); }
bool dp(int s, int x)
{
	if (vis[s][x]) return d[s][x];
	vis[s][x] = 1;
	int &ans = d[s][x];
	if (bitcount(s) == 1) return ans = 1;
	int y = sum[s] / x;
	for (int s0 = s; s0; s0 = (s0 - 1)&s) {
		int s1 = s - s0;
		if (sum[s0] % x == 0 && dp(s0, min(x, sum[s0] / x)) && dp(s1, min(x, sum[s1] / x)))
			return ans=1;
		if (sum[s0] % y == 0 && dp(s0, min(y, sum[s0] / y)) && dp(s1, min(y, sum[s1] / y)))
			return ans=1;
	}

	return ans=0;
}

int main()
{
	int x, y,kase=0;
	while (cin >> n && n != 0)
	{
		cin >> x >> y;
		if (x > y)  swap(x, y);
		for (int i = 0; i < n; i++)
			cin >> q[i];
		memset(sum, 0, sizeof(sum));
		memset(vis, 0, sizeof(vis));		//少了一个判定数组初始化在，报超时;
		for (int s = 0; s < (1 << n); s++){
			for (int i = 0; i < n; i++)
				if ((1 << i)&s)
					sum[s] += q[i];
		}	
		int all = (1 << n) - 1;
		int ans;
		if (x*y != sum[all]) ans = 0;
		else ans = dp(all, x);
		printf("Case %d: %s\n", ++kase, ans ? "Yes" : "No");
	}
	return 0;
}