UOJ#291. 【ZJOI2017】树状数组 树套树
原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ291.html
题解
结论:这个写错的树状数组支持的是后缀加和后缀求和。这里的后缀求和在 x = 0 的时候比较特殊,返回 0 。
于是我们需要查询 v[L-1] 和 v[R] 相同的概率是多少。
我们可以用树套树维护一下左端点在一个区间内,右端点在另一个区间内的修改操作使得对应点发生变化的概率。
注意 L = 1 要特判。
代码
#include <bits/stdc++.h>
#define clr(x) memset(x,0,sizeof (x))
#define For(i,a,b) for (int i=a;i<=b;i++)
#define Fod(i,b,a) for (int i=b;i>=a;i--)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define _SEED_ ('C'+'L'+'Y'+'A'+'K'+'I'+'O'+'I')
#define outval(x) printf(#x" = %d\n",x)
#define outvec(x) printf("vec "#x" = ");for (auto _v : x)printf("%d ",_v);puts("")
#define outtag(x) puts("----------"#x"----------")
using namespace std;
typedef long long LL;
LL read(){
LL x=0,f=0;
char ch=getchar();
while (!isdigit(ch))
f|=ch=='-',ch=getchar();
while (isdigit(ch))
x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
return f?-x:x;
}
const int N=100005,mod=998244353;
int Pow(int x,int y){
int ans=1;
for (;y;y>>=1,x=(LL)x*x%mod)
if (y&1)
ans=(LL)ans*x%mod;
return ans;
}
int Pmul(int x,int y){
return ((LL)x*(1-y)+(LL)y*(1-x))%mod;
}
int n,m;
namespace Seg{
const int S=N*305;
int ls[S],rs[S],cnt=0;
int v1[S],v2[S];
void update(int &rt,int L,int R,int x,int v){
if (!rt)
rt=++cnt,v1[rt]=v2[rt]=0;
v1[rt]=Pmul(v1[rt],v);
v2[rt]=Pmul(v2[rt],v*2%mod);
if (L==R)
return;
int mid=(L+R)>>1;
if (x<=mid)
update(ls[rt],L,mid,x,v);
else
update(rs[rt],mid+1,R,x,v);
}
int Q1(int rt,int L,int R,int xL,int xR){
if (!rt||xR<L||R<xL)
return 0;
if (xL<=L&&R<=xR)
return v1[rt];
int mid=(L+R)>>1;
return Pmul(Q1(ls[rt],L,mid,xL,xR),Q1(rs[rt],mid+1,R,xL,xR));
}
int Q2(int rt,int L,int R,int xL,int xR){
if (!rt||xR<L||R<xL)
return 0;
if (xL<=L&&R<=xR)
return v2[rt];
int mid=(L+R)>>1;
return Pmul(Q2(ls[rt],L,mid,xL,xR),Q2(rs[rt],mid+1,R,xL,xR));
}
}
int root[N<<2];
void upd(int rt,int L,int R,int x,int y,int v){
Seg::update(root[rt],1,n,y,v);
if (L==R)
return;
int mid=(L+R)>>1,ls=rt<<1,rs=ls|1;
if (x<=mid)
upd(ls,L,mid,x,y,v);
else
upd(rs,mid+1,R,x,y,v);
}
int Q1(int rt,int L,int R,int xL,int xR,int yL,int yR){
if (xL>xR||yL>yR||R<xL||L>xR)
return 0;
if (xL<=L&&R<=xR)
return Seg::Q1(root[rt],1,n,yL,yR);
int mid=(L+R)>>1,ls=rt<<1,rs=ls|1;
return Pmul(Q1(ls,L,mid,xL,xR,yL,yR),Q1(rs,mid+1,R,xL,xR,yL,yR));
}
int Q2(int rt,int L,int R,int xL,int xR,int yL,int yR){
if (xL>xR||yL>yR||R<xL||L>xR)
return 0;
if (xL<=L&&R<=xR)
return Seg::Q2(root[rt],1,n,yL,yR);
int mid=(L+R)>>1,ls=rt<<1,rs=ls|1;
return Pmul(Q2(ls,L,mid,xL,xR,yL,yR),Q2(rs,mid+1,R,xL,xR,yL,yR));
}
int main(){
n=read(),m=read();
int cnt=0;
while (m--){
int type=read(),L=read(),R=read();
if (type==1){
cnt^=1;
upd(1,1,n,L,R,Pow(R-L+1,mod-2));
}
else {
L--;
int vL=Q1(1,1,n,1,L,L,R-1);
int vR=Q1(1,1,n,L+1,R,R,n);
int vm=Q2(1,1,n,1,L,R,n);
int ans=Pmul(vm,Pmul(vL,vR));
if (L||!cnt)
ans=(1LL-ans)%mod;
ans=(ans+mod)%mod;
printf("%d\n",ans);
}
}
return 0;
}
UOJ#291. 【ZJOI2017】树状数组 树套树的更多相关文章
- 【BZOJ3196】二逼平衡树(树状数组,线段树)
[BZOJ3196]二逼平衡树(树状数组,线段树) 题面 BZOJ题面 题解 如果不存在区间修改操作: 搞一个权值线段树 区间第K大--->直接在线段树上二分 某个数第几大--->查询一下 ...
- POJ 1195 Mobile phones (二维树状数组或线段树)
偶然发现这题还没A掉............速速解决了............. 树状数组和线段树比较下,线段树是在是太冗余了,以后能用树状数组还是尽量用......... #include < ...
- 树状数组-HDU1541-Stars一维树状数组 POJ1195-Mobile phones-二维树状数组
树状数组,学长很早之前讲过,最近才重视起来,enmmmm... 树状数组(Binary Indexed Tree(B.I.T), Fenwick Tree)是一个查询和修改复杂度都为log(n)的数据 ...
- HDU 5618 Jam's problem again(三维偏序,CDQ分治,树状数组,线段树)
Jam's problem again Time Limit: 5000/2500 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Othe ...
- HDU 5877 2016大连网络赛 Weak Pair(树状数组,线段树,动态开点,启发式合并,可持久化线段树)
Weak Pair Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others) Tota ...
- Codeforces 777E(离散化+dp+树状数组或线段树维护最大值)
E. Hanoi Factory time limit per test 1 second memory limit per test 256 megabytes input standard inp ...
- st表、树状数组与线段树 笔记与思路整理
已更新(2/3):st表.树状数组 st表.树状数组与线段树是三种比较高级的数据结构,大多数操作时间复杂度为O(log n),用来处理一些RMQ问题或类似的数列区间处理问题. 一.ST表(Sparse ...
- Codeforces Round #225 (Div. 1) C. Propagating tree dfs序+ 树状数组或线段树
C. Propagating tree Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/contest/383/p ...
- HYSBZ - 3813 奇数国 欧拉函数+树状数组(线段树)
HYSBZ - 3813奇数国 中文题,巨苟题,巨无敌苟!!首先是关于不相冲数,也就是互质数的处理,欧拉函数是可以求出互质数,但是这里的product非常大,最小都2100000,这是不可能实现的.所 ...
- LOJ2251 [ZJOI2017] 树状数组【线段树】【树套树】
题目分析: 对于一个$add$操作,它的特点是与树状数组的查询相同,会给$1$到它自己产生影响,而$query$操作则会途径所有包含它的树状数组点.现在$add$操作具有前向性(不会影响之后的点).所 ...
随机推荐
- 商务电话思维图(XMind für Geschäftliche Telefongespräche)
在和德国人打交道时,经常会遇到打电话的情景,应该怎么应对呢?不用担心,记住下面这个导图,轻松搞定德语电话的常用句型. 最后,按照惯例,来张美景.人越是上了年纪,就活的越是小心.但无论外界怎么样,请保持 ...
- 贝叶斯推断 && 概率编程初探
1. 写在之前的话 0x1:贝叶斯推断的思想 我们从一个例子开始我们本文的讨论.小明是一个编程老手,但是依然坚信bug仍有可能在代码中存在.于是,在实现了一段特别难的算法之后,他开始决定先来一个简单的 ...
- tensorflow Pipeline 之TextLineReader 和decode_csv多分割替代方案
1.TextLineReader reader = tf.TextLineReader() key, value = reader.read(filename_queue) record_defaul ...
- VMware 安装Linux系统 CentOS
VMware 安装Linux系统 CentOS 1. 下载镜像系统 centos镜像下载地址:https://www.centos.org/download/ 选择DVD下载即可 linux各版本下 ...
- NCE损失(Noise-Constrastive Estimation Loss)
1.算法概述 假设X是从真实的数据(或语料库)中抽取的样本,其服从一个相对可参考的概率密度函数P(d),噪音样本Y服从概率密度函数为P(n),噪音对比估计(NCE)就是通过学习一个分类器把这两类样本区 ...
- Spring cloud gateway
==================================为什么需要API gateway?==================================企业后台微服务互联互通, 因为 ...
- linux oops调试
参考文章: arm 指令定位错误 https://blog.csdn.net/songcdut/article/details/41383483 linux mips指令学习 https://www. ...
- 关于vscode插件 的一些体验
一 vsCode git 首先vscode 继承了 git 在使用git前 先手动创建一个文件夹 用来clone已有项目 然后将 clone下来的项目放入工作区 右上角的小转转就是 pull文件 更 ...
- 【转载】大白话Docker入门(二)
原文:https://yq.aliyun.com/articles/63517?spm=a2c4e.11153940.blogcont63035.15.12011c3fddklk0 上篇的大白话Doc ...
- 读spring源码(二)-XmlBeanDefinitionReader-解析BeanDefinition
上次说到ApplicationContext加载BeanDefinition时会创建一个XmlBeanDefinitionReader,将XML解析.BeanDefinition加载委托给XmlBea ...