题目

我还没疯

发现如果我们将血量抽象成点,一轮操作抽象成图上的一条边,我们如果能求出每一条边的概率,我们就能搞一下这道题

假设我们求出了这个图\(E\),设\(dp_i\)表示从\(i\)点到达\(0\)点的期望路径长度

那么我们可以列出如下的方程

\[dp_u=\sum_{(u,v)\in E}P(u,v)\times(dp_v+1)
\]

发现这个方程可以高斯消元来做

问题变成了如何求出这张图

我们如求出了经过\(k\)次减小的操作,血量\(i\)变成血量\(j\)的概率是多少,我们讨论一下那个增加的操作,就能把这张图求出来了

发现\(k\)很大\(n\)较小,于是觉得可以矩阵优化

优化个鬼啊

设\(dp_{i,j}\)表示经过\(i\)次减小操作血量为\(j\)的概率

显然有

\[dp_{i,j}=\frac{1}{m+1}\times dp_{i-1,j+1}+\frac{m}{m+1}\times dp_{i-1,j}
\]

发现这个柿子确实可以矩阵优化一下,变成\(O(n^3logk)\)

之后套上消元,就有\(40\)分了

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define maxn 205
#define re register
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
const LL mod=1000000007;
inline int read() {
int x=0;char c=getchar();while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
void exgcd(LL a,LL b,LL &x,LL &y) {if(!b) {x=1,y=0;return;}exgcd(b,a%b,y,x);y-=a/b*x;}
inline LL inv(LL t) {LL x,y;exgcd(t,mod,x,y);return (x%mod+mod)%mod;}
int n,p,T;
LL m,k;
LL ans[maxn][maxn],a[maxn][maxn],b[maxn][maxn],E[maxn];
inline void did_ans() {
LL mid[maxn][maxn];
for(re int i=0;i<=n;i++)
for(re int j=0;j<=n;j++) mid[i][j]=ans[i][j],ans[i][j]=0;
for(re int k=0;k<=n;k++)
for(re int i=0;i<=n;i++)
for(re int j=0;j<=n;j++)
ans[i][j]=(ans[i][j]+a[i][k]*mid[k][j]%mod)%mod;
}
inline void did_a() {
LL mid[maxn][maxn];
for(re int i=0;i<=n;i++)
for(re int j=0;j<=n;j++) mid[i][j]=a[i][j],a[i][j]=0;
for(re int k=0;k<=n;k++)
for(re int i=0;i<=n;i++)
for(re int j=0;j<=n;j++)
a[i][j]=(a[i][j]+mid[i][k]*mid[k][j]%mod)%mod;
}
inline void Quick(LL b) {while(b) {if(b&1) did_ans();b>>=1ll;did_a();}}
int main() {
T=read();
while(T--) {
n=read(),p=read(),m=read(),k=read();
memset(b,0,sizeof(b));memset(ans,0,sizeof(ans));
memset(E,0,sizeof(E));memset(a,0,sizeof(a));
if(!p) {puts("0");continue;}
if(!k) {puts("-1");continue;}
if(!m&&k==1) {puts("-1");continue;}
if(!m){
int ans=0;
while(p>0){if(p<n) p++;p-=k;ans++;}
printf("%d\n",ans); continue;
}
LL Inv=inv(m+1ll);
for(re int i=1;i<=n;i++)
a[i][i+1]=Inv,a[i][i]=m*Inv%mod;a[n][n+1]=0;
a[0][0]=1,a[0][1]=Inv;
for(re int i=0;i<=n;i++)
ans[i][i]=1;
Quick(k);
for(re int i=1;i<n;i++) {
for(re int j=1;j<=i;j++)
b[i][j]+=ans[j][i]*Inv%mod*m%mod,b[j][i]%=mod;
for(re int j=1;j<=i+1;j++)
b[i][j]+=ans[j][i+1]*Inv%mod,b[j][i]%=mod;
}
for(re int i=1;i<=n;i++) b[n][i]=ans[i][n];
for(re int i=1;i<=n;i++) b[i][n+1]=mod-1;
for(re int i=1;i<=n;i++) b[i][i]-=1ll,b[i][i]=(b[i][i]+mod)%mod;
for(re int i=1;i<=n;i++) {
LL t=inv(b[i][i]);
for(re int j=n+1;j>=i;--j)
b[i][j]=(b[i][j]*t)%mod;
for(re int j=i+1;j<=n;j++)
for(re int k=n+1;k>=i;--k)
b[j][k]=(b[j][k]-b[j][i]*b[i][k]%mod+mod)%mod;
}
E[n]=b[n][n+1];
for(re int i=n-1;i;--i) {
E[i]=b[i][n+1];
for(re int j=i+1;j<=n;j++)
E[i]=(E[i]-E[j]*b[i][j]%mod+mod)%mod;
}
printf("%lld\n",E[p]);
}
return 0;
}

显然复杂度不对,尤其是矩阵这边

根据高中数学必修三,显然血量减少\(r\)的概率应该是

\[\frac{\binom{k}{r}m^{k-r}}{(m+1)^k}
\]

于是处理一下组合数就不用矩阵了

现在的复杂度就只剩下高斯消元的\(O(n^3)\)了

发现这个矩阵非常特殊,第\(i\)行只有前\(i+1\)列有值,于是我们往下消元的时候一行只需要消三个数就可以了

复杂度\(O(Tn^2)\)

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#define maxn 1505
#define re register
#define LL long long
const LL mod=1000000007;
inline int read() {
int x=0;char c=getchar();while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
int n,p,T;LL m,k;
LL ans[maxn],b[maxn][maxn],E[maxn],C[maxn];
void exgcd(LL a,LL b,LL &x,LL &y) {if(!b) {x=1,y=0;return;}exgcd(b,a%b,y,x);y-=a/b*x;}
inline LL inv(LL t) {LL x,y;exgcd(t,mod,x,y);return (x%mod+mod)%mod;}
inline LL quick(LL a,LL b) {LL S=1;while(b) {if(b&1ll) S=S*a%mod;b>>=1ll;a=a*a%mod;}return S;}
int main() {
T=read();
while(T--) {
n=read(),p=read(),m=read(),k=read();
memset(b,0,sizeof(b));memset(ans,0,sizeof(ans));
memset(E,0,sizeof(E));memset(C,0,sizeof(C));
if(!p) {puts("0");continue;}
if(!k) {puts("-1");continue;}
if(!m&&k==1) {puts("-1");continue;}
if(!m){
int ans=0;
while(p>0){if(p<n) p++;p-=k;ans++;}
printf("%d\n",ans); continue;
}
LL Inv=inv(m+1ll),D=inv(quick(m+1ll,k));
C[0]=1;LL now=k,fac=1,t=k;C[1]=k;
for(re int i=2;i<=n;i++) t--,fac=(fac*i)%mod,now=(now*t)%mod,C[i]=now*inv(fac)%mod;
for(re int i=0;i<=n&&k>=i;i++) ans[i]=C[i]*quick(m,k-i)%mod*D%mod;
for(re int i=1;i<n;i++) {
for(re int j=1;j<=i;j++)
b[i][j]+=ans[i-j]*Inv%mod*m%mod,b[j][i]%=mod;
for(re int j=1;j<=i+1;j++)
b[i][j]+=ans[i+1-j]*Inv%mod,b[j][i]%=mod;
}
for(re int i=1;i<=n;i++) b[n][i]=ans[n-i];
for(re int i=1;i<=n;i++) b[i][n+1]=mod-1;
for(re int i=1;i<=n;i++) b[i][i]-=1ll,b[i][i]=(b[i][i]+mod)%mod;
for(re int i=1;i<=n;i++) {
LL t=inv(b[i][i]);
for(re int j=n+1;j>=i;--j)
b[i][j]=(b[i][j]*t)%mod;
for(re int j=i+1;j<=n;j++) {
b[j][n+1]=(b[j][n+1]-b[j][i]*b[i][n+1]%mod+mod)%mod;
for(re int k=i+1;k>=i;k--)
b[j][k]=(b[j][k]-b[j][i]*b[i][k]%mod+mod)%mod;
}
}
E[n]=b[n][n+1];
for(re int i=n-1;i;--i) {
E[i]=b[i][n+1];
for(re int j=i+1;j<=n;j++)
E[i]=(E[i]-E[j]*b[i][j]%mod+mod)%mod;
}
printf("%lld\n",E[p]);
}
return 0;
}

[BJOI2018]治疗之雨的更多相关文章

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

    [BZOJ5292][BJOI2018]治疗之雨(高斯消元) 题面 BZOJ 洛谷 题解 设\(f[i]\)表示剩余\(i\)点生命时的期望死亡的次数. 考虑打\(k\)次下来脸上被打了\(i\)下的 ...

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

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

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

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

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

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

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

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

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

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

  7. 题解 「BJOI2018 治疗之雨」

    题目传送门 题目大意 有一个初始为 \(p\) 的数,每次操作分为以下两个: 有 \(\frac{1}{m+1}\) 的概率$+1,但是中途 \(p\) 的最大值只能为 \(n\)$ 有 \(k\) ...

  8. 【LOJ】#2513. 「BJOI2018」治疗之雨

    题解 具体就是列一个未知数方程\(dp[i]\)表示有\(i\)滴血的时候期望多少轮 \(dp[i] = 1 + \sum_{j = 1}^{i + 1} a_{i,j}dp[j]\) \(dp[n] ...

  9. 「BJOI2018」治疗之雨

    传送门 Description 有\(m+1\)个数,第一个数为\(p\),每轮:选一个数\(+1\),再依次选\(k\)个数\(-1\) 要求如果第一个数\(=N\),不能选它\(+1\),如果第一 ...

随机推荐

  1. solr 索引库的维护

    一.配置中文分析器:IK-analyzer,在FieldType中指定中文分析器:1 复制IK-analyzer到你的服务器指定目录中.2 在该目录中,我们需要的东西有:IKAnalyzer的jar包 ...

  2. sublime text 3 主题更换

    1.安装colorsublime,里面收藏了大量的主题 2.要选择主题的时候,Ctrl+Shift+P 打开Package Control,输入color,如图 3.Enter,进入选择列表,键盘上下 ...

  3. sublime text 2编辑器中文问题

    Sublime Text 2是一个非常不错的源代码及文本编辑器,但是不支持GB2312和GBK编码在很多情况下会非常麻烦.不过Sublime Package Control所以供的插件可以让Subli ...

  4. Vue指令(四)--v-model

    1.v-model的使用场景 1.v-model的使用,用于表单控件的数据绑定 2.v-model与value共同使用,实现选项框的选中事件,两者相同时,选中 3.v-model 与v-bind:va ...

  5. 很多个java面试题

    1. 为什么说Java是一门平台无关语言? 平台无关实际的含义是“一次编写到处运行”.Java 能够做到是因为它的字节码(byte code)可以运行在任何操作系统上,与底层系统无关. 2. 为什么 ...

  6. springboot项目作为war包运行

    一.首先是pom文件中设置打成war包 < packaging>war< /packaging> 二.然后是修改依赖: <dependency> <group ...

  7. 结构型--代理模式(Proxy)

    一.代理模式是什么? 代理模式属于GOF23设计模式中结构型中的设计模式,通过代理对象来屏蔽(部分或者屏蔽)对真实对象的直接访问,下图为UML图: 在代理模式中组件包括:抽象角色接口.代理角色类.真实 ...

  8. springMVC @Component-@Resource-@Repository-@Service-@Controller的区别和理解

    作用: @Component------------------------泛指组件,当组件不好归类的时候,我们可以使用这个注解进行标注.(Component-------成分; 组分; 零件) @R ...

  9. Codeforces Round #394 (Div. 2)

    前一半是刚刚打完比赛的时候写的……不知为啥手腕有点僵,估计是前一个小时用力过度了吧= = 前四题看着还好,后两题就有点懵逼了……现在还不知道E题的题意到底是啥…… 不管了……还没找着官方题解,贴一下自 ...

  10. Web前端面试指导(十一):样式导入有哪些方式?

    样式导入方式 link import 使用方式 link的使用 <link href="index.css" rel="stylesheet"> i ...