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=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;
}
ll n,m,p[MAXN],rd[MAXN],ans,hsh[MAXN],sl[MAXN<<2],sr[MAXN<<2],l[MAXN<<2],r[MAXN<<2];
int dfn[MAXN<<2],low[MAXN<<2],st[MAXN<<2],bl[MAXN<<2],stp,top,scc;
int q[MAXN<<2],hd=1,tl,ind[MAXN<<2];
struct graph
{
    int cnt,fst[MAXN<<2],nxt[MAXM<<1],to[MAXM<<1],ind[MAXN<<2];
    graph(){memset(fst,0,sizeof(fst));cnt=0;}
    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>>1;
    build(k<<1,l,mid);build(k<<1|1,mid+1,r);
    g1.add(k,k<<1);g1.add(k,k<<1|1);
}
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>>1;
    if(b<=mid) mdf(k<<1,l,mid,a,b,x);
    else if(a>mid) mdf(k<<1|1,mid+1,r,a,b,x);
    else {mdf(k<<1,l,mid,a,mid,x);mdf(k<<1|1,mid+1,r,mid+1,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=0;
        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=1;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(1,1,n);int a,b;
    for(int i=1;i<=n;i++) p[i]=read(),rd[i]=read();
    for(int i=1;i<=n;i++)
    {
        a=lower_bound(p+1,p+n+1,p[i]-rd[i])-p;
        b=upper_bound(p+1,p+n+1,p[i]+rd[i])-p-1;
        mdf(1,1,n,a,b,hsh[i]);
    }
    for(int i=1;i<=m;i++) if(!dfn[i]) tarjan(i);
    build();for(int i=1;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=1;i<=n;i++) (ans+=(i*(r[bl[hsh[i]]]-l[bl[hsh[i]]]+1))%MOD)%=MOD;
    printf("%lld
",ans);
}