查询

题目描述

给出若干条线段，用(x1,y1)，(x2,y2)表示其两端点坐标，现在要求支持两种操作：
0 x1 y1 x2 y2
表示加入一条新的线段，(x1,y1)-(x2,y2)
1 x0
询问所有线段中，x坐标在x0处的最高点的y坐标是什么，如果对应位置没有线段，则输出0。

分治

考虑CDQ分治，问题转化为先做0操作再做1操作。
按照横坐标建线段树，在一个线段树节点上线段可当直线用。
维护凸包即可。
为了保证复杂度，可以按直线斜率排序插入线段树，只涉及队尾操作。
询问按横坐标排序，也只涉及队尾操作。
常数略大，那两个排序是可以归并的。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef double db;
const int maxn=200000+10,maxtot=3400000+10,mx=100000,inf=10000000;
int head[maxn*4],tail[maxn*4],next[maxtot],pre[maxtot],id[maxtot],L[maxn],R[maxn];
bool bz[maxtot];
int qk[maxtot];
db xy[maxtot];
db kk[maxn],bb[maxn],ans[maxn];
int ask[maxn],c[maxn],q[maxn],tc[maxn],tq[maxn],d[maxn],dd[maxn];
int v[maxn];
int i,j,k,l,t,n,m,top,tot,cnt,num,x1,y1,x2,y2;
int read(){
    int x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9'){
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
void ins(int x1,int y1,int x2,int y2,int time){
    if (x2<1||x1>mx) return;
    top++;
    c[top]=top;
    tc[top]=time;
    if (x1==x2){
        y1=max(y1,y2);
        kk[top]=1;
        bb[top]=y1-x1;
        L[top]=R[top]=x1;
        return;
    }
    kk[top]=(db)(y1-y2)/(x1-x2);
    bb[top]=y1-x1*kk[top];
    L[top]=max(x1,1);R[top]=min(x2,mx);
}
db getxy(int a,int b){
    if (kk[a]==kk[b]){
        if (bb[a]<=bb[b]) return 0;
        else return mx;
    }
    return (bb[a]-bb[b])/(kk[b]-kk[a]);
}
void mark(int p,int x){
    if (!bz[p]){
        bz[p]=1;
        qk[++num]=p;
        head[p]=tail[p]=++tot;
        id[tot]=x;
    }
    else{
        while (head[p]!=tail[p]&&getxy(id[tail[p]],x)<=xy[tail[p]]){
            next[pre[tail[p]]]=0;
            tail[p]=pre[tail[p]];
        }
        id[++tot]=x;
        next[tail[p]]=tot;
        pre[tot]=tail[p];
        xy[tot]=getxy(id[tail[p]],x);
        tail[p]=tot;
    }
}
void change(int p,int l,int r,int a,int b,int x){
    if (l==a&&r==b){
        mark(p,x);
        return;
    }
    int mid=(l+r)/2;
    if (b<=mid) change(p*2,l,mid,a,b,x);
    else if (a>mid) change(p*2+1,mid+1,r,a,b,x);
    else{
        change(p*2,l,mid,a,mid,x);
        change(p*2+1,mid+1,r,mid+1,b,x);
    }
}
db work(int p,int x){
    if (!bz[p]) return -inf;
    while (head[p]!=tail[p]&&x<=xy[tail[p]]){
        next[pre[tail[p]]]=0;
        tail[p]=pre[tail[p]];
    }
    int k=id[tail[p]];
    return kk[k]*x+bb[k];
}
db query(int p,int l,int r,int a){
    db t=work(p,a);
    if (l==r) return t;
    int mid=(l+r)/2;
    if (a<=mid) t=max(t,query(p*2,l,mid,a));else t=max(t,query(p*2+1,mid+1,r,a));
    return t;
}
void solve(int l,int r,int u,int v,int x,int y){
    if (l==r) return;
    int mid=(l+r)/2,m1=u-1,m2=x-1;
    int i,j,k;
    fo(i,u,v)
        if (tc[i]<=mid) m1=i;
    fo(i,x,y)
        if (tq[i]<=mid) m2=i;
    solve(l,mid,u,m1,x,m2);
    solve(mid+1,r,m1+1,v,m2+1,y);
    j=u;k=m1+1;
    fo(i,u,v){
        if (j>m1){
            d[i]=c[k];
            dd[i]=tc[k];
            k++;
        }
        else if (k>v){
            d[i]=c[j];
            dd[i]=tc[j];
            j++;
        }
        else if (kk[c[j]]<kk[c[k]]){
            d[i]=c[j];
            dd[i]=tc[j];
            j++;
        }
        else{
            d[i]=c[k];
            dd[i]=tc[k];
            k++;
        }
    }
    fo(i,u,v) c[i]=d[i],tc[i]=dd[i];
    j=x;k=m2+1;
    fo(i,x,y){
        if (j>m2){
            d[i]=q[k];
            dd[i]=tq[k];
            k++;
        }
        else if (k>y){
            d[i]=q[j];
            dd[i]=tq[j];
            j++;
        }
        else if (ask[q[j]]>ask[q[k]]){
            d[i]=q[j];
            dd[i]=tq[j];
            j++;
        }
        else{
            d[i]=q[k];
            dd[i]=tq[k];
            k++;
        }
    }
    fo(i,x,y) q[i]=d[i],tq[i]=dd[i];
    tot=0;num=0;
    if (u<=v&&x<=y){
        fo(i,u,v)
            if (tc[i]<=mid) change(1,1,mx,L[c[i]],R[c[i]],c[i]);
        fo(i,x,y)
            if (tq[i]>mid) ans[q[i]]=max(ans[q[i]],query(1,1,mx,ask[q[i]]));
    }
    fo(i,1,num) bz[qk[i]]=0;
}
int main(){
    freopen("query.in","r",stdin);freopen("query.out","w",stdout);
    n=read();m=read();
    fo(i,1,n){
        x1=read();y1=read();x2=read();y2=read();
        if (x1>x2){
            swap(x1,x2);
            swap(y1,y2);
        }
        ins(x1,y1,x2,y2,i);
    }
    fo(i,1,m){
        t=read();
        if (t==0){
            x1=read();y1=read();x2=read();y2=read();
            if (x1>x2){
                swap(x1,x2);
                swap(y1,y2);
            }
            ins(x1,y1,x2,y2,i+n);
        }
        else{
            ask[++cnt]=read();
            tq[cnt]=i+n;
            q[cnt]=cnt;
            ans[cnt]=-inf;
        }
    }
    solve(1,n+m,1,top,1,cnt);
    fo(i,1,cnt) printf("%.6lf\n",ans[i]==-inf?0:ans[i]);
}