P4097 [HEOI2013]Segment 李超线段树

传送门

文章目录

  • 题意:
  • 思路:

题意:

实现以下两个操作:
(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 1e9n1e5,k,x0,x139989,1y0,y11e9

思路:

复习一下李超线段树。
李超树是用来维护一次函数,将线段划分为若干个段,即线段树的每一个段都维护一个最优线,他在当前段的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;
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315238.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【CF1204D】Kirk and a Binary String【结论题】【LIS】

传送门 题意&#xff1a;给一个01串SSS&#xff0c;求一个等长的01串TTT SSS和TTT所有对应位置的子串最长不下降子序列长度&#xff08;以下简称LIS\text{LIS}LIS&#xff09;相同TTT中0的数量尽量多 ∣S∣≤100000|S| \leq 100000∣S∣≤100000 对于一个01串SSS&#xff0c;…

这周,全球首个IT技术全中文免费学习平台诞生了!太惊艳!

本周三的Developer Wednesday为大家揭秘了Microsoft Learn网站的全新学习模式一时间Microsoft Learn圈粉无数如果你正在为错过直播而暴风哭泣等一等&#xff0c;看这里&#xff01;全程高清回播将为你持续解密这一次你一定不可以再错过了哦&#xff01;什么是Microsoft Learn&…

【CF1182D】Complete Mirror【树的重心】

传送门 题意&#xff1a;给一棵NNN个结点的树&#xff0c;你需要钦定一个根&#xff0c;使得所有深度相同的点的度数相同。 N≤100000N \leq 100000N≤100000 用脑子想一想&#xff0c;就是根节点直接相连的子树都长得一模一样。 如果根节点度数大于1&#xff0c;我们发现它…

P1989 无向图三元环计数 思维 + 建图

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 统计无向图中三元环的个数。 思路&#xff1a; 很明显有一种暴力的方法&#xff0c;就是枚举每条边&#xff0c;让后再跑两个点的所有边&#xff0c;可以卡到复杂度O(m2)O(m^2)O(m2)。 我们可以考虑给点之…

十分钟了解Kubernetes

何为Kubernetes?最简单的一句话来概括Kubernetes。它就是一套成熟的商用服务编排解决方案。Kubernetes定位在Paas层&#xff0c;重点解决了微服务大规模部署时的服务编排问题。Kubernetes组件介绍了解Kubernetes都是从Pod开始的。Pod是Kubernetes最小的调度单元&#xff0c;所…

【CF1020C】Election【贪心】

传送门 题意&#xff1a;NNN个人给MMM个党派投票&#xff0c;开始时每个人选择一个党派投票&#xff08;给定&#xff09;&#xff0c;每个人可以花费一定代价贿赂&#xff08;每个人给定代价&#xff09;&#xff0c;求让1号党派胜出&#xff08;票数严格最大&#xff09;的最…

Codeforces Round #627 (Div. 3) E. Sleeping Schedule dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一天hhh小时&#xff0c;初始时间是000&#xff0c;每天可以使时间aia_iai​或者ai−1a_i-1ai​−1&#xff0c;问最多可以让多少天的时间在[l,r][l,r][l,r]范围内。 思路&#xff1a; 算是个比较简单…

【CF1194E】Count The Rectangles【类扫描线】【单调性】【树状数组】

传送门 题意&#xff1a;给定NNN条与坐标轴平行的线段&#xff0c;保证不垂直的线段没有交点&#xff0c;求一共构成多少个矩形&#xff08;以线段交点为顶点&#xff09;。 1≤N≤50001\leq N\leq50001≤N≤5000 显然是个数据结构乱搞题。 直觉告诉我们先枚举一条线段。 假…

.NET Core 的过去、现在和未来

在最新的微软 Build 大会主题网站上&#xff0c;微软.NET 程序经理 Scott Hunter 发表了一篇文章&#xff0c;指出.NET Core 是.NET 的未来。去年&#xff0c;微软就已经确认 Visual Studio Live 中的.NET 框架将会被.NET Core 取代。将.NET Core 作为开源开发技术栈在 2014 年…

Ozon Tech Challenge 2020 (Div.1 + Div.2) C. Kuroni and Impossible Calcul 抽屉原理

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个数组ana_nan​&#xff0c;求∏1≤i<j≤n∣aj−ai∣modm\begin{matrix} \prod_{1\le i<j\le n} |a_j-a_i| \end{matrix}\bmod m∏1≤i<j≤n​∣aj​−ai​∣​modm。 n≤2e5,m≤1e3n\le2e5…

【HAOI2015】按位或【Min-Max容斥】【FWT】

传送门 题意&#xff1a;开始时你有一个数000,每次选出[0,2n−1][0,2^n-1][0,2n−1]中的一个数进行按位或&#xff0c;每个数选中的概率给定。求得到2n−12^n-12n−1的期望操作次数。 1≤n≤201\leq n\leq 201≤n≤20 神仙题 首先发现每一位都是独立的&#xff0c;可以分开考…

Visual Basic 兴衰记

这是一篇暴露年龄的文章。为什么这么说&#xff1f;因为现在年轻的程序员可能没有接触过当年红极一时的 Visual Basic。28 年前的 1991 年 4 月&#xff0c;Microsoft 发布了 Visual Basic 1.0 for Windows&#xff0c;次年 9 月&#xff0c;发布了 Visual Basic 1.0 for DOS。…

Ozon Tech Challenge 2020 (Div.1 + Div.2, Rated) D. Kuroni and the Celebration 交互 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一颗树&#xff0c;每次可以询问两个点的lcalcalca&#xff0c;询问次数不能超过⌊n2⌋\left \lfloor \frac{n}{2} \right \rfloor⌊2n​⌋&#xff0c;求它的根。 n≤1e3n\le1e3n≤1e3。 思路&#xf…

【CF1047D】Little C Loves 3 II【构造】【赛瓦维斯特定理】

传送门 题意&#xff1a;给一个NMN \times MNM的空棋盘&#xff0c;每次选取两个曼哈顿距离为3的空格子放上棋子&#xff0c;问最多能放多少个。 1≤N,M≤1e91 \leq N,M \leq 1e91≤N,M≤1e9 暴力讨论 假装N≤MN \leq MN≤M ①N1N1N1 容易得到&#xff0c;详见代码 ②N2N…

函数式编程里的Materialization应该翻译成什么?

Materialization是函数式编程里的一个专业术语, 用于特指函数式编程中查询被实际执行并生成结果的这一过程.首先, 搜了一下中文资料, 暂时没有对该词的中文翻译, CSDN\博客园\阿里云上所有关于它的文档都没有做中文翻译, 直接沿用这个英文单词.难道轮到我来创造这个词的中文翻…

Ozon Tech Challenge 2020 (Div.1 + Div.2) E.Kuroni and the Score Distribution 构造

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 不难想到&#xff0c;长度为nnn的数组最多的满足条件的三元组序列是1,2,3....,n1,2,3....,n1,2,3....,n&#xff0c;对于每一个位置贡献为i−12\frac{i-1}{2}2i−1​&#xff0c;那么如果m…

【CF1209E】Rotate Columns【状压DP】【位运算】【贪心】

题意&#xff1a;给一个NMN \times MNM的矩阵&#xff0c;可以进行任意多次操作将一列轮换&#xff0c;求每一行的最大值之和的最大值。多组数据。 Easy VersionN≤4N \leq 4N≤4,M≤100M \leq100M≤100 Hard VersionN≤12N \leq 12N≤12,M≤2000M \leq2000M≤2000 看这数据…

Exceptionless(二) - 使用进阶

作者&#xff1a;markjiang7m2原文地址&#xff1a;https://www.cnblogs.com/markjiang7m2/p/11100563.html官网地址&#xff1a;http://letyouknow.net在上一篇文章Exceptionless - .Net Core开源日志框架中就说到如何对Exceptionless进行本地化部署&#xff0c;不过我也跟大家…

Ozon Tech Challenge 2020 (Div.1 + Div.2) F. Kuroni and the Punishment 随机化

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个数&#xff0c;每次操作可以选择将某个数1,−11,-11,−1&#xff0c;求最少进行多少次操作使得所有数都为正数且gcd>1gcd>1gcd>1。 思路&#xff1a; 考虑gcd2gcd2gcd2的情况&#xff0…

【洛谷P4841】城市规划【指数型生成函数】【麦克劳林级数】【多项式对数】

传送门 题意&#xff1a;求NNN个点的带标号无向连通简单图的个数。 N≤130000N \leq 130000N≤130000 这个问题的主要矛盾在于连通 这个并不好表示&#xff0c;但可以用这个表示出不要求连通的方案数 由于带标号&#xff0c;先构造答案的EGF f(x)∑i0∞aii!xif(x)\sum_{i0}…