description

在N行M列的棋盘上,放若干个炮可以是0个,使得没有任何一个炮可以攻击另一个炮。请问有多少种放置方法?中国象棋中炮的行走方式大家应该很清楚吧.


analysis

  • \(DP\),容易知道每行至多有两个炮,否则会互相打到

  • 设\(f[i][j][k]\)表示到放到第\(i\)行,有\(j\)列放了一个炮,\(k\)列放了两个炮的方案数

  • 该行不放炮,则直接继承上一行的答案

\[f[i][j][k]+=f[i-1][j][k]
\]

  • 一个炮放在没有炮的列上,一个炮的列数\(+1\),且有\(m-k-(j-1)\)个没有炮的列可以放

\[f[i][j][k]+=f[i-1][j-1][k]*[m-k-(j-1)]
\]

  • 一个炮放在一个炮的列上,一个炮的列数\(-1\),两个炮的列数\(+1\),且有\(j+1\)个一个炮的列可以放

\[f[i][j][k]+=f[i-1][j+1][k-1]*(j+1)
\]

  • 一个炮放在一个炮的列上,一个炮放在没有炮的列上,两个炮的列数\(+1\),且分别有\(j\)列、\(m-(k-1)-j\)列可以放

\[f[i][j][k]+=f[i-1][j][k-1]*j*[m-(k-1)-j]
\]

  • 两个炮放在没有炮的列上,一个炮的列数\(+2\),且有\(C^{2}_{m-(j-2)-k}\)种方案

\[f[i][j][k]+=f[i-1][j-2][k]*C^{2}_{m-(j-2)-k}
\]

  • 两个炮放在一个炮的列上,一个炮的列数\(-2\),两个炮的列数\(+2\),且有\(C^{2}_{j+2}\)种方案

\[f[i][j][k]+=f[i-1][j+2][k-2]*C^{2}_{j+2}
\]

  • 如此转移即可,注意判断边界

code

#pragma GCC optimize("O3")
#pragma G++ optimize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define MAXN 105
#define ha 9999973
#define ll long long
#define reg register ll
#define fo(i,a,b) for (reg i=a;i<=b;++i)
#define fd(i,a,b) for (reg i=a;i>=b;--i) using namespace std; ll f[MAXN][MAXN][MAXN];
ll c[MAXN][MAXN];
ll n,m,ans; inline ll read()
{
ll x=0,f=1;char ch=getchar();
while (ch<'0' || '9'<ch){if (ch=='-')f=-1;ch=getchar();}
while ('0'<=ch && ch<='9')x=x*10+ch-'0',ch=getchar();
return x*f;
}
inline ll C(ll n){return n*(n-1)/2;}
int main()
{
n=read(),m=read(),f[0][0][0]=1;
fo(i,1,n)
{
fo(j,0,m)
{
fo(k,0,m-j)
{
f[i][j][k]=f[i-1][j][k];//不填
if (k-1>=0)
(f[i][j][k]+=f[i-1][j+1][k-1]*(j+1)%ha)%=ha;//一颗填一个炮的列
if (j-1>=0)
(f[i][j][k]+=f[i-1][j-1][k]*(m-(j-1)-k))%=ha;//一颗填没有炮的列
if (k-1>=0)
(f[i][j][k]+=f[i-1][j][k-1]*j%ha*(m-j-(k-1)))%=ha;//一颗填一个炮的列,一颗填没有炮的列
if (j-2>=0)
(f[i][j][k]+=f[i-1][j-2][k]*C(m-(j-2)-k))%=ha;//两颗填没有炮的列
if (k-2>=0)
(f[i][j][k]+=f[i-1][j+2][k-2]*C(j+2))%=ha;//两颗填一个炮的列
}
}
}
fo(i,0,m)fo(j,0,m-i)(ans+=f[n][i][j])%=ha;
printf("%lld\n",ans);
return 0;
}

最新文章

  1. Delphi的三目运算 ifthen 和iif
  2. canvas绘图、WebGL、SVG
  3. 使用netty实现的tcp通讯中如何实现同步返回
  4. MYSQL 数据库导入导出命令
  5. iOS开发中的错误整理,pod文件Install失败
  6. Ubuntu 14.10 下CPU实时监控mpstat命令详解
  7. android 下载图片出现SkImageDecoder::Factory returned null,BitmapFactory.Options压缩
  8. Notepad++中常用的插件
  9. 将string转化为char*的方法
  10. FormsAuthentication.GetRedirectUrl 方法
  11. Razor 语法快速参考
  12. SVN库迁移整理方法总结
  13. 【CKEditor ASP.NET】解决360安全浏览器极速模式下不显示
  14. Oracle:对表的CREATE、ALTER、INSERT、RENAME、DELETE操作练习以及主外键约束
  15. linux下mysql修改字符集
  16. POJ1509 Glass Beads(最小表示法 后缀自动机)
  17. Centos7下chkconfig设置MySql自动启动
  18. 洛谷P4570 [BJWC2011]元素 线性基
  19. FreeModbus LINUXTCP Compile ERROR
  20. IOS中@property的属性weak、nonatomic、strong、readonly等介绍

热门文章

  1. STL 队列模板实现
  2. 对象拷贝extend
  3. 内网端口转发[netsh]
  4. ARM TK1 安装ROS- indigo
  5. 两台linux服务器传文件方法
  6. Laravel 项目运行 phpunit 测试结果只显示点号
  7. Maven pom配置(Spring+SpringMvc+mybaties)
  8. Java初识方法
  9. delphi根据不同图片生成不规则窗口的实现(仅限于BMP格式)
  10. 牛客多校第十场 E Hilbert Sort 递归,排序