Luogu5280 ZJOI2019线段树(线段树)
容易发现相当于求2m种操作序列所得的每种线段树tag数量之和。显然考虑每个点的贡献,也即有多少种方案会使该点上有tag。可以将点分为四类:
1.修改时被经过且有儿子被修改的节点
2.修改时被经过且没有儿子被修改的节点
3.修改时未被经过且有兄弟被修改的节点
4.修改时未被经过且没有兄弟被修改的节点
考虑一次修改对这些点的影响。设该次操作之前已进行p次修改。易得:
1.方案数不变
2.方案数+2p
3.方案数+修改前该点到根路径(包括该点自身)上有tag的方案数
4.方案数*2
只要考虑维护3类点所需的东西即可。对此易得:
1.方案数不变
2.方案数+2p
3.方案数*2
4.与其最近的非4类点祖先相同
维护两棵线段树即可。其中第一棵可以转化为维护概率来避免打lazy标记,当然并不重要。
注意维护3类点取出方案数时需要pushdown,并且对此处的pushdown需要小心是否有写丑导致四倍空间不够用。
upd:上一句可以划掉,突然发现自己莫名其妙的把一个log硬是写成两个log了……两棵线段树看成一棵完全没问题,顺便pushdown一下就好。
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define P 998244353
#define N 100010
#define inv 499122177
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
int x=0,f=1;char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
return x*f;
}
int n,m,tree[N<<2],p[N],p2[N],stk[N],ans;
int val[N<<2],taga[N<<2],tagm[N<<2];
void del(int x){ans=(ans+P-x)%P;}
void ins(int x){ans=(ans+x)%P;}
void updatea(int k,int x)
{
val[k]=(val[k]+x)%P;
taga[k]=(taga[k]+x)%P;
}
void downa(int k)
{
updatea(k<<1,taga[k]);
updatea(k<<1|1,taga[k]);
taga[k]=0;
}
void updatem(int k,int x)
{
val[k]=1ll*val[k]*x%P;
tagm[k]=1ll*tagm[k]*x%P;
taga[k]=1ll*taga[k]*x%P;
}
void downm(int k)
{
updatem(k<<1,tagm[k]);
updatem(k<<1|1,tagm[k]);
tagm[k]=1;
}
void add(int k,int l,int r,int x,int y,int p)
{
if (l==x&&r==y) {updatea(k,p);return;}
if (tagm[k]!=1) downm(k);
if (taga[k]!=0) downa(k);
int mid=l+r>>1;
if (y<=mid)
{
updatem(k<<1|1,2);
add(k<<1,l,mid,x,y,p);
}
else if (x>mid)
{
updatem(k<<1,2);
add(k<<1|1,mid+1,r,x,y,p);
}
else add(k<<1,l,mid,x,mid,p),add(k<<1|1,mid+1,r,mid+1,y,p);
}
int query(int k)
{
int top=0;
for (int i=k;i;i>>=1) stk[++top]=i;
for (int i=top;i>1;i--)
{
if (tagm[stk[i]]!=1) downm(stk[i]);
if (taga[stk[i]]!=0) downa(stk[i]);
}
return val[k];
}
void modify(int k,int l,int r,int x,int y,int p)
{
if (l==x&&r==y)
{
del(tree[k]);
tree[k]=1ll*(tree[k]+1)*inv%P;
ins(tree[k]);
return;
}
if (tagm[k]!=1) downm(k);
if (taga[k]!=0) downa(k);
del(tree[k]);
tree[k]=1ll*tree[k]*inv%P;
ins(tree[k]);
int mid=l+r>>1;
if (y<=mid)
{
del(tree[k<<1|1]);
tree[k<<1|1]=1ll*(tree[k<<1|1]+1ll*val[k<<1|1]*p%P)*inv%P;
ins(tree[k<<1|1]);
modify(k<<1,l,mid,x,y,p);
}
else if (x>mid)
{
del(tree[k<<1]);
tree[k<<1]=1ll*(tree[k<<1]+1ll*val[k<<1]*p%P)*inv%P;
ins(tree[k<<1]);
modify(k<<1|1,mid+1,r,x,y,p);
}
else modify(k<<1,l,mid,x,mid,p),modify(k<<1|1,mid+1,r,mid+1,y,p);
}
void fill(int k,int l,int r){tagm[k]=1;if (l==r) return;fill(k<<1,l,l+r>>1);fill(k<<1|1,(l+r>>1)+1,r);}
int main()
{
#ifndef ONLINE_JUDGE
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
const char LL[]="%I64d\n";
#else
const char LL[]="%lld\n";
#endif
n=read(),m=read();
p[0]=1;for (int i=1;i<=m;i++) p[i]=(p[i-1]<<1)%P;
p2[0]=1;for (int i=1;i<=m;i++) p2[i]=1ll*p2[i-1]*inv%P;
int cnt=0;fill(1,1,n);
for (int _=1;_<=m;_++)
{
int op=read();
if (op==1)
{
int l=read(),r=read();
modify(1,1,n,l,r,p2[cnt]);
add(1,1,n,l,r,p[cnt]);
cnt++;
}
else printf("%d\n",1ll*ans*p[cnt]%P);
}
return 0;
}
Luogu5280 ZJOI2019线段树(线段树)的更多相关文章
- HDU 5877     dfs+ 线段树(或+树状树组)
		
1.HDU 5877 Weak Pair 2.总结:有多种做法,这里写了dfs+线段树(或+树状树组),还可用主席树或平衡树,但还不会这两个 3.思路:利用dfs遍历子节点,同时对于每个子节点au, ...
 - 学习笔记--函数式线段树(主席树)(动态维护第K极值(树状数组套主席树))
		
