- 题意:P5024
- 思路:首先\(O(n^2)\)向能炸到的点连边,所以能到达的点的个数就是能到达的点的个数。然后显然要缩点+拓扑排序(我写的记搜)。
然后再写一个线段树优化建图。
然后就WA了,我想了很久才明白我的记搜有问题,一个点可能用来更新另一个点多次。我看了下题解,题解很nice,也很好写,就维护一个联通分量能到的范围[左右端点](此题的性质保证能到的点集是一个子段)。
- 代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2e6+5;
const int M=4e7+5;
struct seg {int l,r;}T[N];
struct node {ll p,r;}a[N];
bool cmp(node u,node v) {return u.p<v.p;}
bool In_s[N];
ll sum[N],mod=1e9+7,mp[N],cnt[N],ans[N];
int L[N],R[N],mx[N],mn[N],r[N],o[N],tim,Ecnt,In[N],Bl[N],mark[N],SCC,mxx,lf[N],Time,low[N],dfn[N],st[N],tp,n,nxt[M],to[M],head[N],ecnt;
void add_edge(int u,int v) {nxt[++ecnt]=head[u];to[ecnt]=v;head[u]=ecnt;}
void Build(int x,int l,int r) {
T[x]=(seg){l,r},mxx=max(mxx,x);
if(l==r) {lf[l]=x,mark[x]=l;return;}
int mid=(l+r)>>1,ls=x<<1,rs=x<<1|1;
add_edge(x,ls),add_edge(x,rs);
Build(ls,l,mid),Build(rs,mid+1,r);
}
void Add(int x,int u,int l,int r) {
if(l<=T[x].l&&T[x].r<=r) {add_edge(lf[u],x);return;}
int mid=(T[x].l+T[x].r)>>1;
if(l<=mid)Add(x<<1,u,l,r);
if(r>mid)Add(x<<1|1,u,l,r);
}
void Tarjan(int u) {
dfn[u]=low[u]=++Time,In_s[u]=1;
st[++tp]=u;
for(int i=head[u];i;i=nxt[i]) {
int v=to[i];
if(!dfn[v]) {
Tarjan(v);
low[u]=min(low[u],low[v]);
}
else if(In_s[v]) low[u]=min(low[u],dfn[v]);
}
if(low[u]==dfn[u]) {
int v;
SCC++;
L[SCC]=1e9;
do {
v=st[tp--],L[SCC]=min(L[SCC],T[v].l),R[SCC]=max(R[SCC],T[v].r),sum[SCC]+=mark[v],Bl[v]=SCC,In_s[v]=0;
}while(u!=v);
}
}
vector<int> G[N];
void dfs(int u) {
if(ans[u]) return;
ans[u]=1;
for(int i=0;i<G[u].size();i++) {
int v=G[u][i];
dfs(v);
L[u]=min(L[v],L[u]),R[u]=max(R[u],R[v]);
}
}
int main() {
scanf("%d",&n);
for(int i=1;i<=n;i++) {
scanf("%lld%lld",&a[i].p,&a[i].r);
}
Build(1,1,n);
sort(a+1,a+1+n,cmp);
for(int i=1;i<=n;i++) mp[i]=a[i].p;
for(int i=1;i<=n;i++) {
int l=lower_bound(mp+1,mp+1+i,a[i].p-a[i].r)-mp;
int r=upper_bound(mp+1+i,mp+1+n,a[i].p+a[i].r)-mp-1;
Add(1,i,l,r);
}
Tarjan(1);
n=mxx;
for(int u=1;u<=n;u++) {
for(int i=head[u];i;i=nxt[i]) {
int v=to[i];
if(Bl[u]!=Bl[v])G[Bl[u]].push_back(Bl[v]),In[Bl[v]]++;
}
}
for(int i=1;i<=SCC;i++) {
if(!ans[i]&&!In[i]) {
dfs(i);
}
}
ll res=0;
for(int i=1;i<=SCC;i++) {
// printf("%d %d\n",L[i],R[i]);
res=(res+1ll*sum[i]%mod*(R[i]-L[i]+1)%mod)%mod;
}
printf("%lld",res);
return 0;
}