传送门·

对于询问 $(a,b)$ ,感觉一维很不好维护,考虑把询问看成平面上的一个点,坐标为 $(a,b)$

每个坐标 $(x,y)$ 的值表示到当前 $x$ 和 $y$ 联通的时间和

考虑一个修改的贡献,它其实就是把左边一段区间 $[l,x]$ 和右边一段区间 $[x+1,r]$ 联通或断开

放到平面上发现其实就是横坐标在 $[l,x]$ ,纵坐标在 $[x+1,r]$ 的矩形里修改,那么矩形左下角为 $[l,x+1]$,右上角为 $[x,r]$

如果每个时间点都把相应矩形 $+1$ 的话显然是不可行的,考虑起点和终点的时间差

如果当前操作是联通,设当前时间为 $t$,则把相应矩形 $-t$,如果是断开则把矩形加 $t$

这样我们询问时直接矩形单点查值即可,但是要注意,如果查询时当前区间仍然联通,那么对应矩形还没加 $t$ ,我们查完值以后答案还要再加 $t$

矩形加具体就是差分成四个前缀修改,查询的时候直接查左下所有位置的和即可

如果把左下角为 $(xa,ya)$,右上角 $(xb,yb)$ 的闭区间都加一个 $v$,那么其实就是把 $(0,0),(xa,ya)$ 加一个 $v$,$(0,0),(xa,yb+1)$ 减 $v$,$(0,0),(xb+1,ya)$ 减 $v$,$(0,0),(xb+1,yb+1)$ 加 $v$

二维平面加一维时间轴, $CDQ$ 分治即可

然后就是具体的维护了,我维护亮灯区间用的是闭区间,好像比较多细节...

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<set>
using namespace std;
typedef long long ll;
inline int read()
{
int x=,f=; char ch=getchar();
while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
return x*f;
}
const int N=6e5+;
char s[N];
int n,Q,tot;
struct dat{//区间
int l,r;
dat (int a=,int b=) { l=a,r=b; }
inline bool operator < (const dat &tmp) const {
return l<tmp.l;
}
};
set <dat> S;//维护开灯的区间
set <dat>::iterator ita,itb;
struct Poi{//CDQ分治的点
int x,y,id,v;//id=0表示是修改,v是修改的值
Poi (int a=,int b=,int c=,int d=) { x=a,y=b,id=c,v=d; }
}d[N<<],tmp[N<<];
int t[N];
inline void add(int x,int v) { while(x<=n+) t[x]+=v,x+=x&-x; }
inline int ask(int x) { int res=; while(x) res+=t[x],x-=x&-x; return res; }
inline bool fc(Poi &a,Poi &b) { if(a.x!=b.x) return a.x<b.x; return a.y!=b.y ? a.y<b.y : a.id<b.id; }
int ans[N];
void solve(int l,int r)
{
if(l==r) return; int mid=l+r>>;
solve(l,mid); solve(mid+,r);
int i=l,j=mid+,p=l-;
while(i<=mid&&j<=r)
{
if(fc(d[i],d[j]))
{
if(!d[i].id) add(d[i].y,d[i].v);
tmp[++p]=d[i++]; continue;
}
if(d[j].id) ans[d[j].id]+=ask(d[j].y);
tmp[++p]=d[j++];
}
while(i<=mid) { if(!d[i].id) add(d[i].y,d[i].v); tmp[++p]=d[i++]; }
while(j<=r) { if(d[j].id) ans[d[j].id]+=ask(d[j].y); tmp[++p]=d[j++]; }
for(int k=l;k<=mid;k++) if(!d[k].id) add(d[k].y,-d[k].v);//撤销
for(int k=l;k<=r;k++) d[k]=tmp[k];
}
int main()
{
n=read(),Q=read();
scanf("%s",s+);
S.insert(dat(,)); S.insert(dat(n+,n+));//初始虚区间防越界
for(int i=;i<=n;i++)
{
if(s[i]=='') continue;
int R=i; while(s[R+]=='') R++;
S.insert(dat(i,R)); i=R;//初始情况
}
char ch[]; int a,b,cnt=;
for(int i=;i<=Q;i++)
{
scanf("%s",ch); a=read();
if(ch[]=='q')
{
b=read();
d[++tot]=Poi(a,b,++cnt,);//询问点
ita=S.lower_bound(dat(a+,)); ita--;
itb=S.lower_bound(dat(b,)); itb--;
if( ita==itb && (*ita).r>=b- ) ans[cnt]+=i;//判断区间是否还联通
continue;
}
if(s[a]=='')//联通
{
ita=S.lower_bound(dat(a+,)); ita--; itb=ita; ita++;//找到左右区间
dat A=(*ita),B=(*itb); int l=a,r=a;
if(B.r==l-&&B.l>=) l=B.l,S.erase(B);//注意判断,注意边界
if(A.l==r+&&A.r<=n) r=A.r,S.erase(A);
S.insert(dat(l,r));
d[++tot]=Poi(l,a+,,-i); d[++tot]=Poi(a+,r+,,-i);
d[++tot]=Poi(l,r+,,i); d[++tot]=Poi(a+,a,,i);
//联通时矩形-i
s[a]=''; continue;
}
//断开
ita=S.lower_bound(dat(a+,)); ita--; dat A=(*ita);
S.erase(A); if(A.l<a) S.insert(dat(A.l,a-)); if(A.r>a) S.insert(dat(a+,A.r));//注意边界
d[++tot]=Poi(A.l,a+,,i); d[++tot]=Poi(a+,A.r+,,i);
d[++tot]=Poi(A.l,A.r+,,-i); d[++tot]=Poi(a+,a,,-i);
//断开时矩形+i
s[a]='';
}
solve(,tot);
for(int i=;i<=cnt;i++) printf("%d\n",ans[i]);
return ;
}

