[NOIP模拟赛] seq
seq
试题分析
介绍一种方法叫做回滚莫队。
- 回滚莫队是一种只加不删的莫队。
首先处理\(l,r\)都在同一个块内的询问,暴力即可。
然后对于\(l,r\)不同在一个块,我们将左端点挂在其所在的块。
将挂在每个块上的右端点排序,这样对于左边我们每次暴力滚回当前块的右端点,这里可以用主席树之类的数据结构维护。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<cmath>
#include<algorithm>
using namespace std;
#define LL long long
inline int read(){
int x=0,f=1; char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
for(;isdigit(c);c=getchar()) x=x*10+c-'0';
return x*f;
}
const int INF = 2147483600;
const int MAXN = 100000;
int N,M; int a[MAXN+1];
int bel[MAXN+1];
struct data{int l,r,id;}q[MAXN+1];
bool cmp(data a,data b){
if(bel[a.l]!=bel[b.l]) return bel[a.l]<bel[b.l];
return a.r<b.r;
}
int L[MAXN+1],R[MAXN+1]; int res;
inline void init_over(){
res=0; memset(L,0,sizeof(L));
memset(R,0,sizeof(R)); return ;
}
struct stak{int pos,x;}staR[MAXN+1],staL[MAXN+1];
int ans[MAXN+1];
int main(){
//freopen(".in","r",stdin);
//freopen(".out","w",stdout);
N=read(),M=read(); int sqr=(int)sqrt(N)+100,top;
for(int i=1;i<=N;i++) a[i]=read();
for(int i=1;i<=M;i++) q[i].l=read(),q[i].r=read(),q[i].id=i;
for(int i=1;i<=N;i++) bel[i]=(i-1)/sqr+1;
sort(q+1,q+M+1,cmp); int l=bel[q[1].l]*sqr,r=bel[q[1].l]*sqr;
for(int i=1;i<=M;i++){
//cout<<q[i].l<<" "<<q[i].r<<":"<<bel[q[i].l]<<" "<<bel[q[i].r]<<endl;
l=bel[q[i].l]*sqr; //l=min(l,N);
if(bel[q[i].l]!=bel[q[i-1].l]) r=bel[q[i].l]*sqr,init_over();
/*if(bel[q[i].l]==bel[q[i].r]){
res=0;
for(int j=q[i].l;j<=q[i].r;j++){
res=max(res,L[a[j]-1]+R[a[j]+1]+1);
R[a[j]-L[a[j]-1]]=L[a[j]+R[a[j]+1]]=L[a[j]-1]+R[a[j]+1]+1;
}
ans[q[i].id]=res;
for(int j=q[i].l;j<=q[i].r;j++) L[a[j]]=R[a[j]]=0;
continue;
}*/
while(r<q[i].r){
++r; res=max(res,L[a[r]-1]+R[a[r]+1]+1);
R[a[r]-L[a[r]-1]]=L[a[r]+R[a[r]+1]]=L[a[r]-1]+R[a[r]+1]+1;
} top=0;
int res2=res;
for(l=min(l,q[i].r);l>=q[i].l;--l){
res2=max(res2,L[a[l]-1]+R[a[l]+1]+1);
int ll=a[l]-L[a[l]-1],rr=a[l]+R[a[l]+1];
staR[++top].pos=ll; staR[top].x=R[ll];
staL[top].pos=rr; staL[top].x=L[rr];
R[ll]=L[rr]=L[a[l]-1]+R[a[l]+1]+1;
} ans[q[i].id]=res2;
for(;top;--top){
R[staR[top].pos]=staR[top].x;
L[staL[top].pos]=staL[top].x;
}
}
for(int i=1;i<=M;i++) printf("%d\n",ans[i]);
return 0;
}
[NOIP模拟赛] seq的更多相关文章
- NOIP模拟赛 6.29
2017-6-29 NOIP模拟赛 Problem 1 机器人(robot.cpp/c/pas) [题目描述] 早苗入手了最新的Gundam模型.最新款自然有着与以往不同的功能,那就是它能够自动行走, ...
- NOIP模拟赛20161022
NOIP模拟赛2016-10-22 题目名 东风谷早苗 西行寺幽幽子 琪露诺 上白泽慧音 源文件 robot.cpp/c/pas spring.cpp/c/pas iceroad.cpp/c/pas ...
- contesthunter暑假NOIP模拟赛第一场题解
contesthunter暑假NOIP模拟赛#1题解: 第一题:杯具大派送 水题.枚举A,B的公约数即可. #include <algorithm> #include <cmath& ...
- NOIP模拟赛 by hzwer
2015年10月04日NOIP模拟赛 by hzwer (这是小奇=> 小奇挖矿2(mining) [题目背景] 小奇飞船的钻头开启了无限耐久+精准采集模式!这次它要将原矿运到泛光之源的矿 ...
- 大家AK杯 灰天飞雁NOIP模拟赛题解/数据/标程
数据 http://files.cnblogs.com/htfy/data.zip 简要题解 桌球碰撞 纯模拟,注意一开始就在袋口和v=0的情况.v和坐标可以是小数.为保险起见最好用extended/ ...
- 队爷的讲学计划 CH Round #59 - OrzCC杯NOIP模拟赛day1
题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的讲学计划 题解:刚开始理解题意理解了好半天,然后发 ...
- 队爷的Au Plan CH Round #59 - OrzCC杯NOIP模拟赛day1
题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的Au%20Plan 题解:看了题之后觉得肯定是DP ...
- 队爷的新书 CH Round #59 - OrzCC杯NOIP模拟赛day1
题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的新书 题解:看到这题就想到了 poetize 的封 ...
- CH Round #58 - OrzCC杯noip模拟赛day2
A:颜色问题 题目:http://ch.ezoj.tk/contest/CH%20Round%20%2358%20-%20OrzCC杯noip模拟赛day2/颜色问题 题解:算一下每个仆人到它的目的地 ...
随机推荐
- HDU 2546 饭卡 (dp)
题目链接 Problem Description 电子科大本部食堂的饭卡有一种很诡异的设计,即在购买之前判断余额.如果购买一个商品之前,卡上的剩余金额大于或等于5元,就一定可以购买成功(即使购买后卡上 ...
- error 0152: No Entity Framework provider found for the ADO.NET provider with invariant name 'System.Data.SqlClient'
error 0152: No Entity Framework provider found for the ADO.NET provider with invariant name 'System. ...
- 更改控件中DrawableLeft图片的大小,图片与文字的距离
Drawable drawable=getResources().getDrawable(R.drawable.xx); //获取图片 drawable.setBounds(left, top, ri ...
- 网络应用框架Netty快速入门
一 初遇Netty Netty是什么? Netty 是一个提供 asynchronous event-driven (异步事件驱动)的网络应用框架,是一个用以快速开发高性能.可扩展协议的服务器和客户端 ...
- 通过实例来学习XML DTD
使用DTD的原因: 注意:由于它自身的一些缺点,DTD终将被淘汰,但是它还是要学习的.学习完DTD后,后面继续学习XML Schema. 1,通过 DTD,您的每一个 XML 文件均可携带一个有关其自 ...
- 最简单的基于FFMPEG的图像编码器(YUV编码为JPEG)(转)
原文转自 https://blog.csdn.net/leixiaohua1020/article/details/25346147/ 伴随着毕业论文的完成,这两天终于腾出了空闲,又有时间搞搞FFMP ...
- LeetCode 19 Valid Parentheses
Given a string containing just the characters '(', ')', '{', '}', '[' and ']', determine if the inpu ...
- win7旗舰版64位缺失tbb.dll文件
win7旗舰版64位缺失tbb.dll文件 https://zhidao.baidu.com/question/688589990330312804.html 到好的电脑中复制一个,黏贴到下同的路径下 ...
- delphi2006语言新特性:Record类型高级用法
delphi语言在传统的Records类型的基础上增加了许多像类一样的高级功能,如:Records可以有属性和方法(包括构造constructors),类属性,类方法,类静态字段和内嵌类型.下面这个示 ...
- Geoserver发布缓存切片(制定Gridsets)
EPSG:4326 Level Pixel Size Scale Name Tiles 0 1: 2 x 1 1 1: 4 x 2 2 1: 8 x 4 3 1: 16 x 8 4 ...