bzoj 5017 炸弹

题目大意：

直线上有n个炸弹有坐标x和半径r

当一个炸弹被引爆时 若有炸弹的坐标在该炸弹坐标+-r范围内则另一个炸弹也被引爆

求先引爆每一个炸弹最终会引爆多少炸弹

思路：

可以想到n平方连边然后tarjan缩点跑拓扑

可以通过线段树来优化建图

对每个点向它能直接引爆的左右范围连边

即用线段树中的线段作为点来建图

 #include<iostream>
 #include<cstdio>
 #include<cmath>
 #include<cstdlib>
 #include<cstring>
 #include<algorithm>
 #include<vector>
 #include<queue>
 #define inf 2139062143
 #define ll long long
 #define MOD 1000000007
 #define MAXN 500100
 #define MAXM 10010000
 #define V1 g1.to[i]
 #define V2 g2.to[i]
 using namespace std;
 inline ll read()
 {
     ll x=,f=;char ch=getchar();
     while(!isdigit(ch)) {if(ch=='-') f=-;ch=getchar();}
     while(isdigit(ch)) {x=x*+ch-'';ch=getchar();}
     return x*f;
 }
 ll n,m,p[MAXN],rd[MAXN],ans,hsh[MAXN],sl[MAXN<<],sr[MAXN<<],l[MAXN<<],r[MAXN<<];
 int dfn[MAXN<<],low[MAXN<<],st[MAXN<<],bl[MAXN<<],stp,top,scc;
 int q[MAXN<<],hd=,tl,ind[MAXN<<];
 struct graph
 {
     int cnt,fst[MAXN<<],nxt[MAXM<<],to[MAXM<<],ind[MAXN<<];
     graph(){memset(fst,,sizeof(fst));cnt=;}
     inline void add(int u,int v) {nxt[++cnt]=fst[u],fst[u]=cnt,to[cnt]=v,ind[v]++;}
 }g1,g2;
 void build(int k,int l,int r)
 {
     sl[k]=l,sr[k]=r;
     if(l==r) {hsh[l]=k,m=k;return ;}
     int mid=l+r>>;
     build(k<<,l,mid);build(k<<|,mid+,r);
     g1.add(k,k<<);g1.add(k,k<<|);
 }
 void mdf(int k,int l,int r,int a,int b,int x)
 {
     if(l==a&&r==b) {if(l!=r||(l==r&&k!=x)) g1.add(x,k);return ;}
     int mid=l+r>>;
     if(b<=mid) mdf(k<<,l,mid,a,b,x);
     else if(a>mid) mdf(k<<|,mid+,r,a,b,x);
     else {mdf(k<<,l,mid,a,mid,x);mdf(k<<|,mid+,r,mid+,b,x);}
 }
 void tarjan(int x)
 {
     dfn[x]=low[x]=++stp,st[++top]=x;
     for(int i=g1.fst[x];i;i=g1.nxt[i])
         if(!dfn[V1]) {tarjan(V1);low[x]=min(low[x],low[V1]);}
         else if(!bl[V1]) low[x]=min(low[x],dfn[V1]);
     if(low[x]==dfn[x])
     {
         l[++scc]=inf;int now=;
         while(now!=x) now=st[top--],bl[now]=scc,l[scc]=min(l[scc],sl[now]),r[scc]=max(r[scc],sr[now]);
     }
 }
 void build()
 {
     for(int x=;x<=m;x++)
         for(int i=g1.fst[x];i;i=g1.nxt[i])
             if(bl[x]!=bl[V1]) g2.add(bl[x],bl[V1]);
 }
 int main()
 {
     n=read();build(,,n);int a,b;
     for(int i=;i<=n;i++) p[i]=read(),rd[i]=read();
     for(int i=;i<=n;i++)
     {
         a=lower_bound(p+,p+n+,p[i]-rd[i])-p;
         b=upper_bound(p+,p+n+,p[i]+rd[i])-p-;
         mdf(,,n,a,b,hsh[i]);
     }
     for(int i=;i<=m;i++) if(!dfn[i]) tarjan(i);
     build();for(int i=;i<=scc;i++) if(!g2.ind[i]) q[++tl]=i;
     while(hd<=tl)
     {
         a=q[hd++];
         for(int i=g2.fst[a];i;i=g2.nxt[i])
             {g2.ind[V2]--;if(!g2.ind[V2]) q[++tl]=V2;}
     }
     for(int x=tl;x;x--)
         for(int i=g2.fst[q[x]];i;i=g2.nxt[i])
             l[q[x]]=min(l[q[x]],l[V2]),r[q[x]]=max(r[q[x]],r[V2]);
     for(int i=;i<=n;i++) (ans+=(i*(r[bl[hsh[i]]]-l[bl[hsh[i]]]+))%MOD)%=MOD;
     printf("%lld\n",ans);
 }