【BZOJ5292】[BJOI2018]治疗之雨(高斯消元)

题面

BZOJ

洛谷

题解

设\(f[i]\)表示剩余\(i\)点生命时的期望死亡的次数。

考虑打\(k\)次下来脸上被打了\(i\)下的概率:\(\displaystyle \frac{{k\choose i}m^{k-i}}{(m+1)^k}\)。

\(m=0\)时全部打脸上了,直接判掉。

设\(P[i][j]\)表示\(i\)点血量奶完后再被打一轮下来变成\(j\)点血的概率,这个很容易算出来。

那么我们可以列出和\(f[i]\)相关的\(n+1\)个方程。

比如说:

\[\begin{cases}
f[0]=0\\
f[1]=P[1][0]*f[0]+P[1][1]*f[1]+P[1][2]*f[2]+1\\
...\\
f[n]=P[n][0]*f[0]+P[n][1]*f[1]+P[n][2]*f[2]+...+P[n][n]*f[n]+1
\end{cases}\]

考虑怎么解这个玩意,首先肯定可以高斯消元三方解决。

这个矩阵观察后发现类似于一个下三角矩阵,那么每次用第\(i\)行消掉\(i-1\)行,这样子就是一个下三角了,然后从上往下就可以直接求解。

还有一种方法就是移项后不难发现\(f[1]\)只和\(f[2]\)相关,因此可以设\(f[1]=x\),这样子其他所有数都可以通过\(x\)给表示出来,最后带到\(f[n]\)的方程里就可以解出\(x\)。

#include<iostream>
#include<cstdio>
using namespace std;
#define ll long long
#define MAX 1600
#define MOD 1000000007
void add(int &x,int y){x+=y;if(x>=MOD)x-=MOD;}
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
int fpow(int a,int b){int s=1;while(b){if(b&1)s=1ll*s*a%MOD;a=1ll*a*a%MOD;b>>=1;}return s;}
int n,p,m,k;
int P[MAX],a[MAX][MAX];
int A[MAX],B[MAX];
int main()
{
int T=read();
while(T--)
{
n=read();p=read();m=read();k=read();
if(!k){puts("-1");continue;}
if(m==0)
{
int cnt=0;if(k==1){puts("-1");continue;}
while(p>0)++cnt,p=min(n,p+1)-k;
printf("%d\n",cnt);continue;
}
for(int i=0;i<=n;++i)P[i]=0;
int invm1=fpow(m+1,MOD-2);P[0]=1;
for(int i=1;i<n&&i<=k;++i)P[i]=1ll*P[i-1]*fpow(i,MOD-2)%MOD*(k-i+1)%MOD;
for(int i=0;i<n&&i<=k;++i)P[i]=1ll*P[i]*fpow(m,k-i)%MOD*fpow(fpow(m+1,k),MOD-2)%MOD;
if(k>=n){P[n]=1;for(int i=0;i<n;++i)P[n]=(P[n]+MOD-P[i])%MOD;}
for(int i=0;i<=n;++i)
for(int j=0;j<=n;++j)a[i][j]=0;
for(int i=1;i<n;++i)
for(int j=1;j<=i+1;++j)
add(a[i][j],(1ll*invm1*P[i-j+1]+1ll*m*invm1%MOD*P[i-j])%MOD);
for(int i=1;i<=n;++i)add(a[n][i],P[n-i]);
for(int i=0;i<=n;++i)A[i]=B[i]=0;
A[1]=1;B[1]=0;
for(int i=1;i<n;++i)
{
int sa=A[i],sb=(MOD-1+B[i])%MOD;
for(int j=1;j<=i;++j)sa=(sa+MOD-1ll*a[i][j]*A[j]%MOD)%MOD;
for(int j=1;j<=i;++j)sb=(sb+MOD-1ll*a[i][j]*B[j]%MOD)%MOD;
int v=fpow(a[i][i+1],MOD-2);
A[i+1]=1ll*sa*v%MOD;B[i+1]=1ll*sb*v%MOD;
}
int sa=A[n],sb=(MOD-B[n]+1)%MOD;
for(int i=1;i<=n;++i)sa=(sa+MOD-1ll*a[n][i]*A[i]%MOD)%MOD;
for(int i=1;i<=n;++i)sb=(sb+1ll*a[n][i]*B[i])%MOD;
int x=1ll*sb*fpow(sa,MOD-2)%MOD;
int ans=(1ll*A[p]*x+B[p])%MOD;
printf("%d\n",ans);
continue;
}
}

