Jungle Outpost HDU

Jungle Outpost

HDU - 3761

题意：

半平面交～

用了动态分配内存但是忘记释放了一直ＭＬＥ．．．．．．

#include <bits/stdc++.h>
using namespace std;
const int maxn = 50010;
const int inf = 0x3f3f3f3f;
const double eps = 1e-8;
struct Point {
    double x,y;
    Point (double x = 0, double y = 0) : x(x), y(y) {}
};
typedef Point Vector;
Vector operator + (Vector a, Vector b) {
    return Vector (a.x + b.x, a.y + b.y);
}
Vector operator * (Vector a, double s) {
    return Vector (a.x * s, a.y * s);
}
Vector operator / (Vector a, double p) {
    return Vector (a.x / p, a.y / p);
}
Vector operator - (Point a, Point b) {
    return Vector (a.x - b.x, a.y - b.y);
}
bool operator < (Point a, Point b) {
    return a.x < b.x || (a.x == b.x && a.y < b.y);
}
int dcmp (double x) {
    if(fabs(x) < eps) return 0;
    return x < 0 ? -1 : 1;
}
bool operator == (const Point &a, const Point &b) {
    return dcmp(a.x - b.x) == 0 && dcmp(a.y - b.y) == 0;
}
double Dot(Vector a, Vector b) {
    return a.x * b.x + a.y * b.y;
}
double Cross (Vector a, Vector b) {
    return a.x * b.y - a.y * b.x;
}
//半平面交
struct Line{
    Point p;
    Vector v;
    double rad;
    Line () {}
    Line (Point p, Vector v) : p(p), v(v) {
        rad = atan2(v.y,v.x);
    }
    bool operator < (const Line &L) const {
        return rad < L.rad;
    }
};
bool OnLeft(Line L, Point p) {
    return Cross(L.v, p - L.p) > 0;
}
Point GetLineIntersection (Line a, Line b) {
    Vector u = a.p - b.p;
    double t = Cross(b.v, u) / Cross(a.v, b.v);
    return a.p + a.v*t;
}
int HalfplaneIntersection(Line *L, int n) {
    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(dcmp(Cross(q[last].v, q[last-1].v)) == 0) {
            last--;
            if(OnLeft(q[last], L[i].p)) q[last] = L[i];
        }
        if(first < last) p[last-1] = GetLineIntersection (q[last-1],q[last]);
    }
    while(first < last && !OnLeft(q[first], p[last-1])) last--;  //删除无用平面
    delete []p;
    delete []q;
    if(last - first <= 1) return 0;  //空集
    return 1;
}
Point pp[maxn];
Line line[maxn];
int n;
bool ck(int m){
    int cnt = 0;
    for(int i = 0; i < n; i++){
        line[cnt++] = Line(pp[i], pp[(i+m+1)%n] - pp[i]);
    }
    int ans = HalfplaneIntersection(line, cnt);
    return ans==0;
}
int main(){
    int t;
    //freopen("in.txt", "r", stdin);
    scanf("%d", &t);
    while(t--){
        scanf("%d", &n);
        for(int i = n-1; i >= 0; i--){
            scanf("%lf %lf", &pp[i].x, &pp[i].y);
        }
        int L = 1, R = n-3;
        while(L <= R) {
            int m = (L+R)/2;
            if(ck(m)) R = m-1;
            else L = m+1;
        }
        printf("%d
", L);
    }
}