[BJOI2018]治疗之雨
我还没疯
发现如果我们将血量抽象成点,一轮操作抽象成图上的一条边,我们如果能求出每一条边的概率,我们就能搞一下这道题
假设我们求出了这个图\(E\),设\(dp_i\)表示从\(i\)点到达\(0\)点的期望路径长度
那么我们可以列出如下的方程
\]
发现这个方程可以高斯消元来做
问题变成了如何求出这张图
我们如求出了经过\(k\)次减小的操作,血量\(i\)变成血量\(j\)的概率是多少,我们讨论一下那个增加的操作,就能把这张图求出来了
发现\(k\)很大\(n\)较小,于是觉得可以矩阵优化
优化个鬼啊
设\(dp_{i,j}\)表示经过\(i\)次减小操作血量为\(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\)的概率应该是
\]
于是处理一下组合数就不用矩阵了
现在的复杂度就只剩下高斯消元的\(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]治疗之雨的更多相关文章
- 【BZOJ5292】[BJOI2018]治疗之雨(高斯消元)
[BZOJ5292][BJOI2018]治疗之雨(高斯消元) 题面 BZOJ 洛谷 题解 设\(f[i]\)表示剩余\(i\)点生命时的期望死亡的次数. 考虑打\(k\)次下来脸上被打了\(i\)下的 ...
- [BZOJ5292][BJOI2018]治疗之雨(概率DP+高斯消元)
https://blog.csdn.net/xyz32768/article/details/83217209 不难找到DP方程与辅助DP方程,发现DP方程具有后效性,于是高斯消元即可. 但朴素消元显 ...
- luoguP4457 [BJOI2018]治疗之雨 概率期望 + 高斯消元
应该是最后一道紫色的概率了....然而颜色啥也代表不了.... 首先看懂题意: 你现在有$p$点体力,你的体力上限为$n$ 在一轮中, 1.如果你的体力没有满,你有$\frac{1}{m + 1}$的 ...
- [BZOJ5292] [BJOI2018]治疗之雨
题目链接 BZOJ:https://lydsy.com/JudgeOnline/problem.php?id=5292 洛谷:https://www.luogu.org/problemnew/show ...
- BZOJ5292 & 洛谷4457 & LOJ2513:[BJOI2018]治疗之雨——题解
https://www.lydsy.com/JudgeOnline/problem.php?id=5292 https://www.luogu.org/problemnew/show/P4457 ht ...
- 洛谷P4457/loj#2513 [BJOI2018]治疗之雨(高斯消元+概率期望)
题面 传送门(loj) 传送门(洛谷) 题解 模拟赛的时候只想出了高斯消元然后死活不知道怎么继续--结果正解居然就是高斯消元卡常? 首先有个比较难受的地方是它一个回合可能不止扣一滴血--我们得算出\( ...
- 题解 「BJOI2018 治疗之雨」
题目传送门 题目大意 有一个初始为 \(p\) 的数,每次操作分为以下两个: 有 \(\frac{1}{m+1}\) 的概率$+1,但是中途 \(p\) 的最大值只能为 \(n\)$ 有 \(k\) ...
- 【LOJ】#2513. 「BJOI2018」治疗之雨
题解 具体就是列一个未知数方程\(dp[i]\)表示有\(i\)滴血的时候期望多少轮 \(dp[i] = 1 + \sum_{j = 1}^{i + 1} a_{i,j}dp[j]\) \(dp[n] ...
- 「BJOI2018」治疗之雨
传送门 Description 有\(m+1\)个数,第一个数为\(p\),每轮:选一个数\(+1\),再依次选\(k\)个数\(-1\) 要求如果第一个数\(=N\),不能选它\(+1\),如果第一 ...
随机推荐
- bzoj 4540: [Hnoi2016]序列
Description 给定长度为n的序列:a1,a2,-,an,记为a[1:n].类似地,a[l:r](1≤l≤r≤N)是指序列:al,al+1,-,ar- 1,ar.若1≤l≤s≤t≤r≤n,则称 ...
- 深入理解JavaScript系列(48):对象创建模式(下篇)
介绍 本篇主要是介绍创建对象方面的模式的下篇,利用各种技巧可以极大地避免了错误或者可以编写出非常精简的代码. 模式6:函数语法糖 函数语法糖是为一个对象快速添加方法(函数)的扩展,这个主要是利用pro ...
- [转].Net Core Web应用发布至IIS后报“An error occurred while starting the application”错误
本文转自:http://www.cnblogs.com/TomGui/p/6438686.html An error occurred while starting the application. ...
- js扩展
http://www.css88.com/doc/underscore/#findWhere
- JavaScript - 表单提交前预览图片
其实这东西网上到处都是,但并不完整. 正好我也遇到了这个问题,不仅仅是预览,还需要得到图片的属性. 于是东凑西凑整理出一个完整的版本,并根据个人的理解加上了一点点说明. 首先做一些准备工作,HTML方 ...
- Backbone之温故而知新1-MVC
在忙碌了一段时间之后,又有了空余时间来学习新的东西,自从上次研究了backbone之后,一直不得入门,今天有时间有温故了一次,有了些许进步在此记录下, 在开始之前,不得不提一下我的朋友给了我“豆瓣音乐 ...
- [javaEE] 三层架构案例-用户模块(二)
使用junit测试框架,测试查找用户和添加用户功能 com.tsh.test.xmlUserDaoTest package com.tsh.test; import org.junit.Test; i ...
- Javaee的Dao层的抽取
有时候我们在实现不同功能的时候回看到很多的Dao层的增加.修改.删除.查找都很相似,修改我们将他们提取BaseDao 一.提取前 1. 提取前的LinkDao层: public interface L ...
- mac os下载安装jmeter
一.简介 jmeter是属于apache的一个开源产品,纯Java应用.最初用来进行功能测试,而后又扩展了更多的测试功能. 二.下载 进入apache的jmeter下载页:http://jmeter. ...
- 基于bootstrap的内容折叠功能
加入js及css支持: <link rel="stylesheet" href="css/bootstrap.min.css"/> <scri ...