[线段树] 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)和((x1+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;}