【洛谷 P1240】诸侯安置

题目描述
很久以前，有一个强大的帝国，它的国土成正方形状，如图所示。
这个国家有若干诸侯。由于这些诸侯都曾立下赫赫战功，国王准备给他们每人一块封地(正方形中的一格)。但是，这些诸侯又非常好战，当两个诸侯位于同一行或同一列时，他们就会开战。如下图2—3为n＝3时的国土，阴影部分表示诸侯所处的位置。前两幅图中的诸侯可以互相攻击，第三幅则不可以。

国王自然不愿意看到他的诸侯们互相开战，致使国家动荡不安。 因此，他希望通过合理的安排诸侯所处的位置，使他们两两之间都不能攻击。现在，给出正方形的边长n，以及需要封地的诸侯数量k，要求你求出所有可能的安置方案数。(n≤l00，k≤2n2-2n+1)
由于方案数可能很多，你只需要输出方案数除以504的余数即可。

输入格式：
仅一行，两个整数n和k，中间用一空格隔开。

输出格式：
一个整数，表示方案数除以504的余数

输入样例#1：
2 2
输出样例#1：
4

解题思路

本题最重要的东西——图的转化：

从左到右Dp,设f[i][j]表示前i列第i列选了j个且第i列必须选的方案数，于是就有了状态转移方程：
f[i][j]=∑f[k][j-1]*(len[i]-(j-1))
因为前面已经选了j-1个，那么在第i列就有len[i]-（j-1）个位置可选。
代码：

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<string>
#include<iostream>
#include<cstdlib>
using namespace std;
int f[205][40005];
int ans,n,pk;
const int p=504;
int main()
{
    scanf("%d%d",&n,&pk);
    f[0][0]=1;          
    for(int i=1;i<n*2;i++)
    for(int j=1;j<=pk && j<=i;j++)
    for(int k=i-1;k>=0;k--)
    {
       if(i&1) f[i][j]=(f[i][j]+f[k][j-1]*(i-j+1))%p;
       else f[i][j]=(f[i][j]+f[k][j-1]*(i-j))%p;
    }

    for(int i=0;i<n*2;i++)
    ans=(ans+f[i][pk])%p;
    printf("%d",ans);
    return 0;
}