首页 > 代码库 > 【BZOJ】3339: Rmq Problem & 3585: mex(线段树+特殊的技巧)
【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, rcconst 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
2 1 0 2 1
3 3
2 3
2 4
1 2
3 5
Sample Output
1
2
3
0
3
2
3
0
3
HINT
数据规模和约定
对于100%的数据:
1<=n,m<=200000
0<=ai<=109
1<=l<=r<=n
对于30%的数据:
1<=n,m<=1000
Source
By 佚名提供
【BZOJ】3339: Rmq Problem & 3585: mex(线段树+特殊的技巧)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。