【BZOJ】3339: Rmq Problem & 3585: mex(线段树+特殊的技巧)

http://www.lydsy.com/JudgeOnline/problem.php?id=3585

好神的题。

但是!!!!!!!!!!!!!!我线段树现在要开8倍空间才能过!!!!!!!!!!这什么梗。。。。。。。。。。。。。。。。。。。。。。

我思考了很久空间的问题,因为我在pushdown的时候可能会越界,或许是pushup?

不管了。然后看了zyf的写法。看来以后得注意下。。。pushdown弄成先放了。。。

本题的做法:

好神orz

首先像莫队一样离线区间,左端点在前。

考虑如何从[l, r]转移到[l+1, r]:

显然,a[l]此时是去掉了,l+1~next[l]是空着了a[l]这个自然数的。next[l]是下一个a[l]的位置。

所以我们直接线段树更新l+1~next[i]的sg值即可!即如果这个区间内有询问的sg值>a[l],那么更新!

好神。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define rdm(x, i) for(int i=ihead[x]; i; i=e[i].next)
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; }

#define lc x<<1
#define rc x<<1|1
#define MID (l+r)>>1
#define lson l, mid, lc
#define rson mid+1, r, rc

const int N=200005, oo=~0u>>1;
struct dat { int l, r, id; }q[N];
int sg[N], mn[N<<2];
void upd(int x, int k) { mn[x]=min(mn[x], k); }
void pushdown(int x) {
	if(mn[x]!=oo) {
		upd(lc, mn[x]);
		upd(rc, mn[x]);
		mn[x]=oo;
	}
}
void build(int l, int r, int x) {
	if(l==r) { mn[x]=sg[l]; return; }
	mn[x]=oo;
	int mid=MID;
	build(lson); build(rson);
}
void update(int l, int r, int x, int L, int R, int k) {
	if(L<=l && r<=R) { mn[x]=min(mn[x], k); return; }
	pushdown(x);
	int mid=MID;
	if(L<=mid) update(lson, L, R, k);
	if(mid<R) update(rson, L, R, k);
}
int query(int l, int r, int x, int k) {
	if(l==r) return mn[x];
	pushdown(x);
	int mid=MID;
	if(k<=mid) return query(lson, k);
	return query(rson, k);
}
int a[N], b[N], n, m, tot, vis[N], ans[N], inext[N];

void getnext() {
	for1(i, 1, n) b[i]=a[i];
	sort(b+1, b+1+n); tot=unique(b+1, b+1+n)-b-1;
	for1(i, 1, n) a[i]=lower_bound(b+1, b+1+tot, a[i])-b;
	for1(i, 1, tot) vis[i]=0;
	for3(i, n, 1) inext[i]=vis[a[i]], vis[a[i]]=i;
}
bool cmp(const dat &a, const dat &b) { return a.l<b.l; }
void init() {
	int k=0;
	for1(i, 1, n) {
		if(a[i]<N) vis[a[i]]=1;
		while(vis[k]) ++k;
		sg[i]=k;
	}
	build(1, n, 1);
	sort(q+1, q+1+m, cmp);
	getnext();
	int now=1;
	for1(i, 1, m) {
		int l=q[i].l;
		while(now<l) {
			int nxt=inext[now];
			if(nxt==0) nxt=n;
			else --nxt; // printf("now:%d nxt:%d a[now]:%d\n", now, nxt, b[a[now]]);
			if(now+1<=nxt) update(1, n, 1, now+1, nxt, b[a[now]]);
			++now;
		}
		ans[q[i].id]=query(1, n, 1, q[i].r);
	}
	for1(i, 1, m) printf("%d\n", ans[i]);
}

int main() {
	read(n); read(m);
	for1(i, 1, n) read(a[i]);
	for1(i, 1, m) read(q[i].l), read(q[i].r), q[i].id=i;
	init();
	return 0;
}

  

 


 

 

Description

  有一个长度为n的数组{a1,a2,...,an}。m次询问,每次询问一个区间内最小没有出现过的自然数。

Input

  第一行n,m。
  第二行为n个数。
  从第三行开始,每行一个询问l,r。

Output

  一行一个数,表示每个询问的答案。

Sample Input

5 5
2 1 0 2 1
3 3
2 3
2 4
1 2
3 5

Sample Output

1
2
3
0
3

HINT

 

数据规模和约定

  对于100%的数据:

  1<=n,m<=200000

  0<=ai<=109

  1<=l<=r<=n


  对于30%的数据:


  1<=n,m<=1000

 

Source

转载于:https://www.cnblogs.com/iwtwiioi/p/4149394.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值