题意:有一个 $n * n$ 的网格,其中 $m$ 个格子上涂了色。每次随机选择一个格子涂色,允许重复涂,求让网格每一行每一列都至少有一个格子涂了色的操作次数期望。
题解:,,这种一般都要倒推才行。
设$f[i][j]$表示还有$i$行,$j$列未满足的情况下的期望次数。
因为每次选择都是完全随机,不受其他东西的影响。
所以对于题中给出的$m$,实际上就是告诉了我们要求什么东西,假设在已经有那$m$个涂色方块的情况下,我们还有$t1$行,$t2$列未满足,那么我们要求的就是$f[t1][t2]$.
那么我们可以列出转移方程:(一行写不下,分2行写)
$$f[i][j] = 1 + \frac{ij}{n ^ 2} f[i - 1][j - 1] + \frac{(n - i)j}{n ^ 2} f[i][j - 1] $$
$$f[i][j] += \frac{i(n - j)}{n ^ 2} f[i - 1][j] + \frac{(n - i)(n - j)}{n^2} f[i][j]$$
$1$是每次选择的代价,后面的就是进入每一种状态的概率,对于任意后继状态,它对当前状态的贡献就是它的期望 * 进入这个状态的概率(全期望公式)

然后移项化简,对于后面这一堆东西提出一个$\frac{1}{n ^ 2}$,然后把$f[i][j]$放到等式左边,这样就只需要在最后面除一次,可以降低一点精度误差?
$$f[i][j] = \frac{n ^ 2 + ijf[i - 1][j - 1] + (n - i)jf[i][j - 1] + i(n - j)f[i - 1][j]}{[n ^ 2 - (n - i) (n - j)]}$$
然后因为计算的时候,可能会出现为满足行或列的数量为0的情况,这种时候还放在一起计算就不太方便了(需要特判),因此考虑把这些情况单独拿出来看。
那么因为$f[i][0]$这个状态之受行的影响,所以可以看做一个优惠券收集问题,即:
$$f[i][0] = \sum_{j = 1}^{i} \frac{n}{j}$$
那么由于$$f[i - 1][0] = \sum_{j = 1}^{i - 1} \frac{n}{j}$$
可以得到关于$f[i][0]$和$f[i - 1][0]$的一个递推式,即:
$$f[i][0] = f[i - 1][0] + \frac{n}{i}$$
解释一下那个和式:如果我们现在还有$i$行未满足,共$n$行,那么选一次可以导致一行新的被满足的概率就是$\frac{i}{n}$,那么期望就为$\frac{n}{i}$.
这个东西的感性理解大概是:如果一个事件发生的概率是$\frac{1}{5}$,那么显然期望$5$天这个事件就会发生,所以期望是概率的倒数。。。
不过这个也是可以证明的,只只要将计算式列出,做一个错位相减,利用一下极限的思想,最后可以算出期望确实是概率的倒数。

 #include<bits/stdc++.h>
using namespace std;
#define R register int
#define AC 2200
#define db double int n, m, l, r;
double f[AC][AC];
bool zl[AC], zr[AC];//表示第i行or第i列有没有被标记 inline int read()
{
int x = ;char c = getchar();
while(c > '' || c < '') c = getchar();
while(c >= '' && c <= '') x = x * + c - '', c = getchar();
return x;
} void pre()
{
n = l = r = read(), m = read();
for(R i = ; i <= m; i ++)
{
int x = read(), y = read();
if(!zl[x]) zl[x] = true, -- l;
if(!zr[y]) zr[y] = true, -- r;
}
} double cal(double x, double y){
return x / y;
} void work()
{//因为把i or j为0的状态放在下面一起枚举不太方便(要特判),所以在前面单独求
for(R i = ; i <= n; i ++)//
{
f[i][] = f[i - ][] + (double)n / i;//因为这个时候只有行的影响,所以只需要考虑行,那么就相当于一个购物券收集问题
f[][i] = f[][i - ] + (double)n / i;
}
for(R i = ; i <= l; i ++)
for(R j = ; j <= r; j ++)
{
/*f[i][j] = 1;
f[i][j] += cal(i * j, n * n) * f[i - 1][j - 1];
f[i][j] += cal((n - i) * j, n * n) * f[i][j - 1];
f[i][j] += cal(i * (n - j), n * n) * f[i - 1][j];
f[i][j] /= 1 - cal((n - i) * (n - j), n * n);*/
f[i][j] = n * n;//把除法放在最后以降低精度误差
f[i][j] += i * j * f[i - ][j - ];
f[i][j] += (n - i) * j * f[i][j - ];
f[i][j] += i * (n - j) * f[i - ][j];
f[i][j] /= n * n - (n - i) * (n - j);
}
printf("%.10lf\n", f[l][r]);
} int main()
{
freopen("in.in", "r", stdin);
pre();
work();
fclose(stdin);
return ;
}

