Description
2333年的某一天，临冬突降大雪，主干道已经被雪覆盖不能使用。城
主 囧·雪 决定要对主干道进行一次清扫。
临冬城的主干道可以看为一条数轴。囧·雪 一共找来了n个清理工，第
i个清理工的工作范围为[li,ri]，也就是说这个清理工会把[li,ri]这一
段主干道清理干净(当然已经被清理过的部分就被忽略了)。当然有可能主
干道不能全部被清理干净，不过这也没啥关系。
虽然 囧·雪 啥都不知道，但是他还是保证了不会出现某一个清理工的
工作范围被另一个清理工完全包含的情况(不然就太蠢了)。
作为临冬城主，囧·雪 给出了如下的清扫方案：
在每一天开始的时候，每一个还没有工作过的清理工会观察自己工作
范围内的道路，并且记下工作范围内此时还没有被清理的道路的长度(称
为这个清理工的工作长度)。然后 囧·雪 会从中选择一个工作长度最小的
清理工(如果两个清理工工作长度相同，那么就选择编号小的清理工)。然
后被选择的这个清理工会清理自己的工作范围内的道路。为了方便检查工
作质量，囧·雪 希望每一天只有一个清理工在工作。
你要注意，清理工的工作长度是可能改变的，甚至有可能变成0。尽管
如此，这个清理工也还是会在某一天工作。
现在，囧·雪 想要知道每一天都是哪个清理工在工作？
Input
第一行两个整数t,n。分别表示主干道的长度(也就是说，主干道是数
轴上[1,t]的这一段)以及清理工的人数。
接下来n行，每行两个整数li,ri。意义如题。
n<=3*10^5, 1<=li<ri<=t<=10^9，保证输入的li严格递增
Output
输出n行，第i行表示第i天工作的清理工的编号。
区间互不包含说明左端点升序时右端点严格单调，区间端点离散化，用线段树维护按左端点排序后工作长度的最小值，每次查出最小值后将其删去并枚举这次清扫的区间，对应在线段树中是连续的一段，可以标记实现区间减。
#include<cstdio>
#include<algorithm>
const int N=300077;
char buf[10000007],*ptr=buf-1;
int _(){
int x=0,c=*++ptr;
while(c<48)c=*++ptr;
while(c>47)x=x*10+c-48,c=*++ptr;
return x;
}
int min(int a,int b){return a<b?a:b;}
int mx,n,ls[N],rs[N],xs[N*2],xp=0,ans,f[N*2],_l,_r,_a;
struct node{
node*lc,*rc;
int L,R,mn,a;
void up(){mn=min(lc->mn,rc->mn);}
void dn(){
if(a){
lc->mn+=a;
lc->a+=a;
rc->mn+=a;
rc->a+=a;
a=0;
}
}
void dec(){
if(_l<=L&&R<=_r){
mn+=_a;
a+=_a;
return;
}
int M=L+R>>1;
dn();
if(_l<=M)lc->dec();
if(_r>M)rc->dec();
up();
}
void cal(){
if(L==R){
ans=L;
mn=0x7fffffff;
return;
}
dn();
(lc->mn<=rc->mn?lc:rc)->cal();
up();
}
}ns[N*2],*np=ns,*rt;
node*build(int L,int R){
node*w=np++;
w->L=L;w->R=R;
if(L==R){
w->mn=rs[L]-ls[L];
}else{
int M=L+R>>1;
w->lc=build(L,M);
w->rc=build(M+1,R);
w->up();
}
return w;
}
int get(int x){
int a=x,c;
while(x!=f[x])x=f[x];
while(x!=f[a])c=f[a],f[a]=x,a=c;
return x;
}
int main(){
fread(buf,1,sizeof(buf),stdin);
mx=_();n=_();
for(int i=1;i<=n;++i){
ls[i]=_();
rs[i]=_();
}
rt=build(1,n);
int p1=1,p2=1;
while(p1<=n&&p2<=n){
if(ls[p1]<rs[p2])xs[++xp]=ls[p1],ls[p1++]=xp;
else if(ls[p1]>rs[p2])xs[++xp]=rs[p2],rs[p2++]=xp;
else xs[++xp]=ls[p1],ls[p1++]=rs[p2++]=xp;
}
while(p1<=n)xs[++xp]=ls[p1],ls[p1++]=xp;
while(p2<=n)xs[++xp]=rs[p2],rs[p2++]=xp;
for(int i=1;i<=xp;++i)f[i]=i;
for(int i=1;i<=n;++i){
rt->cal();
printf("%d\n",ans);
for(int p=get(ls[ans]);p<rs[ans];p=f[p]=get(p+1)){
_l=std::upper_bound(rs+1,rs+n+1,p)-rs;
_r=std::upper_bound(ls+1,ls+n+1,p)-ls-1;
_a=xs[p]-xs[p+1];
rt->dec();
}
}
return 0;
}