题意

一个\(n \times m\)的矩形空间,起点是\((1,1)\),终点是\((n,m)\)。

假设当前位于\((x,y)\):

  • 如果当前位于最后一行,那么下一步只能走向\((x,y+1)\)
  • 如果当前位于最后一列,那么下一步只能走向\((x+1,y)\)
  • 否则,以相等的概率走向\((x,y+1)\)和\((x+1,y)\)中的一个。

矩形空间中有一个小的矩形黑洞,用左上角和右下角的坐标表示,走进黑洞视为游戏失败,走到\((n,m)\)视为游戏成功,问游戏成功的概率。

解题思路

先不考虑黑洞。

对于\(1 \leq x < n, 1 \leq y < m\),走到\((x,y)\)一共有\(C_{x+y-2}^{x-1}\)种可能的路径,然后走每一条可能的路径的概率为\(\frac{1}{2^{x+y-2}}\),所以走到\((x,y)\)的概率是\(\frac{C_{x+y-2}^{x-1}}{2^{x+y-2}}\)。

但是\(n,m\)的取值最高可以到\(1e5\),如果不加处理会炸精度,如果使用\(\log\)处理就可以把数值压在能接受的范围内。再加上

\[2^{\log(\frac{C_{x+y-2}^{x-1}}{2^{x+y-2}})} = 2 ^ {\log((x + y - 2)!)-\log((y-1)!) - -\log((x-1)!) - x -y +2}
\]

\(O(n)\)预处理出\(\log(i!)\)后就可以方便的计算\(P(x,y)\)了。

然后根据规则,最后一行和最后一列的概率要另外算。这里仅以最后一行为例,最后一列也是用同样的方法。

记走到\((x,y)\)的概率为\(P(x,y)\),那么如果\((x,y)\)位于最后一行,则\(P(x,y)=P(x,y-1)+\frac{1}{2}P(x-1,y)\)。然后\(P(1,n)\)很容易推出等于\(\frac{1}{2^{n-1}}\),所以这一行的概率就可以\(O(n)\)的递推出来。

假设上图中,黑色部分为黑洞,\((1,1)\)位于坐上角,那么很容易就可以得出游戏成功的概率为走到灰色格子的概率之和。将之前推导出的公式代入即可。

总结

早上七点的比赛差点错过了。前三题几乎都是直接秒,但是卡在了第四题。公式都推对了,但是没有想到炸精度怎么处理。脑海中出现了之前llg想在学校月赛搞用log处理大数的出题思路,但是看精度一直觉得不太行,然后就下班了,哪想到这就是正解。果然有时候就应该xjbg。

然后补题的时候又是因为担心精度问题用了long double,然后各种TLE,换成double就过了。

AC代码

#include<bits/stdc++.h>
using namespace std;
const int N = 2e5 + 5; double lg[N], lstr[N], lstc[N];
int main()
{
#ifdef BACKLIGHT
freopen("in.txt", "r", stdin);
#endif lg[0] = 0;
for (int i =1; i <= 2e5; ++i) lg[i] = lg[i-1] + log2(i); int T;
scanf("%d", &T);
for (int Case = 1; Case <= T; Case ++) {
int n, m, l, r, u, d;
scanf("%d %d %d %d %d %d", &n, &m, &l, &u, &r, &d); lstr[1] = pow(2, lg[1 + n - 2] - lg[1 - 1] - lg[n - 1] - 1 - n + 2);
for (int i = 2; i <= m; i++) {
lstr[i] = lstr[i-1] + 0.5 * pow(2, lg[i + (n-1)- 2] - lg[i - 1] - lg[(n-1) - 1] - i - (n-1) + 2);
} lstc[1] = pow(2, lg[1 + m - 2] - lg[1 - 1] - lg[m - 1] - 1 - m + 2);
for (int i = 2; i <= n; i++) {
lstc[i] = lstc[i-1] + 0.5 * pow(2, lg[i + (m-1)- 2] - lg[i - 1] - lg[(m-1) - 1] - i - (m-1) + 2);
} double ans = 0, tmp;
for (int i = 1; i <= l - 1; ++i) {
int D = l + d - i;
if(D > m) continue;
else if(D == m) {
ans += lstc[i];
}
else {
tmp = lg[i + D - 2] - lg[i - 1] - lg[D - 1] - i - D + 2;
ans += pow(2, tmp);
}
} for (int i = 1; i <= u - 1; ++i) {
int R = r + u - i;
if(R > n) continue;
else if(R == n) {
ans += lstr[i];
}
else {
tmp = lg[i + R - 2] - lg[i - 1] - lg[R - 1] - i - R + 2;
ans += pow(2, tmp);
}
} printf("Case #%d: %.12lf\n", Case, ans);
}
return 0;
}

