[线段树] BZOJ 3165: [Heoi2013]Segment

Description
要求在平面直角坐标系下维护两个操作：
1.在平面上加入一条线段。记第i条被插入的线段的标号为i。
2.给定一个数k,询问与直线 x = k相交的线段中，交点最靠上的线段的编号。
Input
第一行一个整数n，表示共n 个操作。
接下来n行，每行第一个数为0或1。
若该数为 0，则后面跟着一个正整数 k，表示询问与直线
x = ((k +lastans–1)%39989+1)相交的线段中交点（包括在端点相交的情形）最靠上的线段的编号，其中%表示取余。若某条线段为直线的一部分，则视作直线与线段交于该线段y坐标最大处。若有多条线段符合要求，输出编号最小的线段的编号。
若该数为 1，则后面跟着四个正整数 x0, y0, x 1, y 1，表示插入一条两个端点为
((x0+lastans-1)%39989+1,(y0+lastans-1)%10^9+1)和((x
1+lastans-1)%39989+1,(y1+lastans-1)%10^9+1) 的线段。
其中lastans为上一次询问的答案。初始时lastans=0。
Output
对于每个 0操作，输出一行，包含一个正整数，表示交点最靠上的线段的编号。若不存在与直线相交的线段，答案为0。

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
#define ll long long
#define sl (s<<1)
#define sr (s<<1|1)
inline char tc(void){
static char fl[10000],*A=fl,*B=fl;
return A==B&&(B=(A=fl)+fread(fl,1,10000,stdin),A==B)?EOF:*A++;
}
inline int read(void){
int a=0;char c;
while((c=tc())<'0'||c>'9');
while(c>='0'&&c<='9')a=a*10+c-'0',c=tc();
return a;
}
int n,cnt,lastans,t[131072],x_1[100001],y_1[100001],x_2[100001],y_2[100005],s,ans;
long double o[100002];
inline int pd(const int &a,const int &b){
if(o[a]==o[b])return -1;
return o[a]>o[b];
}
inline int mx_1(const int &a,const int &b,const int &w){
if(fabs(y_1[a]+o[a]*(w-x_1[a])-y_1[b]-o[b]*(w-x_1[b]))<1e-9)return -1;
return y_1[a]+o[a]*(w-x_1[a])>y_1[b]+o[b]*(w-x_1[b]);
}
void change(const int &s,const int &l,const int &r,const int &x){
if(l>x_2[x]||r<x_1[x])return ;
if(l==r){
if(mx_1(x,t[s],l)==1)t[s]=x;
return ;
}
if(l>=x_1[x]&&r<=x_2[x]){
int mid=l+r>>1;
int w1=pd(x,t[s]),w2=mx_1(x,t[s],mid);
if(w1==1){
if(w2==1)change(sl,l,mid,t[s]),t[s]=x;
else change(sr,mid+1,r,x);
}else if(w1==0){
if(w2==1)change(sr,mid+1,r,t[s]),t[s]=x;
else change(sl,l,mid,x);
}
}else{
int mid=l+r>>1;
change(sl,l,mid,x),
change(sr,mid+1,r,x);
}
return ;
}
void query(const int &s,const int &l,const int &r,const int &x){
if(r<x||l>x)return ;
if(x<=r&&x>=l){
int w=mx_1(t[s],ans,x);
if(w==-1)
ans=min(ans,t[s]);
else if(w)ans=t[s];
}
if(l==r)return ;
int mid=l+r>>1;
query(sl,l,mid,x),query(sr,mid+1,r,x);
return ;
}
int main(void){
register int i,q;
n=read();
while(n--){
s=read();
if(s==1){
x_1[++cnt]=(read()+lastans-1)%39989+1,y_1[cnt]=(read()+lastans-1)%1000000000+1,
x_2[cnt]=(read()+lastans-1)%39989+1,y_2[cnt]=(read()+lastans-1)%1000000000+1,
x_1[cnt]>x_2[cnt]?swap(x_1[cnt],x_2[cnt]),swap(y_1[cnt],y_2[cnt]),1:0;
if(x_2[cnt]!=x_1[cnt])o[cnt]=(y_2[cnt]-y_1[cnt])/(long double)(x_2[cnt]-x_1[cnt]);
else y_1[cnt]=y_2[cnt]=max(y_1[cnt],y_2[cnt]);
change(1,1,40000,cnt);
}else{
q=(read()+lastans-1)%39989+1,ans=0,query(1,1,40000,q),printf("%d\n",lastans=ans);
}
}
return 0;
}