[CQOI2011]放棋子 题解(dp+组合数学)
Description
Input
Output
输出仅一行,即方案总数除以 1,000,000,009的余数。
Sample Input
3 1
Sample Output
$Solution$
20%:爆搜,没甚么技术含量虽然我考场上还是没打对只骗到10分Orz
100%:
考虑dp
设$f[i][j][k]$为前k种颜色的棋子占任意i行j列的方案数
那么这个值肯定是前面一系列值的$\sum$
显然需要枚举两层$0<=l<i\ ,\ 0<=r<j$
之后就可以得到$f[l][r][k-1]$并将其累加
但因为我们设的状态是任意行列
需要在剩下的$n-l$行中选$i-l$行,列的话同理
所以要$*C_{n-l}^{i-l}*C_{m-r}^{j-r}$,
而且如果要转移过去还必须乘上某一种颜色占任意i行j列的方案数
这时设$g[i][j][k]$表示k枚同色棋子占任意i行j列的方案数
可得:
$f[i][j][k] = \sum _ {l = 0} ^ {i - 1} \sum _ {r = 0} ^ {j - 1} f[l][r][k - 1] * g[i - l][j - r][a[k]] * C_{n - l} ^ {i - l} * C_{m - r} ^ {j - r}$
正向求g比较困难,我们可以逆向思维,用所有方案数-不合法方案数之和
$g[i][j][k] = C_{i j} ^ {k} - \sum _ {l = 1} ^ {i} \sum _ {r = 1} ^ {j} g[l][r][k] * C_{i}^{l} * C_{j} ^ {r}$
最后统计$ans=\sum _ {i = 1} ^ {n} \sum _ {j = 1} ^ {m} f[i][j][c]$
收获:如果觉得状态设计得当,而缺少转移方程的某一部分时,不妨设一个辅助数组单独考虑。
#include<cstdio>
#include<iostream>
using namespace std;
typedef long long ll;
int n,m,c,a[];
const ll mod=1e9+;
ll f[][][],g[][][],ans=,C[][];
int main()
{
scanf("%d%d%d",&n,&m,&c);
for(int i=;i<=c;i++)
scanf("%d",&a[i]);
if(c>min(n,m))
{
puts("");
return ;
}
f[][][]=;C[][]=;
for(int i=;i<=n*m;i++)
{
C[i][]=;
for(int j=;j<=i;j++)
C[i][j]=(C[i-][j]+C[i-][j-])%mod;
}
for(int k=;k<=c;k++)
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
{
if(a[k]>i*j)continue;
ll res=;
g[i][j][a[k]]=C[i*j][a[k]];
for(int l=;l<=i;l++)
for(int r=;r<=j;r++)
if(l<i||r<j)
(res+=C[i][l]*C[j][r]%mod*g[l][r][a[k]]%mod)%=mod;
g[i][j][a[k]]=(g[i][j][a[k]]-res+mod)%mod;
}
for(int k=;k<=c;k++)
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
for(int l=;l<i;l++)
for(int r=;r<j;r++)
(f[i][j][k]+=C[n-l][i-l]*C[m-r][j-r]%mod*f[l][r][k-]%mod*g[i-l][j-r][a[k]]%mod)%=mod;
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
(ans+=f[i][j][c])%=mod;
cout<<ans<<endl;
return ;
}
[CQOI2011]放棋子 题解(dp+组合数学)的更多相关文章
- 【BZOJ 3294】 3294: [Cqoi2011]放棋子 (DP+组合数学+容斥原理)
3294: [Cqoi2011]放棋子 Description Input 输入第一行为两个整数n, m, c,即行数.列数和棋子的颜色数.第二行包含c个正整数,即每个颜色的棋子数.所有颜色的棋子总数 ...
- [CQOI2011]放棋子 (DP,数论)
[CQOI2011]放棋子 \(solution:\) 看到这道题我们首先就应该想到有可能是DP和数论,因为题目已经很有特性了(首先题面是放棋子)(然后这一题方案数很多要取模)(而且这一题的数据范围很 ...
- P3158 [CQOI2011]放棋子(dp+组合数)
P3158 [CQOI2011]放棋子 放棋子的顺序和方案数无关,所以可以从按颜色递推 设$f[u][p][k]$为放到第$u$种颜色,所剩空间$p*k$的方案数 $g[u][i][j]$表示第$u$ ...
- BZOJ 3294: [Cqoi2011]放棋子(计数dp)
传送门 解题思路 设\(f[i][j][k]\)表示前\(k\)个颜色的棋子占领了\(i\)行\(j\)列的方案数,那么转移时可以枚举上一个颜色时占领的位置,\(f[i][j][k]=\sum\lim ...
- bzoj3294[Cqoi2011]放棋子 dp+组合+容斥
3294: [Cqoi2011]放棋子 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 755 Solved: 294[Submit][Status] ...
- [洛谷P3158] [CQOI2011]放棋子
洛谷题目链接:[CQOI2011]放棋子 题目描述 在一个m行n列的棋盘里放一些彩色的棋子,使得每个格子最多放一个棋子,且不同 颜色的棋子不能在同一行或者同一列.有多少祌方法?例如,n=m=3,有两个 ...
- BZOJ 3294: [Cqoi2011]放棋子
3294: [Cqoi2011]放棋子 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 628 Solved: 238[Submit][Status] ...
- bzoj千题计划261:bzoj3294: [Cqoi2011]放棋子
http://www.lydsy.com/JudgeOnline/problem.php?id=3294 如果一个颜色的棋子放在了第i行第j列,那这种颜色就会占据第i行第j列,其他颜色不能往这儿放 设 ...
- [CQOI2011]放棋子--DP
题目描述: 输入格式 输入第一行为两个整数n, m, c,即行数.列数和棋子的颜色数.第二行包含c个正整数,即每个颜色的棋子数.所有颜色的棋子总数保证不超过nm.N,M<=30 C<=10 ...
随机推荐
- 继续我们的学习。这次鸟哥讲的是LVM。。。磁盘管理 最后链接文章没有看
LVM...让我理解就是一个将好多分区磁盘帮到一起的玩意,类似于烙大饼...然后再切 新建了一个虚拟机,然后又挂了一个5G的硬盘,然后分出了5块空间,挂载到了虚拟机上.这些步骤很简单 fdisk ...
- Excel表格 函数
1.四则运算(加.减.乘.除).求和.平均.计数.最值. 2. 逻辑函数 (IF函数.NOT函数.等) 3.时间和日期 ( NOW 返回当前日期和时间.等) 4.数学与三角函数 5.文本 ( LOWE ...
- Use Git Credential Managers to Authenticate to Azure Repos
https://docs.microsoft.com/en-us/azure/devops/repos/git/set-up-credential-managers?view=azure-devops ...
- php环境搭建以及优化
WampServer 配置伪静态 httpd.conf文件 搜索找到“LoadModule rewrite_module modules/mod_rewrite.so”这一行,去掉前面的“#”: 搜索 ...
- kafka ProducerConfig 配置
kafka-clients : 1.0.1
- 73、salesforce通过JAVA来Call在salesforce中已经写好的Restful处理接口
/** *使用salesforce通过REST方式作为webservice,需要以下几点 *1.类和方法需要global,方法需要静态 *2.类需要通过RestResource(UrlMapping= ...
- python 使用yaml模块
python:yaml模块一.yaml文件介绍YAML是一种简洁的非标记语言.其以数据为中心,使用空白,缩进,分行组织数据,从而使得表示更加简洁.1. yaml文件规则基本规则: 大小写敏感 ...
- spring mvc 使用freemark配置
直接在配置中加入 <mvc:annotation-driven/> <!--freemark--> <bean id="freemarkerConfig&quo ...
- 高水线 High water mark(HWM)
所有的Oracle表都有一个容纳数据的上限(很像一个水库历史最高的水位),我们把这个上限称为“High water mark"或HWM.这个HWM是一个标记(专门有一个数据块来记录高水标记等 ...
- python调用tushare的pro_bar通用行情接口
接口名称:pro_bar 更新时间:股票和指数通常在15点-17点之间,数字货币实时更新,具体请参考各接口文档明细. 描述:目前整合了股票(未复权.前复权.后复权).指数.数字货币.ETF基金.期货. ...