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)%109+1)和((x
1+lastans-1)%39989+1,(y1+lastans-1)%109+1) 的线段。
其中lastans为上一次询问的答案。初始时lastans=0。
Output
对于每个 0操作,输出一行,包含一个正整数,表示交点最靠上的线段的编
号。若不存在与直线相交的线段,答案为0。
Sample Input
1 8 5 10 8
1 6 7 2 6
0 2
0 9
1 4 7 6 7
0 5
Sample Output
0 3
HINT
对于100%的数据,1 ≤ n ≤ 10^5 , 1 ≤ k, x0, x1 ≤ 39989, 1 ≤ y0 ≤ y1 ≤ 10^9。
线段树。每次修改的时候看与当前区间的最高线段有没有交点。有的话继续往下查找。没有的话判断能否替代
询问的时候自底向上,单次复杂度logn
#include<cstdio>
#include<algorithm>
using namespace std;
struct tree
{
int l,r;
double ll,rr;
int s;
}tr[200001];
int pl[40001];
inline void build(int p,int l,int r)
{
tr[p].l=l;
tr[p].r=r;
if(l!=r)
{
int mid=(l+r)/2;
build(p*2,l,mid);
build(p*2+1,mid+1,r);
}
else
pl[l]=p;
}
inline void change(int p,int x0,int y0,int x1,int y1,double k,int d)
{
if(x0<=tr[p].l&&tr[p].r<=x1)
{
double yy0=double(y0)+double(tr[p].l-x0)*k,yy1=yy0+double(tr[p].r-tr[p].l)*k;
if(yy0>=tr[p].ll&&yy1>=tr[p].rr)
{
tr[p].ll=yy0;
tr[p].rr=yy1;
tr[p].s=d;
return ;
}
if(yy0<=tr[p].ll&&yy1<=tr[p].rr)
return ;
if(tr[p].l==tr[p].r)
{
int height=max(tr[p].ll,tr[p].rr);
if(yy0>height||yy1>height)
{
tr[p].ll=yy0;
tr[p].rr=yy1;
tr[p].s=d;
}
}
else
{
change(p*2,x0,y0,x1,y1,k,d);
change(p*2+1,x0,y0,x1,y1,k,d);
}
}
else
{
int mid=(tr[p].l+tr[p].r)/2;
if(x0<=mid)
change(p*2,x0,y0,x1,y1,k,d);
if(x1>mid)
change(p*2+1,x0,y0,x1,y1,k,d);
}
}
struct ans
{
double h;
int x;
};
ans xt,xc;
inline void ask(int p,int x)
{
if(p==0)
return ;
if(tr[p].l==tr[p].r)
{
xt.x=tr[p].s;
xt.h=max(tr[p].ll,tr[p].rr);
ask(p/2,x);
}
double k=(tr[p].rr-tr[p].ll)/double(tr[p].r-tr[p].l);
double height=tr[p].ll+k*double(x-tr[p].l);
if(height>xt.h)
{
xt.h=height;
xt.x=tr[p].s;
}
ask(p/2,x);
}
int main()
{
// freopen("segment.in","r",stdin);
// freopen("segment.out","w",stdout);
int n;
scanf("%d",&n);
build(1,1,39989);
int i;
int x;
int x0,y0,x1,y1;
int lastans=0;
int d=0;
for(i=1;i<=n;i++)
{
scanf("%d",&x);
if(x==0)
{
scanf("%d",&x0);
x0=((x0+lastans-1)%39989+1);
xt=xc;
ask(pl[x0],x0);
printf("%d\n",xt.x);
lastans=xt.x;
}
else
{
d++;
scanf("%d%d%d%d",&x0,&y0,&x1,&y1);
x0=((x0+lastans-1)%39989+1);
y0=((y0+lastans-1)%1000000000+1);
x1=((x1+lastans-1)%39989+1);
y1=((y1+lastans-1)%1000000000+1);
if(x0>x1)
{
int t=x0;
x0=x1;
x1=t;
t=y0;
y0=y1;
y1=t;
}
double k=double(y1-y0)/double(x1-x0);
change(1,x0,y0,x1,y1,k,d);
}
}
return 0;
}