【LOJ #3273】「JOISC 2020 Day1」扫除（线段树分治 / 线段树）

传送门

考虑对于所有被扫除过的点
有若$x_i\le x_j,则y_i\ge y_j$
即所有灰尘位置是单调的

先考虑没有插入的情况
对于一次修改，有可能其中一部分地方被之前的一次操作修改过
而因此现在这一块没有灰尘，这一块就是这次修改影响不到的地方
、
考虑用权值线段树维护每个位置可以修改得到的地方
比如如果之前有一次$l=l_0$的$V$
那么对于这一次$l=l_1$的$H$
实际上只会将$x\in [l_0+1,n-l_1]$内的灰尘的$x\to n-l_1$
相当于区间赋值

对于$x,y$分别维护一个权值线段树支持区间取$max$，单点查询即可得到每次影响的范围
再对所有灰尘求一次答案即可

考虑如果有中途加入的灰尘
用线段树分治转成没有插入的情况即可

复杂度$O(nlo{g}^{2}n)$

#include<bits/stdc++.h>
using namespace std;
#define cs const
#define re register
#define pb push_back
#define pii pair<int,int>
#define ll long long
#define fi first
#define se second
#define bg begin
cs int RLEN=1<<20|1;
inline char gc(){
    static char ibuf[RLEN],*ib,*ob;
    (ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
    return (ib==ob)?EOF:*ib++;
}
inline int read(){
    char ch=gc();
    int res=0;bool f=1;
    while(!isdigit(ch))f^=ch=='-',ch=gc();
    while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
    return f?res:-res;
}
inline ll readll(){
    char ch=gc();
    ll res=0;bool f=1;
    while(!isdigit(ch))f^=ch=='-',ch=gc();
    while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
    return f?res:-res;
}
inline int readstring(char *s){
	int top=0;char ch=gc();
	while(isspace(ch))ch=gc();
	while(!isspace(ch)&&ch!=EOF)s[++top]=ch,ch=gc();
	return top;
}
template<typename tp>inline void chemx(tp &a,tp b){a<b?a=b:0;}
template<typename tp>inline void chemn(tp &a,tp b){a>b?a=b:0;}
cs int N=1500005;
char xxx;
struct dust{int x,y,t;}a[N];
struct node{int x,y;}ans[N];
struct opt{int op,len;}b[N];
struct qry{int id,l,r;}c[N];
vector<pii> upd[2];
int n,m,q,tn,qn,rt0,rt1;
namespace Seg{
	int tot;
	cs int N=::N*40;
	int lc[N],rc[N],mx[N];
	#define mid ((l+r)>>1)
	inline void clear(){
		rt0=rt1=tot=0;
	}
	inline int newnode(){
		int u=++tot;lc[u]=rc[u]=mx[u]=0;
		return u;
	}
	void update(int &u,int l,int r,int st,int des,int k){
		if(l>r||r<st||l>des)return;
		if(!u)u=newnode();
		if(st<=l&&r<=des){chemx(mx[u],k);return;}
		if(st<=mid)update(lc[u],l,mid,st,des,k);
		if(mid<des)update(rc[u],mid+1,r,st,des,k);
	}
	int query(int u,int l,int r,int p){
		if(!u||l>r)return 0;
		if(l==r)return mx[u];
		if(p<=mid)return max(mx[u],query(lc[u],l,mid,p));
		else return max(mx[u],query(rc[u],mid+1,r,p));
	}
	#undef mid
}
char yyy;
void update(int op,int l){
	int p=Seg::query(op?rt0:rt1,0,n,l);
	upd[op].pb(pii(l,p));//cout<<op<<" "<<p<<" "<<l<<'\n';
	Seg::update(op?rt1:rt0,0,n,p,n-l-1,l+1);
}
void solve(int l,int r,vector<int> qi){
	if(!qi.size())return;
	vector<int> tp,lq,rq;
	int mid=(l+r)>>1;
	for(int &i:qi){
		int st=c[i].l,des=c[i].r;
		if(st<=l&&r<=des)tp.pb(i);
		else{if(st<=mid)lq.pb(i);
			 if(mid<des)rq.pb(i);}
	}
	for(int i=l;i<=r;i++)update(b[i].op,b[i].len);
	sort(tp.bg(),tp.end(),[](cs int &i,cs int &j){return ans[i].y>ans[j].y;});
	sort(upd[0].bg(),upd[0].end(),greater<pii>()),Seg::clear();
	int j=0;
	for(int &i:tp){
		int x=ans[i].x,y=ans[i].y;
		for(;j<upd[0].size()&&upd[0][j].fi>=y;j++)
		Seg::update(rt0,0,n,upd[0][j].se,n-upd[0][j].fi,n-upd[0][j].fi);
		chemx(ans[i].x,Seg::query(rt0,0,n,x));
	}
	sort(tp.bg(),tp.end(),[](cs int &i,cs int &j){return ans[i].x>ans[j].x;});
	sort(upd[1].bg(),upd[1].end(),greater<pii>()),j=0;
	for(int &i:tp){
		int x=ans[i].x,y=ans[i].y;
		for(;j<upd[1].size()&&upd[1][j].fi>=x;j++)
		Seg::update(rt1,0,n,upd[1][j].se,n-upd[1][j].fi,n-upd[1][j].fi);
		chemx(ans[i].y,Seg::query(rt1,0,n,y));
	}
//	cout<<l<<" "<<r<<" "<<ans[1].x<<" "<<ans[1].y<<'\n';
	Seg::clear(),upd[0].clear(),upd[1].clear();
	if(l==r)return;
	solve(l,mid,lq),solve(mid+1,r,rq);
}
int main(){
	#ifdef Stargazer
	freopen("lx.in","r",stdin);
	cerr<<(&yyy-&xxx)/1024/1024<<'\n';
	#endif
	n=read(),m=read(),q=read();
	for(int i=1;i<=m;i++)
		a[i].x=read(),a[i].y=read(),a[i].t=1;
	for(int i=1;i<=q;i++){
		int op=read(),x;
		switch(op){
			case 1:c[++qn].id=x=read(),c[qn].l=a[x].t,c[qn].r=tn;break;
			case 2:b[++tn].op=0,b[tn].len=read();break;
			case 3:b[++tn].op=1,b[tn].len=read();break;
			case 4:a[++m].x=read(),a[m].y=read(),a[m].t=tn+1;
		}
	}vector<int> vec;
	for(int i=1,x;i<=qn;i++){
		x=c[i].id,ans[i].x=a[x].x,ans[i].y=a[x].y;
		if(c[i].l<=c[i].r)vec.pb(i);
	}
	solve(1,tn,vec);
	for(int i=1;i<=qn;i++)cout<<ans[i].x<<" "<<ans[i].y<<'\n';
	return 0;
}