题面

传送门

题解

调得咱自闭了……

不难发现这就是个李超线段树,不过因为这里加入的是线段而不是直线,所以得把线段在线段树上对应区间内拆开之后再执行李超线段树的操作,那么复杂度就是\(O(n\log^2n)\)

以上是题解,以下是吐(zang)槽(hua)

为什么我插入竖直线段的时候会出现\(l=r\)且\(l\times k+b\neq r\times k+b\)的情况呢……为什么我本地测会突然\(RE\)加了个\(puts("qwq")\)就能输出了呢……

//minamoto
#include<bits/stdc++.h>
#define R register
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
R int res,f=1;R char ch;
while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
char sr[1<<21],z[20];int C=-1,Z=0;
inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
void print(R int x){
if(C>1<<20)Ot();if(x<0)sr[++C]='-',x=-x;
while(z[++Z]=x%10+48,x/=10);
while(sr[++C]=z[Z],--Z);sr[++C]='\n';
}
const int M=39989,L=1e9,N=1e5+5;
struct node{
node *lc,*rc;double b,k;bool flag;int id;
inline void ins(R double bb,R double kk,R int ii){b=bb,k=kk,id=ii,flag=1;}
inline double calc(const R int &x){return k*x+b;}
}pool[N<<2],*rt;int tot;
int n,q,id,ii,ql,qr,res,x;double b,k,bb,kk,mx;
inline node *newnode(){return &pool[tot++];}
void build(node* &p,int l,int r){
p=newnode();if(l==r)return;
int mid=(l+r)>>1;
build(p->lc,l,mid),build(p->rc,mid+1,r);
}
void update(node *p,int l,int r,double b,double k,int id){
if(ql<=l&&qr>=r){
if(!p->flag)return p->ins(b,k,id),void();
double lv1=l*k+b,rv1=r*k+b,lv2=p->calc(l),rv2=p->calc(r);
if(lv1<=lv2&&rv1<=rv2)return;
if(lv1>lv2&&rv1>rv2)return p->ins(b,k,id),void();
int mid=(l+r)>>1;
double x=(b-p->b)/(p->k-k);
if(lv1<=lv2){
if(x<=mid)bb=p->b,kk=p->k,ii=p->id,p->ins(b,k,id),update(p->lc,l,mid,bb,kk,ii);
else update(p->rc,mid+1,r,b,k,id);
}else{
if(x<=mid)update(p->lc,l,mid,b,k,id);
else bb=p->b,kk=p->k,ii=p->id,p->ins(b,k,id),update(p->rc,mid+1,r,bb,kk,ii);
}
return;
}
int mid=(l+r)>>1;
if(ql<=mid)update(p->lc,l,mid,b,k,id);
if(qr>mid)update(p->rc,mid+1,r,b,k,id);
}
void query(node *p,int l,int r){
if(cmax(mx,p->calc(x)))res=p->id;
if(l==r)return;
int mid=(l+r)>>1;
x<=mid?query(p->lc,l,mid):query(p->rc,mid+1,r);
}
int cnt,op,val[N],ID[N];
int main(){
q=read(),n=39989;
build(rt,1,n);
while(q--){
op=read();
if(op==1){
int x0=(read()+res-1)%n+1,y0=(read()+res-1)%L+1;
int x1=(read()+res-1)%n+1,y1=(read()+res-1)%L+1;
++cnt;
if(x0==x1&&cmax(val[x0],max(y0,y1))){ID[x0]=cnt;continue;}
if(x0>x1)swap(x0,x1),swap(y0,y1);
k=1.0*(y1-y0)/(x1-x0),b=y0-k*x0,id=cnt,ql=x0,qr=x1;
update(rt,1,n,b,k,id);
}else{
x=(read()+res-1)%n+1,mx=res=0;
query(rt,1,n);
cmax(mx,1.0*val[x])?res=ID[x]:0;
print(res);
}
}
return Ot(),0;
}

