「ZJOI2017」树状数组

以下均基于模2意义下,默认\(n,m\)同阶。

熟悉树状数组的应该可以发现,这题其实是求\(l-1\)和\(r\)位置值相同的概率。

显然\(l=1\)的情况需要特盘。

大暴力

对于\(l=1\)的情况,可以发现一个操作不会产生影响当且仅当增加\(r\)的值,而其他情况会改变\(l-1\)或\(r\)。

对于\(l!=1\)的情况:

​ 针对一次修改区间\([ql,qr]\)。

  1. \([ql,qr]\)包含\(l-1,r\),那么有\(\displaystyle 2 \over qr-ql+1\)概率使\(l-1\)或\(r\)的值改变。
  2. \([ql,qr]\)不包含\(l-1,r\),不会发生变化。
  3. \([ql,qr]\)包含\(l-1,r\)中一个,那么有\(\displaystyle 1 \over qr-ql+1\)使\(l-1\)或\(r\)的值改变。

把每次修改记下来,就可以写出\(o(n^2)\)的暴力了。

代码见namespace fc​

正解

显然可以只记相等的概率。

可以发现,对于一次询问,改变前面修改的顺序并不会改变该询问的答案。

也就是说它满足交换律。

对于\(l=1\)的情况,显然可以一棵线段树维护。

对于\(l!=1\)的情况,把\([l,r ]\)区间当做一个二维的点\((l,r)\),那么每一次修改都会对二维区间产生贡献。

具体的:对于暴力分的第1类,即\(x \in [ql,qr],y \in [ql,qr]\)的点,对于暴力的第3类,即\(x \in [1,ql -1],y \in [ql,qr]\),\(x \in [ql,qr],y \in [qr+1,n]\)。

只需二维线段树区间标记即可,为了方便,标记表示的是 对应的二维区间的点 发生变化的概率。

