题目:https://loj.ac/problem/2552

70 分就是 f[i][j] 表示第 i 个人血量为 j 的概率。这部分是 O( n*Q ) 的;g[i][j][0/1] 表示询问的人中,前/后 i 个人,存活 j 个人的概率。做 g[ ][ ] 是 n^2 的,算答案是 n3 的。

考虑 g[ i ] 表示询问的人中有 i 个存活的概率。因为每个人加入 g[ ] 的顺序无关,所以可以每次 O(n) 地从g[ ] 里剔除第 i 个人的贡献。

令第 i 个人不存活的概率是 u ,存活的概率是 v 。

当初的转移是 g[ i ] * u -> g'[ i ] , g[ i ] * v -> g'[ i+1 ] ,所以现在可以倒着做:g[ i ] = g'[ i+1 ] / v , g'[ i ] -= g[ i ] ,或者正着做:g[ i ] = g'[ i ] / u , g'[ i+1 ] -= g[ i ] * v 。

发现如果 u=0 ,那么 g[ i ] 的值只体现在了 g'[ i+1 ] 里;所以倒着做。如果 v=0 就正着做。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
int rdn()
{
int ret=;bool fx=;char ch=getchar();
while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
return fx?ret:-ret;
}
const int N=,M=,mod=;
int upt(int x){while(x>=mod)x-=mod;while(x<)x+=mod;return x;}
int pw(int x,int k)
{int ret=;while(k){if(k&)ret=(ll)ret*x%mod;x=(ll)x*x%mod;k>>=;}return ret;} int n,a[N],f[N][M],g[N],tp[N],h[N],q[N],inv[N];
int main()
{
n=rdn(); inv[]=;
for(int i=;i<=n;i++)inv[i]=(ll)upt(-mod/i)*inv[mod%i]%mod;
for(int i=,d;i<=n;i++)
{ a[i]=rdn();f[i][a[i]]=;}
int Q=rdn(),op,x,u,v;
while(Q--)
{
op=rdn();
if(!op)
{
x=rdn();u=rdn();v=rdn();
u=(ll)u*pw(v,mod-)%mod;
for(int j=;j<=a[x];j++)
{
f[x][j-]=(f[x][j-]+(ll)f[x][j]*u)%mod;
f[x][j]=(ll)f[x][j]*upt(-u)%mod;
}
if(!f[x][a[x]])a[x]--;
}
else
{
x=rdn(); int tl=;
for(int i=;i<=x;i++)
q[++tl]=rdn();
g[]=;for(int i=;i<=tl;i++)g[i]=;
for(int i=;i<=tl;i++)
{
u=f[q[i]][]; v=upt(-u);
for(int j=i;j>=;j--)
{
g[j]=(ll)g[j]*u%mod;
if(j)g[j]=(g[j]+(ll)g[j-]*v)%mod;
}
}
for(int i=;i<=tl;i++)
{
u=f[q[i]][]; v=upt(-u);
memcpy(tp,g,sizeof g);
if(!u)
{
int iv=pw(v,mod-);
for(int j=tl-;j>=;j--)
{
h[j]=(ll)tp[j+]*iv%mod;
tp[j]=upt((tp[j]-(ll)u*h[j])%mod);
}
}
else
{
int iu=pw(u,mod-);
for(int j=;j<tl;j++)
{
h[j]=(ll)tp[j]*iu%mod;
tp[j+]=upt((tp[j+]-(ll)v*h[j])%mod);
}
}
int ans=;
for(int j=;j<tl;j++)
ans=(ans+(ll)h[j]*v%mod*inv[j+])%mod;
printf("%d ",ans);
}
puts("");
}
}
for(int i=;i<=n;i++)
{
int ans=;
for(int j=;j<=a[i];j++)
ans=(ans+(ll)j*f[i][j])%mod;
printf("%d ",ans);
}
puts(""); return ;
}

