• BZOJ4049 : [Cerc2014] Mountainous landscape


    对于区间[l,r],若询问的直线与该区间的凸壳有交,则有解。

    在线段树每个区间上维护凸壳,然后查询时在相应凸壳上二分斜率即可,时间复杂度$O(nlog^2n)$。

    #include<cstdio>
    #define N 262145
    typedef long long ll;
    int T,n,i,c,l[N],r[N],t;
    struct P{
      int x,y;
      P(){}
      P(int _x,int _y){x=_x,y=_y;}
      inline P operator-(P b){return P(x-b.x,y-b.y);}
      inline ll operator*(P b){return(ll)x*b.y-(ll)y*b.x;}
    }p[100010],A,B,q[1900000];
    void build(int x,int a,int b){
      q[l[x]=++t]=p[a];
      for(int i=a+1;i<=b+1;q[++t]=p[i++])while(t>l[x]&&(p[i]-q[t])*(q[t]-q[t-1])<=0)t--;
      r[x]=t;
      if(a==b)return;
      int mid=(a+b)>>1;
      build(x<<1,a,mid),build(x<<1|1,mid+1,b);
    }
    inline bool ex(int l,int r){
      int mid,t=l++;
      while(l<=r){
        mid=(l+r)>>1;
        if((q[mid]-A)*(B-q[mid])<(q[mid-1]-A)*(B-q[mid-1]))l=(t=mid)+1;else r=mid-1;
      }
      return(q[t]-A)*(B-q[t])<0;
    }
    int ask(int x,int a,int b){
      if(c<=a){
        if(!ex(l[x],r[x]))return 0;
        if(a==b)return a;
      }
      int mid=(a+b)>>1;
      if(c<=mid){
        int t=ask(x<<1,a,mid);
        if(t)return t;
      }
      return ask(x<<1|1,mid+1,b);
    }
    inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
    int main(){
      for(read(T);T--;puts("0")){
        for(read(n),i=1;i<=n;i++)read(p[i].x),read(p[i].y);
        for(t=0,build(1,1,n-1),i=1;i<n-1;i++)A=p[i],B=p[i+1],c=i+1,printf("%d ",ask(1,1,n-1));
      }
      return 0;
    }
    

      

  • 相关阅读:
    5月16日 AJAX
    5月13日 Jquery基础
    5月13日 PDO封装
    会话控制练习
    文件上传
    会话控制
    注册审核
    练习
    三级联动
    AJAX之XML
  • 原文地址:https://www.cnblogs.com/clrs97/p/4700870.html
Copyright © 2020-2023  润新知