[CTSC2018]假面
先来考虑一下第一问,血量有\(P\)的概率减\(1\)
由于我们最后需要求每一个人的期望血量,于是考虑维护出每个人处于不同血量时候的概率
一个简单\(dp\)即可
\]
\(dp_{i,j}\)表示第\(i\)个人血量为\(j\)的概率
第二问,发现概率的选中的人中的存活人数有关
于是每一个人只有两种状态,生或者是死,我们考虑求出每种存活人数对应的概率是多少
显然这里我们有一个这样的\(dp\)
\]
\(f_{i,j}\)表示前\(i\)个人里存活\(j\)个人的概率
发现做一次这个\(dp\)是\(O(n^2)\)的,对每一个人做一遍复杂度就是\(O(n^3)\)了
考虑消除某一个人的影响
设\(f_{j}=f_{n,j},g_{j}=g_{n-1,j}\),我们要消除第\(n\)的人的影响
显然边界条件有
\]
于是
\]
求得了\(g_0\),我们就可以递推了
\]
于是就有
\]
但是需要特别注意\(p_n=1\)的时候我们需要特殊处理一波,也非常简单就是
\]
代码
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define re register
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
const int maxn=205;
const LL mod=998244353;
inline int read() {
char c=getchar();int x=0;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,h[maxn],m,opt,a[maxn];
LL dp[maxn][105],t[maxn],f[maxn],inv[maxn],g[maxn];
inline LL ksm(LL a,LL b) {
LL S=1;
while(b) {if(b&1) S=S*a%mod;b>>=1;a=a*a%mod;}
return S;
}
int main() {
n=read();
for(re int i=1;i<=n;i++) h[i]=read(),dp[i][h[i]]=1;
m=read();
inv[1]=1;
for(re int i=2;i<=n;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
while(m--) {
opt=read();
if(opt==0) {
int pos=read();
LL u=read(),v=read();
LL P=ksm(v,mod-2)*u%mod;
for(re int i=0;i<=h[pos];i++)
dp[pos][i]=(dp[pos][i]*(1-P+mod)%mod+dp[pos][i+1]*P%mod)%mod;
}
if(opt==1) {
int k=read();
for(re int i=1;i<=k;i++) a[i]=read();
for(re int i=1;i<=k;i++) {
t[a[i]]=0;
for(re int j=1;j<=h[a[i]];j++)
t[a[i]]=(t[a[i]]+dp[a[i]][j])%mod;
}
memset(f,0,sizeof(f));
f[0]=1;
for(re int i=1;i<=k;i++)
for(re int j=i;j>=0;--j) {
f[j]=(f[j]*(1-t[a[i]]+mod)%mod)%mod;
if(j-1>=0) f[j]=(f[j]+f[j-1]*t[a[i]]%mod)%mod;
}
for(re int i=1;i<=k;i++) {
memset(g,0,sizeof(g));
if(t[a[i]]!=1) {
LL d=ksm((1-t[a[i]]+mod)%mod,mod-2);
g[0]=f[0]*d%mod;
for(re int j=1;j<k;j++)
g[j]=(f[j]-g[j-1]*t[a[i]]%mod+mod)*d%mod;
}
else for(re int j=1;j<=k;j++) g[j-1]=f[j];
LL ans=0;
for(re int j=0;j<k;j++)
ans=(ans+g[j]*inv[j+1]%mod)%mod;
printf("%lld ",ans*t[a[i]]%mod);
}
putchar(10);
}
}
for(re int i=1;i<=n;i++) {
LL now=0;
for(re int j=1;j<=h[i];j++)
now=(now+dp[i][j]*(LL)j%mod)%mod;
printf("%lld ",now);
}
return 0;
}
[CTSC2018]假面的更多相关文章
- BZOJ5340: [Ctsc2018]假面
BZOJ5340: [Ctsc2018]假面 https://lydsy.com/JudgeOnline/problem.php?id=5340 分析: 背包,只需要求\(g_{i,j}\)表示强制活 ...
- [CTSC2018] 假面 | 期望 DP
题目链接 LOJ 2552 Luogu P4564 考场上这道题我先是写了个70分暴力,然后发现似乎可以NTT,然鹅问题是--我没学过NTT,遂脑补之,脑补出来了,下午出成绩一看,卡成暴力分(70)- ...
- BZOJ5340: [Ctsc2018]假面【概率+期望】【思维】
LINK 思路 首先考虑减血,直接一个dp做过去,这个部分分不难拿 然后是\(op=1\)的部分 首先因为要知道每个人被打的概率,所以需要算出这个人活着的时候有多少个人活着时概率是什么 那么用\(g_ ...
- BZOJ5340 & 洛谷4564 & LOJ2552:[CTSC2018]假面——题解
https://www.lydsy.com/JudgeOnline/problem.php?id=5340 https://www.luogu.org/problemnew/show/P4564 ht ...
- bzoj 5340: [Ctsc2018]假面
Description 题面 Solution 生命值范围比较小,首先维护每一个人在每个血量的概率,从而算出生存的概率,设为 \(a[i]\) 询问时,只需要考虑生存的人数,可以 \(DP\) 设 \ ...
- BZOJ5340 [Ctsc2018]假面 【概率dp】
题目链接 BZOJ5340 题解 我们能很容易维护每个人当前各种血量的概率 设\(p[u][i]\)表示\(u\)号人血量为\(i\)的概率 每次攻击的时候,讨论一下击中不击中即可转移 是\(O(Qm ...
- [CTSC2018]假面(概率DP)
考场上以为CTSC的概率期望题都不可做,连暴力都没写直接爆零. 结果出来发现全场70以上,大部分AC,少于70的好像极少,感觉血亏. 设a[i][j]表示到当前为止第i个人的血量为j的概率(注意特判血 ...
- 并不对劲的bzoj5340:loj2552:uoj399:p4564: [Ctsc2018]假面
题目大意 有\(n\)(\(n\leq200\))个非负整数\(m_1,m_2,...,m_n\)(\(\forall i\in[1,n],m_i\leq100\)),有\(q\)(\(q\leq2* ...
- UOJ399 CTSC2018 假面 期望、DP
传送门 \(Q \leq 200000 , C \leq 1000 , m_i \leq 100\)-- 先考虑如何维护最后一次操作时所有人的血量期望.不难发现我们需要的复杂度是\(O(Qm_i)\) ...
随机推荐
- DataGrridView 当前行显示不同颜色
如果想让选中DataGridview的行显示不同颜色,就要通过DataGridview控件RowPerpaint事件中重新设置所选行的DefauleCellStyle属性来实现 private voi ...
- forms身份认证仍然能访问html页面解决办法
asp.net的forms身份认证保护是一个非常棒的东西,用VS2010创建一个Web应用程序即可看到范例 在web.config中配置 <authentication mode="F ...
- Java - 谨慎实现Comparable接口
类实现了Comparable接口就表明类的实例本身具有内在的排序关系(natural ordering). 因此,该类可以与很多泛型算法和集合实现进行协作. 而我们之需要实现Comparable接口唯 ...
- [C语言] 数据结构-预备知识指针
所有的伟大源于一个勇敢的开始 数据结构预备知识 指针 1.指针:是C语言的灵魂,指针=地址 地址:内存单元的编号 指针变量:存放内存单元地址的变量 int *p;//p是指针变量,int *表示该p变 ...
- linux下elasticsearch集成mongodb详细教程
由于公司业务需要,要用elasticsearch做索引库实现搜索功能,历尽千辛万苦,最后总算把mongodb和elasticsearch集成成功 1.搭建mongodb集群 参考https://www ...
- 1.Java设计模式-工厂模式
1.简单工厂模式(Factory Method) 常用的工厂模式是静态工厂模式,利用static修饰方法,作为一种类似于常见的工具类Utils等辅助效果,一般情况下工厂类不需要实例化. //1.定义一 ...
- 《JavaWeb从入门到改行》多重外键关系在java中的处理方案
目录:(点击红色方框展开子目录) 问题描述 无 项目案例说明 业务描述 数据库说明 项目源码及下载 无 问题描述 如上两图,数据库中各个表之间有很多的外键关系,其中业务关系是一个用户下有该用户的订单, ...
- [HNOI2006]最短母串问题
题目大意:给定一个字符串集,求一个最短字串,使得该集合内的串都是该串的一个子串 算法:AC自动机+最短路+状压DP 注意空间限制 #include"cstdio" #include ...
- [SD2015]序列统计——solution
http://www.lydsy.com/JudgeOnline/problem.php?id=3992 很容易得出DP方程: f[i][c]=f[i-1][a]*f[1][b]① 其中a*b%M=c ...
- 转:Redis和Memcache的区别分析
Redis和Memcache的区别分析 原文链接:http://blog.csdn.net/u013474436/article/details/48632665 简单区别: 1. Redis中,并 ...