题解-PKUWC2018 猎人杀
Problem
题意概要:给定 \(n\) 个人的倒霉度 \(\{w_i\}\),每回合会有一个人死亡,每个人这回合死亡的概率为 自己的倒霉度
/目前所有存活玩家的倒霉度之和
,求第 \(1\) 个人最后一个死亡的概率
Solution
设 \(B = \sum_{i=2}^nw_i\)
要求 \(1\) 号最后一个被选中有点不好做,但是求 \(1\) 号第一个被选中还是比较好做的(\(\frac {w_1}{\sum_{i=1}^nw_i}\))
至于这两者怎么联系起来,使用 \(\mathrm {min-max}\) 容斥(和 HAOI2015按位或 有点像:前者是求 \(1\) 号最后一个被选中的概率;后者是求集合内最后一个被选中的期望次数)
\(\mathrm{min-max}\) 容斥的式子为:
\]
由于 \(1\) 为需要求的点,将其从 \(S,T\) 的定义中刨除,即
\]
而同时
\]
发现 \(\min\{T\}\) 最多有 \(B\) 种,可以考虑计算出每一项的容斥系数再 \(O(B)\) 计算
至于计算容斥系数,可使用母函数求得,即求出下面式子的每一项系数:
\]
分治 ntt 求解,每一层母函数度数和为 \(B\),复杂度 \(O(B\log B)\),总复杂度 \(O(B\log^2B)\)
Code
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
inline void read(int&x){
char ch=getchar();x=0;while(!isdigit(ch))ch=getchar();
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
}
const int N = 200200, T = 37, p = 998244353;
int stk[T], tp;
int brr[T][N], L[T];
int w[N], rev[N];
inline int qm(const int&x) {return x < p ? x : x - p;}
inline int qpow(int A, int B) {
int res = 1; while(B) {
if(B&1) res = (ll)res * A%p;
A = (ll)A * A%p, B >>= 1;
} return res;
}
void dft(int*a,int n,int sgn) {
for(int i=1;i<n;++i) if(i < rev[i]) swap(a[i], a[rev[i]]);
for(int i=1;i<n;i<<=1) {
int gn = qpow(3, (p-1)/(i<<1));
for(int j=0;j<n;j+=(i<<1)) {
int g = 1;
for(int k=0;k<i;++k,g=(ll)g*gn%p) {
int x = a[j+k], y = (ll)g * a[j+k+i]%p;
a[j+k] = qm(x + y), a[j+k+i] = qm(x - y+p);
}
}
}
if(!sgn) {
int iv = qpow(n, p-2); reverse(a+1,a+n);
for(int i=0;i<n;++i) a[i] = (ll)a[i] * iv%p;
}
}
void mul(int ai, int bi, int ci) {
int*ar = brr[ai], *br = brr[bi], *cr = brr[ci];
int n = L[ai], m = L[bi], &t = L[ci];
t = n + m; int nn = 1, l = 0;
while(nn < t) nn <<= 1, ++ l;
for(int i=n;i<nn;++i) ar[i] = 0;
for(int i=m;i<nn;++i) br[i] = 0;
for(int i=1;i<nn;++i) rev[i] = (rev[i>>1]>>1) | ((i&1) << l-1);
dft(ar, nn, 1), dft(br, nn, 1);
for(int i=0;i<nn;++i) cr[i] = (ll)ar[i] * br[i]%p;
dft(cr, nn, 0);
}
int binary(int l, int r) {
if(l == r) {
int id = stk[--tp], *arr = brr[id];
for(int i=w[l]<<1;i;--i) arr[i] = 0;
arr[0] = 1, arr[w[l]] = p-1, L[id] = w[l]+1;
return id;
}
int mid = l + r >> 1;
int ls = binary(l, mid), rs = binary(mid+1, r);
int id = stk[--tp]; mul(ls, rs, id);
stk[tp++] = ls, stk[tp++] = rs;
return id;
}
int main() {
for(int i=0;i<T;++i) stk[tp++] = T-i-1;
int n; read(n);
for(int i=1;i<=n;++i) read(w[i]);
int id = binary(2, n);
int ans = 0, *arr = brr[id];
for(int i=0;i<L[id];++i)
ans = (ans + (ll)arr[i] * qpow(w[1] + i, p-2))%p;
ans = (ll)ans * w[1]%p;
printf("%d\n",ans);
return 0;
}
题解-PKUWC2018 猎人杀的更多相关文章
- [PKUWC2018]猎人杀
题解 感觉是一道神题,想不出来 问最后\(1\)号猎人存活的概率 发现根本没法记录状态 每次转移的分母也都不一样 可以考虑这样一件事情: 如果一个人被打中了 那么不急于从所有人中将ta删除,而是给ta ...
- LOJ2541 PKUWC2018 猎人杀 期望、容斥、生成函数、分治
传送门 首先,每一次有一个猎人死亡之后\(\sum w\)会变化,计算起来很麻烦,所以考虑在某一个猎人死亡之后给其打上标记,仍然计算他的\(w\),只是如果打中了一个打上了标记的人就重新选择.这样对应 ...
- LOJ2541 PKUWC2018猎人杀(概率期望+容斥原理+生成函数+分治NTT)
考虑容斥,枚举一个子集S在1号猎人之后死.显然这个概率是w1/(Σwi+w1) (i∈S).于是我们统计出各种子集和的系数即可,造出一堆形如(-xwi+1)的生成函数,分治NTT卷起来就可以了. #i ...
- 洛谷 P5644 - [PKUWC2018]猎人杀(分治+NTT)
题面传送门 很久之前(2020 年)就听说过这题了,这么经典的题怎么能只听说而亲自做一遍呢 首先注意到每次开枪打死一个猎人之后,打死其他猎人概率的分母就会发生变化,这将使我们维护起来非常棘手,因此我们 ...
- [LOJ2541][PKUWC2018]猎人杀(容斥+分治+FFT)
https://blog.csdn.net/Maxwei_wzj/article/details/80714129 n个二项式相乘可以用分治+FFT的方法,使用空间回收可以只开log个数组. #inc ...
- 【洛谷5644】[PKUWC2018] 猎人杀(容斥+生成函数+分治NTT)
点此看题面 大致题意: 有\(n\)个人相互开枪,每个人有一个仇恨度\(a_i\),每个人死后会开枪再打死另一个还活着的人,且第一枪由你打响.设当前剩余人仇恨度总和为\(k\),则每个人被打中的概率为 ...
- [LOJ2541] [PKUWC2018] 猎人杀
题目链接 LOJ:https://loj.ac/problem/2541 Solution 很巧妙的思路. 注意到运行的过程中概率的分母在不停的变化,这样会让我们很不好算,我们考虑这样转化:假设所有人 ...
- 【LOJ2541】【PKUWC2018】猎人杀(容斥,FFT)
[LOJ2541][PKUWC2018]猎人杀(容斥,FFT) 题面 LOJ 题解 这题好神仙啊. 直接考虑概率很麻烦,因为分母总是在变化. 但是,如果一个人死亡之后,我们不让他离场,假装给他打一个标 ...
- 「PKUWC2018」猎人杀
「PKUWC2018」猎人杀 解题思路 首先有一个很妙的结论是问题可以转化为已经死掉的猎人继续算在概率里面,每一轮一直开枪直到射死一个之前没死的猎人为止. 证明,设所有猎人的概率之和为 \(W\) , ...
随机推荐
- Java核心复习 —— J.U.C 并发工具类
一.CountDownLatch 文档描述 A synchronization aid that allows one or more threads to wait until* a set of ...
- 在不切换分支的情况下,如何在all branches中快速查看指定分支,相对其他分支的状态
假设需要查看的分支为temp 1.git show temp 通过这个可以拿到commit id 2.查看TortoiseGit的日志,左下角勾选所有分支 3.在日志界面搜索commit id,然后右 ...
- 等待 Redis 应答
https://zhuanlan.zhihu.com/p/58608323 mq消息合并:由于mq请求发出到响应的时间,即往返时间, RTT(Round Time Trip),每次提交都要消耗RTT, ...
- 使用docker运行mysql
以前开发的时候都是用本地的sqlite开发,但是极少数情况下,sqlite支持的语法发布到服务器上链接mysql会报错. 为了避免这种现象,还是链接本地mysql开发还是更稳定的, 可是开发的项目多了 ...
- Windows 操作系统 端口转发
在Windows 下可以使用netsh interface portproxy 命令实现端口转发功能. 例:netsh interface portproxy add v4tov4 listenpor ...
- MauiMETA工具的使用(一)
MauiMETA工具的使用(一) 摘自:https://www.jianshu.com/p/a377119947f8 tianxiaoMCU 关注 2018.12.21 14:15 字数 267 ...
- Java中将一个反斜杠转换成两个反斜杠
代码示例: s = s.replaceAll("\\\\", "\\\\\\\\");
- AWS 核心服务概述(二)
目录 AWS网络服务 VPC Direct Connect Route53 AWS 计算服务 EC2 EMR(Elastic MapReduce) AWS Lambda Auto Scaling El ...
- Halcon 学习2 金属雕刻字识别
*HALCON 里面的例程名称:engraved_cnn.hdevread_image (Image, 'engraved')get_image_size (Image, Width, Height) ...
- docker安装并持久化postgresql数据库
安装docker步骤略过 1.拉取postgresql镜像 docker pull postgresql 2.创建本地卷,数据卷可以在容器之间共享和重用, 默认会一直存在,即使容器被删除(docker ...