description

  偶然间,chnlich 发现了他小时候玩过的一个游戏“魂斗罗”,于是决定怀旧。但是这是一个奇怪的魂斗罗 MOD。

  有 N 个关卡,初始有 Q 条命。

  每通过一个关卡,会得到 u 分和1条命,生命上限为 Q。其中 u=min(最近一次连续通过的关数,R)。

  若没有通过这个关卡,将会失去1条命,并进入下一个关卡。

  当没有生命或没有未挑战过的关卡时,游戏结束,得到的分数为每关得到的分数的总和。

  由于 chnlich 好久不玩这个游戏了,每条命通过每个关卡的概率均为p(0<=p<=1),原先 chnlich 的最高分纪录是 S。

  现在 chnlich 想要知道,当 p 至少为多少时,chnlich 期望获得的总分数能够超过原先的最高分。

analysis

  • 很好的二分\(+\)矩乘\(+DP\)题

  • 由题解可得\(50pts\)的\(DP\)方程

\[g[i+1][min(j+1,R)][min(k+1),Q]+=g[i][j][k]*p
\]

\[g[i+1][0][k-1]+=g[i][j][k]*(1-p)
\]

\[Answer+=g[i][j][k]*p*(j+1)
\]

  • 二分一个最小可能的\(p\),用矩乘判断是否可行

  • 对于每一种积分\(x\)命\(y\)条的状态,用不同的数字表示来压状态

  • 可以知道当积分\(≥q\)时(连赢至少\(q\)次),命一定为\(q\)条,这样可以省去很多状态

  • 压完状态之后,就只有\(30\)种左右的状态,可以矩乘快速幂

  • 具体就是把按照方程每种状态转移到其他状态的矩阵位置加上贡献\(p,1-p\)或\(p*(j+1)\)

  • 注意那个\(Answer\)由所有的情况累加到最终状态,即为积分\(0\)命\(0\)条的情况

  • 判断最终矩阵第起始状态行(积分\(0\)命\(q\)条)第终止状态列(积分\(0\)命\(0\)条)的值是否大于\(s\)即可

  • 下次遇到矩乘题目,要认真想想\(DP\),不要被矩阵怎么推给卡住


code

#pragma GCC optimize("O3")
#pragma G++ optimize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define MAXN 40
#define ll long long
#define reg register ll
#define fo(i,a,b) for (reg i=a;i<=b;++i)
#define fd(i,a,b) for (reg i=a;i>=b;--i) using namespace std; ll statu[10][50];
ll n,m,r,q,s; struct matrix
{
double f[MAXN][MAXN];
matrix(){memset(f,0,sizeof(f));}
};
inline ll read()
{
ll x=0,f=1;char ch=getchar();
while (ch<'0' || '9'<ch){if (ch=='-')f=-1;ch=getchar();}
while ('0'<=ch && ch<='9')x=x*10+ch-'0',ch=getchar();
return x*f;
}
inline ll max(ll x,ll y){return x>y?x:y;}
inline ll min(ll x,ll y){return x<y?x:y;}
inline matrix operator*(matrix a,matrix b)
{
matrix c;
fo(i,1,n)fo(j,1,n)fo(k,1,n)c.f[i][j]+=a.f[i][k]*b.f[k][j];
return c;
}
inline matrix pow(matrix x,ll y)
{
matrix z;
fo(i,1,n)z.f[i][i]=1;
if (y==0)return z;
while (y)
{
if (y&1)z=z*x;
x=x*x,y>>=1;
}
return z;
}
inline bool judge(double p)
{
matrix a;n=0;
fo(i,0,q)//statu[i][j]表示剩i条命有j积分的情况
{
ll tmp=i==q?r:min(i-1,r);
fo(j,0,max(tmp,0))statu[i][j]=++n;
}
a.f[1][1]=1;
fo(i,1,q)
{
ll tmp=i==q?r:min(i-1,r);
fo(j,0,max(tmp,0))
{
ll x=statu[i][j],y=statu[i-1][0],z=statu[min(i+1,q)][min(j+1,r)];
if (i>1)a.f[x][y]=1-p;a.f[x][z]=p,a.f[x][statu[0][0]]=p*min(j+1,r);
}
}
matrix b=pow(a,m);
return b.f[statu[q][0]][statu[0][0]]>s;
}
int main()
{
freopen("T2.in","r",stdin);
m=read(),r=read(),q=read(),s=read();
if (!judge(1)){printf("Impossible.\n");return 0;}
ll l=0,r=10000000,tmp=r,mid=(l+r)>>1;
while (l<=r)
{
mid=(l+r)>>1;
if (judge(1.0*mid/tmp))r=mid-1;else l=mid+1;
}
printf("%.6lf\n",1.0*l/tmp);
return 0;
}

