[CF895E] Eyes Closed(线段树,期望)
Desctiption
传送门:Portal
大致题意: 给你一个序列, 支持两种操作:
1 l1 r1 l2 y2
在\([l1, r1]\)随机选择一个数a
, \([l2, r2]\) 内随机选择一个数b
, 交换a, b.2 l r
询问一个区间的期望.
\]
Solution
根据期望线性性,只需要维护出每个位置元素值的期望就可以.
区间操作期望问题的经典套路是维护出每个位置的期望.
考虑一个数字\(a_i\)
他有\(\frac{len - 1}{len}\) 的概率保持原数, 否则,根据全期望公式, 它会变成\(E(b_i) = \sum\frac{b_i}{lenB}\)
那么有:
\]
\(b_i\)也同理.
所以我们只需要维护一个区间加/乘, 区间求和的线段树就可以了.
Code
#include<bits/stdc++.h>
using namespace std;
#define rep(i, a, b) for(int i = (a), i##_end_ = (b); i <= i##_end_; ++i)
#define drep(i, a, b) for(int i = (a), i##_end_ = (b); i >= i##_end_; --i)
#define clar(a, b) memset((a), (b), sizeof(a))
#define debug(...) fprintf(stderr, __VA_ARGS__)
typedef long long LL;
typedef long double LD;
int read() {
char ch = getchar();
int x = 0, flag = 1;
for (;!isdigit(ch); ch = getchar()) if (ch == '-') flag *= -1;
for (;isdigit(ch); ch = getchar()) x = x * 10 + ch - 48;
return x * flag;
}
void write(int x) {
if (x < 0) putchar('-'), x = -x;
if (x >= 10) write(x / 10);
putchar(x % 10 + 48);
}
const int Maxn = 200009;
int a[Maxn];
namespace SGMTtree {
double tree[Maxn << 3], addTag[Maxn << 3], mulTag[Maxn << 3];
#define lc(x) ((x) << 1)
#define rc(x) (((x) << 1) | 1)
#define ls rt << 1, l, mid
#define rs (rt << 1) | 1, mid + 1, r
void setAdd(int root, int l, int r, double v) {
tree[root] += v * (r - l + 1.0); addTag[root] += v;
}
void setMul(int root, int l, int r, double v) {
tree[root] *= v;
mulTag[root] *= v; addTag[root] *= v;
}
void (*setTag)(int, int, int, double);
void set(int p) {
if(p == 1) setTag = setAdd;
if(p == 2) setTag = setMul;
}
void pushdown(int rt, int l, int r) {
int mid = (l + r) >> 1;
if (mulTag[rt] != 1.0) {
setMul(ls, mulTag[rt]), setMul(rs, mulTag[rt]);
mulTag[rt] = 1;
}
if (addTag[rt] != 0.0) {
setAdd(ls, addTag[rt]); setAdd(rs, addTag[rt]);
addTag[rt] = 0;
}
}
void pushup(int rt) { tree[rt] = tree[lc(rt)] + tree[rc(rt)]; }
void build(int rt, int l, int r) {
addTag[rt] = 0, mulTag[rt] = 1;
if (l == r) {
tree[rt] = a[l];
return ;
}
int mid = (l + r) >> 1;
build(ls), build(rs);
pushup(rt);
}
void modify(int rt, int l, int r, int p, int q, double v) {
if (p > q) return ;
if (p <= l && r <= q) {
setTag(rt, l, r, v);
return ;
}
int mid = (l + r) >> 1; pushdown(rt, l, r);
if (q <= mid) modify(ls, p, q, v);
else if (p > mid) modify(rs, p, q, v);
else modify(ls, p, q, v), modify(rs, p, q, v);
pushup(rt);
}
double query(int rt, int l, int r, int p, int q) {
if (p <= l && r <= q) return tree[rt];
int mid = (l + r) >> 1; pushdown(rt, l, r);
if (q <= mid) return query(ls, p, q);
else if (p > mid) return query(rs, p, q);
else return query(ls, p, q) + query(rs, p, q);
}
#undef lc
#undef rc
#undef ls
#undef rs
};
int n, q;
void init() {
n = read(); q = read();
rep (i, 1, n) a[i] = read();
SGMTtree :: build(1, 1, n);
}
void solve() {
rep (i, 1, q) {
int opt = read();
if (opt == 1) {
int l1 = read(), r1 = read(), l2 = read(), r2 = read();
double c = r1 - l1 + 1.0, d = r2 - l2 + 1.0, tmp1 = SGMTtree :: query(1, 1, n, l1, r1), tmp2 = SGMTtree :: query(1, 1, n, l2, r2);
SGMTtree :: set(2), SGMTtree :: modify(1, 1, n, l1, r1, c - 1.0);
SGMTtree :: set(2), SGMTtree :: modify(1, 1, n, l2, r2, d - 1.0);
SGMTtree :: set(1), SGMTtree :: modify(1, 1, n, l1, r1, tmp2 / d);
SGMTtree :: set(1), SGMTtree :: modify(1, 1, n, l2, r2, tmp1 / c);
SGMTtree :: set(2), SGMTtree :: modify(1, 1, n, l1, r1, 1 / c);
SGMTtree :: set(2), SGMTtree :: modify(1, 1, n, l2, r2, 1 / d);
}
if (opt == 2) {
int l = read(), r = read();
printf("%.7lf\n", SGMTtree :: query(1, 1, n, l, r));
}
}
}
int main() {
// freopen("CF895E.in", "r", stdin);
// freopen("CF895E.out", "w", stdout);
init();
solve();
#ifdef Qrsikno
debug("\nRunning time: %.3lf(s)\n", clock() * 1.0 / CLOCKS_PER_SEC);
#endif
return 0;
}
调试:
线段树函数指针必须要用namespace.
[CF895E] Eyes Closed(线段树,期望)的更多相关文章
- CF895 E. Eyes Closed(线段树 期望)
题意 Sol 今天考试的T3,,我本来留了一个小时去写.但是T2一刚就刚了两个小时 最后也没来的及写.. 然后考完 开始写,,25min就A了.. 感觉自己太高估自己的思维,太低估自己的码力了... ...
- BZOJ 2752: [HAOI2012]高速公路(road) [线段树 期望]
2752: [HAOI2012]高速公路(road) Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 1219 Solved: 446[Submit] ...
- JZYZOJ1527 [haoi2012]高速公路 线段树 期望
http://172.20.6.3/Problem_Show.asp?id=1527 日常线段树的pushdown写挂,果然每次写都想得不全面,以后要注意啊……求期望部分也不熟练,和平均数搞混也是or ...
- BZOJ2752: [HAOI2012]高速公路(road)(线段树 期望)
Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 1820 Solved: 736[Submit][Status][Discuss] Descripti ...
- BZOJ 4262 线段树+期望
思路: 把询问离线下来,查询max和查询min相似,现在只考虑查询max 令sum[l,r,x]表示l到r内的数为左端点,x为右端点的区间询问的答案 那么询问就是sun[l1,r1,r2]-sum[l ...
- 洛谷P3924 康娜的线段树(期望 前缀和)
题意 题目链接 Sol 思路就是根据期望的线性性直接拿前缀和算贡献.. 这题输出的时候是不需要约分的qwq 如果你和我一样为了AC不追求效率的话直接#define int __int128就行了.. ...
- Codeforces 895E Eyes Closed(线段树)
题目链接 Eyes Closed 题意 两个人玩一个游戏,现在有两种操作: 1.两个人格子挑选一个区间,保证两个的区间不相交.在这两个区间里面各选出一个数,交换这两个数. 2.挑选一个区间,求这个 ...
- P3924 康娜的线段树(期望)
P3924 康娜的线段树 看起来$O(nlogn)$可过其实由于巨大常数是无法通过的 $O(nlogn)$:70pts 我们手玩样例发现 线段树上某个节点的期望值$f[o]=(f[lc]+f[rc]) ...
- hdu-5805 NanoApe Loves Sequence(线段树+概率期望)
题目链接: NanoApe Loves Sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 262144/131072 ...
随机推荐
- Android关于SurfaceView,SurfaceHolder,SurfaceHolder.CallBack详解
官方的定义: 1.SurfaceView SurfaceView是视图(View)的继承类,这个视图里内嵌了一个专门用于绘制的Surface.你可以控制这个Surface的格式和尺寸.Surfacev ...
- ECharts 图表导出
Echarts图形是由Javascript亲自在前端网页上绘制的 1.用ECharts配置项手册中的toolbox.feature.saveAsImage toolbox: { show: true, ...
- Saying goodbye to Flash in Chrome一代人的回忆FLASH
一早打开chorme就推送了这条FLASH将在2020年推出CHORME 想起了当年风靡全球的flash热潮,游戏视频动画,都由flash运行,最熟悉的童年游戏4399,小时候的天堂. 说起这个不得不 ...
- git 常用命令与上传步骤
git 上传步骤: git init 初始化Git仓库 git add . 提交你修改的文件 git status 查看项目当中的状态(红色代表的是 未add 的文件 绿色的是已 ...
- redis的hmset乐观锁的实现
1.lua脚本(集成实现了乐观锁,hmset ,expire等) local key=KEYS[1]; local oldVerion=tonumber(ARGV[1]); local seconds ...
- Codeforces Round #393 (Div. 2) - B
题目链接:http://codeforces.com/contest/760/problem/B 题意:给定n张床,m个枕头,然后给定某个特定的人(n个人中的其中一个)他睡第k张床,问这个人最多可以拿 ...
- Ubuntu18.04 安装 Idea 2018.2
https://blog.csdn.net/weixx3/article/details/81136822 Ubuntu18.04 安装 Idea 2018.2环境信息:OS:Ubuntu18.04J ...
- 关系型数据库---MYSQL---优化
1.为什么要进行SQL优化? 1.1 查询性能低 1.2 执行时间过长 1.3 等待时间过长 1.4 SQL写的太差(尤其是多表关联查询) 1.5 索引失效 1.6 服务器参数(缓存.线程数)设置不合 ...
- Ubuntu安装openmpi
Ubuntu 下安装 openmpi 需要同时安装下面三个包: sudo apt-get install openmpi-bin openmpi-doc libopenmpi-dev 建议更新系统源之 ...
- django 添加分页功能的包
Django pagination based upon the core pagination module