UVA-12371 Guard

动态规划的一道题目，有一定的难度。首先要分清楚情况的转移，一共有两种情况，第一种是从（i-1）*（i-1）的矩阵进行转移，这一种情况又可以分为两种小的情况，第一种是对于后续产生的矩阵i*i，该矩阵的最右下角（也就是坐标（i，i））是没有棋子的，那么我们可以从i-1列中选择一列x，将这一列的内容填充到第i列中，然后在前面的i-1行中选择一行t，将这一行最右的棋子放到第x列中，同时在i行x列中放置一个棋子，在t行i列中补上一个棋子就满足要求了，那么这种情况的总数为dp[i-1][j]*(i-1)*(i-1)，dp的第一位代表的是（i-1）*（i-1）的矩阵，第二位代表的是这个矩阵中一共存在多少个圈；那么在第二中情况，也就是最终的i*i矩阵的最右下的位置放置好一颗棋子，那么这个时候可以从i-1行中选择某一行，假设为x，同时假设该行对应的最右边的棋子在y列，那么可以将这一行的棋子拿掉，同时在x行i列以及i行y列中均放置棋子就满足条件了，那么这种情况对应的总数是dp[i-1][j]*(i-1)。对于第二种情况，我们可以发现矩阵的行列数仅仅增加1是不会增加圈子的总数的，所以第二种情况是从dp[i-2][j-1]来考虑的，那么也就是我们从第i行中选择某两个位置，同时从i-1行中选择某一行，这样正好组成新的圈，让原有的圈的个数增加1就行了，那么这种情况的总数也就是dp[i-2][j-1]*C(i,2)*(i-1)。综合上面的这些情况，列出最终的状态转移方程，打表即可。具体实现见如下代码：

#include<iostream>
#include<vector>
#include<string>
#include<set>
#include<stack>
#include<queue>
#include<map>
#include<algorithm>
#include<cmath>
#include<iomanip>
#include<cstring>
#include<sstream>
#include<cstdio>
#include<deque>
#include<functional>
using namespace std;

typedef long long LL;
const LL mod = 1000000007;
const LL maxn = 100000;
LL dp[maxn+10][55];

void Init(){
	memset(dp, 0, sizeof(dp));
	dp[2][1] = 1;
	for (LL i = 3; i <= maxn; i++){
		for (LL j = 1; j <= 50&&j<=i; j++){
			dp[i][j] = (dp[i - 1][j] * (i-1) % mod*(i-1)%mod+
				dp[i - 1][j] * (i - 1) % mod + 
				dp[i - 2][j - 1] * ((LL)i*(i - 1) / 2) % mod*(i - 1) % mod) % mod;
		}
	}
}

int main(){
	int T;
	cin >> T;
	Init();
	for (int i = 1; i <= T; i++){
		LL N, K;
		cin >> N >> K;
		cout << "Case " << i << ": " << dp[N][K] << endl;
	}
	return 0;
}