李超线段树 [Heoi2013]Segment

问题 D: [Heoi2013]Segment
时间限制: 4 Sec 内存限制: 256 MB
题目描述
要求在平面直角坐标系下维护两个操作：
1.在平面上加入一条线段。记第i条被插入的线段的标号为i。
2.给定一个数k,询问与直线 x = k相交的线段中，交点最靠上的线段的编号。

输入

第一行一个整数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。

输出
对于每个 0操作，输出一行，包含一个正整数，表示交点最靠上的线段的编号。若不存在与直线相交的线段，答案为0。

样例输入
6
1 8 5 10 8
1 6 7 2 6
0 2
0 9
1 4 7 6 7
0 5
样例输出
2
0 3
提示
对于100%的数据，1 ≤ n ≤ 10^5 , 1 ≤ k, x0, x1 ≤ 39989, 1 ≤ y0 ≤ y1 ≤ 10^9。

其实这道题就是个板子，李超线段树。（这个板子就是为了解决这个问题，只不过板子太难了而已，汗颜）

就题论算法把，李超线段树用来处理向一个区间加有斜率的线段，之后判断某位置权值最大的线段是哪条之类的问题。而最重要的是多了的insert函数。当分到的区间已属于这条线段覆盖的区间是，就要进行更神奇的操作了，去比较当前节点所存的最优线段。先附上代码。

void insert(int l,int r,int x,int k)
{if(!t[x].h)t[x].h=k;if(cmp(t[x].h,k,l))swap(t[x].h,k);if(l==r||a[t[x].h].k==a[k].k)return;int mid=(l+r)/2;double g=(double)(a[t[x].h].b-a[k].b)/(a[k].k-a[t[x].h].k);if(g<l||g>r)return;if(g<=mid)insert(l,mid,x*2,t[x].h),t[x].h=k;else insert(mid+1,r,x*2+1,k);
}

因为比较分为两种情况，完全压制和在区间内有交点。完全压制不动或直接修改后返回就好了。而对于相交的就必须求出交点，根据交点的位置（其实也就是判断那条线段在交点那一边处于优势）去修改子区间。
对于这段代码理解还是有点困难（本蒟蒻太水，自己口胡。。）要把本区间右侧较大的线段置成本区间最优（我也不是太懂为什么。。望神犇来解释），若不懂，手动模拟一下过程即可。
还要注意一些细节：直线斜率是否存在，以及两条线是否平行。。（RE到死）

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define N 101000
#define ll long long
using namespace std;
struct node
{int l,r,id;double k,b;
}a[N];
struct tree{int l,r,h;}t[50000*4];
int n,m;
void build(int l,int r,int x)
{t[x].l=l;t[x].r=r;if(l==r){t[x].h=0;return;}int mid=(l+r)/2;build(l,mid,x*2);build(mid+1,r,x*2+1);
}
bool cmp(int x,int y,double l)
{if(!x)return 1;double l1=a[x].k*l+a[x].b,l2=a[y].k*l+a[y].b;return l1!=l2?l1<l2:x<y;
}
int q(int k,int x)
{if(t[x].l==t[x].r)return t[x].h;int mid=(t[x].l+t[x].r)/2,tmp;if(k<=mid)tmp=q(k,x*2);else tmp=q(k,x*2+1);return cmp(t[x].h,tmp,k)?tmp:t[x].h;
}
void insert(int l,int r,int x,int k)
{if(!t[x].h)t[x].h=k;if(cmp(t[x].h,k,l))swap(t[x].h,k);if(l==r||a[t[x].h].k==a[k].k)return;int mid=(l+r)/2;double g=(double)(a[t[x].h].b-a[k].b)/(a[k].k-a[t[x].h].k);if(g<l||g>r)return;if(g<=mid)insert(l,mid,x*2,t[x].h),t[x].h=k;else insert(mid+1,r,x*2+1,k);
}
void c(int k,int x)
{if(t[x].l>=a[k].l&&t[x].r<=a[k].r){insert(t[x].l,t[x].r,x,k);return;}int mid=(t[x].l+t[x].r)/2;if(a[k].l<=mid)c(k,x*2);if(a[k].r>mid)c(k,x*2+1);
}
int main()
{scanf("%d",&n);int x1,x2,y1,y2,z,ans=0;build(1,50000,1);while(n--){scanf("%d",&z);if(z==0){scanf("%d",&x1);x1=(x1+ans-1)%39989+1;//out<<x1<<endl;ans=q(x1,1);printf("%d\n",ans);}else{m++;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)%1000000000+1;y2=(y2+ans-1)%1000000000+1;//cout<<x1<<" "<<x2<<" "<<y1<<" "<<y2<<endl;if(x1>x2)swap(x1,x2),swap(y1,y2);if(x1==x2)a[m].k=0,a[m].b=max(y1,y2);else{a[m].k=(double)(y1-y2)/(x1-x2);a[m].b=(double)y1-a[m].k*x1;}a[m].l=x1;a[m].r=x2;a[m].id=m;c(m,1);}}
}