LOJ#2552. 「CTSC2018」假面(期望 背包)
题意
Sol
多年以后,我终于把这题的暴力打出来了qwq 好感动啊。。
刚开始的时候想的是:
设\(f[i][j]\)表示第\(i\)轮, 第\(j\)个人血量的期望值
转移的时候若要淦这个人,那么\(f[i][j] = (f[i - 1][j] + 1) * p + (f[i - 1][j]) * (1 - p)\)
然后发现自己傻逼了。。因为期望不能正着推。
考虑直接推概率,设\(t[k][i][j]\)表示第\(k\)轮,第\(i\)个人,血量为\(j\)的概率
这玩意儿是可以转移的,就是判一下这次打中了没有
第二问可以对每个点分别算答案,设\(g[i][j]\)表示除必须活着的人外,前\(i\)个人中,有\(j\)个活着的概率,背包转移一下
这样复杂度是\(O(qn + n^3)\)的
显然第二问看起来非常暴力,
标算的做法好像叫“退背包”,也就是从背包中删除一个元素
先不考虑某个元素必须存活,推一遍得到\(g[i][j]\)表示前\(i\)个人中,有\(j\)个存活的概率
考虑转移的式子,设\(ali[i]\)表示第\(i\)个人活着的概率
\(g[i][j] = g[i - 1][j - 1] * ali[i] + g[i - 1][j] * (1 - ali[i])\)
而我们要得到的实际上就是\(g[i-1][j]\)这一项
那么\(g[i - 1][j] = \frac{g[i][j] - g[i - 1][j - 1] * ali[i]}{1 - ali[i]}\)
倒着推一遍即可,注意当\(1 - ali[i] = 0\)的时候需要特判,此时\(g[i - 1][j] = g[i][j + 1]\)
70分
#include<bits/stdc++.h>
//#define int long long
using namespace std;
const int MAXN = 201, mod = 998244353;
int f[2][MAXN], g[MAXN][MAXN], t[2][MAXN][MAXN];
// f: expect
inline int read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
int N, a[MAXN], Q, em[MAXN];
int fp(int a, int p) {
int base = 1;
while(p) {
if(p & 1) base = 1ll * base * a % mod;
a = 1ll * a * a % mod; p >>= 1;
}
return base;
}
int inv(int a) {
return fp(a, mod - 2);
}
int add(int x, int y) {
if(x + y < 0) return x + y + mod;
else return x + y >= mod ? x + y - mod : x + y;
}
int mul(int x, int y) {
x = (x + mod) % mod; y = (y + mod) % mod;
return 1ll * x * y % mod;
}
int solve(int id, int o, int N) {//这里dp的时候不能直接表示有j个活着,必须表示除i之外有j个活着。。
memset(g, 0, sizeof(g));
g[0][0] = 1;
for(int i = 1; i <= N; i++) {
for(int j = 0; j <= N; j++) {
if(em[i] ^ id) {
g[i][j] = mul(g[i - 1][j], t[o][em[i]][0]);
if(j) g[i][j] = add(g[i][j], mul(g[i - 1][j - 1], 1 - t[o][em[i]][0]));
}
else g[i][j] = g[i - 1][j];
}
}
int ans = 0;
for(int i = 0; i < N; i++)
ans = add(ans, mul(mul(1 - t[o][id][0], g[N][i]), inv(i + 1)));
return ans;
}
signed main() {
// freopen("a.in", "r", stdin);
// freopen("b.out", "w", stdout);
N = read();
for(int i = 1; i <= N; i++) a[i] = read(), t[0][i][a[i]] = 1, f[0][i] = a[i];
Q = read();
int o = 1;
for(int i = 1; i <= Q; i++, o ^= 1) {
int opt = read();
memcpy(t[o], t[o ^ 1], sizeof(t[o]));
if(opt == 0) {//
int id = read(), u = read(), v = read(), p = 1ll * u * inv(v) % mod;
t[o][id][0] = add(t[o][id][0], mul(p, t[o][id][1]));
for(int j = 1; j <= a[id]; j++) t[o][id][j] = add(mul(p, t[o ^ 1][id][j + 1]), mul(1 - p, t[o ^ 1][id][j]));
} else if(opt == 1) {
int k = read(), cnt = 0;
for(int i = 1; i <= k; i++) em[++cnt] = read();
for(int i = 1; i <= k; i++) printf("%d ", solve(em[i], o, cnt)); puts("");
}
}
for(int i = 1; i <= N; i++) {
int ans = 0;
for(int j = 1; j <= a[i]; j++)
ans = add(ans, mul(j, t[o ^ 1][i][j]));
printf("%d ", ans);
}
return 0;
}
/*
*/
100分
#include<bits/stdc++.h>
//#define int long long
using namespace std;
const int MAXN = 201, mod = 998244353;
int f[2][MAXN], g[MAXN][MAXN], t[2][MAXN][MAXN];
// f: expect
inline int read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
int N, a[MAXN], Q, em[MAXN], ans[MAXN], ali[MAXN], tp[MAXN], Inv[MAXN];
int add(int x, int y) {
if(x + y < 0) return x + y + mod;
else return x + y >= mod ? x + y - mod : x + y;
}
int mul(int x, int y) {
x = (x + mod) % mod; y = (y + mod) % mod;
return 1ll * x * y % mod;
}
int fp(int a, int p) {
int base = 1;
while(p) {
if(p & 1) base = 1ll * base * a % mod;
a = 1ll * a * a % mod; p >>= 1;
}
return base;
}
int inv(int a) {
a = add(a, mod);
return fp(a, mod - 2);
}
void Pre(int o, int N) {
// memset(g, 0, sizeof(g));
g[0][0] = 1;
for(int i = 1; i <= N; i++) {
ali[i] = (1 - t[o][em[i]][0] + mod) % mod;//alive
for(int j = 0; j <= i; j++) {
g[i][j] = mul(g[i - 1][j], t[o][em[i]][0]);
if(j) g[i][j] = add(g[i][j], mul(g[i - 1][j - 1], ali[i]));
}
}
}
int solve(int id, int o, int N) {
//memset(tp, 0, sizeof(tp));
if(!ali[id]) return 0;
if(ali[id] == 1) {
for(int i = 1; i <= N; i++) tp[i - 1] = g[N][i];
} else {
int down = inv(1 - ali[id]);
tp[0] = mul(g[N][0], down);
for(int i = 1; i <= N; i++)
tp[i] = mul(g[N][i] - mul(tp[i - 1], ali[id]), down);
}
int ans = 0;
for(int i = 1; i <= N; i++)
ans = add(ans, mul(mul(ali[id], tp[i - 1]), Inv[i]));
return ans;
}
signed main() {
//freopen("faceless10.in", "r", stdin);
// freopen("b.out", "w", stdout);
N = read();
for(int i = 1; i <= N; i++) a[i] = read(), t[0][i][a[i]] = 1, f[0][i] = a[i], Inv[i] = inv(i);
Q = read();
int o = 1;
for(int i = 1; i <= Q; i++, o ^= 1) {
int opt = read();
memcpy(t[o], t[o ^ 1], sizeof(t[o]));
if(opt == 0) {//
int id = read(), u = read(), v = read(), p = 1ll * u * inv(v) % mod;
t[o][id][0] = add(t[o][id][0], mul(p, t[o][id][1]));
for(int j = 1; j <= a[id]; j++) t[o][id][j] = add(mul(p, t[o ^ 1][id][j + 1]), mul(1 - p, t[o ^ 1][id][j]));
} else if(opt == 1) {
int k = read();
for(int i = 1; i <= k; i++) em[i] = read();
Pre(o, k);
for(int i = k; i >= 1; i--) ans[i] = solve(i, o, k);
for(int i = 1; i <= k; i++) printf("%d ", ans[i]); puts("");
}
}
for(int i = 1; i <= N; i++) {
int ans = 0;
for(int j = 1; j <= a[i]; j++)
ans = add(ans, mul(j, t[o ^ 1][i][j]));
printf("%d ", ans);
}
return 0;
}
/*
*/
LOJ#2552. 「CTSC2018」假面(期望 背包)的更多相关文章
- loj#2552. 「CTSC2018」假面
题目链接 loj#2552. 「CTSC2018」假面 题解 本题严谨的证明了我菜的本质 对于砍人的操作好做找龙哥就好了,blood很少,每次暴力维护一下 对于操作1 设\(a_i\)为第i个人存活的 ...
- LOJ 2552 「CTSC2018」假面——DP
题目:https://loj.ac/problem/2552 70 分就是 f[i][j] 表示第 i 个人血量为 j 的概率.这部分是 O( n*Q ) 的:g[i][j][0/1] 表示询问的人中 ...
- Loj #2554. 「CTSC2018」青蕈领主
Loj #2554. 「CTSC2018」青蕈领主 题目描述 "也许,我的生命也已经如同风中残烛了吧."小绿如是说. 小绿同学因为微积分这门课,对"连续"这一概 ...
- Loj #2553. 「CTSC2018」暴力写挂
Loj #2553. 「CTSC2018」暴力写挂 题目描述 temporaryDO 是一个很菜的 OIer .在 4 月,他在省队选拔赛的考场上见到了<林克卡特树>一题,其中 \(k = ...
- LOJ 2553 「CTSC2018」暴力写挂——边分治+虚树
题目:https://loj.ac/problem/2553 第一棵树上的贡献就是链并,转化成 ( dep[ x ] + dep[ y ] + dis( x, y ) ) / 2 ,就可以在第一棵树上 ...
- 「CTSC2018」假面
真~签到题qwq 昨天在考场上先写了个70分暴力dp,然后发现好像可以优化.因为结界技能的模型相当于要求出 对于每个物品,仅仅不选它的背包是什么.... 于是当场脑补出两种做法: 前缀和后缀背包卷积 ...
- LOJ 2557 「CTSC2018」组合数问题 (46分)
题目:https://loj.ac/problem/2557 第一个点可以暴搜. 第三个点无依赖关系,k=3,可以 DP .dp[ cr ][ i ][ j ] 表示前 cr 个任务.第一台机器最晚完 ...
- LOJ 2555 「CTSC2018」混合果汁——主席树
题目:https://loj.ac/problem/2555 二分答案,在可以选的果汁中,从价格最小的开始选. 按价格排序,每次可以选的就是一个前缀.对序列建主席树,以价格为角标,维护体积和.体积*价 ...
- LOJ 2554 「CTSC2018」青蕈领主——结论(思路)+分治FFT
题目:https://loj.ac/problem/2554 一个“连续”的区间必然是一个排列.所有 r 不同的.len 最长的“连续”区间只有包含.相离,不会相交,不然整个是一个“连续”区间. 只有 ...
随机推荐
- mongoengine使用示例
#!/usr/bin/env python # -*- coding: utf-8 -*- ''' Created on 2017年10月19日 @author: zzy ''' from mongo ...
- Python3之os模块
一:简介 os模块主要用于提供系统高级别的操作. 二:常用方法 os.access(path, mode) # 检验权限模式 os.chdir(path) # 改变当前工作目录 os.chflags( ...
- c语言数据结构学习心得——栈
栈(Stack) 只允许在一端进行插入或删除操作的线性表 栈顶(Top):栈中允许进行插入和删除操作的那一端: 栈底(Bottom):固定的,不允许进行插入或删除的另一端 1.栈是受限的线性表,所以自 ...
- CentOS7更换yum源为阿里云镜像源
1. 备份原来的yum源 cp /etc/yum.repos.d/CentOS-Base.repo /etc/yum.repos.d/CentOS-Base.repo.bak 2.设置aliyun的y ...
- Exadata 18.1新特性--计算节点升级增强
新的计算节点升级增强能大幅缩短升级花费的时间,最快能减少40%的时间,具体的增强主要体现在以下方面: 1.计算节点操作系统的备份工作. 以前计算节点操作系统的备份工作是在正式升级之前执行,在滚动升级过 ...
- 开发环境,eclipse编辑器java代码自动提示
Eclipse+ADT+Android SDK 搭建安卓开发环境 eclipse编辑器java代码自动提示 window-->Preferences-->JAva-->Content ...
- ZQUOJ 22854 (优先队列+剪枝)
题目:给出K , N , M : N为顶点数 , M为边数 : 求K个从1到N的不重复的最短边 , 可以来回的走: 分析:很自然的就可以想到用个优先队列广收下K次终点嘛 , 但是.0.0 爆了 ...
- 【研究】ms17-010永恒之蓝漏洞复现
1 永恒之蓝漏洞复现(ms17-010) 1.1 漏洞描述: Eternalblue通过TCP端口445和139来利用SMBv1和NBT中的远程代码执行漏洞,恶意代码会扫描开放44 ...
- 转: centos系统home下的中文目录改为英文目录
转自h t t p : / /xugang-1017-126-com.iteye.com/blog/2081845 如果安装了中文版的Cent OS之后,root目录和home目录下会出现中文的路径名 ...
- JS 判断节点类型
节点类型的分类 节点类型 说明 值 元素节点 每一个HTML标签都是一个元素节点,如 <div> . <p>.<ul>等 1 属性节点 元素节点(HTML标签)的属 ...