[cf 997 E] Good Subsegments
(这是石神找到的一道hiao题。)
题意:
你有一个长度为n的排列,有Q组询问$[l,r]$,每次询问$[l,r]$的子区间中有多少是好的。
一个区间是好的区间当且仅当该区间中的元素在排序后是连续的。
$n,Q\leq 120000$。
题解:
为什么说这是一道hiao题呢?
像我一样的数据结构弱智选手看到之后可能会想一些树套树的操作。
但这个题有一种具有代表性的经典方法。
注意到一段区间是好的当且仅当$(max-min)-(r-l)=0$。
我们可以将将询问离线下来,每次处理以$r$为右端点的所有询问。
考虑从左往右枚举右端点$r$,并用一棵线段树维护每个点$l$的$(max-min)-(r-l)$的值。
我们不妨将以$r$为右端点时的这颗线段树称作版本$r$。
那么,以$r$为右端点的好区间个数就是版本$r$中0的个数。
以$r$为右端点,左端点在$[l,r]$间的好区间个数就是版本$r$中$[l,r]$中0的个数之和,记为$siz[r][l,r]$。
但是注意到,我们要求的答案其实是$\sum_{i=l}^{r}{siz[i][l,i]}$。
那么需要对线段树上的每个节点额外记录一个$sum$,代表这个区间所有版本中0的个数之和。
也就是${sum[l,r]}=\sum_{i=1}^{now}{siz[i][l,r]}$。
显然答案就是$sum[l,r]$,那么如何维护这个$sum$呢?
不妨维护一个类似$lazy$标记的$tag$标记,表示该区间现在的$siz$要对$sum$产生几次的贡献。
这个标记实际上就是把$min=0$时的$siz$一直继承下来,直到访问到它的时候再更新$siz$。
注意无论该区间现在的$min$值是不是等于0(可能被上面区间的$pushdown$改掉了),在访问到它之前都需要保留该标记,因为那个0的版本始终是有贡献的。
在$pushdown$这个$tag$的时候要注意:假如该区间当前的$min$值是由一个而不是两个子区间转移上来的,那我们只应该更新转移上来的那个子区间的$tag$。
如何判断哪个是转移上来的子区间?显然当且仅当$min_{son}=min_{now}$(注意$min_{now}$未必等于0)。
时间复杂度$O((n+Q)logn)$。
积累套路:对于询问区间子序列的信息,可以离线移动右指针,类似扫描线一样计算贡献。
代码:
#include<bits/stdc++.h>
#define maxn 120005
#define maxm 500005
#define inf 0x7fffffff
#define ll long long using namespace std;
ll N,Q,A[maxn],mn[maxn<<],ans[maxn];
ll sum[maxn<<],nd[maxn<<],lz[maxn<<];
ll siz[maxn<<],tim[maxn<<];
vector<pair<ll,ll> > ques[maxn];
struct node{ll v,l,r;};
stack<node> smx,smn; inline ll read(){
ll x=,f=; char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-;
for(;isdigit(c);c=getchar()) x=x*+c-'';
return x*f;
} inline void pushup(ll k){
if(mn[k<<]==mn[k<<|]) mn[k]=mn[k<<],siz[k]=siz[k<<]+siz[k<<|];
else if(mn[k<<]<mn[k<<|]) mn[k]=mn[k<<],siz[k]=siz[k<<];
else mn[k]=mn[k<<|],siz[k]=siz[k<<|];
}
inline void pushdown(ll k){
if(lz[k]){
lz[k<<]+=lz[k],lz[k<<|]+=lz[k];
mn[k<<]+=lz[k],mn[k<<|]+=lz[k];
lz[k]=;
}
if(nd[k]){
if(mn[k<<]==mn[k]) sum[k<<]+=siz[k<<]*nd[k],nd[k<<]+=nd[k];
if(mn[k<<|]==mn[k]) sum[k<<|]+=siz[k<<|]*nd[k],nd[k<<|]+=nd[k];
nd[k]=;
}
}
inline void build(ll l,ll r,ll k){
if(l==r){mn[k]=inf,siz[k]=;return;}
ll mid=l+r>>;
build(l,mid,k<<),build(mid+,r,k<<|);
pushup(k); return;
}
inline void add(ll l,ll r,ll x,ll y,ll v,ll t,ll k){
if(x<=l && r<=y){mn[k]+=v,lz[k]+=v;return;}
pushdown(k); ll mid=l+r>>;
if(x<=mid) add(l,mid,x,y,v,t,k<<);
if(y>mid) add(mid+,r,x,y,v,t,k<<|);
pushup(k); return;
}
inline void ins(ll l,ll r,ll p,ll v,ll t,ll k){
if(l==r){mn[k]=v;return;}
pushdown(k); ll mid=l+r>>;
if(p<=mid) ins(l,mid,p,v,t,k<<);
else ins(mid+,r,p,v,t,k<<|);
pushup(k); return;
}
inline ll query(ll l,ll r,ll x,ll y,ll t,ll k){
if(x<=l && r<=y) return sum[k];
pushdown(k); ll mid=l+r>>,res=;
if(x<=mid) res+=query(l,mid,x,y,t,k<<);
if(y>mid) res+=query(mid+,r,x,y,t,k<<|);
return res;
} int main(){
N=read();
for(ll i=;i<=N;i++) A[i]=read();
Q=read();
for(ll i=;i<=Q;i++){
ll l=read(),r=read();
ques[r].push_back(make_pair(i,l));
}
build(,N,);
for(ll i=;i<=N;i++){
ll x=A[i],nl1=i,nl2=i;
while(!smx.empty()){
if(smx.top().v<x){
ll l=smx.top().l,r=smx.top().r;
add(,N,l,r,x-smx.top().v,i,);
nl1=l,smx.pop();
}
else break;
}
while(!smn.empty()){
if(smn.top().v>x){
ll l=smn.top().l,r=smn.top().r;
add(,N,l,r,smn.top().v-x,i,);
nl2=l,smn.pop();
}
else break;
}
node tp1; tp1.v=x,tp1.l=nl1,tp1.r=i;
node tp2; tp2.v=x,tp2.l=nl2,tp2.r=i;
smx.push(tp1),smn.push(tp2);
if(i-) add(,N,,i-,-,i,);
ins(,N,i,,i,);
sum[]+=siz[],nd[]+=;
for(ll j=;j<ques[i].size();j++){
ll l=ques[i][j].second,r=i,id=ques[i][j].first;
ans[id]=query(,N,l,r,i,);
}
}
for(ll i=;i<=Q;i++)
printf("%I64d\n",ans[i]);
return ;
}
cf 997 E
[cf 997 E] Good Subsegments的更多相关文章
- #C++初学记录(set进阶#acm cf 190802 B. Subsegments)
B. Subsegments#set进阶 Programmer Sasha has recently begun to study data structures. His coach Stas to ...
- CF 372B Counting Rectangles is Fun [dp+数据维护]
题意,给出一个n行m列的矩阵 里面元素是0或者1 给出q个询问 a,b,c,d 求(a,b)到(c,d)有多少个由0组成的矩形 我们定义 watermark/2/text/aHR0cDovL2Jsb2 ...
- ORA-00494: enqueue [CF] held for too long (more than 900 seconds) by 'inst 1, osid 5166'
凌晨收到同事电话,反馈应用程序访问Oracle数据库时报错,当时现场现象确认: 1. 应用程序访问不了数据库,使用SQL Developer测试发现访问不了数据库.报ORA-12570 TNS:pac ...
- cf之路,1,Codeforces Round #345 (Div. 2)
cf之路,1,Codeforces Round #345 (Div. 2) ps:昨天第一次参加cf比赛,比赛之前为了熟悉下cf比赛题目的难度.所以做了round#345连试试水的深浅..... ...
- cf Round 613
A.Peter and Snow Blower(计算几何) 给定一个点和一个多边形,求出这个多边形绕这个点旋转一圈后形成的面积.保证这个点不在多边形内. 画个图能明白 这个图形是一个圆环,那么就是这个 ...
- ARC下OC对象和CF对象之间的桥接(bridge)
在开发iOS应用程序时我们有时会用到Core Foundation对象简称CF,例如Core Graphics.Core Text,并且我们可能需要将CF对象和OC对象进行互相转化,我们知道,ARC环 ...
- [Recommendation System] 推荐系统之协同过滤(CF)算法详解和实现
1 集体智慧和协同过滤 1.1 什么是集体智慧(社会计算)? 集体智慧 (Collective Intelligence) 并不是 Web2.0 时代特有的,只是在 Web2.0 时代,大家在 Web ...
- CF memsql Start[c]UP 2.0 A
CF memsql Start[c]UP 2.0 A A. Golden System time limit per test 1 second memory limit per test 256 m ...
- CF memsql Start[c]UP 2.0 B
CF memsql Start[c]UP 2.0 B B. Distributed Join time limit per test 1 second memory limit per test 25 ...
随机推荐
- yarn和npm的对比,以及项目中使用方式
相比npm 的优点 1.npm安装是串行,而yarn是并行,速度大大提升 2.已经下载过的包会被缓存,无需重复下载,更关键的是,支持离线安装 3.精准的版本控制,加上验证每个包的完整性,保证每次安装的 ...
- 【转】Impala 中的 Invalidate Metadata 和 Refresh
前言Impala采用了比较奇葩的多个impalad同时提供服务的方式,并且它会由catalogd缓存全部元数据,再通过statestored完成每一次的元数据的更新到impalad节点上,Impala ...
- linux高性能服务器编程 (四) --TCP/IP通信案例
第四章 TCP/IP通信案例 HTTP代理服务器的大致工作原理 在HTTP通信链上,客户端和服务器之间通常存在某些中转代理服务器.它们提供对目标资源的中转访问.一个HTTP请求可能被多个 ...
- python 路径拼接
>>> import os>>> os.path.join('/hello/','good/boy/','doiido')>>>'/hello/g ...
- java项目里classpath具体指哪儿个路径
一.classpath路径指什么 只知道把配置文件如:mybatis.xml.spring-web.xml.applicationContext.xml等放到src目录(就是存放代码.java文件的目 ...
- 【CSP膜你赛】柠檬的密码(manacher 二分 单调性 st表)
题目描述 Lemon觉得他需要一个复杂的密码来保证他的帐号的安全.他经过多日思考,决定使用一个长度为奇数的回文串来作为他的密码. 但是这个回文串太长了,Lemon记不住,于是Lemon决定把它记在本 ...
- answer
https://www.cnblogs.com/549294286/p/10451394.html 基于BIO实现的Server端,当建立了100个连接时,会有多少个线程?如果基于NIO,又会是多少个 ...
- Parallel.For循环与普通的for循环
前两天看书发现了一个新的循环Parallel.For,这个循环在循环期间可以创建多个线程并行循环,就是说循环的内容是无序的.这让我想到了我前面的牛牛模拟计算是可以用到这个循环的,我前面的牛牛模拟计算是 ...
- vs2017+qt5.x编译32位应用<转>
原文地址:https://www.cnblogs.com/woniu201/p/10862170.html 概述 最近有同学私信我,问如何使用vs2017+qt5.10编译出32位的应用,需要使用ms ...
- Python + Selenium 自动发布文章(一):开源中国
https://blog.csdn.net/qq_28804275/article/details/80891949 https://blog.csdn.net/qq_28804275/article ...