P5445 [APIO2019]路灯的更多相关文章

  1. P5445 [APIO2019]路灯(树套树)

    P5445 [APIO2019]路灯 转化为平面上的坐标(x,y),set维护连续区间. 用树套树维护矩阵加法,单点查询. 注意维护矩阵差分的时候, $(x,y,v)$是对$(x,y)(n+1,n+1 ...

  2. 【LOJ#3146】[APIO2019]路灯(树套树)

    [LOJ#3146][APIO2019]路灯(树套树) 题面 LOJ 题解 考场上因为\(\text{bridge}\)某个\(\text{subtask}\)没有判\(n=1\)的情况导致我卡了\( ...

  3. LOJ3146 APIO2019路灯(cdq分治+树状数组)

    每个时刻都形成若干段满足段内任意两点可达.将其视为若干正方形.则查询相当于求历史上某点被正方形包含的时刻数量.并且注意到每个时刻只有O(1)个正方形出现或消失,那么求出每个矩形的出现时间和消失时间,就 ...

  4. 题解-APIO2019路灯

    problem \(\mathtt {loj-3146}\) 题意概要:一条直线上有 \(n+1\) 个点和 \(n\) 条道路,每条道路连通相邻两个点.在 \(q\) 个时刻内,每个时刻有如下两种操 ...

  5. APIO2019简要题解

    Luogu P5444 [APIO2019]奇怪装置 看到这种题,我们肯定会想到\((x,y)\)一定有循环 我们要找到循环节的长度 推一下发现\(x\)的循环节长为\(\frac{AB}{B+1}\ ...

  6. 2021record

    2021-10-14 P2577 [ZJOI2004]午餐 2021-10-13 CF815C Karen and Supermarket(小小紫题,可笑可笑) P6748 『MdOI R3』Fall ...

  7. P3722 [AH2017/HNOI2017]影魔(单调栈+扫描线+线段树)

    题面传送门 首先我们把这两个贡献翻译成人话: 区间 \([l,r]\) 产生 \(p_1\) 的贡献当且仅当 \(a_l,a_r\) 分别为区间 \([l,r]\) 的最大值和次大值. 区间 \([l ...

  8. 【APIO2019】路灯(ODT & (树套树 | CDQ分治))

    Description 一条 \(n\) 条边,\(n+1\) 个点的链,边有黑有白.若结点 \(a\) 可以到达 \(b\),需要满足 \(a\to b\) 的路径上的边不能有黑的.现给出 \(0\ ...

  9. [APIO2019] [LOJ 3146] 路灯 (cdq分治或树状数组套线段树)

    [APIO2019] [LOJ 3146] 路灯 (cdq分治或树状数组套线段树) 题面 略 分析 首先把一组询问(x,y)看成二维平面上的一个点,我们想办法用数据结构维护这个二维平面(注意根据题意这 ...

随机推荐

  1. 缓存算法LRU笔记

      LRU原理与分析 LRU是Least Recently Used 的缩写,翻译过来就是“最近最少使用”,也就是说,LRU缓存把最近最少使用的数据移除,让给最新读取的数据.而往往最常读取的,也是读取 ...

  2. 【机器人M号】题解

    题目 题目描述 3030年,Macsy正在火星部署一批机器人. 第1秒,他把机器人1号运到了火星,机器人1号可以制造其他的机器人. 第2秒,机器人1号造出了第一个机器人--机器人2号. 第3秒,机器人 ...

  3. vue引入css和js的方法

    1.引入js的话 require('地址') imgport xxx from ‘src’ import xxx 2.引入css @import ‘  ’

  4. jetSonNano darknet ubdefined reference to 'pow',undefined reference to 'sqrtf'....

    我在用CMakelist编译工程时,遇到了这个一连串基础数学函数找不到的问题,如下图所示: 我当时在工程中明明引用了 #include "math.h"头函数,这是因为你的工程在预 ...

  5. node的cropto加密

    本文转自https://blog.csdn.net/sinat_35670989/article/details/78224214 'use strict' //crypto(kri:pto)意为加密 ...

  6. os.environ.get()的用法

    os.environ.get()是python中os模块获取环境变量的一个方法 import os JS_ADDRESS = os.environ.get("PALM_JS_ADDRESS& ...

  7. element的隐藏组件滚动条el-scrollbar使用

    elementui中有个隐藏的组件,就是element官网使用的滚动条,tree 左右滑动滚动条 ①首先全局引入element,import ElementUI from 'element-ui'; ...

  8. grep的用法,小技巧,模板中含有\t时:grep -P "^\t" file

    linux中grep和find的用法区别 本文章详细的介绍了关于在linux中的grep和find两个命令的用法介绍,以及后面总结了它们两年用法区别哦. 先我们来介绍一下关于grep用法和一些小注意事 ...

  9. Apache2.4常用编译参数

    转载文章,亲试 Apache2.4以后的版本编译依赖apr,所以,编译之前需要先安装apr及apr-util. 编译参数只是一个参考作用,这个参数是我平常使用的,具体工作中的需求还是要区别对待的 一些 ...

  10. qbzt day3 上午

    内容提要 堆 lca(最近公共祖先) st表 hash 并查集 树状数组 线段树 数据结构 1.堆 Priority_queue 他滋兹:插入删除查询最大值(最小值) 分为大根堆小根堆 2.LCA 首 ...