数据结构中的判定转状态+扫描线: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;
}
相关推荐
刘马想放假3 天前
Modbus 全栈技术解析:TCP、RTU、ASCII、RTU over TCP
数据结构·网络协议
北域码匠3 天前
冒泡排序太慢?鸡尾酒排序双向优化,原生 C# 零第三方库完整代码
数据结构·排序算法·泛型·c# 算法·鸡尾酒排序·原生 c# 开发·冒泡排序优化·嵌入式算法
Darling噜啦啦10 天前
列表转树算法深度解析:从 Map 到 Reduce 的两种实现,面试高频考点
数据结构·算法·面试
小小工匠11 天前
Redis - 事务机制:能实现 ACID 属性吗
数据结构·redis·性能优化·并发·持久化
玖玥拾11 天前
C/C++ 数据结构(七)栈、容器适配器
c语言·数据结构·c++··容器适配器
Qres82111 天前
算法复键——树状数组
数据结构·算法
牛油果子哥q11 天前
并查集(DSU)超精讲,路径压缩、按秩合并、万能模板、连通性判定、最小生成树与刷题实战全解
数据结构·c++·最小生成树·并查集
凌波粒12 天前
LeetCode--491.递增子序列(回溯算法)
数据结构·算法·leetcode
WL学习笔记12 天前
单项不带头不循环链表
数据结构·链表
小糯米60112 天前
JS 数组
数据结构·算法·排序算法