hdu 5155(DP+排列组合)

问题描述
一个盒子由n*m个格子组成，有一些格子里会有闪闪发光的宝石。
现在有求盒子从左边看过去，每一行都闪烁着光芒，从前面看过去，每一列也都闪烁着光芒。
问：盒子里的宝石有多少种分布情况。

答案有可能很大，所以输出答案对1000000007取模。

解题思路：状态方程很容易想到，dp[i][j]表示到达(i,j)时的方案数。关键是这里的递推关系怎么出来。

这里参考了别人的思路：http://blog.csdn.net/helloworld10086/article/details/42456775

这里最难理解的就是第二种情况：这里实际上是把dp[i-k][j-1]的情况插空了，即前i-k行应该是连在一起的，现在把它们分开，但不会影响我的方案数。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;

typedef __int64 ll;
const int MOD = 1000000007;
const int N = 51;
ll c[N][N],dp[N][N],t[N];
int n,m;

void init() {
	memset(c,0,sizeof(c));
	c[0][0] = 1;
	for (int i = 1;i <= N; i++) {
		c[i][0] = 1;
		for (int j = 1;j <= i; j++)
			c[i][j] = (c[i-1][j-1] + c[i-1][j]) % MOD;
	}
	t[0] = 1;
	for(int i = 1; i <= N; i++) {
		t[i] = t[i-1] *2 % MOD;
	}
}

int main() {
	init();
	while(scanf("%d%d",&n,&m) != EOF) {
		for(int i = 0; i <= n; i++) 
			dp[i][1] = 1;
		for(int j = 0; j <= m; j++)
			dp[1][j] = 1;

		for(int i = 2; i <= n; i++) {
			for(int j = 2; j <= m; j++) {
				dp[i][j] = dp[i][j-1] * (t[i] - 1) % MOD;
				for(int k = 1; k < i; k++) {
					dp[i][j] = (dp[i][j] + dp[i-k][j-1] * c[i][k] % MOD * t[i-k] % MOD) % MOD;
				}
			}	
		}
		printf("%I64d\n",dp[n][m]);
	}
	return 0;
}