loj 6034 线段游戏

题目大意：

给出若干条线段，用 (x1,y2),(x2,y2) 表示其两端点坐标，现在要求支持两种操作：

• 0 x1 y1 x2 y2 表示加入一条新的线段 (x1,y2),(x2,y2)
• 1 x0 询问所有线段中，x坐标在 x0 ​​ 处的最高点的 y 坐标是什么，如果对应位置没有线段，则输出

思路：

使用李超树维护

每一个加入的线段 在它能覆盖的线段树区间中修改

修改的时候判断它是否在中点的位置更优

然后将这个区间原来维护的线段与新线段比较

不优的那个考虑下放到左右儿子区间中继续更新

查询的时候答案为经过的所有线段树的线段中最大的

 #include<iostream>
 #include<cstdio>
 #include<cmath>
 #include<cstdlib>
 #include<cstring>
 #include<algorithm>
 #include<vector>
 #include<queue>
 #include<map>
 #define inf 2139062143
 #define ll long long
 #define MAXN 100100
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(!isdigit(ch)) {if(ch=='-') f=-;ch=getchar();}
     while(isdigit(ch)) {x=x*+ch-'';ch=getchar();}
     return x*f;
 }
 int n,m,mx[MAXN<<];
 struct seg{double k,b;}tag[MAXN<<];
 inline void build(int k,int l,int r)
 {
     tag[k].k=,tag[k].b=-inf;
     if(l==r) return ;
     int mid=(l+r)>>;
     build(k<<,l,mid);build(k<<|,mid+,r);
 }
 inline void real_mdf(int k,int l,int r,seg x)
 {
     //cout<<k<<" "<<l<<" "<<r<<" "<<x.k<<" "<<x.b<<endl;
     if(l==r) {if(tag[k].k*l+tag[k].b<x.k*l+x.b) tag[k]=x;return ;}
     int mid=(l+r)>>;
     if(tag[k].k*mid+tag[k].b<x.k*mid+x.b) swap(tag[k],x);
     //cout<<"i: "<<k<<" "<<l<<" "<<r<<" "<<x.k<<" "<<x.b<<endl;
     if(tag[k].k*l+tag[k].b<x.k*l+x.b) real_mdf(k<<,l,mid,x);
     if(tag[k].k*r+tag[k].b<x.k*r+x.b) real_mdf(k<<|,mid+,r,x);
 }
 inline void mdf(int k,int l,int r,int a,int b,seg x)
 {
     if(a>b) return ;
     //cout<<k<<" "<<l<<" "<<r<<" "<<a<<" "<<b<<" "<<x.k<<" "<<x.b<<endl;
     if(l==a&&r==b) {real_mdf(k,l,r,x);return ;}
     int mid=(l+r)>>;
     if(b<=mid) mdf(k<<,l,mid,a,b,x);
     else if(a>mid) mdf(k<<|,mid+,r,a,b,x);
     else {mdf(k<<,l,mid,a,mid,x);mdf(k<<|,mid+,r,mid+,b,x);}
 }
 inline double query(int k,int l,int r,int x)
 {
     //cout<<x<<" "<<tag[k].k*x+tag[k].b<<endl;
     if(l==r) return tag[k].k*x+tag[k].b;
     int mid=(l+r)>>;double res;
     //cout<<k<<" "<<res<<endl;
     if(x<=mid) res=query(k<<,l,mid,x);
     else res=query(k<<|,mid+,r,x);
     return max(tag[k].k*x+tag[k].b,res);
 }
 int main()
 {
     //freopen("C2.in","r",stdin);
     n=read(),m=read();int a,b,x,y,z;double k,B;
     build(,,MAXN-);
     for(int i=;i<=n;i++)
     {
         x=read(),y=read(),a=read(),b=read();
         //cout<<6<<endl;
         if(x==a) k=,B=max(y,b);
         else k=(double)(y-b)/(x-a),B=y-k*x;
         //cout<<k<<" "<<b<<endl;
         mdf(,,MAXN-,max(,min(a,x)),min(MAXN-,max(x,a)),(seg){k,B});
     }
     while(m--)
     {
         z=read();
         if(z) {a=read(),k=query(,,MAXN-,a);printf("%.3lf\n",k<=-inf?:k);}
         else
         {
             x=read(),y=read(),a=read(),b=read();
             //cout<<6<<endl;
             if(x==a) k=,B=max(y,b);
             else k=(double)(y-b)/(x-a),B=y-k*x;
             //cout<<k<<" "<<b<<endl;
             mdf(,,MAXN-,max(,min(a,x)),min(MAXN-,max(x,a)),(seg){k,B});
         }
         //system("pause");
     }
 }