题意：

实现以下两个操作：
(1)(1)(1)在平面上加入一条线段。记第iii条被插入的线段的标号为iii
(2)(2)(2)给定一个数kkk，询问与直线x=kx=kx=k相交的线段中，交点纵坐标最大的线段的编号。
n≤1e5,k,x0,x1≤39989,1≤y0,y1≤1e9n\le 1e5,k,x_0,x_1\le39989,1\le y_0,y_1 \le 1e9n≤1e5,k,x0,x1≤39989,1≤y0,y1≤1e9。

思路：

复习一下李超线段树。
李超树是用来维护一次函数，将线段划分为若干个段，即线段树的每一个段都维护一个最优线，他在当前段的midmidmid处的位置是最大值，我注意李超树线段不能向上pushuppushuppushup。当然我们如果要查询区间最小值或者最大值怎么办呢？我们总不能从[x1,x2][x_1,x_2][x1,x2]跑一遍，每个点都查询一遍吧？这样显然是不科学的，我们可以发现一次函数是个单调函数，我们在插入的时候顺便维护以下最小值就好了，最小值一定在当前区间的端点处。

下面我们考虑这个题，这个题不是维护单点最大值是多少，而是问最大值的线段编号，我们只需要再加一个变量维护以下编号即可，注意修改一下查询函数以及判断当前区间是否存在最优段。

//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
struct Node {int l,r,flag,id; double k,b;double calc(const int pos) const {return k*pos+b;}int cross(const Node &x) const {return floor((b-x.b)/(x.k-k));}
}tr[N<<2];
struct Query {int id;double val;
};void build(int u,int l,int r) {tr[u]={l,r,0,0,0,0};if(l==r) return;int mid=(l+r)>>1;build(L,l,mid); build(R,mid+1,r);
}void modify(int u,int l,int r,Node k) {if(k.l<=l&&k.r>=r) {if(!tr[u].flag) { tr[u]=k; }else if(k.calc(l)-tr[u].calc(l)>eps&&k.calc(r)-tr[u].calc(r)>eps) tr[u]=k;   else if(k.calc(l)-tr[u].calc(l)>eps||k.calc(r)-tr[u].calc(r)>eps) {int mid=(l+r)>>1;if(k.calc(mid)-tr[u].calc(mid)>eps) swap(tr[u],k);else if(fabs(k.calc(mid)-tr[u].calc(mid))<eps) {if(tr[u].id>k.id) swap(tr[u],k);}if(k.cross(tr[u])-mid<-eps) modify(L,l,mid,k);else modify(R,mid+1,r,k);}}else {int mid=(l+r)>>1;if(k.l<=mid) modify(L,l,mid,k);if(k.r>mid) modify(R,mid+1,r,k);}
}Query query(int u,int l,int r,int x,int flag) {//if(!tr[u].flag) return {0,0}; if(l==r) return {tr[u].id,tr[u].calc(x)};int mid=(l+r)>>1;double ans=tr[u].calc(x);int id=tr[u].id;if(x<=mid) {Query cmp=query(L,l,mid,x,flag);if(cmp.id==0) return {id,ans};if(ans>cmp.val) return {id,ans};else if(ans==cmp.val) return {min(id,cmp.id),ans};else return cmp;  }else {Query cmp=query(R,mid+1,r,x,flag);if(cmp.id==0) return {id,ans};if(ans>cmp.val) return {id,ans};else if(ans==cmp.val) return {min(id,cmp.id),ans};else return cmp;  }
}/*
double query(int u,int l,int r,int x) {//if(l==1&&r==1) printf("%.2f %.2f\n",tr[u].k,tr[u].b);if(l==r) return tr[u].calc(x);int mid=(l+r)>>1;double ans=tr[u].calc(x);if(x<=mid) return max(ans,query(L,l,mid,x));else return max(ans,query(R,mid+1,r,x));
}
*/int main() {//  ios::sync_with_stdio(false);
//  cin.tie(0);scanf("%d",&n);build(1,1,50000);int ans=0,id=0;while(n--) {int op; scanf("%d",&op);if(op==0) {int x; scanf("%d",&x);x=(x+ans-1)%39989+1;Query now=query(1,1,50000,x,1);printf("%d\n",ans=now.id);}else {id++;int x1,x2,y1,y2;scanf("%d%d%d%d",&x1,&y1,&x2,&y2);x1=(x1+ans-1)%39989+1; x2=(x2+ans-1)%39989+1;y1=(y1+ans-1)%mod+1;   y2=(y2+ans-1)%mod+1;if(x1>x2) swap(x1,x2),swap(y1,y2);if(x1==x2) modify(1,1,50000,{x1,x2,1,id,0,(double)max(y1,y2)});else {double k=double(y1-y2)/(x1-x2);double b=y1-k*x1;modify(1,1,50000,{x1,x2,1,id,k,b});}}}return 0;
}
/**/

P4097 [HEOI2013]Segment 李超线段树相关推荐

P4097 [HEOI2013]Segment 李超线段树
$ \color{#0066ff}{ 题目描述 }$ 要求在平面直角坐标系下维护两个操作: 在平面上加入一条线段.记第 i 条被插入的线段的标号为 i 给定一个数 k,询问与直线 x = k 相交的线 ...
【BZOJ 3165】 [Heoi2013]Segment 李超线段树
所谓李超线段树就是解决此题一类的问题(线段覆盖查询点最大(小)),把原本计算几何的题目变成了简单的线段树,巧妙地结合了线段树的标记永久化与标记下传,在不考虑精度误差的影响下,打法应该是这样的. #in ...
BZOJ.3165.[HEOI2013]Segment(李超线段树)
BZOJ 洛谷对于线段,依旧是存斜率即可. 表示精度误差一点都不需要管啊/托腮就我一个人看成了mod(10^9+1)吗.. //4248kb 892ms #include <cstdio&g ...
李超线段树 [Heoi2013]Segment
问题 D: [Heoi2013]Segment 时间限制: 4 Sec 内存限制: 256 MB 题目描述要求在平面直角坐标系下维护两个操作: 1.在平面上加入一条线段.记第i条被插入的线段的标号为 ...
【李超线段树】BZOJ3165 [Heoi2013]Segment
题面在这里李超线段树的裸题,不解释示例程序: #include<cstdio> #include<cmath> #include<algorithm> usin ...
【BZOJ3165】Segment（李超线段树）
题目来源:BZOJ3165 考虑以横坐标为下标维护线段树. 在每个结点维护一个标记,表示覆盖整个结点的最高线段,注意这个标记对整个区间都有作用,无需下传. 因为只有单点询问,所以可以不用维护区间的最高 ...
李超线段树(Li-Chao Segment Tree)
李超线段树李超线段树是一种用于维护平面直角坐标系内线段关系的数据结构.它常被用来处理这样一种形式的问题:给定一个平面直角坐标系,支持动态插入一条线段,询问从某一个位置 (X,+∞)(X,+\inft ...
BZOJ 1568 李超线段树
思路: 李超线段树裸题 //By SiriusRen #include <cmath> #include <cstdio> #include <cstring> # ...
【牛客NOIP模拟】牛牛的RPG游戏【二维偏序】【任意坐标斜率优化】【CDQ 分治】【李超线段树】
题意: n×mn\times mn×m 的网格图,每个点有两个权值 vali,j,bufi,jval_{i,j},buf_{i,j}vali,j,bufi,j,从 (1,1)(1,1)(1,1) ...

P4097 [HEOI2013]Segment 李超线段树

文章目录

题意：

思路：

P4097 [HEOI2013]Segment 李超线段树相关推荐

最新文章

热门文章