Codeforces 1237F. Balanced Domino Placements
很妙的题
首先先考虑一个简化的问题,现在有一行格子让你填
你要么填一格 要么填两格 有的格子不让你填 问你填了 $a$ 个一格和填了 $b$ 个两格有多少种方案
那么显然先只考虑放两格的方案,这个可以简单 $dp$ 得到,设 $f[i][j]$ 表示前 $i$ 个格子放了 $j$ 个两格的方案数
那么如果 $i,i-1$ 都没障碍,那么 $f[i][j]=f[i-1][j]+f[i-2][j-1]$ ,否则 $f[i][j]=f[i-1][j]$
然后再来考虑填一格的,显然剩下的 $tot-2j$ 个位置都可以随便填,那么方案数为 $C[tot-2j][i]$ ,直接乘法原理乘起来即可
接下来可以考虑怎么把这道题简化到这个情况,假设放了 $a$ 个水平的多米诺,$b$ 个垂直的多米诺
对于行来说,相当于放 $a$ 个 $1$ ,$b$ 个 $2$,对于列就相当于放 $a$ 个 $2$ ,$b$ 个 $1$
注意到每个多米诺可以根据在第几行和第几列来唯一确定,所以我们对行列分别求一下之前那个东西然后乘起来再乘上 $a!b!$ 即可
乘上 $a!b!$ 就相当于把骨牌不同的放置顺序看成不同的放置方案,意思是强制行的第 $i$ 个放置和列的第 $i$ 个放置配对
就原本一维的方案我们让它放置有顺序,然后强制行和列两两匹配,这样才能确定二维平面上的具体位置
因为如果只是行列乘起来,那么没法确定某个骨牌的具体位置,考虑对于 $a$ ,行放了位置 $1,3$ ,列放了位置 $(1,2),(3,4)$
那么多米诺骨牌可以是 $((1,1),(1,2))$ 和 $((3,3),(3,4))$ ,但是也有可能是 $((1,3),(1,4))$ 和 $((3,1),(3,2))$
自己画画图就很容易理解了
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
inline int read()
{
int x=,f=; char ch=getchar();
while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
return x*f;
}
const int N=,mo=;
inline int fk(int x) { return x>=mo ? x-mo : x; }
int n,m,K,Ans,fac[N],C[N][N];
int f[N][N],g[N][N];
bool px[N],py[N];
int main()
{
n=read(),m=read(),K=read(); int mx=max(n,m);
fac[]=; for(int i=;i<=mx;i++) fac[i]=1ll*fac[i-]*i%mo;
for(int i=;i<=mx;i++)
{
C[i][]=;
for(int j=;j<=i;j++)
C[i][j]=fk(C[i-][j-]+C[i-][j]);
}
for(int i=;i<=K;i++)
{
int a=read(),b=read(),c=read(),d=read();
px[a]=px[c]=; py[b]=py[d]=;
}
int cntx=,cnty=;
for(int i=;i<=n;i++) cntx+=px[i];
for(int i=;i<=m;i++) cnty+=py[i];
for(int i=;i<=n;i++) f[i][]=;
for(int i=;i<=n;i++)
for(int j=;j<=mx;j++)
if(px[i]||px[i-]) f[i][j]=f[i-][j];
else f[i][j]=fk(f[i-][j]+f[i-][j-]);
for(int i=;i<=m;i++) g[i][]=;
for(int i=;i<=m;i++)
for(int j=;j<=mx;j++)
if(py[i]||py[i-]) g[i][j]=g[i-][j];
else g[i][j]=fk(g[i-][j]+g[i-][j-]);
for(int i=;i<=mx;i++)
for(int j=;j<=mx;j++)
if(i+j*<=n-cntx&&i*+j<=m-cnty)
Ans=fk(Ans + 1ll*f[n][j]*C[n-cntx-j*][i]%mo *g[m][i]%mo *C[m-cnty-i*][j]%mo *fac[i]%mo *fac[j]%mo );
printf("%d\n",Ans);
return ;
}
Codeforces 1237F. Balanced Domino Placements的更多相关文章
- CF1237F Balanced Domino Placements
题意 给定一个 \(h\) 行 \(w\) 列的方格图,上面已经放置了一些 \(1\times 2\) 的多米诺骨牌. 我们称一个放置多米诺骨牌的方案是好的,当且仅当任何两个多米诺骨牌不占用相同的行与 ...
- [Codeforces 873B]Balanced Substring
Description You are given a string s consisting only of characters 0 and 1. A substring [l, r] of s ...
- CodeForces - 873B Balanced Substring(思维)
inputstandard input outputstandard output You are given a string s consisting only of characters 0 a ...
- Codeforces 873B - Balanced Substring(思维)
题目链接:http://codeforces.com/problemset/problem/873/B 题目大意:一个字符串全部由‘0’和‘1’组成,当一段区间[l,r]内的‘0’和‘1’个数相等,则 ...
- Codeforces 1237E. Balanced Binary Search Trees
传送门 这一题是真的坑人,时间空间都在鼓励你用 $NTT$ 优化 $dp$...(但是我并不会 $NTT$) 看到题目然后考虑树形 $dp$ ,设 $f[i][0/1]$ 表示 $i$ 个节点的树,根 ...
- Codeforces 1237D. Balanced Playlist
传送门 首先显然的,如果一个位置开始播放了两圈还没结束,那么就永远不会结束 先考虑位置 $1$ 开始播放,用一个 $multisetset$ 维护一下当前听的所有歌,直到某一首歌 $r$ 不合法了就停 ...
- Codeforces 1237C2. Balanced Removals (Harder)
传送门 先来考虑一下二维时的情况,那么对于 $x$ 相同的点,我们按 $y$ 排序,然后相邻的一对对消除 最后 $x$ 坐标相同的点最多剩下一个,那么此时所有点的 $x$ 坐标都不一样 再按 $x$ ...
- Codeforces 1237B. Balanced Tunnel
传送门 这一题有点意思 首先预处理出 $pos[x]$ 表示编号 $x$ 的车是第几个出隧道的 然后按进入隧道的顺序枚举每辆车 $x$ 考虑有哪些车比 $x$ 晚进入隧道却比 $x$ 早出隧道 显然是 ...
- Codeforces Good Bye 2015 C. New Year and Domino 前缀和
C. New Year and Domino 题目连接: http://www.codeforces.com/contest/611/problem/C Description They say &q ...
随机推荐
- mysql: error while loading shared libraries: libnuma.so
安装mysql后,执行初始化配置脚本,创建系统自带的数据库和表时报异常: [root@VM_0_12_centos mysql]# scripts/mysql_install_db --basedir ...
- $\LaTeX$数学公式大全
$\LaTeX$数学公式大全$1\ Geek\ and\ Hebrew\ letters$ $\LaTeX$数学公式大全$2\ Math\ Constructs$ $\LaTeX$数学公式大全$3\ ...
- Android webView加载图片显示过大的问题
webview的基本使用流程这里我就不重复说明了,本篇针对的是文章详情加载完成后出现的情况,这里我们使用的方法是:通过js脚本,重置img标签中图片的宽度和高度. 使用步骤: 1.此方法需要使用js, ...
- postgresql获取表最后更新时间(通过触发器将时间写入另外一张表)
通过触发器方式获取表最后更新时间,并将时间信息写入到另外一张表 一.创建测试表和表记录更新时间表 CREATE TABLE weather( city varchar(80), temp_lo int ...
- Qt控制台输出QString
有时候想在控制台输出我们想要的QString变量. 1.qDebug可以实现在控制台终端打印,但我们还是想使用C++中的std::cout<<variable This function ...
- js实现动态显示时间
思路: *得到当前时间 var date = new Date(); //格式化为本地时间 var d1 = date.toLocaleString(); *使页面每秒显示一次时间 setInterv ...
- php支持多个地址跨域访问
//跨域访问的时候才会存在此字段 $origin = isset($_SERVER['HTTP_ORIGIN']) ? $_SERVER['HTTP_ORIGIN'] : ''; $allow_ori ...
- Trie代码学习
感觉不把这个Trie理解一下,AC自动机的代码看起来有点费劲. 这里代码的学习仿照训练指南209页. 这里如果只是查询单词,感觉用map更好,但是如果查前缀,还是用Trie. 1.Trie查询前缀字符 ...
- gdb调试caffe工程
环境:UBUNTU 16.04 CMake caffe 1.0.0-rc3 1.首先编译caffe为debug版本: 在caffe根目录下编辑CMakeFileList.txt,加入如下几行: SET ...
- curl指令的坑
今天使用curl指令构造一个docker api访问,一直得不到预期的结果.调试了半天,发现是网址没加引号. token=$(curl -v -XGET -H >& 由于网址跟了一串参数 ...