[线段树] 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;
}