洛谷 P2051 [AHOI2009]中国象棋

题目描述

这次小可可想解决的难题和中国象棋有关，在一个N行M列的棋盘上，让你放若干个炮（可以是0个），使得没有一个炮可以攻击到另一个炮，请问有多少种放置方法。大家肯定很清楚，在中国象棋中炮的行走方式是：一个炮攻击到另一个炮，当且仅当它们在同一行或同一列中，且它们之间恰好 有一个棋子。你也来和小可可一起锻炼一下思维吧！

输入输出格式

输入格式：
一行包含两个整数N，M，之间由一个空格隔开。

输出格式：
总共的方案数，由于该值可能很大，只需给出方案数模9999973的结果。

输入输出样例

输入样例#1：
1 3
输出样例#1：
7
说明

样例说明

除了3个格子里都塞满了炮以外，其它方案都是可行的，所以一共有2*2*2-1=7种方案。

数据范围

100%的数据中N和M均不超过100

50%的数据中N和M至少有一个数不超过8

30%的数据中N和M均不超过6

---

【分析】
特别神奇的状态转移
唉又是一道省选DP神题 还带组合数= =

我手玩几组数据之后发现行列貌似可以交换的

然后还是不会做 只好求助于各神犇 题解比较少= =凑个热闹

dp[i][j][k]表示在前i行中有j列已经有2个 有k列已经有1个 那么0的个数就是m-j-k个

转移方程看代码

这一行可以不放

那么直接加上dp[i-1][j][k]

放一个的话

可以把一列0变成1

可以把一列1变成2

放两个的话

可以把两列0变成1

可以把两列1变成2

还有 可以把一列0变成1 再把另外一列1变成2

关于组合数要说一下。。

放一个好说 直接乘1或0的个数完事（加在任意一个上面都是可以的）

放两个的话 前两个也好弄 直接C（n，2） （n为可选个数）

最后一个呢 先乘上原来0的个数 再乘上原来1的个数 乘法原理嘛= =

---

【代码】

//P2051 [AHOI2009]中国象棋
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mod=9999973;
ll f[105][105][105];
int n,m,ans;
int main()
{
    int i,j,k;
    scanf("%d%d",&n,&m);
    f[0][0][0]=1;
    fo(i,1,n)
      fo(j,0,m)
        fo(k,0,m-j)
        {
            f[i][j][k]=(f[i][j][k]+f[i-1][j][k])%mod;
            if(j) f[i][j][k]=(f[i][j][k]+f[i-1][j-1][k]*(m-j-k+1))%mod;  //0->1
            if(k) f[i][j][k]=(f[i][j][k]+f[i-1][j+1][k-1]*(j+1))%mod;  //1->2
            if(k&&j) f[i][j][k]=(f[i][j][k]+f[i-1][j][k-1]*j*(m-j-k+1))%mod;  //0->1 and 1->2
            if(j>=2) f[i][j][k]=(f[i][j][k]+f[i-1][j-2][k]*(m-j-k+2)*(m-j-k+1)/2)%mod; //0->1 (*2)
            if(k>=2) f[i][j][k]=(f[i][j][k]+f[i-1][j+2][k-2]*(j+2)*(j+1)/2)%mod; //1->2 (*2)
        }
    fo(j,0,m)
      fo(k,0,m-j)
        ans=(ans+f[n][j][k])%mod;
    printf("%d\n",ans);
    return 0;
}