Google Kick Start 2020 Round B T4 Wandering Robot的更多相关文章

  1. Google Kick Start 2020 Round B T1-3

    这场题目除了最后一题稍微难了点,其他都是1眼题. T1 Bike Tour 没啥好说的,一个循环解决. T2 Bus Routes 没啥好说的,从第\(n\)的车站开始贪心取最晚的. T3 Robot ...

  2. Google Kick Start 2020 Round C

    ac代码 A. Countdown for循环跑一跑,没啥好说的. B. Stable Wall 如果\(s_{i,j} \ne s_{i+1,j}\),那么说明\(s_{i+1,j}\)必须在\(s ...

  3. Google Kick Start Round G 2019

    Google Kick Start Round G 2019 Book Reading 暴力,没啥好说的 #include<bits/stdc++.h> using namespace s ...

  4. Google Kick Start 2019 C轮 第一题 Wiggle Walk 题解

    Google Kick Start 2019 C轮 第一题 Wiggle Walk 题解 题目地址:https://codingcompetitions.withgoogle.com/kickstar ...

  5. [C++]Store Credit——Google Code Jam Qualification Round Africa 2010

    Google Code Jam Qualification Round Africa 2010 的第一题,很简单. Problem You receive a credit C at a local ...

  6. Google Code Jam 2010 Round 1C Problem A. Rope Intranet

    Google Code Jam 2010 Round 1C Problem A. Rope Intranet https://code.google.com/codejam/contest/61910 ...

  7. kick start 2019 round D T3题解

    ---恢复内容开始--- 题目大意:共有N个房子,每个房子都有各自的坐标X[i],占据每个房子需要一定花费C[i].现在需要选择K个房子作为仓库,1个房子作为商店(与题目不同,概念一样),由于仓库到房 ...

  8. dp - Google Code jam Qualification Round 2015 --- Problem B. Infinite House of Pancakes

    Problem B. Infinite House of Pancakes Problem's Link:   https://code.google.com/codejam/contest/6224 ...

  9. Google Code jam Qualification Round 2015 --- Problem A. Standing Ovation

    Problem A. Standing Ovation Problem's Link:   https://code.google.com/codejam/contest/6224486/dashbo ...

随机推荐

  1. 你不是说你会Aop吗?

    一大早,小王就急匆匆的跑过来找我,说:周哥,那个记录日志的功能我想请教一下. 因为公司某个项目要跟别的平台做对接,我们这边需要给他们提供一套接口.昨天,我就将记录接口日志的工作安排给了小王. 下面是我 ...

  2. Chrome太占内存?试试这个

    " The Great Suspender" 是一个免费的开源 Google Chrome 扩展程序,适用于那些发现chrome占用过多系统资源或经常遭受chrome崩溃的人. 一 ...

  3. CI4框架应用一 - 环境搭建

    CI框架 (codeigniter)算是一个古老的框架了,由于在工作中一直在使用这个框架,还是比较有感情的.我对CI的感觉就是,简单易用,学习曲线平滑,对于新手友好. 目前CI框架已经更新到CI4了, ...

  4. 001_centos7配置网络动态获取IP地址

    笔者今天刚装完centos7的虚拟机,发现无法获取IP地址,经过网上查询资料,发现centos7是默认没有网络配置的,需要手工配置. 而centos7与centos6不同,没有了config命令,所以 ...

  5. 13、Visitor 访问者模式 访问数据结构并处理数据 行为型设计模式

    1.模式的定义与特点 访问者(Visitor)模式的定义:将作用于某种数据结构中的各元素的操作分离出来封装成独立的类,使其在不改变数据结构的前提下可以添加作用于这些元素的新的操作,为数据结构中的每个元 ...

  6. 将vscode打造成强大的C/C++ IDE

    一.安装 你可以直接从微软官网下载,如果你想要一个纯净的vscode(微软官方的有一项商标.一个插件库.一个 C# 调试器以及遥测),可以手动编译https://github.com/microsof ...

  7. JavaScript基础-02

    1. 六种数据类型: string字符串:number数值:boolean布尔值:null空值:undefined 未定义:object对象 基本数据类型(值类型): string字符串:number ...

  8. next()与nextLine()的区别

    abc def ghij kl mno pqr st uvw xyz 你用next(),第一次取的是abc,第二次取的是def,第三次取的是ghij 你用nextLine(),第一次取的是abc de ...

  9. Core + Vue 后台管理基础框架9——统一日志

    1.背景 前阵子有园友留言,提到日志相关的东西,同时,最近圈子里也有提到日志这个东西.一个充分.集中的统一日志平台还是很有必要的,否则系统出问题了只能靠猜或者干瞪眼.何谓充分,日志记录满足最低要求.出 ...

  10. 了解JS压缩图片,这一篇就够了

    前言 公司的移动端业务需要在用户上传图片是由前端压缩图片大小,再上传到服务器,这样可以减少移动端上行流量,减少用户上传等待时长,优化用户体验. 插播一下,本文案例已整理成插件,已上传npm ,可通过 ...