Google Kick Start 2020 Round B T4 Wandering Robot
题意
一个\(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\)处理就可以把数值压在能接受的范围内。再加上
\]
\(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的更多相关文章
- Google Kick Start 2020 Round B T1-3
这场题目除了最后一题稍微难了点,其他都是1眼题. T1 Bike Tour 没啥好说的,一个循环解决. T2 Bus Routes 没啥好说的,从第\(n\)的车站开始贪心取最晚的. T3 Robot ...
- 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 ...
- Google Kick Start Round G 2019
Google Kick Start Round G 2019 Book Reading 暴力,没啥好说的 #include<bits/stdc++.h> using namespace s ...
- Google Kick Start 2019 C轮 第一题 Wiggle Walk 题解
Google Kick Start 2019 C轮 第一题 Wiggle Walk 题解 题目地址:https://codingcompetitions.withgoogle.com/kickstar ...
- [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 ...
- 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 ...
- kick start 2019 round D T3题解
---恢复内容开始--- 题目大意:共有N个房子,每个房子都有各自的坐标X[i],占据每个房子需要一定花费C[i].现在需要选择K个房子作为仓库,1个房子作为商店(与题目不同,概念一样),由于仓库到房 ...
- 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 ...
- 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 ...
随机推荐
- 关于Exceptionless日志收集框架如何关闭磁盘缓存
问题:在使用Exceptionless的时候,Exception在收集到日志时会默认在appdata里面缓存当条日志的缓存文件,时间久了之后,如果收集到的日志越多磁盘的空间就会不足. 我使用的环境是 ...
- three.js 着色器材质之变量(二)
上一篇郭先生在例子中用到了着色器变量中的uniform和varying.这篇继续结合例子将一下attribute变量,在使用过程中也发现由于three.js的版本迭代,之前的一些属性和参数已经发生了改 ...
- UIAutomator环境Android8.0 环境异常解决
个人PC环境 ANDROID_HOME:F:\1Study\Andriod\51zxw_2018-0102\Sdk ANT_HOME:D:\ant\apache-ant-1.10.5\ CLASSPA ...
- Vue + ccropper.js裁切图片(vue-cropper)
按原比例裁剪图片并且不失真. 安装: cnpm install vue-cropper --save-dev 使用: <template> <div style="disp ...
- JDK8 String类知识总结
一.概述 java的String类可以说是日常实用的最多的类,但是大多数时候都只是简单的拼接或者调用API,今天决定深入点了解一下String类. 要第一时间了解一个类,没有什么比官方的javaDoc ...
- Ubuntu无法安装 英伟达显卡
安装Ubuntu无法正常驱动英伟达,这时需要在启动参数中添加nomodset 如果不会添加参数可以参考这篇文章:安装ubuntu时黑屏三种解决办法
- [C#] (原创)一步一步教你自定义控件——02,ScrollBar(滚动条)
一.前言 技术没有先进与落后,只有合适与不合适. 本篇的自定义控件是:滚动条(ScollBar). 我们可以在网上看到很多自定义的滚动条控件,它们大都是使用UserControl去做,即至少使用一个P ...
- Flutter简介
Flutter简介 Flutter 是 Google推出并开源的移动应用开发框架,主打跨平台.高保真.高性能.开发者可以通过 Dart语言开发 App,一套代码同时运行在 iOS 和 Android平 ...
- Dubbo系列之 (四)服务订阅(1)
辅助链接 Dubbo系列之 (一)SPI扩展 Dubbo系列之 (二)Registry注册中心-注册(1) Dubbo系列之 (三)Registry注册中心-注册(2) Dubbo系列之 (四)服务订 ...
- python 10 else EasyGui(转载)
else语句 if else 要么怎么样,要么不怎么样 while else 干完了能怎样,干不完就不怎样 (异常处理) else 没有问题,就干吧 try: int('a') except Valu ...