LOJ 2552 「CTSC2018」假面——DP的更多相关文章

  1. loj#2552. 「CTSC2018」假面

    题目链接 loj#2552. 「CTSC2018」假面 题解 本题严谨的证明了我菜的本质 对于砍人的操作好做找龙哥就好了,blood很少,每次暴力维护一下 对于操作1 设\(a_i\)为第i个人存活的 ...

  2. LOJ#2552. 「CTSC2018」假面(期望 背包)

    题意 题目链接 Sol 多年以后,我终于把这题的暴力打出来了qwq 好感动啊.. 刚开始的时候想的是: 设\(f[i][j]\)表示第\(i\)轮, 第\(j\)个人血量的期望值 转移的时候若要淦这个 ...

  3. Loj #2554. 「CTSC2018」青蕈领主

    Loj #2554. 「CTSC2018」青蕈领主 题目描述 "也许,我的生命也已经如同风中残烛了吧."小绿如是说. 小绿同学因为微积分这门课,对"连续"这一概 ...

  4. Loj #2553. 「CTSC2018」暴力写挂

    Loj #2553. 「CTSC2018」暴力写挂 题目描述 temporaryDO 是一个很菜的 OIer .在 4 月,他在省队选拔赛的考场上见到了<林克卡特树>一题,其中 \(k = ...

  5. LOJ 2553 「CTSC2018」暴力写挂——边分治+虚树

    题目:https://loj.ac/problem/2553 第一棵树上的贡献就是链并,转化成 ( dep[ x ] + dep[ y ] + dis( x, y ) ) / 2 ,就可以在第一棵树上 ...

  6. LOJ 2557 「CTSC2018」组合数问题 (46分)

    题目:https://loj.ac/problem/2557 第一个点可以暴搜. 第三个点无依赖关系,k=3,可以 DP .dp[ cr ][ i ][ j ] 表示前 cr 个任务.第一台机器最晚完 ...

  7. 「CTSC2018」假面

    真~签到题qwq 昨天在考场上先写了个70分暴力dp,然后发现好像可以优化.因为结界技能的模型相当于要求出 对于每个物品,仅仅不选它的背包是什么....  于是当场脑补出两种做法: 前缀和后缀背包卷积 ...

  8. LOJ 2555 「CTSC2018」混合果汁——主席树

    题目:https://loj.ac/problem/2555 二分答案,在可以选的果汁中,从价格最小的开始选. 按价格排序,每次可以选的就是一个前缀.对序列建主席树,以价格为角标,维护体积和.体积*价 ...

  9. LOJ 2554 「CTSC2018」青蕈领主——结论(思路)+分治FFT

    题目:https://loj.ac/problem/2554 一个“连续”的区间必然是一个排列.所有 r 不同的.len 最长的“连续”区间只有包含.相离,不会相交,不然整个是一个“连续”区间. 只有 ...

随机推荐

  1. 用闭包解决 js 循环中函数变量暂存问题

    需求:有一个数组,根据数组的值渲染对应的数字div,单击对应的div 在控制台打印对应的数字.如点击1,控制台打印1. 问题: 不管点击哪个值 打出来都是4 代码如下 <!DOCTYPE htm ...

  2. (转)Java 原子性引用 AtomicReference

    链接:https://www.jianshu.com/p/882d0e2c3ea6 來源:简书  作者:专职跑龙套 AtomicReference An object reference that m ...

  3. python 字符串str和json格式转换

    最近在写一个脚本,需要处理从excel中读取的数据,发现读取的json格式数据进行转换时报错 ValueError: Expecting property name enclosed in doubl ...

  4. Oracle 简单统计示例

    有数据如下: eg1:现在需要统计所有男性人员数量,所有女性人员数量,sclassno=10000的男性人员的总年龄,年龄大于20的女性人员数量 ----sign( number )/*If numb ...

  5. 测开之路四十九:用Django实现扑克牌游戏

    用Django实现和之前flask一样的扑克牌游戏 项目结构 html <!DOCTYPE html><html lang="en"><head> ...

  6. COALESCE 函数作用

    用途. 将空值替换成其他值 返回第一个非空值. 任意一个不为空的值.比较有用.

  7. PHP+JS的信息提示弹窗

    基于PHP函数的Msg信息提示框 1.可以设置弹出信息,跳转地址,跳转的时间,跳转的信息标题提示: 2.代码实例: <?php function ShowMsg($msg, $gourl,$ti ...

  8. 转:父类私有变量是否被子类继承详细解说(答案:内存中存在,但sun公司定义为不继承)

    应作者要求,本处提供一个连接,表示对原作者版权尊重. https://blog.csdn.net/mr_duantao/article/details/50966471

  9. Loadrunner实现多台测试机负载时可能遇到的难题

    压测机安装的loadrunner版本为11.0 压测机器的操作系统是windows server 2016 所有的负载机器操作系统和压侧机器是一样的,都是windows server 2016 Err ...

  10. 【Python—待解惑的问题】

    疑问1: map()函数第一个参数接收函数,所以 dict.keys() 是函数?这个不是方法吗? map(dict.keys,s) s1.keys() == dict.keys(s1) 输入: Tr ...