UVA LA2797 二分加半平面交，很巧妙的二分

二分真是个神奇的东西，又是一个二分题目，而且这题分析的半平面交也比较巧妙，本人智商不够，二分和半平面交均想不到，

并且这题二分为什么会对，本渣渣也不太明白，留着慢慢的明白吧。

分析：

如果敌人只有一颗炸弹，你会把总部建在哪里呢？对于每个点，炸掉它以后不会暴露在外面的区域是一条有向直线的”左边“。这让我们想到了什么？

没错，半平面交！综合考虑所有点，所以建总部的范围就是所有这些半平面的交。

如果敌人有俩颗炸弹，总部应该建在哪里呢？分析后发现，敌人最聪明的做法是炸掉俩个连续的顶点，而不是分散火力（想一想，为什么）。这样，

每两个连续顶点对应一个新的半平面，可以建总部的范围仍然是所有半平面的交。

有了上面的分析，整个问题迎刃而解：二分答案，用上述方法判断大难是否可行（也就是半平面交是否为非空）。二分需要O（logn)时间，半平面交需要

O（nlogn）时间，总时间复杂度为O（n^2logn).

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cctype>
#include<string>
#include<set>
#include<map>
#include<queue>
#include<stack>
using namespace std;
const double eps=1e-16;
const int maxn=50000+10;
int dcmp(double x)
{
    if(fabs(x)<eps) return 0;
    return x>0?1:-1;
    //return fabs(x) < eps ? 0 : (x > 0 ? 1 : -1);
}

struct point
{
  double x;
  double y;
  point(){}
  point(double x,double  y):x(x),y(y){}
  void in()
  {
      cin>>x>>y;
  }
  void out()
  {
      cout<<x<<' '<<y<<endl;
  }
  point operator + (const point &t) const
  {
      return point(x+t.x,y+t.y);
  }
  point operator - (const point &t) const
  {
      return point(x-t.x,y-t.y);
  }
  point operator * (const double &t) const
  {
      return point(x*t,y*t);
  }
};
double cross(point a,point b)
{
    return a.x*b.y-a.y*b.x;
}
double dot(point a,point b)
{
    return a.x*b.x+a.y*b.y;
}
double length(point a)
{
    return sqrt(dot(a,a));
}
point nomal(point t)
{
    double l=length(t);
    return  point(-t.y/l,t.x/l);
}
struct line
{
    point p;
    point v;
    double ang;
    line() {}
    line(point p,point v):p(p),v(v){
        ang=atan2(v.y,v.x);
    }
    bool operator < (const line &l) const
    {
        return ang<l.ang;
    }
};

bool onleft(line l,point p)
{
    return cross(l.v,p-l.p)>0;
}
point getintersection(line a,line b)
{
    point u=a.p-b.p;
    double t=cross(b.v,u)/cross(a.v,b.v);
    return a.p+a.v*t;
}
point p[maxn];
line q[maxn];
point poly[maxn],pt[maxn];
line l[maxn];
int halfplaneintersection(line *l,int n,point *poly)
{
    sort(l,l+n);

    int first,last;
  //  point *p=new point[n];
 //   line *q=new  line[n];
    q[first=last=0]=l[0];
    for(int i=1;i<n;i++)
    {
        while(first<last && !onleft(l[i],p[last-1])) last--;
        while(first<last && !onleft(l[i],p[first])) first++;
        q[++last]=l[i];
        if(fabs(cross(q[last].v,q[last-1].v))<eps)
        {
            last--;
            if(onleft(q[last],l[i].p)) q[last]=l[i];
        }
        if(first<last) p[last-1]=getintersection(q[last-1],q[last]);
    }
    while(first<last && !onleft(q[first],p[last-1])) last--;
    if(last-first<=1) return 0;
    p[last]=getintersection(q[last],q[first]);
    int m=0;
    for(int i=first;i<=last;i++) poly[m++]=p[i];
    return m;
}


int slove(int n)
{
    if(n==3)
        return 1;
    int ll=1,r=n-3;
    while(ll<=r)
    {
        int c=0;
        int m=ll+r>>1;
        for(int i=0;i<n;i++)
            l[c++]=line(pt[(i+m+1)%n],pt[i]-pt[(i+m+1)%n]);
        if(halfplaneintersection(l,c,poly))
            ll=m+1;
        else
            r=m-1;
    }
    return ll;
}
int main()
{
    int n;
    while(cin>>n)
    {
        for(int i=0;i<n;i++)
            pt[i].in();
        cout<<slove(n)<<endl;
    }
    return 0;
}