求助!!!线段树维护最远左右端点为什么86分,WA#1#2?
查看原帖
求助!!!线段树维护最远左右端点为什么86分,WA#1#2?
326639
zyc070419楼主2022/7/15 15:46

变量名说明:

1、结构体 aa 存炸弹信息,结构体中变量 aa 是位置, bb 是半径, l,rl,r 是炸弹能引爆的炸弹的范围,idid 是炸弹编号。

2、 f[i]f[i] 表示引爆第 ii 个炸弹后所有能引爆到的炸弹的最左边的炸弹的编号, g[i]g[i] 是最右边。

3、线段树内存区间 [l, r][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;
}

蒟蒻实在不知道哪里错了

2022/7/15 15:46
加载中...