CF398B Painting The Wall 概率期望的更多相关文章

  1. 【CF398B】B. Painting The Wall(期望)

    B. Painting The Wall time limit per test 1 second memory limit per test 256 megabytes input standard ...

  2. Codeforces Round #233 (Div. 2)D. Painting The Wall 概率DP

                                                                                   D. Painting The Wall ...

  3. [Codefoeces398B]Painting The Wall(概率DP)

    题目大意:一个$n\times n$的棋盘,其中有$m$个格子已经被染色,执行一次染色操作(无论选择的格子是否已被染色)消耗一个单位时间,染色时选中每个格子的概率均等,求使每一行.每一列都存在被染色的 ...

  4. Painting The Wall 期望DP Codeforces 398_B

    B. Painting The Wall time limit per test 1 second memory limit per test 256 megabytes input standard ...

  5. 【BZOJ-1419】Red is good 概率期望DP

    1419: Red is good Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 660  Solved: 257[Submit][Status][Di ...

  6. uvalive 7331 Hovering Hornet 半平面交+概率期望

    题意:一个骰子在一个人正方形内,蜜蜂在任意一个位置可以出现,问看到点数的期望. 思路:半平面交+概率期望 #include<cstdio> #include<cstring> ...

  7. OI队内测试一【数论概率期望】

    版权声明:未经本人允许,擅自转载,一旦发现将严肃处理,情节严重者,将追究法律责任! 序:代码部分待更[因为在家写博客,代码保存在机房] 测试分数:110 本应分数:160 改完分数:200 T1: 题 ...

  8. 2016 多校联赛7 Balls and Boxes(概率期望)

    Mr. Chopsticks is interested in random phenomena, and he conducts an experiment to study randomness. ...

  9. 牛客网多校赛第9场 E-Music Game【概率期望】【逆元】

    链接:https://www.nowcoder.com/acm/contest/147/E 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524 ...

随机推荐

  1. Mybaits: MyBaits的xml文件中大于号和小于号的转义

    < 小于号  <     > 大于号  & & 和 & ' 单引号 &apos; " 双引号  "

  2. (二)SpringBoot2.0基础篇- 静态资源的访问及Thymeleaf模板引擎的使用

    一.描述 在应用系统开发的过程中,不可避免的需要使用静态资源(浏览器看的懂,他可以有变量,例:HTML页面,css样式文件,文本,属性文件,图片等): 并且SpringBoot内置了Thymeleaf ...

  3. 高能福利 |"荐"者有份,有"福"同享

    WeTest 导读 越来越多的开发者加入WeTest大家庭了,感谢大家一直以来的支持,WeTest又有新一步“大福利”赠送了! 即日起,参加WeTest用户推荐有礼活动,推荐者和被推荐者皆可获得福利. ...

  4. CentOS 6.5关闭防火墙

    关闭命令:  service iptables stop 永久关闭防火墙:chkconfig iptables off 两个命令同时运行,运行完成后查看防火墙关闭状态 service iptables ...

  5. 模拟IDE上的run过程

    看了一下老陈写的模仿JDK动态代理,从中取一部分单独扩展,模拟一下IDE上的run过程(不愧是老陈,去年写的东西我要现在才能理解) 对run过程的猜想 在点击run的过程中应该做了不少事.先编译运行r ...

  6. 如何在DCS管理控制台将两个Redis主备实例建立全球灾备。

    华为云分布式缓存服务DCS,具有强大的功能,现在小编教大家如何在DCS管理控制台将两个Redis主备实例建立全球灾备. 建立全球灾备,会对主实例和备实例进行升级,实例进程会重启,连接会中断.同时备实例 ...

  7. 观察者模式——Java实例

    一.定义 观察者模式(有时又被称为模型-视图(View)模式.源-收听者(Listener)模式或从属者模式)是软件设计模式的一种.观察者模式定义了一种一对多的依赖关系,让多个观察者对象同时监听某一个 ...

  8. Python:默认参数

    Python是个人最喜欢的语言,刚开始接触Python时,总觉得有很多槽点,不太喜欢.后来,不知不觉中,就用的多了.习惯了.喜欢上了.Python的功能真的很强大,自己当初学习这门语言的时候,也记录过 ...

  9. 哈夫曼(Huffman)树+哈夫曼编码

    前天acm实验课,老师教了几种排序,抓的一套题上有一个哈夫曼树的题,正好之前离散数学也讲过哈夫曼树,这里我就结合课本,整理一篇关于哈夫曼树的博客. 主要摘自https://www.cnblogs.co ...

  10. 《Linux内核与分析》第六周

    20135130王川东 1.操作系统的三大管理功能包括:进程管理,内存管理,文件系统. 2. Linux内核通过唯一的进程标识PID来区别每个进程.为了管理进程,内核必须对每个进程进行清晰的描述,进程 ...