调了一早上80pt不知道那里错了。。。结果发现我看错题意了。。。一个清洁工扫的是区间而不是点。。。
先离散化
首先对于一个人它扫了地，他影响的其他人是一个区间
而且每个区间都只会被扫一次，所以我们可以暴力枚举
那么算法的架构已经出来了：线段树维护清洁工清理区间长度的最小值，然后每次选出一个清洁工，就把它扫的区间枚举一次，修改会被影响的其他清洁工
并查集瞎维护一下就可以了。。。今年NOIP才考
#include<cstdio> #include<iostream> #include<cstring> #include<cstdlib> #include<algorithm> #include<cmath> #define mid (ql+qr)/2 #define lc now<<1 #define rc now<<1|1 using namespace std; typedef pair<int,int> pa; const int _=100; const int maxn=3*1e5+_; const int fbin=(1<<19)+_; int read() { int x=0,f=1;char ch=getchar(); while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();} while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();} return x*f; } void write(int x) { if(x<0)putchar('-'); if(x>=10)write(x/10); putchar(x%10+'0'); } struct node{int l,r,id;}p[maxn]; bool cmp(node n1,node n2){return n1.l<n2.l;} int lslen,ls[maxn*4],sd[maxn*4]; //----------------------------------def-------------------------------------------------- struct trnode { int mn,id,la; }tr[fbin*2];int z; void update(int now) { if(tr[lc].mn<0){tr[now].mn=tr[rc].mn,tr[now].id=tr[rc].id;return ;} if(tr[rc].mn<0){tr[now].mn=tr[lc].mn,tr[now].id=tr[lc].id;return ;} if((tr[lc].mn<tr[rc].mn||(tr[lc].mn==tr[rc].mn&&tr[lc].id<tr[rc].id))) tr[now].mn=tr[lc].mn,tr[now].id=tr[lc].id; else tr[now].mn=tr[rc].mn,tr[now].id=tr[rc].id; } void pushdown(int now) { if(tr[now].la==0)return ; if(tr[lc].mn!=-1)tr[lc].mn+=tr[now].la,tr[lc].la+=tr[now].la; if(tr[rc].mn!=-1)tr[rc].mn+=tr[now].la,tr[rc].la+=tr[now].la; tr[now].la=0; } //~~~~~~~~in~~~~~~~~~~~ void bt(int now,int ql,int qr) { if(ql==qr){ tr[now].mn=sd[p[ql].r]-sd[p[ql].l],tr[now].id=++z; return ; } else{ bt(lc,ql,mid),bt(rc,mid+1,qr); update(now); } } void change(int now,int ql,int qr,int l,int r,int d) { if(ql==l&&qr==r){tr[now].mn+=d,tr[now].la+=d;return ;} pushdown(now); if(r<=mid) change(lc,ql,mid,l,r,d); else if(mid+1<=l)change(rc,mid+1,qr,l,r,d); else change(lc,ql,mid,l,mid,d),change(rc,mid+1,qr,mid+1,r,d); update(now); } int findmn(int now,int ql,int qr,int p) { if(ql==qr){return tr[now].mn;} pushdown(now); if(p<=mid)return findmn(lc,ql,mid,p); else return findmn(rc,mid+1,qr,p); } //~~~~~~~~~~~out~~~~~~~~~~~~~~~ //-----------------------------------seg---------------------------------------------- int R[maxn*4]; int findR(int x) { if(R[x]==x)return x; R[x]=findR(R[x]);return R[x]; } int main() { freopen("a.in","r",stdin); freopen("a.out","w",stdout); int T,n; T=read(),n=read(); for(int i=1;i<=n;i++) p[i].l=read(),p[i].r=read(),p[i].id=i; sort(p+1,p+n+1,cmp); for(int i=1;i<=n;i++) ls[++lslen]=p[i].l,ls[++lslen]=p[i].r; sort(ls+1,ls+lslen+1); lslen=unique(ls+1,ls+lslen+1)-ls-1; for(int i=2;i<=lslen;i++) sd[i]=sd[i-1]+ls[i]-ls[i-1]; for(int i=1;i<=n;i++) { p[i].l=lower_bound(ls+1,ls+lslen+1,p[i].l)-ls, p[i].r=lower_bound(ls+1,ls+lslen+1,p[i].r)-ls; } //-------------------------LSH------------------------- bt(1,1,n); for(int i=0;i<=lslen+1;i++)R[i]=i; int x; for(int t=1;t<=n;t++) { x=tr[1].id; write(p[x].id);puts(""); int l,r,st,ed; for(int k=findR(p[x].l+1);k<=p[x].r;R[k]=findR(k+1),k=R[k]) { l=1,r=n,st=n+1; while(l<=r) { int m=(l+r)/2; if(k<=p[m].r)r=m-1,st=m; else l=m+1; } l=1,r=n,ed=0; while(l<=r) { int m=(l+r)/2; if(p[m].l<k)l=m+1,ed=m; else r=m-1; } if(st<=ed) change(1,1,n,st,ed,sd[k-1]-sd[k]); } change(1,1,n,x,x,-1); } return 0; }