函数式线段树..资瓷 区间第K极值查询 似乎不过似乎划分树的效率更优于它,但是如果主席树套树状数组后,可以处理动态的第K极值.即资瓷插入删除,划分树则不同- 那么原理也比较易懂: 建造一棵线段树(权值 ...
 - BZOJ_3196_二逼平衡树_(树套树,线段树+Treap)
		
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=3196 可以处理区间问题的平衡树. 3196: Tyvj 1730 二逼平衡树 Time Lim ...
 - [BZOJ 1901] Dynamic Rankings 【树状数组套线段树 || 线段树套线段树】
		
题目链接:BZOJ - 1901 题目分析 树状数组套线段树或线段树套线段树都可以解决这道题. 第一层是区间,第二层是权值. 空间复杂度和时间复杂度均为 O(n log^2 n). 线段树比树状数组麻 ...
 - BZOJ 3110 ZJOI 2013 K大数查询 树套树(权值线段树套区间线段树)
		
题目大意:有一些位置.这些位置上能够放若干个数字. 如今有两种操作. 1.在区间l到r上加入一个数字x 2.求出l到r上的第k大的数字是什么 思路:这样的题一看就是树套树,关键是怎么套,怎么写.(话说 ...
 - 归并树  划分树  可持久化线段树(主席树) 入门题 hdu 2665
		
如果题目给出1e5的数据范围,,以前只会用n*log(n)的方法去想 今天学了一下两三种n*n*log(n)的数据结构 他们就是大名鼎鼎的 归并树 划分树 主席树,,,, 首先来说两个问题,,区间第k ...
 - HDOJ 4417 - Super Mario 线段树or树状数组离线处理..
		
题意: 同上 题解: 抓着这题作死的搞~~是因为今天练习赛的一道题.SPOJ KQUERY.直到我用最后一种树状数组通过了HDOJ这题后..交SPOJ的才没超时..看排名...时间能排到11名了..有 ...
 - hdu 4836 The Query on the Tree(线段树or树状数组)
		
The Query on the Tree Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Ot ...
 - 【BZOJ3295】动态逆序对(线段树,树状数组)
		
[BZOJ3295]动态逆序对(线段树,树状数组) 题面 Description 对于序列A,它的逆序对数定义为满足iAj的数对(i,j)的个数.给1到n的一个排列,按照某种顺序依次删除m个元素,你的 ...
 - 线段树(单标记+离散化+扫描线+双标记)+zkw线段树+权值线段树+主席树及一些例题
		
“队列进出图上的方向 线段树区间修改求出总量 可持久留下的迹象 我们 俯身欣赏” ----<膜你抄> 线段树很早就会写了,但一直没有总结,所以偶尔重写又会懵逼,所以还是要总结一下. ...
 
随机推荐
- Nodejs mongoose 详解
			
前言 Mongoose 是在nodejs环境下,对mongodb进行便捷操作的对象模型工具.本文介绍解(翻)密(译)Mongoose插件. Schema 开始我们就要讲到Schema,一个Schema ...
 - JVM参数配置  java内存区域
			
java内存区域 一些基本概念 http://www.importnew.com/18694.html https://www.cnblogs.com/wangyayun/p/6557851.html ...
 - python3 pip 安装Scrapy在win10 安装报错error: Microsoft Visual C++ 14.0 is required. Get it with "Microsoft Visual C++ Build Tools": http://landinghub.visualstudio.com/visual-cpp-build-tools
			
问题描述 当前环境win10,python_3.6.1,64位. 在windows下,在dos中运行pip install Scrapy报错: building 'twisted.test.raise ...
 - H5 audio标签
			
37-audio标签 注意点: audio标签的使用和video标签的使用基本一样, video中能够使用的属性在audio标签中大部分都能够使用, 并且功能都一样 只不过有3个属性不能用, heig ...
 - 闽江学院软件学院2016级JAVA构建之法-学生自学兴趣小组招募通知
			
为提升我2016级学生提升JAVA软件开发学习氛围,鼓励更多同学通过自学.团队学习.在线(社区)学习等方式学习并掌握JAVA课程,尤其是鼓励同学们通过微软中国邹欣老师所倡导的"构建之法&qu ...
 - JS典记
			
var href = ""; //遍历a标签 $ ( "a"). each (function () { href = ...
 - org.apache.ibatis.exceptions.PersistenceException:  ### Error querying database.  Cause: com.mysql.jdbc.exceptions.jdbc4.MySQLSyntaxErrorException: You have an error in your SQL syntax; check the manu
			
这个是sql 语句 错误 仔细检查 SQL语句是否写错了 org.apache.ibatis.exceptions.PersistenceException: ### Error queryi ...
 - asp.net mvc 自定义全局过滤器 验证用户是否登录
			
一般具有用户模块的系统都需要对用户是否登录进行验证,如果用户登录了就可以继续操作,否则退回用户的登录页面 对于这样的需求我们可以通过自定义一个独立的方法来完成验证的操作,但是这样代码的重复率就大大提高 ...
 - PV、TPS、QPS计算公式(转)
			
英文解释: PV=page viewTPS=transactions per secondQPS=queries per secondRPS=requests per second RPS=并发数/平 ...
 - C#复习笔记(3)--C#2:解决C#1的问题(泛型)
			
这一章会描述在C#2中所做的主要的变化 泛型 泛型的概念中包含类型参数和类型实参,类型参数相当于类型实参的蓝图. 泛型类型分为未绑定泛型类型和已构造泛型类型.已构造泛型类型又分为开放的泛型类型和封闭的 ...