【JZOJ5039】【NOI2017模拟4.2】查询

Description

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

Data Constraint

对于10%的数据：n,m<=1000
另外20%的数据：所有的1操作都在0操作之后
另外20%的数据：所有线段的两端的x坐标都包含所有的询问的x坐标，你可以将每条线段当做直线处理
对于100%的数据：n<=50000,m<=150000，x1,x2,y1,y2,x0均为整数，0< x0<=1e5, -1e6<=x1,x2,y1,y2<=1e6。

Code

学了个新思路，原来线段树还可以这样用……
我们在线段树的每个节点存下一条线段，表示在该区间内该线段相对较优。下面讨论插入问题。
1、当你插入一条线段时，假如发现该线段在区间内线段的上方，直接覆盖就好。下方就退出。
2、当你插如的线段与区间内线段有交时，我们判断交点在mid的左边还是右边。假如另一边是插入线段较优，我们就把区间内的线段修改成插入线段，并把原有线段往相交的一边下传。否则就直接让插入线段往相交的一边下传。

Code

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=1e6+1;
struct code{
    long double k;
    long double b;
}f[maxn*8];
int n,m,i;
long double p,q,t,j,k,l,x,y,z,xx,yy,bz;
double ans;
long double pan(long double k,long double b,long double x){
    return k*x+b;
}
void insert(int l,int r,int v,int x,int y){
    int mid=(l+r)/2;long double t,k,z;
    if (l>=x && r<=y){
        if (bz==f[v].b) f[v].k=q,f[v].b=p;
        else if (l==r){
            t=pan(f[v].k,f[v].b,l)-pan(q,p,l);
            if (t<0) f[v].k=q,f[v].b=p;
        }else if (l!=r){
            t=pan(f[v].k,f[v].b,l)-pan(q,p,l);k=pan(f[v].k,f[v].b,r)-pan(q,p,r);
            if (t*k>=0){
                if (t<=0 && k<=0) f[v].k=q,f[v].b=p;
            }else{
                z=pan(f[v].k,f[v].b,mid)-pan(q,p,mid);
                if (t*z>0){
                    if (t<0) swap(f[v].k,q),swap(f[v].b,p);
                    insert(mid+1,r,v*2+1,x,y);
                    if (t<0) q=f[v].k,p=f[v].b;
                }
                else if (t*z<0){
                    if (z<0) swap(f[v].k,q),swap(f[v].b,p);
                    insert(l,mid,v*2,x,y);  
                    if (z<0) q=f[v].k,p=f[v].b;
                }
            }
        }
        return;
    }
    if (l<=y && mid>=x) insert(l,mid,v*2,x,y);
    if (mid<y && r>=x) insert(mid+1,r,v*2+1,x,y);
}
void find(int l,int r,int v,int x){
    int mid=(l+r)/2;
    if (pan(f[v].k,f[v].b,x)>pan(p,q,x)) p=f[v].k,q=f[v].b;
    if (l==r) return;
    if (mid>=x) find(l,mid,v*2,x);
    else find(mid+1,r,v*2+1,x);
}
int main(){
    freopen("query.in","r",stdin);freopen("query.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (i=1;i<maxn*8;i++)
        f[i].b=-1e10;
    bz=f[1].b;t=-1e10;
    for (i=1;i<=n;i++){
        scanf("%Lf%Lf%Lf%Lf",&x,&y,&xx,&yy);
        if (xx<x) swap(x,xx),swap(y,yy);
        if (x==xx) y=yy=max(y,yy);
        x+=maxn;xx+=maxn;
        if (xx!=x) q=(yy-y)/(xx-x);
        else q=0;
        p=y-q*x;
        insert(1,maxn*2,1,x,xx);
    }
    for (i=1;i<=m;i++){
        scanf("%Lf%Lf",&z,&x);
        if (z){
            p=q=-1e9;x+=maxn;
            find(1,maxn*2,1,x);
            p=pan(p,q,x);ans=p;
            if (p!=-1e10) printf("%.6lf\n",ans);
            else printf("0.0000000\n");
        }else{
            scanf("%Lf%Lf%Lf",&y,&xx,&yy);  
            if (xx<x) swap(x,xx),swap(y,yy);
            if (x==xx) y=yy=max(y,yy);
            x+=maxn;xx+=maxn;
            if (xx!=x) q=(yy-y)/(xx-x);
            else q=0;
            p=y-q*x;
            insert(1,maxn*2,1,x,xx);
        }
    }
}