洛谷P4097 [HEOI2013]Segment(李超线段树)的更多相关文章

  1. 【洛谷P4097】Segment 李超线段树

    题目大意:维护一个二维平面,给定若干条线段,支持询问任意整数横坐标处对应的纵坐标最靠上的线段的 id,相同高度取 id 值较小的,强制在线. 题解:初步学习了李超线段树.李超线段树的核心思想在于通过标 ...

  2. Luogu P4097 [HEOI2013]Segment 李超线段树

    题目链接 \(Click\) \(Here\) 李超线段树的模板.但是因为我实在太\(Naive\)了,想象不到实现方法. 看代码就能懂的东西,放在这里用于复习. #include <bits/ ...

  3. P4097 [HEOI2013]Segment 李超线段树

    $ \color{#0066ff}{ 题目描述 }$ 要求在平面直角坐标系下维护两个操作: 在平面上加入一条线段.记第 i 条被插入的线段的标号为 i 给定一个数 k,询问与直线 x = k 相交的线 ...

  4. 2018.07.23 洛谷P4097 [HEOI2013]Segment(李超线段树)

    传送门 给出一个二维平面,给出若干根线段,求出x" role="presentation" style="position: relative;"&g ...

  5. 洛谷 P4097 [HEOI2013]Segment 解题报告

    P4097 [HEOI2013]Segment 题目描述 要求在平面直角坐标系下维护两个操作: 在平面上加入一条线段.记第 \(i\) 条被插入的线段的标号为 \(i\) 给定一个数 \(k\),询问 ...

  6. 洛谷P4069 [SDOI2016]游戏(李超线段树)

    题意 题目链接 Sol 这题细节好多啊qwq..稍不留神写出一个小bug就要调1h+.. 思路就不多说了,把询问区间拆成两段就是李超线段树板子题了. 关于dis的问题可以直接维护. // luogu- ...

  7. BZOJ3165: [Heoi2013]Segment(李超线段树)

    题意 题目链接 Sol 李超线段树板子题.具体原理就不讲了. 一开始自己yy着写差点写自闭都快把叉积搬出来了... 后来看了下litble的写法才发现原来可以写的这么清晰简洁Orz #include& ...

  8. 【BZOJ 3165】 [Heoi2013]Segment 李超线段树

    所谓李超线段树就是解决此题一类的问题(线段覆盖查询点最大(小)),把原本计算几何的题目变成了简单的线段树,巧妙地结合了线段树的标记永久化与标记下传,在不考虑精度误差的影响下,打法应该是这样的. #in ...

  9. [洛谷P4097] [HEOI2013] Segment

    Description 要求在平面直角坐标系下维护两个操作: 1.在平面上加入一条线段.记第 \(i\) 条被插入的线段的标号为 \(i\) 2.给定一个数 \(k\) ,询问与直线 \(x = k\ ...

随机推荐

  1. Python实现七牛云视频播放

    这篇文章是使用Python的Web框架Django Rest Framework来提供视频相关的api接口,主要功能包括视频上传.视频转码.视频访问授权.删除视频文件.视频截图功能. 七牛云上的基本概 ...

  2. java算法 第七届 蓝桥杯B组(题+答案) 10.压缩变换

    10.压缩变换  (程序设计) 小明最近在研究压缩算法.他知道,压缩的时候如果能够使得数值很小,就能通过熵编码得到较高的压缩比.然而,要使数值很小是一个挑战. 最近,小明需要压缩一些正整数的序列,这些 ...

  3. 不要在遍历子结点时修改parent

    [不要在遍历子结点时修改parent] 在用for/foreach遍历子结点时,如果在这过程中有改变子结点的parent,会导致不可预料的结果.我所遇到的问题是,在此种情况下,并非所有的子结点都能遍历 ...

  4. Python基础:字符串的常见操作

    # 切片 # 切片 获取对象中一部分数据 [起始位置:结束位置(不包含):步长] qpstr = "山东张学友" result = qpstr[1: 3: 1] # 东张 prin ...

  5. 游戏动作师使用Unity3D遇到过的所有问题

    http://blog.csdn.net/onafioo/article/details/50865169 http://www.gameres.com/thread_480489.html 文/拉撒 ...

  6. js setTimeout 和 setInterval 区别

    setTimeout和setInterval的语法相同.它们都有两个参数,一个是将要执行的代码字符串,还有一个是以毫秒为单位的时间间隔,当过了那个时间段之后就将执行那段代码.不过这两个函数还是有区别的 ...

  7. Python代码规范利器Flake8

    写代码其实是需要规范的,团队中更是如此:不然 Google 也不会发布各种编码规范,耳熟能详的有Google C++ 风格指南,Google Python 风格指南,等等. 这些规范有用吗?有用也没用 ...

  8. itunes win10应用商店安装 升级固件路径

    win + r 在运行中输入: %localappdata%\Packages\AppleInc.iTunes_nzyj5cx40ttqa\LocalCache\Roaming\Apple Compu ...

  9. 3.1.5 倒计时器:CountDownLatch

    package 第三章.倒计时器CountDownLatch; import java.util.concurrent.CountDownLatch; /** * Created by zzq on ...

  10. MSSQL数据库高版本迁移到低版本

    起因是因为客户要把系统从阿里云迁移到本地服务器,阿里云上的数据库版本是MSSQL2016,客户提供的服务器是Server2008R2的,问题就来了,Server2008不支持2016版本,最后只能装的 ...