[bzoj5017][Snoi2017]炸弹 tarjan缩点+线段树优化建图+拓扑

5017: [Snoi2017]炸弹

Time Limit: 30 Sec  Memory Limit: 512 MB
Submit: 608  Solved: 190
[Submit][Status][Discuss]

Description

在一条直线上有 N 个炸弹，每个炸弹的坐标是 Xi，爆炸半径是 Ri，当一个炸弹爆炸时，如果另一个炸弹所在位置 Xj 满足： 

Xi−Ri≤Xj≤Xi+Ri,那么，该炸弹也会被引爆。 

现在，请你帮忙计算一下，先把第 i 个炸弹引爆，将引爆多少个炸弹呢？ 

 

Input

第一行，一个数字 N，表示炸弹个数。 

第 2∼N+1行，每行 2 个数字，表示 Xi，Ri，保证 Xi 严格递增。 

N≤500000

−10^18≤Xi≤10^18

0≤Ri≤2×10^18

 

Output

一个数字，表示Sigma(i*炸弹i能引爆的炸弹个数),1<=i<=N mod10^9+7。 

 

Sample Input

4
1 1
5 1
6 5
15 15

Sample Output

32

HINT

 

Source

 

 

显然一个点可以引爆的炸弹是一个连续的区间，对于每一个点我们向他可以引爆的最左边的炸弹到最右边的炸弹连边。可以使用线段树优化建图。

我们对建出来的图tarjan缩点，每个点维护minl和maxr。

对于这个图按拓扑序倒序dp，求出每个点的minl和maxr，统计答案即可。

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<cstdio>
#include<algorithm>
#define maxn 500005
#define ll long long
#define mod 1000000007
using namespace std;
inline ll read() {
    char ch=getchar();ll x=0,f=1;
    for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-1;
    for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';
    return x*f;
}
int n,rt;
ll s[maxn*5],a[maxn*5],id[maxn*5],tsz;
struct seg {int s[2];ll mn,mx;}t[maxn*5];
struct Edge {int fr,to,nxt;}e[maxn*20];
int head[maxn*5],sz;
void addedge(int u,int v) {e[sz].fr=u;e[sz].nxt=head[u];e[sz].to=v;head[u]=sz++;}
void build(int l,int r,int &x) {
    x=++tsz;
    if(l==r) {id[l]=x;t[x].mn=t[x].mx=l;return;}
    int mid=l+r>>1;
    build(l,mid,t[x].s[0]);build(mid+1,r,t[x].s[1]);
    if(t[x].s[0]) addedge(x,t[x].s[0]);
    if(t[x].s[1]) addedge(x,t[x].s[1]);
    t[x].mn=l;t[x].mx=r;
    return ;
}
void add(int l,int r,int x,int L,int R,int p) {
    if(L<=l&&R>=r) {if(id[p]==x) return;addedge(id[p],x);return;}
    int mid=l+r>>1;
    if(L<=mid) add(l,mid,t[x].s[0],L,R,p);
    if(R>mid) add(mid+1,r,t[x].s[1],L,R,p);
    return;
}
bool inq[maxn*5];
ll dfn[maxn*5],low[maxn*5],tim,q[maxn*5],top,lm[maxn*5],rm[maxn*5];
int bel[maxn*5],scc;
void tarjan(int x) {
    dfn[x]=low[x]=++tim;q[++top]=x;inq[x]=1;
    for(int i=head[x];i>=0;i=e[i].nxt) {
        int to=e[i].to;
        if(!dfn[to]){
            tarjan(to);low[x]=min(low[x],low[to]);
        }else if(inq[to]) low[x]=min(low[x],dfn[to]);
    }
    if(dfn[x]==low[x]) {
        scc++;
        lm[scc]=214748364700000000ll;
        ll now;
        do {
            now=q[top--];inq[now]=0;bel[now]=scc;
            lm[scc]=min(lm[scc],t[now].mn);rm[scc]=max(rm[scc],t[now].mx);
        }while(now!=x);
    }
}
ll rd[maxn*5];
int main() {
    //freopen("bomb9.in","r",stdin);
    memset(head,-1,sizeof(head));
    n=read();
    for(int i=1;i<=n;i++) {s[i]=read();a[i]=read();}
    build(1,n,rt);
    for(int i=1;i<=n;i++) {
        int now=lower_bound(s+1,s+n+1,s[i])-s,L=lower_bound(s+1,s+n+1,s[i]-a[i])-s,R=upper_bound(s+1,s+n+1,s[i]+a[i])-s-1;
        if(L==R) continue;add(1,n,1,L,R,i);
    }
    for(int i=1;i<=tsz;i++) if(!dfn[i]) {tarjan(i);}
    int tmp=sz;sz=0;memset(head,-1,sizeof(head));
    for(int i=0;i<tmp;i++) {
        int u=bel[e[i].fr],v=bel[e[i].to];
        if(u==v) continue;
        rd[v]++;addedge(u,v);
    }
    int hd=0,tl=0;
    for(int i=1;i<=scc;i++) if(!rd[i]) q[tl++]=i;
    while(hd!=tl) {
        int now=q[hd++];
        for(int i=head[now];i>=0;i=e[i].nxt) {
            int to=e[i].to;rd[to]--;
            if(!rd[to]) q[tl++]=to;
        }
    }
    for(int i=scc;i>=1;i--) {
        int now=q[i];
        for(int j=head[now];j>=0;j=e[j].nxt) {
            int to=e[j].to;
            lm[now]=min(lm[now],lm[to]);rm[now]=max(rm[now],rm[to]);
        }
    }
    ll ans=0;
    for(int i=1;i<=n;i++) {
        ans+=(ll)i*(rm[bel[id[i]]]-lm[bel[id[i]]]+1);ans%=mod;
    }
    printf("%lld\n",ans);
}