BZOJ 3165: [Heoi2013]Segment 标记永久化

3165: [Heoi2013]Segment

Time Limit: 40 Sec  Memory Limit: 256 MB
Submit: 581  Solved: 231
[Submit][Status][Discuss]

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。 

Sample Input

6
1 8 5 10 8
1 6 7 2 6
0 2
0 9
1 4 7 6 7
0 5

Sample Output

2
0

3

HINT

对于100%的数据，1 ≤ n ≤ 10^5 , 1 ≤  k, x0, x1 ≤ 39989, 1 ≤ y0 ≤ y1 ≤ 10^9。

标记永久化 （李超线段树）

和传送门基本一样

插入的是线段 刚开始竟然GG的看成直线

#include<ctime>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>

using namespace std;

typedef double db;
typedef long long ll;


inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=100100,mod1=39989,mod2=int(1e9),xlim=41000;
const ll inf=0X3f3f3f3f3f3f3f3fll;

double k[N],b[N];
int cnt,n,ans;
db mxy;

int tr[xlim<<2];

inline bool check(int x,int y,int p)
{return k[x]*p+b[x]>k[y]*p+b[y];}

inline db getres(int x,int p)
{return k[x]*p+b[x];}

void modify(int now,int l,int r,int x,int x1,int x2)
{
	if(l==r){if(l>=x1&&l<=x2)if(check(x,tr[now],l))tr[now]=x;return ;}
	int mid=l+r>>1;
	if(l>=x1&&r<=x2)
	{
		if(k[x]>=k[tr[now]])
			if(check(x,tr[now],mid))
				modify(now<<1,l,mid,tr[now],x1,x2),tr[now]=x;
			else modify(now<<1|1,mid+1,r,x,x1,x2);
		else
			if(check(x,tr[now],mid))
				modify(now<<1|1,mid+1,r,tr[now],x1,x2),tr[now]=x;
			else modify(now<<1,l,mid,x,x1,x2);
	}
	else if(mid<x1)
	modify(now<<1|1,mid+1,r,x,x1,x2);
	else if(mid>x2)
	modify(now<<1,l,mid,x,x1,x2);
	else modify(now<<1|1,mid+1,r,x,x1,x2),modify(now<<1,l,mid,x,x1,x2);
}

void query(int now,int l,int r,int p)
{
	int mid=l+r>>1;db tmp=getres(tr[now],p);
	if(mxy<tmp)mxy=tmp,ans=tr[now];
	else if(mxy==tmp)ans=min(ans,tr[now]);
	if(l==r)return ;
	if(mid<p)query(now<<1|1,mid+1,r,p);
	else query(now<<1,l,mid,p);
}

int main()
{
	n=read();
	register int i,opt,x1,y1,x2,y2;
	for(i=1;i<=n;++i)
	{
		opt=read();
		if(opt)
		{
			x1=(read()+ans-1)%mod1+1;y1=(read()+ans-1)%mod2+1;
			x2=(read()+ans-1)%mod1+1;y2=(read()+ans-1)%mod2+1;
			if(x1<x2)swap(x1,x2),swap(y1,y2);
			k[++cnt]=1.0*(y1-y2)/(x1-x2);
			b[cnt]=1.0*y1-k[cnt]*x1;
			modify(1,1,xlim,cnt,x2,x1);
		}
		else 
		{
			mxy=-inf;
			x1=(read()+ans-1)%mod1+1;
			query(1,1,xlim,x1);if(mxy==-inf)ans=0;
			print(ans);puts("");
		}
	}
	return 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
*/