~~~~~ P5025 [SNOI2017] 炸弹 ~~~~~ 总题单链接
思路
~~~~~ 把每个点与其爆炸范围内的点连边,用线段树优化建图可以做到 O ( N ∗ l o g ( N ) ) O(N*log(N)) O(N∗log(N))。
~~~~~ 发现在同一个连通块里的点可以互相到达,所以缩点。
~~~~~ 所点后 D P DP DP 求得每个连通块能到达的左右端点(也可以用拓扑排序或搜索,但 D P DP DP 码量最少),这里给出 D P DP DP 过程的代码:
for(ll u=1;u<=tr.tot;u++)
for(ll v:eg[u])
if(scc[u]!=scc[v])
ng[scc[v]].push_back(scc[u]);
for(ll i=1;i<=tr.tot;i++)dp[i]={INF,-INF};
for(ll i=1;i<=n;i++)dp[scc[i]]={min(dp[scc[i]].fir,i),max(dp[scc[i]].sec,i)};
for(ll u=1;u<=cnt;u++){
for(ll v:ng[u])
dp[v]={min(dp[v].fir,dp[u].fir),max(dp[v].sec,dp[u].sec)};
}
代码
#include<bits/stdc++.h>
#define ll long long
#define fir first
#define sec second
#define MOD 1000000007
#define INF 0x3f3f3f3f3f3f3f
using namespace std;
vector<ll>eg[25000005];
ll n,g[500005],r[500005];
class{
public:
ll tot=n;
struct Point{
ll L,R,id;
}po[2000005];
void build(ll x=1,ll y=n,ll p=1){
po[p]={x,y,++tot};
if(x==y){
eg[po[p].id].push_back(x);
return;
}
ll mid=(x+y)>>1;
build(x,mid,p<<1);
build(mid+1,y,p<<1|1);
eg[po[p].id].push_back(po[p<<1].id);
eg[po[p].id].push_back(po[p<<1|1].id);
}
void query(ll wc,ll x,ll y,ll p=1){
if(po[p].R<x||po[p].L>y)return;
if(po[p].L>=x&&po[p].R<=y){
eg[wc].push_back(po[p].id);
return;
}
query(wc,x,y,p<<1);
query(wc,x,y,p<<1|1);
}
}tr;
class{
public:
ll dfn[2500005],low[2500005],tot;
ll scc[2500005],siz[2500005],cnt;
ll stk[2500005],ins[2500005],top;
vector<ll>ng[2500005];pair<ll,ll>dp[200005];
void Tarjan(ll p){
dfn[p]=low[p]=++tot;
stk[++top]=p;ins[p]=1;
for(ll v:eg[p]){
if(!dfn[v]){
Tarjan(v);
low[p]=min(low[p],low[v]);
}
else if(ins[v])low[p]=min(low[p],dfn[v]);
}
if(dfn[p]==low[p]){
cnt++;
while(1){
ll v=stk[top--];
ins[v]=0;
scc[v]=cnt;
if(v<=n)siz[cnt]++;
if(v==p)break;
}
}
}
void work(){
for(ll i=1;i<=tr.tot;i++)
if(!dfn[i])Tarjan(i);
for(ll u=1;u<=tr.tot;u++)
for(ll v:eg[u])
if(scc[u]!=scc[v])
ng[scc[v]].push_back(scc[u]);
for(ll i=1;i<=tr.tot;i++)dp[i]={INF,-INF};
for(ll i=1;i<=n;i++)dp[scc[i]]={min(dp[scc[i]].fir,i),max(dp[scc[i]].sec,i)};
for(ll u=1;u<=cnt;u++){
for(ll v:ng[u])
dp[v]={min(dp[v].fir,dp[u].fir),max(dp[v].sec,dp[u].sec)};
}
ll ans=0;
for(ll i=1;i<=n;i++){
(ans+=i*(dp[scc[i]].sec-dp[scc[i]].fir+1)%MOD)%=MOD;
}
cout<<ans;
}
}TJ;
signed main(){
ios::sync_with_stdio(false);
cin>>n;
tr.tot=n;tr.build();
for(ll i=1;i<=n;i++)
cin>>g[i]>>r[i];
for(ll i=1;i<=n;i++){
ll L=lower_bound(g+1,g+1+n,g[i]-r[i])-g;
ll R=upper_bound(g+1,g+1+n,g[i]+r[i])-g-1;
if(L<=i-1)tr.query(i,L,i-1);
if(i+1<=R)tr.query(i,i+1,R);
}
TJ.work();
return 0;
}