数据结构中的判定转状态+扫描线:P1502

https://www.luogu.com.cn/problem/P1502

发现正常扫描线很难维护恰好大小为 W W W 的区间

反过来,对于每个星星维护合法的左下角下标

把原先的判定转成了和点有关的状态,把点变成矩形后求并即可

cpp 复制代码
#include <algorithm>
#include<iostream>
#include <stdio.h>
using namespace std;
#define int long long
inline int read(){int x=0,f=1;char ch=getchar(); while(ch<'0'||
ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){
x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}return x*f;}
#define Z(x) (x)*(x)
#define pb push_back
//mt19937 rand(time(0));
//mt19937_64 rand(time(0));
//srand(time(0));
#define N 20010
//#define M
//#define mo
struct node {
	int l, r, x, w; 
}a[N];
int n, m, i, j, k, T;
int ans, sum, l, r, x, y, w, Lx, Rx; 
int W, H, rt, b[N]; 

bool cmp(node x, node y) {
	if(x.x==y.x) return x.w<y.w; 
	return x.x<y.x; 
}

struct Segment_tree {
	int tot, ls[N<<2], rs[N<<2]; 
	int mx[N<<2], tag[N<<2]; 
	void build(int &k, int l, int r) {
		if(!k) k=++tot, mx[k]=tag[k]=ls[k]=rs[k]=0; 
		if(l==r) return ; 
		int mid=(l+r)>>1; 
		build(ls[k], l, mid); 
		build(rs[k], mid+1, r); 
	}
	void push_down(int k) {
		tag[ls[k]]+=tag[k]; mx[ls[k]]+=tag[k]; 
		tag[rs[k]]+=tag[k]; mx[rs[k]]+=tag[k]; 
		tag[k]=0; 
	}
	void add(int k, int l, int r, int x, int y, int z) {
		if(l>=x && r<=y) return mx[k]+=z, tag[k]+=z, void(); 
		int mid=(l+r)>>1; 
		push_down(k); 
		if(x<=mid) add(ls[k], l, mid, x, y, z); 
		if(y>=mid+1) add(rs[k], mid+1, r, x, y, z); 
		mx[k]=max(mx[ls[k]], mx[rs[k]]); 
	}
}Seg;

signed main()
{
//	freopen("in.txt", "r", stdin);
//	freopen("out.txt", "w", stdout);
	T=read(); 
	while(T--) {
//	while(~scanf("%lld%lld%lld", &n, &W, &H)) {
		n=read(); W=read(); H=read(); 
	
		ans=k=j=0; 
		for(i=1; i<=n; ++i) {
			x=read(); y=read(); w=read();
			Lx=x; Rx=x+W-1; //[Lx, Rx]
			b[++k]=Lx; b[++k]=Rx; 
			a[++j].l=Lx; a[j].r=Rx; a[j].x=y; a[j].w=w; 
			a[++j].l=Lx; a[j].r=Rx; a[j].x=y+H; a[j].w=-w; 
		}
		sort(b+1, b+2*n+1); 
		for(i=1; i<=2*n; ++i) a[i].l=lower_bound(b+1, b+2*n+1, a[i].l)-b; 
		for(i=1; i<=2*n; ++i) a[i].r=lower_bound(b+1, b+2*n+1, a[i].r)-b; 
		Seg.tot=rt=ans=0; Seg.build(rt, 1, 2*n); 
		
		sort(a+1, a+2*n+1, cmp); 
		
		for(i=1; i<=2*n; ++i) {
			Seg.add(1, 1, 2*n, a[i].l, a[i].r, a[i].w); 
			ans=max(ans, Seg.mx[1]); 
		}
		printf("%lld\n", ans); 
	}
	return 0;
}
相关推荐
cpp_25011 小时前
P1024 [NOIP 2001 提高组] 一元三次方程求解
数据结构·c++·算法·题解·二分答案·洛谷·csp
AI人工智能+电脑小能手9 小时前
【大白话说Java面试题】【Java基础篇】第15题:JDK1.7中HashMap扩容为什么会发生死循环?如何解决
java·开发语言·数据结构·后端·面试·哈希算法
Mrlxl.cn11 小时前
计算机网络——网络层
c语言·数据结构·计算机网络·考研
寒秋花开曾相惜13 小时前
(学习笔记)4.2 逻辑设计和硬件控制语言HCL(4.2.1 逻辑门&4.2.2 组合电路和HCL布尔表达式)
linux·网络·数据结构·笔记·学习·fpga开发
码完就睡15 小时前
数据结构——哈希表原理与C语言实现总结
数据结构·散列表
山甫aa17 小时前
二叉树遍历----从零开始的数据结构
数据结构·c++·二叉树
北顾笙98017 小时前
day35-数据结构力扣
数据结构·算法·leetcode
cpp_250118 小时前
P2249 【深基13.例1】查找
数据结构·c++·算法·题解·二分·洛谷
山甫aa18 小时前
二叉树算法-----从零开始的算法
数据结构·算法
睡觉就不困鸭19 小时前
第十七天 翻转字符串里的单词
数据结构·算法·哈希算法·散列表