[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 ...
随机推荐
- POJ 1426 Find The Multiple (dfs??!!)
Description Given a positive integer n, write a program to find out a nonzero multiple m of n whose ...
- [jzoj5840]Miner 题解(欧拉路)
首先考虑第一问.每个联通块的情况是相对独立的,所以可以分别求每个联通块的答案.无向图中存在欧拉路的条件是奇点数为0或2,那么合法方案肯定是tp到一个奇点,通过一条欧拉路到另一个奇点,再tp到另一个奇点 ...
- [NOIP模拟25]题解
A.字符串 Catalan数不能再裸了 #include<cstdio> #include<iostream> #include<cstring> using na ...
- 在 MacBook 上安装 Ubuntu
建立ubuntu 安裝U盤,加入usb port,按住option key 開機,正常安裝. 完成後,會缺少wifi drive sudo apt-get install bcmwl-kernel-s ...
- (转)使用OpenGL ES显示图像
编写:jdneo - 原文:http://developer.android.com/training/graphics/opengl/index.html 转:http://hukai.me/and ...
- Linux随笔 - DNS搭建
域名系统(英文:Domain Name System,缩写:DNS)是因特网的一项服务.它作为将域名和IP地址相互映射的一个分布式数据库,能够使人更方便的访问互联网.DNS 使用TCP和UDP端口53 ...
- [Java Performance] 线程及同步的性能之线程池/ThreadPoolExecutors/ForkJoinPool
线程池和ThreadPoolExecutors 虽然在程序中可以直接使用Thread类型来进行线程操作,但是更多的情况是使用线程池,尤其是在Java EE应用服务器中,一般会使用若干个线程池来处理 ...
- [轉]Linux Data Structures
Table of Contents, Show Frames, No Frames Chapter 15 Linux Data Structures This appendix lists the m ...
- JS浏览器事件循环机制
文章来自我的 github 博客,包括技术输出和学习笔记,欢迎star. 先来明白些概念性内容. 进程.线程 进程是系统分配的独立资源,是 CPU 资源分配的基本单位,进程是由一个或者多个线程组成的. ...
- Spring中AOP的基于xml开发和配置
pom文件: <?xml version="1.0" encoding="UTF-8"?> <project xmlns="http ...