BZOJ3165 : [Heoi2013]Segment

建立线段树，每个节点维护该区间内的最优线段。

插入线段时，在线段树上分裂成$O(\log n)$棵子树，若与当前点的最优线段不相交，那么取较优的，否则暴力递归子树。

查询时在叶子到根路径上所有点的最优线段中取个最优的即可。

时间复杂度$O(n\log^2n)$。

#include<cstdio>
#include<cmath>
#include<algorithm>
#define N 39989
using namespace std;
struct Seg{
double k,b;
Seg(){}
Seg(int x0,int y0,int x1,int y1){if(x0==x1)k=0,b=max(y0,y1);else k=1.0*(y0-y1)/(x0-x1),b=-k*x0+y0;}
double gety(int x){return k*x+b;}
}s[100010];
int m,op,cnt,X0,Y0,X1,Y1,ans,v[131000];
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
inline int sig(double x){return fabs(x)<1e-8?0:(x>0?1:-1);}
void ins(int x,int a,int b,int c,int d,int p){
if(c<=a&&b<=d){
if(sig(s[p].gety(a)-s[v[x]].gety(a))>0&&sig(s[p].gety(b)-s[v[x]].gety(b))>0){v[x]=p;return;}
if(sig(s[p].gety(a)-s[v[x]].gety(a))<=0&&sig(s[p].gety(b)-s[v[x]].gety(b))<=0)return;
if(a==b)return;
}
int mid=(a+b)>>1;
if(c<=mid)ins(x<<1,a,mid,c,d,p);
if(d>mid)ins(x<<1|1,mid+1,b,c,d,p);
}
void ask(int x,int a,int b,int c){
if(sig(s[ans].gety(c)-s[v[x]].gety(c))<0)ans=v[x];
else if(!sig(s[ans].gety(c)-s[v[x]].gety(c))&&ans>v[x])ans=v[x];
if(a==b)return;
int mid=(a+b)>>1;
c<=mid?ask(x<<1,a,mid,c):ask(x<<1|1,mid+1,b,c);
}
int main(){
s[0].b=-1;
read(m);
while(m--){
read(op);
if(!op){
read(X0),X0=(X0+ans-1)%39989+1;
ans=0,ask(1,1,N,X0);
printf("%d\n",ans);
}else{
read(X0),read(Y0),read(X1),read(Y1);
X0=(X0+ans-1)%39989+1,Y0=(Y0+ans-1)%1000000000+1;
X1=(X1+ans-1)%39989+1,Y1=(Y1+ans-1)%1000000000+1;
s[++cnt]=Seg(X0,Y0,X1,Y1);
if(X0>X1)swap(X0,X1);
ins(1,1,N,X0,X1,cnt);
}
}
return 0;
}