【BZOJ5292】[BJOI2018]治疗之雨(高斯消元)的更多相关文章

  1. [BZOJ5292][BJOI2018]治疗之雨(概率DP+高斯消元)

    https://blog.csdn.net/xyz32768/article/details/83217209 不难找到DP方程与辅助DP方程,发现DP方程具有后效性,于是高斯消元即可. 但朴素消元显 ...

  2. [BZOJ5292] [BJOI2018]治疗之雨

    题目链接 BZOJ:https://lydsy.com/JudgeOnline/problem.php?id=5292 洛谷:https://www.luogu.org/problemnew/show ...

  3. luoguP4457 [BJOI2018]治疗之雨 概率期望 + 高斯消元

    应该是最后一道紫色的概率了....然而颜色啥也代表不了.... 首先看懂题意: 你现在有$p$点体力,你的体力上限为$n$ 在一轮中, 1.如果你的体力没有满,你有$\frac{1}{m + 1}$的 ...

  4. 洛谷P4457/loj#2513 [BJOI2018]治疗之雨(高斯消元+概率期望)

    题面 传送门(loj) 传送门(洛谷) 题解 模拟赛的时候只想出了高斯消元然后死活不知道怎么继续--结果正解居然就是高斯消元卡常? 首先有个比较难受的地方是它一个回合可能不止扣一滴血--我们得算出\( ...

  5. P4457-[BJOI2018]治疗之雨【期望dp,高斯消元】

    正题 题目链接:https://www.luogu.com.cn/problem/P4457 题目大意 开始一个人最大生命值为\(n\),剩余\(hp\)点生命,然后每个时刻如果生命值没有满那么有\( ...

  6. BZOJ5292 & 洛谷4457 & LOJ2513:[BJOI2018]治疗之雨——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=5292 https://www.luogu.org/problemnew/show/P4457 ht ...

  7. 【BZOJ-3143】游走 高斯消元 + 概率期望

    3143: [Hnoi2013]游走 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2264  Solved: 987[Submit][Status] ...

  8. 【BZOJ-3270】博物馆 高斯消元 + 概率期望

    3270: 博物馆 Time Limit: 30 Sec  Memory Limit: 128 MBSubmit: 292  Solved: 158[Submit][Status][Discuss] ...

  9. *POJ 1222 高斯消元

    EXTENDED LIGHTS OUT Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 9612   Accepted: 62 ...

随机推荐

  1. ImageProcessor组件

    ImageProcessor组件 开源免费的.NET图像即时处理的组件ImageProcessor   承接以前的组件系列,这个组件系列旨在介绍.NET相关的组件,让大家可以在项目中有一个更好的选择, ...

  2. RabbitMQ防止消息丢失

    转载请注明出处 0.目录 RabbitMQ-从基础到实战(1)— Hello RabbitMQ RabbitMQ-从基础到实战(3)— 消息的交换 1.简介 RabbitMQ中,消息丢失可以简单的分为 ...

  3. ES5中文分词(IK)

    ElasticSearch5中文分词(IK) ElasticSearch安装 官网:https://www.elastic.co 1.ElasticSearch安装 1.1.下载安装公共密钥 rpm ...

  4. [2017BUAA软工助教]团队建议

    关于团队项目的个人建议 (以下排名不分先后) 一.hotcode5 你们组要做一个"课件-心得"共享平台 目前最大的竞争对手其实不是北航课程中心网站,而是每个系自己的大班群. 热心 ...

  5. use redis instance in docker hub

    redis - Docker Hubhttps://hub.docker.com/_/redis

  6. Python3练习题 035:Project Euler 007:第10001个素数

    import time def f(x): #判断 x 是否为素数,返回bool值 if x == 2: return True elif x <= 1: return False else: ...

  7. 数据处理 array json 格式 转换成 数组形式

    处理这种数据应该使用的方式是 this.cities= res.data.data.cities.sort((a,b)=>{ //排序 进行字母排序 return a.pinyin[0].cha ...

  8. 下拉框插件select2的使用

    它的优点有: 样式还算好看,支持多选,支持索搜 下面来介绍下select2的用法 1.最简单的用法 只需要加载css和js即可使用 <select name="" id=&q ...

  9. 转《service worker在移动端H5项目的应用》

    1. PWA和Service Worker的关系 PWA (Progressive Web Apps) 不是一项技术,也不是一个框架,我们可以把她理解为一种模式,一种通过应用一些技术将 Web App ...

  10. 腾讯机试题 AcWing 603 打怪兽

    题目链接:https://www.acwing.com/problem/content/605/ 题目大意: 略 分析: 用dp[i][j]表示用j元钱能在前i只怪兽上所能贿赂到的最大武力值. 有一种 ...