当前位置: 首页 > news >正文

【Luogu】 P5445 [APIO2019] 路灯

题目链接

点击打开链接

题目解法

转化很妙

考虑关路灯 x x x 的操作
令左边第一个未关的路灯为 L L L,右边第一个未关的路灯为 R R R,那么这一次会影响的区间即为 l ∈ [ L + 1 , x ] , r ∈ [ x , R − 1 ] l\in[L+1,x],\;r\in[x,R-1] l[L+1,x],r[x,R1],既然有 2 个限制,那么我们把限制变成二维平面上的矩阵,即把左上为 ( L + 1 , x ) (L+1,x) (L+1,x),右下为 ( x , R − 1 ) (x,R-1) (x,R1) 的矩阵变为 0
开路灯的操作类似

考虑询问的是什么?
1 1 1 q − 1 q-1 q1 次每个版本 ( x , y ) (x,y) (x,y) 位置的和
考虑给矩阵定义一个增加(或减少)值,使查询变成只要查询第 i i i 个版本的 ( x , y ) (x,y) (x,y)
这里直接给出:对于第 i i i 个操作,加上或减去值为 i − q i-q iq
考虑先开路灯,再关路灯的操作,即为 T − i − ( T − j ) = j − i T-i-(T-j)=j-i Ti(Tj)=ji,恰好为答案,这可以拓展到多个开关路灯的操作
这里需要一个特判,如果查询第 i i i 个版本从 x x x 可以到 y y y,那么答案需要 − ( q − i ) -(q-i) (qi)

可以把矩形变成 4 个区域的加减,然后不难发现,这就是一个三维偏序,可以考虑离线 c d q cdq cdq
时间复杂度 O ( n l o g 2 n ) O(nlog^2n) O(nlog2n)

#include <bits/stdc++.h>
#define lowbit(x) x&-x
using namespace std;
const int N=300100;
struct Node{ int x1,x2,v,op,id;}q[N*10],t[N*10];
int idx,n,m,state[N],ans[N],tr[N];
char str[N];
set<int> se;
inline int read(){int FF=0,RR=1;char ch=getchar();for(;!isdigit(ch);ch=getchar()) if(ch=='-') RR=-1;for(;isdigit(ch);ch=getchar()) FF=(FF<<1)+(FF<<3)+ch-48;return FF*RR;
}
void ADD(int x,int y,int v){ q[++idx]={x,y,v,0,0};}
void add_square(int x1,int x2,int y1,int y2,int v){ADD(x2,y2,v);if(x1>1) ADD(x1-1,y2,-v);if(y1>1) ADD(x2,y1-1,-v);if(x1>1&&y1>1) ADD(x1-1,y1-1,v);
}
void work(int x,int t){if(!state[x]){auto it=se.upper_bound(x);int R=*it,L=*(--it);se.insert(x);add_square(L+1,x,x,R-1,-(m-t));}else{se.erase(x);auto it=se.upper_bound(x);int R=*it,L=*(--it);add_square(L+1,x,x,R-1,m-t);}
}
void add(int x,int y){ for(;x<=n;x+=lowbit(x)) tr[x]+=y;}
int ask(int x){if(!x) return 0;int res=0;for(;x;x-=lowbit(x)) res+=tr[x];return res;
}
void cdq(int l,int r){if(l==r) return;int mid=(l+r)>>1;cdq(l,mid),cdq(mid+1,r);int i=l,j=mid+1,k=0;while(i<=mid||j<=r){if(j>r||(i<=mid&&q[i].x1>=q[j].x1)){if(!q[i].op) add(q[i].x2,q[i].v);t[++k]=q[i],i++;}else{if(q[j].op) ans[q[j].id]+=ask(n)-ask(q[j].x2-1);t[++k]=q[j],j++;}}for(i=l;i<=mid;i++) if(!q[i].op) add(q[i].x2,-q[i].v);for(int i=1,j=l;i<=k;i++,j++) q[j]=t[i];
}
bool con(int x,int y){ return ask(y)-ask(x-1)==y-x+1;}
int main(){n=read(),m=read();scanf("%s",str+1);for(int i=1;i<=n;i++) state[i]=str[i]-48;add_square(1,n,1,n,m);se.insert(0),se.insert(n+1);for(int i=1;i<=n;i++) if(state[i]) add(i,1);for(int i=1;i<=n;i++) if(!state[i]) work(i,0);for(int i=1;i<=m;i++){char op[10];scanf("%s",op);if(op[0]=='t'){int x=read();state[x]^=1,work(x,i);ans[i]=-1e9;if(state[x]) add(x,1);else add(x,-1);}else{int x=read(),y=read();y--;if(con(x,y)) ans[i]=-(m-i);q[++idx]={x,y,0,1,i};}}memset(tr,0,sizeof(tr));cdq(1,idx);for(int i=1;i<=m;i++) if(ans[i]!=-1e9) printf("%d\n",ans[i]);return 0;
}
http://www.lryc.cn/news/158206.html

相关文章:

  • Kafka3.0.0版本——消费者(独立消费者消费某一个主题中某个分区数据案例__订阅分区)
  • 基于Simulink的用于电力系统动态分析
  • 日200亿次调用,喜马拉雅网关的架构设计
  • 构造函数和析构函数(个人学习笔记黑马学习)
  • GPT引领前沿与应用突破之GPT4科研实践技术与AI绘图教程
  • Git上传新项目
  • C语言文件操作总结
  • 原生js之dom如何进行事件监听(事件捕获/冒泡)
  • 使用SimPowerSystems并网光伏阵列研究(Simulink实现)
  • BUUCTF-WEB-[ACTF2020 新生赛]Includel
  • 算法通关村十四关:白银挑战-堆能高效解决的经典问题
  • 跨站请求伪造(CSRF)攻击与防御原理
  • 从0到1实现播放控制器
  • 【Vue-Element-Admin】导出el-table全部数据
  • MFC 更改控件的大小和位置
  • 【向量数据库】相似向量检索Faiss数据库的安装及余弦相似度计算(C++)
  • 教育培训小程序的设计与功能解析
  • 【ES】illegal_argument_exception“,“reason“:“Result window is too large
  • SpringBoot实现登录拦截
  • 浅谈泛在电力物联网、能源互联网与虚拟电厂
  • 深度学习框架安装与配置指南:PyTorch和TensorFlow详细教程
  • vue中属性执行顺序
  • 【代码随想录】Day 50 动态规划11 (买卖股票Ⅲ、Ⅳ)
  • PHP反序列化漏洞
  • 容器编排学习(一)k8s集群管理
  • js去除字符串空格的几种方式
  • Spring 自带工具——URI 工具UriComponentsBuilder
  • 优化案例5:视图目标列改写优化
  • Origin绘制彩色光谱图
  • 项目复盘:从实践中学习