POJ2441 Arrange the Bulls 状态压缩DP

        N头牛M个牲口棚，每头牛不希望与其它牛共享牲口棚，并且每头牛只希望待在自己喜欢的某几个牲口棚中，给定N头牛以及各自可以待的牲口棚编号，求一共有多少种满足条件的分配方法。 N,M不超过20.

        根据前k-1头牛以及它们各自占有牲口棚的情况，我们可以推算出第k头牛可行的分配方案。一共有M个牲口棚，所以总共的状态有2^M个，我们可以将这些状态压缩成一个整数，第i个二进制位为1则代表第i个牲口棚已经被占有，我们设dp(n, k)为前n头牛分配后各个牲口棚占有状态为k的总共方案数，则对于第n+1头牛的每一个可以占有的牲口棚s，有：

dp(n + 1, k | (1 << s)) = dp(n + 1, k | (1 << s))  + dp(n, k)，并且k的第s二进制位上为0，否则就会有重复占用的情况了，初始时，dp(0,0) = 1。

#include <stdio.h>
#include <vector>
#include <math.h>
#include <string.h>
#include <string>
#include <iostream>
#include <queue>
#include <list>
#include <algorithm>
#include <stack>
#include <map>
#include <time.h>
using namespace std;

int dp[2][2000000];

int main()
{
#ifdef _DEBUG
	freopen("e:\\in.txt", "r", stdin);
#endif
	int n, m;
	memset(dp, 0, sizeof(dp));
	vector<vector<int> >barns;
	scanf("%d %d", &n, &m);
	barns.resize(n);
	for (int i = 0; i < n; i++)
	{	
		int size;
		scanf("%d", &size);
		barns[i].resize(size);
		for (int j = 0; j < size; j++)
		{
			int tmp;
			scanf("%d", &tmp);
			barns[i][j] = tmp - 1;
		}
	}
	dp[0][0] = 1;
	for (int i = 1; i <= n; i++)
	{
		memset(dp[i & 1], 0, sizeof(dp[i & 1]));
		for (int j = 0; j < (1 << m); j++)
		{
			if (dp[(i - 1) & 1][j] == 0)
			{
				continue;
			}
			for (int k = 0; k < barns[i - 1].size(); k++)
			{
				if (((j >> barns[i - 1][k]) & 1) == 0)
				{
					dp[i & 1][j | (1 << barns[i - 1][k])] += dp[(i - 1) & 1][j] ;
				}
			}
		}
	}
	int count = 0;
	for (int i = 0; i < (1<< m); i++)
	{
		count += dp[n & 1][i];
	}
	printf("%d\n", count);
	return 0;
}