【JZOJ2867】Contra的更多相关文章

  1. Python高手之路【六】python基础之字符串格式化

    Python的字符串格式化有两种方式: 百分号方式.format方式 百分号的方式相对来说比较老,而format方式则是比较先进的方式,企图替换古老的方式,目前两者并存.[PEP-3101] This ...

  2. 【原】谈谈对Objective-C中代理模式的误解

    [原]谈谈对Objective-C中代理模式的误解 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 这篇文章主要是对代理模式和委托模式进行了对比,个人认为Objective ...

  3. 【原】FMDB源码阅读(三)

    [原]FMDB源码阅读(三) 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 FMDB比较优秀的地方就在于对多线程的处理.所以这一篇主要是研究FMDB的多线程处理的实现.而 ...

  4. 【原】Android热更新开源项目Tinker源码解析系列之一:Dex热更新

    [原]Android热更新开源项目Tinker源码解析系列之一:Dex热更新 Tinker是微信的第一个开源项目,主要用于安卓应用bug的热修复和功能的迭代. Tinker github地址:http ...

  5. 【调侃】IOC前世今生

    前些天,参与了公司内部小组的一次技术交流,主要是针对<IOC与AOP>,本着学而时习之的态度及积极分享的精神,我就结合一个小故事来初浅地剖析一下我眼中的“IOC前世今生”,以方便初学者能更 ...

  6. Python高手之路【三】python基础之函数

    基本数据类型补充: set 是一个无序且不重复的元素集合 class set(object): """ set() -> new empty set object ...

  7. Python高手之路【一】初识python

    Python简介 1:Python的创始人 Python (英国发音:/ˈpaɪθən/ 美国发音:/ˈpaɪθɑːn/), 是一种解释型.面向对象.动态数据类型的高级程序设计语言,由荷兰人Guido ...

  8. 【开源】简单4步搞定QQ登录,无需什么代码功底【无语言界限】

    说17号发超简单的教程就17号,qq核审通过后就封装了这个,现在放出来~~ 这个是我封装的一个开源项目:https://github.com/dunitian/LoTQQLogin ————————— ...

  9. 【原】FMDB源码阅读(二)

    [原]FMDB源码阅读(二) 本文转载请注明出处 -- polobymulberry-博客园 1. 前言 上一篇只是简单地过了一下FMDB一个简单例子的基本流程,并没有涉及到FMDB的所有方方面面,比 ...

随机推荐

  1. hdu6311 /// 欧拉路径 无向图最小路径覆盖 输出正反路径

    题目大意: 给定n m 为图的点数和边数 接下来m行 u v 为u到v有一条边 要求最少几笔能画完图的所有边 输出每笔画过的路径编号 正数编号正向 负数编号反向 题解:https://www.cnbl ...

  2. IPC$渗透使用

    简介 IPC$(Internet Process Connection)是共享"命名管道"的资源,它是为了让进程间通信而开放的命名管道,通过提供可信任的用户名和口令,连接双方可以建 ...

  3. JS切换

    这个效果就跟换肤效果差不多 需要准备两套或两套以上的css <!DOCTYPE html> <html lang="en"> <head> &l ...

  4. Vue双向数据绑定原理深度解析

    首先,什么是双向数据绑定?Vue是三大MVVM框架之一,数据绑定简单来说,就是当数据发生变化时,相应的视图会进行更新,当视图更新时,数据也会跟着变化. 在分析其原理和代码的时候,大家首先了解如下几个j ...

  5. vt-is-UTF8 - check whether current VT is in UTF8- or byte-mode. 检查当前VT是否处于VTF8模式或是字节模式.

    总览 vt-is-UTF8 [-h|--help] [-V|--version] [-q|--quiet] 描述 vt-is-UTF8 checks whether the current VT is ...

  6. 安装php 在阿里云yum源的环境

    yum -y install httpd mysql mysql-server php php-mysql postgresql postgresql-server php-postgresql ph ...

  7. DEV插件下的控件Grid和Gridlookupedit控件的结合使用

    创建GridlookupEtid控件 设置其对应属性: 设置属性 this.gridLookUpEdit1.Properties.TextEditStyle = DevExpress.XtraEdit ...

  8. echats问题

    echats 横轴显示不下datazoom配置,加入滚动条 实例博客  https://blog.csdn.net/Zheng_xiao_xin/article/details/80882113 常用 ...

  9. CSS3(@media)判断手机横竖屏

    @media all and (orientation : landscape) { h2{color:red;}/*横屏时字体红色*/ } @media all and (orientation : ...

  10. java常用类——比较器

    Comparable和Comparator接口都是为了对类进行比较,众所周知,诸如Integer,double等基本数据类型,java可以对他们进行比较,而对于类的比较,需要人工定义比较用到的字段比较 ...