bzoj5017 炸弹 (线段树优化建图+tarjan+拓扑序dp)

直接建图边数太多，用线段树优化一下

然后缩点，记下来每个点里有多少个炸弹

然后按拓扑序反向dp一下就行了

#include<bits/stdc++.h>
#define pa pair<ll,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=5e5+10,maxp=maxn*4,maxl=maxn*30,P=1e9+7;
const ll inf=1e18+1;

inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}

int N;
ll pos[maxn],r[maxn];
int rg[maxn][2],id[maxn],ch[maxp][2],pct,root;
int eg[maxl][2],egh[maxp],ect;
int eg2[maxl][2],egh2[maxp],ect2,ine[maxp];
int dfn[maxp],low[maxp],tot,stk[maxp],sh;
int bel[maxp],nct,hav[maxp],hh[maxp];
int val[maxp],rk[maxp];
bool instk[maxp],islef[maxp],to[maxp];
queue<int> q;

inline void adeg(int a,int b){
    eg[++ect][0]=b;eg[ect][1]=egh[a];egh[a]=ect;
}
inline void adeg2(int a,int b){
    ine[b]++;
    eg2[++ect2][0]=b;eg2[ect2][1]=egh2[a];egh2[a]=ect2;
}

void build(int &p,int l,int r){
    p=++pct;
    if(l==r) id[l]=p,islef[p]=1;
    if(l<r){
        int m=l+r>>1;
        build(ch[p][0],l,m);
        build(ch[p][1],m+1,r);
        adeg(p,ch[p][0]),adeg(p,ch[p][1]);
    }
}

void conn(int p,int l,int r,int x,int y,int z){
    if(x>y) return;
    if(x<=l&&r<=y){
        adeg(z,p);
    }else{
        int m=l+r>>1;
        if(x<=m) conn(ch[p][0],l,m,x,y,z);
        if(y>=m+1) conn(ch[p][1],m+1,r,x,y,z);
    }
}

void tarjan(int x){
    dfn[x]=low[x]=++tot;
    stk[++sh]=x;instk[x]=1;
    for(int i=egh[x];i;i=eg[i][1]){
        int b=eg[i][0];
        if(instk[b]) low[x]=min(low[x],dfn[b]);
        else if(!dfn[b]){
            tarjan(b);
            low[x]=min(low[x],low[b]);
        }
    }
    if(dfn[x]==low[x]){
        ++nct;
        while(sh){
            bel[stk[sh]]=nct;
            hav[stk[sh]]=hh[nct],hh[nct]=stk[sh];
            val[nct]+=islef[stk[sh]];
            instk[stk[sh]]=0;
            sh--;
            if(stk[sh+1]==x) break;
        }
    }
}

int main(){
    //freopen(".in","r",stdin);
    int i,j,k;
    N=rd();
    for(i=1;i<=N;i++) pos[i]=rd(),r[i]=rd();
    for(i=1;i<=N;i++){
        rg[i][0]=lower_bound(pos+1,pos+N+1,pos[i]-r[i])-pos;
        rg[i][1]=upper_bound(pos+1,pos+N+1,pos[i]+r[i])-pos-1;
        // printf("%d %d %d
",i,rg[i][0],rg[i][1]);
    }
    build(root,1,N);
    for(i=1;i<=N;i++){
        conn(root,1,N,rg[i][0],i-1,id[i]);
        conn(root,1,N,i+1,rg[i][1],id[i]);
    }
    for(i=1;i<=pct;i++){
        if(!dfn[i]) tarjan(i);
    }
    for(i=1;i<=nct;i++){
        for(j=hh[i];j;j=hav[j]){
            for(k=egh[j];k;k=eg[k][1]){
                int bb=bel[eg[k][0]];if(bb==i) continue;
                if(!to[bb]) adeg2(i,bb),to[bb]=1;
            }
        }
        
        for(j=hh[i];j;j=hav[j]){
            for(k=egh[j];k;k=eg[k][1]){
                int bb=bel[eg[k][0]];
                to[bb]=0;
            }
        }
    }
    int nn=0;
    for(i=1;i<=nct;i++){
        if(!ine[i]) q.push(i);
    }
    while(!q.empty()){
        int p=q.front();q.pop();
        rk[++nn]=p;
        for(i=egh2[p];i;i=eg2[i][1]){
            int b=eg2[i][0];
            if(--ine[b]==0) q.push(b);
        }
        
    }
    for(i=nct;i;i--){
        int p=rk[i];
        for(j=egh2[p];j;j=eg2[j][1]){
            int b=eg2[j][0];
            val[p]+=val[b];
        }
    }
    ll ans=0;
    for(i=1;i<=N;i++){
        ans=(ans+1ll*i*val[bel[id[i]]])%P;
    }
    printf("%lld
",ans);
    return 0;
}