洛谷P3122 [USACO15FEB]圈住牛Fencing the Herd(计算几何+CDQ分治)

题面

传送门

题解

题目转化一下就是所有点都在直线(Ax+By-C=0)的同一侧,也就可以看做所有点代入(Ax+By-C)之后的值符号相同,我们只要维护每一个点代入直线之后的最大值和最小值,看看每条直线的最大最小值符号是否相同就好了

以最大值为例,我们强制(Bgeq 0),那么能令这条直线取到最大值的点一定在所有点的上凸壳上,我们把凸壳搞出来,然后把所有直线按斜率从小到大排序,一遍扫过去就可以更新答案了

然而对于每条直线把前面所有点做个凸包?那有个吉儿用我还不如暴力枚举来得快

我们考虑(CDQ)分治,每次把([l,mid])之间的点做一个凸包,然后用来更新([mid+1,r])之间的答案就行了

//minamoto
#include<bits/stdc++.h>
#define R register
#define ll long long
#define inf 1e18
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
ll read(){
    R ll res,f=1;R char ch;
    while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
    for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
    return res*f;
}
const int N=2e5+5;
struct node{
    int x,y;
    node(){}
    node(R int xx,R int yy):x(xx),y(yy){}
    inline node operator -(const node &b)const{return node(x-b.x,y-b.y);}
    inline ll operator *(const node &b)const{return 1ll*x*b.y-1ll*y*b.x;}
    inline bool operator <(const node &b)const{return x==b.x?y<b.y:x<b.x;}
}p[N],st[N];
struct query{
    ll a,b,c;node p;int id;
    query(){}
    query(R ll a,R ll b,R ll c,R node P,R int Id):a(a),b(b),c(c),p(P),id(Id){}
}c[N],las[N];
struct qwq{ll a,b,c;int op,id;}q[N];
inline bool cmp1(const query &x,const query &y){return x.p*y.p>=0;}
inline bool cmp2(const query &x,const query &y){return x.p*y.p<=0;}
int ans[N];ll mn[N],mx[N],vc[N],maxn,minn;int op[N],n,m,tot,top,cnt,t;
inline void upd(R int id,R int x,R int y){
//	printf("%d %d %d
",id,x,y);
    R ll val=las[id].a*x+las[id].b*y-las[id].c;
    cmax(mx[id],val),cmin(mn[id],val);
}
void solve(int l,int r){
    if(l==r)return;
    int mid=(l+r)>>1;
    solve(l,mid),solve(mid+1,r),top=cnt=tot=0;
    fp(i,l,mid)if(q[i].op&1)p[++tot]=node(q[i].a,q[i].b);
    fp(i,mid+1,r)if(q[i].op&2)c[++cnt]=query(q[i].a,q[i].b,q[i].c,node(q[i].b,-q[i].a),q[i].id);
    if(!tot||!cnt)return;
    sort(p+1,p+1+tot),st[top=1]=p[1];
    fp(i,2,tot){
        while(top>1&&(p[i]-st[top-1])*(st[top]-st[top-1])>=0)--top;
        st[++top]=p[i];
    }
    sort(c+1,c+1+cnt,cmp1);
//    printf("%d %d %d
",l,r,cnt);
//    fp(i,1,cnt)printf("%d
",c[i].id);
    for(R int i=1,j=1;i<=top;++i){
        while(j<=cnt&&(i+1>top||(c[j].p*(st[i+1]-st[i]))>=0))
            upd(c[j].id,st[i].x,st[i].y),++j;
    }
    st[top=1]=p[1];
    fp(i,2,tot){
        while(top>1&&(p[i]-st[top-1])*(st[top]-st[top-1])<=0)--top;
        st[++top]=p[i];
    }
    sort(c+1,c+1+cnt,cmp2);
    for(R int i=1,j=1;i<=top;++i){
        while(j<=cnt&&(i+1>top||c[j].p*(st[i+1]-st[i])<=0))
            upd(c[j].id,st[i].x,st[i].y),++j;
    }
}
int main(){
//	freopen("testdata.in","r",stdin);
    n=read(),m=read(),maxn=-inf,minn=inf;
    fp(i,1,n)q[i].op=1,q[i].a=read(),q[i].b=read(),cmax(maxn,q[i].a),cmin(minn,q[i].a);
    t=n;
    fp(i,1,m)ans[i]=-1;
    fp(i,1,m){
        ++t,op[i]=q[t].op=read();
        if(q[t].op&1)q[t].a=read(),q[t].b=read(),cmax(maxn,q[t].a),cmin(minn,q[t].a);
        else{
            q[t].a=read(),q[t].b=read(),q[t].c=read();
            if(!q[t].b)ans[i]=(-1.0*q[t].c/q[t].a<minn||-1.0*q[t].c/q[t].a>maxn),--t;
            else{
                if(q[t].b<0)q[t].a=-q[t].a,q[t].b=-q[t].b,q[t].c=-q[t].c;
                mn[i]=inf,mx[i]=-inf,q[t].id=i;
                las[i].a=q[t].a,las[i].b=q[t].b,las[i].c=q[t].c;
            }
        }
    }
//    fp(i,1,t)printf("%d %d %lld %lld %lld
",q[i].op,q[i].id,q[i].a,q[i].b,q[i].c);
    solve(1,t);
    fp(i,1,t)if(q[i].op&2)ans[q[i].id]=(mx[q[i].id]<0||mn[q[i].id]>0);
    fp(i,1,m)if(op[i]&2)puts(ans[i]?"YES":"NO");
    return 0;
}
原文地址:https://www.cnblogs.com/bztMinamoto/p/10520710.html