【BZOJ】【P3165】【Heoi2013】【Segment】【题解】【线段树】

最新推荐文章于 2022-09-23 17:02:26 发布

iamzky

最新推荐文章于 2022-09-23 17:02:26 发布

阅读量976

点赞数

CC 4.0 BY-SA版权

分类专栏： OI 文章标签： bzoj

本文链接：https://blog.youkuaiyun.com/iamzky/article/details/43238021

OI 专栏收录该内容

324 篇文章

订阅专栏

本文介绍了一种使用线段树进行区间更新和查询的高效算法实现，通过具体实例展示了如何构建线段树并实现区间最大值的快速查找。特别针对斜率优化问题，提供了完整的代码示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3165

题解

Code:

#include<bits/stdc++.h>
using namespace std;
const int maxn=40001;
int n=39989,m,ind;
struct seg{
	double k,b;
	int id;
	seg(int x0=0,int y0=0,int x1=0,int y1=0,int _id=0){
		id=_id;
		if(x0==x1){k=0;b=max(y0,y1);}
		else{k=(double)(y0-y1)/(x0-x1);b=y0-k*x0;}
	}
	double get(double x){return k*x+b;}
};
bool les(seg A,seg B,double x){
	if(!A.id)return 1;
	return A.get(x)!=B.get(x)?A.get(x)<B.get(x):A.id<B.id;
}
struct sgt{
	seg t[maxn<<2];
	seg Qmax(int i,int l,int r,int ps){
		if(l==r)return t[i];seg tmp;
		if(ps<=(l+r)/2)tmp=Qmax(i<<1,l,(l+r)/2,ps);
		else tmp=Qmax(i<<1|1,(l+r)/2+1,r,ps);
		return les(t[i],tmp,ps)?tmp:t[i];
	}
	void insert(int i,int l,int r,seg se){
		if(!t[i].id)t[i]=se;
		if(les(t[i],se,l))swap(t[i],se);
		if(l==r||t[i].k==se.k)return;
		double x=(t[i].b-se.b)/(se.k-t[i].k);
		if(x<l||x>r)return;
		if(x<=(l+r)/2)insert(i<<1,l,(l+r)/2,t[i]),t[i]=se;
		else insert(i<<1|1,(l+r)/2+1,r,se);
	}
	void insert(int i,int l,int r,int l0,int r0,seg se){
		if(l0<=l&&r0>=r){insert(i,l,r,se);return;}
		if(l0<=(l+r)/2)insert(i<<1,l,(l+r)/2,l0,r0,se);
		if(r0>(l+r)/2)insert(i<<1|1,(l+r)/2+1,r,l0,r0,se);
	}	
}T;
int la;
int main(){
	scanf("%d",&m);
	while(m--){
		int op;scanf("%d",&op);
		if(op==0){
			int k;scanf("%d",&k);
			k=(k+la-1)%39989+1;
			printf("%d\n",la=T.Qmax(1,1,n,k).id);
		}else{
			int x0,y0,x1,y1;scanf("%d%d%d%d",&x0,&y0,&x1,&y1);
			x0=(x0+la-1)%39989+1;x1=(x1+la-1)%39989+1;
			y0=(y0+la-1)%int(1e9)+1;y1=(y1+la-1)%int(1e9)+1;
			if(x0>x1)swap(x0,x1),swap(y0,y1);
			T.insert(1,1,n,x0,x1,seg(x0,y0,x1,y1,++ind));
		}
	}
	return 0;
}