#include<bits/stdc++.h>
#define rep(q,a,b) for(int q=a,q##_end_=b;q<=q##_end_;++q)
#define dep(q,a,b) for(int q=a,q##_end_=b;q>=q##_end_;--q)
#define mem(a,b) memset(a,b,sizeof a )
#define debug(a) cerr<<#a<<' '<<a<<"___"<<endl
using namespace std;
template<typename T>
void in(T &r) {
static char c;
r=0;
while(c=getchar(),!isdigit(c));
do r=(r<<1)+(r<<3)+(c^48);
while(c=getchar(),isdigit(c));
}
bool cur1;
int n,m;
const int mn=100005;
const int mod=998244353;
int inv[mn];
namespace fc{
int l[3005],r[3005];
void solve(){
int ty,a,b;
int ct=0;
while(m--){
in(ty),in(a),in(b);
if(ty==1)++ct,l[ct]=a,r[ct]=b;
else{
--a;
int wi=1,wo=0;
if(!a){
rep(q,1,ct){
int len=r[q]-l[q]+1;
if(l[q]<=b&&b<=r[q]){
int mid1=wi,mid2=wo;
wi=(1LL*mid2*(1-inv[len])+1LL*mid1*inv[len])%mod;
wo=(1LL*mid1*(1-inv[len])+1LL*mid2*inv[len])%mod;
}else swap(wi,wo);
}
}else{
rep(q,1,ct){
int len=r[q]-l[q]+1;
int mid1=wi,mid2=wo;
if(l[q]<=a&&b<=r[q]){//[l[q],r[q]]->[l[q],r[q]]
wi=(1LL*mid1*(1-2*inv[len])+1LL*mid2*2*inv[len])%mod;
wo=(1LL*mid2*(1-2*inv[len])+1LL*mid1*2*inv[len])%mod;
}else if(l[q]<=b&&b<=r[q]||l[q]<=a&&a<=r[q]){
//[1,l[q]-1]->[l[q],r[q]]
//[l[q],r[q]]->[r[q]+1,n]
wi=(1LL*mid1*(1-inv[len])+1LL*mid2*inv[len])%mod;
wo=(1LL*mid2*(1-inv[len])+1LL*mid1*inv[len])%mod;
}
}
}
printf("%d\n",(wi+mod)%mod);
}
}
}
}
namespace something_just_for_fun{
struct two_dimensional_segment_tree{
int tot,lson[mn*400],rson[mn*400],addv[mn*400],rt[mn<<2];
two_dimensional_segment_tree(){
tot=0;
}
int y_1,y_2,y_3,y_4,ad_v;
void se_add(int &o,int l,int r){
if(!o)o=++tot;
if(y_3<=l&&r<=y_4){
addv[o]=(1LL*(1-addv[o])*ad_v+1LL*addv[o]*(1-ad_v))%mod;
}else{
int mid=l+r>>1;
if(y_3<=mid)se_add(lson[o],l,mid);
if(y_4>mid)se_add(rson[o],mid+1,r);
}
}
void fi_add(int o,int l,int r){
if(y_1<=l&&r<=y_2)se_add(rt[o],1,n);
else{
int mid=l+r>>1;
if(y_1<=mid)fi_add(o<<1,l,mid);
if(y_2>mid)fi_add(o<<1|1,mid+1,r);
}
}
void add(int l,int r,int l1,int r1,int v){
if(l>r||l1>r1)return;
y_1=l,y_2=r,y_3=l1,y_4=r1,ad_v=v;
fi_add(1,1,n);
}
int v;
void se_ask(int &o,int l,int r){
if(!o)return;
v=(1LL*(1-addv[o])*v+1LL*addv[o]*(1-v))%mod;
int mid=l+r>>1;
if(y_2<=mid)se_ask(lson[o],l,mid);
else se_ask(rson[o],mid+1,r);
}
void fi_ask(int o,int l,int r){
se_ask(rt[o],1,n);
if(l==r)return;
else{
int mid=l+r>>1;
if(y_1<=mid)fi_ask(o<<1,l,mid);
else fi_ask(o<<1|1,mid+1,r);
}
}
int ask(int l,int r){
y_1=l,y_2=r,v=1;
fi_ask(1,1,n);
return v;
}
}an;
struct segment_tree{
int addv[mn<<2];
int y_1,y_2,ad_v;
void fi_add(int o,int l,int r){
if(y_1<=l&&r<=y_2)addv[o]=(1LL*(1-addv[o])*ad_v+1LL*addv[o]*(1-ad_v))%mod;
else{
int mid=l+r>>1;
if(y_1<=mid)fi_add(o<<1,l,mid);
if(y_2>mid)fi_add(o<<1|1,mid+1,r);
}
}
void add(int l,int r,int v){
if(l>r)return;
y_1=l,y_2=r,ad_v=v;
fi_add(1,1,n);
}
int v;
void fi_ask(int o,int l,int r){
v=(1LL*(1-addv[o])*v+1LL*addv[o]*(1-v))%mod;
if(l==r)return;
else{
int mid=l+r>>1;
if(y_1<=mid)fi_ask(o<<1,l,mid);
else fi_ask(o<<1|1,mid+1,r);
}
}
int ask(int x){
y_1=x,v=1;
fi_ask(1,1,n);
return v;
}
}at;
int l[mn],r[mn];
void solve(){
int ty,a,b;
while(m--){
in(ty),in(a),in(b);
if(ty==1){
at.add(a,b,1-inv[b-a+1]),at.add(1,a-1,1),at.add(b+1,n,1);
an.add(a,b,a,b,2*inv[b-a+1]%mod),an.add(1,a-1,a,b,inv[b-a+1]),an.add(a,b,b+1,n,inv[b-a+1]);
}
else --a,printf("%d\n",((!a?at.ask(b):an.ask(a,b))+mod)%mod);
}
}
}
bool cur2;
int main(){
// cerr<<(&cur2-&cur1)/1024.0/1024<<endl;
freopen("bit.in","r",stdin);
freopen("bit.out","w",stdout);
in(n),in(m);
inv[0]=inv[1]=1;
rep(q,2,n)inv[q]=1LL*(mod-mod/q)*inv[mod%q]%mod;
something_just_for_fun::solve();
return 0;
}

