变量名说明:
1、结构体 a 存炸弹信息,结构体中变量 a 是位置, b 是半径, l,r 是炸弹能引爆的炸弹的范围,id 是炸弹编号。
2、 f[i] 表示引爆第 i 个炸弹后所有能引爆到的炸弹的最左边的炸弹的编号, g[i] 是最右边。
3、线段树内存区间 [l, r] 之间的炸弹的最左、最右爆炸范围。
#include<bits/stdc++.h>
#define ls (rt << 1)
#define rs ((rt << 1) | 1)
#define INF 1000000000000000000
#define int long long
using namespace std;
const int N = 1e6 + 30;
const int mod = 1e9 + 7;
struct node{
int id, a, b, l, r;
}a[N];
struct segtree{
int mn, mx;
}t[N << 2];
int n, p[N], f[N], g[N], ans;
deque< pair<int, int> > q;
bool cmpl(node x, node y){
return x.l < y.l;
}
bool cmpr(node x, node y){
return x.r > y.r;
}
void build(int rt, int l, int r){
t[rt].mn = INF; t[rt].mx = -INF;
if(l == r) return;
int mid = (l + r) >> 1;
build(ls, l, mid); build(rs, mid + 1, r);
}
int query(int rt, int l, int r, int ql, int qr, int opt){
if(ql <= l && r <= qr){
if(!opt) return t[rt].mn;
else return t[rt].mx;
}
int mid = (l + r) >> 1, res;
if(!opt){
res = INF;
if(ql <= mid) res = min(res, query(ls, l, mid, ql, qr, opt));
if(qr > mid) res = min(res, query(rs, mid + 1, r, ql, qr, opt));
}else{
res = -INF;
if(ql <= mid) res = max(res, query(ls, l, mid, ql, qr, opt));
if(qr > mid) res = max(res, query(rs, mid + 1, r, ql, qr, opt));
}
return res;
}
void update(int rt, int l, int r, int p, int val, int opt){
if(l == r){
if(!opt) t[rt].mn = min(t[rt].mn, val);
else t[rt].mx = max(t[rt].mx, val);
return;
}
int mid = (l + r) >> 1;
if(p <= mid) update(ls, l, mid, p, val, opt);
else update(rs, mid + 1, r, p, val, opt);
t[rt].mn = min(t[ls].mn, t[rs].mn);
t[rt].mx = max(t[ls].mx, t[rs].mx);
}
int Mod(int x){
return x >= mod ? x - mod : x;
}
signed main(){
int tmp;
scanf("%lld", &n);
for(int i = 1; i <= n; ++i){
scanf("%lld%lld", &a[i].a, &a[i].b);
a[i].id = i;
p[i] = a[i].a;
}
build(1, 1, n);
for(int i = 1; i <= n; ++i){
tmp = a[i].a - a[i].b;
a[i].l = lower_bound(p + 1, p + 1 + n, tmp) - p;
tmp = a[i].a + a[i].b;
a[i].r = upper_bound(p + 1, p + 1 + n, tmp) - p - 1;
}
sort(a + 1, a + 1 + n, cmpl);
for(int i = 1; i <= n; ++i){
f[a[i].id] = min(a[i].l, query(1, 1, n, a[i].l, a[i].r, 0));
update(1, 1, n, a[i].id, f[a[i].id], 0);
}
sort(a + 1, a + 1 + n, cmpr);
for(int i = 1; i <= n; ++i){
g[a[i].id] = max(a[i].r, query(1, 1, n, a[i].l, a[i].r, 1));
update(1, 1, n, a[i].id, g[a[i].id], 1);
}
for(int i = 1; i <= n; ++i) ans = Mod(ans + ((Mod(g[i] - f[i] + 1) * i) % mod));
printf("%lld\n", ans);
return 0;
}
蒟蒻实在不知道哪里错了