大学生程序设计邀请赛（华东师范大学）-D-线段树

链接：点击打开链接

题解：二分预处理出来每个点单次传染到的左右范围，再用线段树维护每个点最终的范围最小值、最大值，利用线段树不断扩大范围并更新。随机一个
1
~
n
的序列利用线段树更新，这样速度会变得很快，或其他优化方法均可（这里采用sin函数排序差不多就无序了= =）。还有就是注意题目是按原来的i输出，我因为这个被卡了一下午，MDZZ。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define inf 0x3f3f3f3f
const int mx=1e5+10;
typedef long long ll;
int n,sum[mx<<2][2],a[mx],ans[mx],b[mx];
using namespace std;
struct sec{
int l,r;
bool operator < (sec A)const{
return l<A.l;
}
sec() {}
sec(int ll,int rr){ l=ll,r=rr; }
}s[mx],po[mx];
bool cmp(int a,int b){
return sin(1.0*a)<sin(1.0*b);
}
void build(int l,int r,int rt){
if(l==r){
sum[rt][0]=s[l].l;
sum[rt][1]=s[l].r;
return ;
}
int m=(l+r)>>1;
build(lson);
build(rson);
sum[rt][0]=min(sum[rt<<1][0],sum[rt<<1|1][0]);
sum[rt][1]=max(sum[rt<<1][1],sum[rt<<1|1][1]);
}
sec query(int L,int R,int l,int r,int rt){
if(L<=l&&R>=r){
return sec(sum[rt][0],sum[rt][1]);
}
int m=(l+r)>>1;
sec num[2]={{inf,0},{inf,0}};
if(m>=L) num[0]=query(L,R,lson);
if(R>m) num[1]=query(L,R,rson);
return sec( min(num[0].l,num[1].l) , max(num[0].r,num[1].r) );
}
void update(int L,int R,int l,int r,int rt,int vl,int vr){
if(l==r){
sum[rt][0]=min(vl,sum[rt][0]);
sum[rt][1]=max(vr,sum[rt][1]);
return ;
}
int m=(l+r)>>1;
if(m>=L) update(L,R,lson,vl,vr);
if(R>m) update(L,R,rson,vl,vr);
sum[rt][0]=min(sum[rt<<1][0],sum[rt<<1|1][0]);
sum[rt][1]=max(sum[rt<<1][1],sum[rt<<1|1][1]);
}
int main(){
while(~scanf("%d",&n)){
for(int i=1;i<=n;i++){
scanf("%d%d",&s[i].l,&s[i].r);
po[i].l=s[i].l ,po[i].r=i;
}a[n+1]=-inf;
sort(s+1,s+1+n);
sort(po+1,po+1+n);
for(int i=1;i<=n;i++) a[i]=s[i].l;
for(int i=1;i<=n;i++){
int l=s[i].l,r=s[i].r;
s[i].l=lower_bound(a+1,a+1+n,l-r)-a;
s[i].r=lower_bound(a+1,a+1+n,l+r)-a;
if(a[s[i].r]!=l+r) s[i].r--;
}
build(1,n,1);
for(int i=1;i<=n;i++) a[i]=i;
sort(a+1,a+1+n,cmp);
for(int i=1;i<=n;i++){
int id=0,t=a[i];
sec num[2]={{t,t},{0,0}};
while(num[0].l!=num[1].l || num[0].r!=num[1].r){
num[id^1]=query(num[id].l,num[id].r,1,n,1);
id^=1;
}
update(t,t,1,n,1,num[id].l,num[id].r);//更新t点感染范围
ans[po[t].r]=num[id].r-num[id].l+1;//找回原来的位置赋值
}printf("%d",ans[1]);
for(int i=2;i<=n;i++) printf(" %d",ans[i]);
puts("");
}
return 0;
}