「ZJOI2017」树状数组的更多相关文章

  1. 「ZJOI2017」树状数组(二维线段树)

    「ZJOI2017」树状数组(二维线段树) 吉老师的题目真是难想... 代码中求的是 \(\sum_{i=l-1}^{r-1}a_i\),而实际求的是 \(\sum_{i=l}^{r}a_i\),所以 ...

  2. UOJ#291. 【ZJOI2017】树状数组 树套树

    原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ291.html 题解 结论:这个写错的树状数组支持的是后缀加和后缀求和.这里的后缀求和在 x = 0 的时 ...

  3. 【ZJOI2017】树状数组

    题目描述 漆黑的晚上,九条可怜躺在床上辗转反侧.难以入眠的她想起了若干年前她的一次悲惨的 OI 比赛经历.那是一道基础的树状数组题. 给出一个长度为 $n$ 的数组 $A$,初始值都为 $0$,接下来 ...

  4. loj #535. 「LibreOJ Round #6」花火 树状数组求逆序对+主席树二维数点+整体二分

    $ \color{#0066ff}{ 题目描述 }$ 「Hanabi, hanabi--」 一听说祭典上没有烟火,Karen 一脸沮丧. 「有的哦-- 虽然比不上大型烟花就是了.」 还好 Shinob ...

  5. 「CodePlus 2017 11 月赛」Yazid 的新生舞会(树状数组/线段树)

    学习了新姿势..(一直看不懂大爷的代码卡了好久T T 首先数字范围那么小可以考虑枚举众数来计算答案,设当前枚举到$x$,$s_i$为前$i$个数中$x$的出现次数,则满足$2*s_r-r > 2 ...

  6. 「CodePlus 2017 12 月赛」火锅盛宴(模拟+树状数组)

    1A,拿来练手的好题 用一个优先队列按煮熟时间从小到大排序,被煮熟了就弹出来. 用n个vector维护每种食物的煮熟时间,显然是有序的. 用树状数组维护每种煮熟食物的数量. 每次操作前把优先队列里煮熟 ...

  7. 【loj6041】「雅礼集训 2017 Day7」事情的相似度 后缀自动机+STL-set+启发式合并+离线+扫描线+树状数组

    题目描述 给你一个长度为 $n$ 的01串,$m$ 次询问,每次询问给出 $l$ .$r$ ,求从 $[l,r]$ 中选出两个不同的前缀的最长公共后缀长度的最大值. $n,m\le 10^5$ 题解 ...

  8. 「模拟赛20180306」回忆树 memory LCA+KMP+AC自动机+树状数组

    题目描述 回忆树是一棵树,树边上有小写字母. 一次回忆是这样的:你想起过往,触及心底--唔,不对,我们要说题目. 这题中我们认为回忆是这样的:给定 \(2\) 个点 \(u,v\) (\(u\) 可能 ...

  9. 「模拟赛20180307」三元组 exclaim 枚举+树状数组

    题目描述 给定 \(n,k\) ,求有多少个三元组 \((a,b,c)\) 满足 \(1≤a≤b≤c≤n\)且\(a + b^2 ≡ c^3\ (mod\ k)\). 输入 多组数据,第一行数据组数\ ...

随机推荐

  1. Go语言核心36讲(新年彩蛋)--学习笔记

    新年彩蛋 | 完整版思考题答案 基础概念篇 Go 语言在多个工作区中查找依赖包的时候是以怎样的顺序进行的? 答:你设置的环境变量GOPATH的值决定了这个顺序.如果你在GOPATH中设置了多个工作区, ...

  2. JVM 内存布局

    JVM 内存布局规定了 Java 在运行过程中内存申请.分配.管理的策略,保证了 JVM 的高效稳定运行. 线程是否共享 Heap (堆区) 堆是 OOM 故障最主要的发生区域.它是内存区域中最大的一 ...

  3. BL8810|USB 2.0单芯片解决方案闪存读卡器|BL8810替代GL823K

    创惟GL823K是一款USB 2.0单LUN读卡器控制器,可支持SD/MMC/MSPRO闪存卡.它支持USB 2.0高速传输,将Digital TM(SD).SDHC.SDXC.Mini DTM.Mi ...

  4. Capstone CS5267|CS5267参数|CS5267规格书

    CS5267 USB Type-C to HDMI2.0b 4k@60Hz Converter with PD3.0 Support 1.CS5267概述 Capstone CS5267是一款高性能T ...

  5. 初识python之词组截取及翻译

    d = {} k = [] v = [] with open('dir','r',encoding='utf-8') as f: for i in f.readlines(): j = i.strip ...

  6. scp 文件传输

    1.推送 scp -r imageAPP/ root@ip:/data/soft/ 本地当前文件下的 imageAPP文件 推送到ip服务器 /data/soft/ 目录下 2.拉取 scp -r r ...

  7. nginx - win系统启动一闪而过 ,服务没有启动成功

    这种现象是因为配置文件里配置的服务监听端口被占了

  8. js几种常见排序的实现

    1. 冒泡排序 定义: 比较相邻的前后二个数据,如果前面数据大于后面的数据,就将二个 数据交换. 这样对数组的第0个数据到N-1个数据进行一次遍历后,最大的一个数据就"沉"到数组第 ...

  9. SQL高级优化系列

    目录 SQL高级优化系列(一)之MySQL优化 SQL高级优化系列(二)之MySQL架构 SQL高级优化系列(三)之存储引擎 SQL高级优化系列(四)之SQL优化 SQL高级优化系列(五)之执行计划 ...

  10. Static 静态+this

    (一):静态 1.Static修饰的都是静态的,都是类相关的,不需要new对象,直接采用类名.的方式访问 2.当一个属性是类级别的,所有对象的这个属性都是一样的,直接定义为静态 类=属性+方法 属性描 ...