[LnOI2019]加特林轮盘赌(DP,概率期望)
[LnOI2019]加特林轮盘赌(DP,概率期望)
题解:
首先特判掉\(p=0/1\)的情况...
先考虑如果\(k=1\)怎么做到\(n^2\)的时间复杂度
设\(f[i]\)表示有\(i\)个人,\(k=1\)的时候幸存的概率
设\(g[i][j]\)表示\(i\)个人每个人挨一下恰好死\(j\)个人的概率
我们就可以列出转移方程:
\]
- 含义:枚举打了一圈后剩下多少人,乘\(1-p\)是因为要保证自己不死
把含有\(f[i]\)的项全部移到左边,化简得:
\]
于是我们就可以愉快地\(dp\)啦
数组\(g\)存不下来,\(dp\)的时候动态更新就好啦
如果\(k!=1\)怎么办呢?
我们可以把\(1\sim k-1\)的这些人先打一遍,对于每种剩下的人数分别计算答案,然后按照概率加起来就好了
\(double\)运算常数很大,比赛时一直\(TLE\)在\(65\)分,关于这个常数的问题有两种解决方法:
- 可以把所有数全部乘上一个\(2e9\)转化为\(long long\)进行运算,算完在除回来,缺点是精度较低,容易写错
- 直接把\(double\)换成\(longdouble\)就行了,不过稍微慢点
看了看排行榜,好像有\(O(n)\)的做法?
代码:
#include<map>
#include<set>
#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define qmax(x,y) (x=max(x,y))
#define qmin(x,y) (x=min(x,y))
#define mp(x,y) make_pair(x,y)
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
inline int read(){
int ans=0,fh=1;
char ch=getchar();
while(ch<'0'||ch>'9'){
if(ch=='-') fh=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
ans=ans*10+ch-'0',ch=getchar();
return ans*fh;
}
const ll maxn=1e4+100,w=2e9;
int n,k;
ll p,g[maxn],f[maxn],pp[maxn];
double php;
int main(){
// freopen("nh.in","r",stdin);
// freopen("zhy.out","w",stdout);
scanf("%lf%d%d",&php,&n,&k);
if(php==0){
printf(n==1?"1\n":"0\n");
return 0;
}
p=php*w;
int x=k-1;
f[1]=1*w,g[0]=w-p,g[1]=p;
if(x==1)
for(int j=0;j<=n;j++)
pp[j]=g[j];
for(int i=2;i<=n;i++){
int j;
for(j=1;j<i-3;j+=4){
f[i]+=f[j]*g[i-j]/w;
f[i]+=f[j+1]*g[i-j-1]/w;
f[i]+=f[j+2]*g[i-j-2]/w;
f[i]+=f[j+3]*g[i-j-3]/w;
}
for(;j<i;j++) f[i]+=f[j]*g[i-j]/w;
(f[i]*=w-p)/=w;
for(j=i;j>=4;j-=4){
g[j]=(g[j-1]*p+g[j]*(w-p))/w;
g[j-1]=(g[j-2]*p+g[j-1]*(w-p))/w;
g[j-2]=(g[j-3]*p+g[j-2]*(w-p))/w;
g[j-3]=(g[j-4]*p+g[j-3]*(w-p))/w;
}
for(;j;j--) g[j]=(g[j-1]*p+g[j]*(w-p))/w;
(g[0]*=w-p)/=w;
if(i==x)
for(int j=0;j<=n;j++)
pp[j]=g[j];
(f[i]*=w)/=w-g[0];
}
if(x==0){
printf("%.12lf\n",f[n]*1.0/w);
return 0;
}
double ans=0;
for(int i=0;i<=x;i++){
ans+=pp[i]*f[n-i]/w;
}
printf("%.12lf\n",ans*1.0/w);
return 0;
}
[LnOI2019]加特林轮盘赌(DP,概率期望)的更多相关文章
- 洛谷 P5249 [LnOI2019]加特林轮盘赌 题解【概率期望】【DP】
很有意思的题目. 题目背景 加特林轮盘赌是一个养生游戏. 题目描述 与俄罗斯轮盘赌等手枪的赌博不同的是,加特林轮盘赌的赌具是加特林. 加特林轮盘赌的规则很简单:在加特林的部分弹夹中填充子弹.游戏的参加 ...
- 洛谷 P5249 - [LnOI2019]加特林轮盘赌(期望 dp+高斯消元)
题面传送门 期望真 nm 有意思,所以蒟蒻又来颓期望辣 先特判掉 \(P_0=0\) 的情况,下面假设 \(P_0\ne 0\). 首先注意到我们每次将加特林对准一个人,如果这个人被毙掉了,那么相当于 ...
- Codeforces1097D. Makoto and a Blackboard(数论+dp+概率期望)
题目链接:传送门 题目大意: 给出一个整数n写在黑板上,每次操作会将黑板上的数(初始值为n)等概率随机替换成它的因子. 问k次操作之后,留在黑板上的数的期望. 要求结果对109+7取模,若结果不是整数 ...
- hdu-5781 ATM Mechine(dp+概率期望)
题目链接: ATM Mechine Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Other ...
- BZOJ3566:[SHOI2014]概率充电器(树形DP,概率期望)
Description 著名的电子产品品牌 SHOI 刚刚发布了引领世界潮流的下一代电子产品——概率充电器: “采用全新纳米级加工技术,实现元件与导线能否通电完全由真随机数决定!SHOI 概率充电器, ...
- hdu-5816 Hearthstone(状压dp+概率期望)
题目链接: Hearthstone Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Other ...
- BZOJ4008:[HNOI2015]亚瑟王(DP,概率期望)
Description 小 K 不慎被 LL 邪教洗脑了,洗脑程度深到他甚至想要从亚瑟王邪教中脱坑. 他决定,在脱坑之前,最后再来打一盘亚瑟王.既然是最后一战,就一定要打得漂亮.众所周知,亚瑟王是一个 ...
- [LnOI2019]加特林轮盘赌
Luogu5249 轮流开枪打一个环上的人 , 每次\(p\)的概率打死 , \(p\)始终相同 , 从第\(1\)个人开始 , 求第\(k\)个人成为唯一幸存者的概率 \(19.3.30\) 官方题 ...
- 【BZOJ-1419】Red is good 概率期望DP
1419: Red is good Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 660 Solved: 257[Submit][Status][Di ...
随机推荐
- Integer自动装拆箱
public static void main(String[] args) { Integer a1 = 1; Integer a2 = 1; Integer b1 = 127; Integer b ...
- [HAOI2008]排名系统& [Zjoi2006]GameZ游戏排名系统
1056: [HAOI2008]排名系统 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2487 Solved: 711[Submit][Statu ...
- JDK源码分析之concurrent包(二) -- 线程池ThreadPoolExecutor
上一篇我们简单描述了Executor框架的结构,本篇正式开始并发包中部分源码的解读. 我们知道,目前主流的商用虚拟机在线程的实现上可能会有所差别.但不管如何实现,在开启和关闭线程时一定会耗费很多CPU ...
- IE11 Enterprise Mode Template missing from GPMC
IE11 Enterprise Mode Template missing from GPMC Reason:You have not copied the new IE11 Enterpri ...
- C# WinForm 中进行UrlEncode
public static string ToUrlEncode(string strCode) { StringBuilder sb = new StringBuilder(); byte[] by ...
- Java基础 - 面向对象 - 构造方法
在类中除了成员方法之外,还存在一种特殊类型的方法,那就是构造方法.构造方法是一个与类同名的方法,对象的创建就是通过构造方法完成的.每当类实例化一个对象时,类都会自动调用构造方法. 构造方法的特点: 构 ...
- 转!!java中File的delete()方法删除文件失败的原因
一般来说 java file.delete失败 有以下几个原因 1.看看是否被别的进程引用,手工删除试试(删除不了就是被别的进程占用)2.file是文件夹 并且不为空,有别的文件夹或文件, 3.极有可 ...
- 阻塞IO 非阻塞IO 异步IO
阻塞IO 一般表现为 进程/线程 调用IO操作后就一直死循环等待,直至IO操作结束,返回IO结果 非阻塞IO 一般表现为 进程/线程 调用IO操作后,可以先去干别的事情,但是每隔一段时间,回去询问一下 ...
- git学习------>如何用git log命令来查看某个指定文件的提交历史记录
有时候接手一份新代码时,看到某些文件的改动,但不清楚这个改动的作者和原因,想查看该文件的具体提交历史记录. 今天一个同事是这样做的,直接敲git log命令,然后再使用vim命令的搜索关键字的方法来查 ...
- sql server性能调优
转自:https://www.cnblogs.com/woodytu/tag/%E6%80%A7%E8%83%BD%E8%B0%83%E4%BC%98%E5%9F%B9%E8%AE%AD/defaul ...