hdu 4353 Finding Mine 计算几何

找到一个三角形，使得A/B的值最小，可证明，比这个三角形大的多边形，A/B的值都会比这个三角形的比值大。

暴力肯定会T，需要加一个优化。

传送门：http://acm.hdu.edu.cn/showproblem.php?pid=4353

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <iostream>
#include <algorithm>

#define eps 1e-8

using namespace std;

struct Node {
    int x,y; 
};

Node a[205],b[505];
int num[505][505];
int n,m,t;

int cross(Node a1,Node a2,Node a3){
    return (a2.x-a1.x)*(a3.y-a2.y) - (a2.y-a1.y)*(a3.x-a2.x);
}

void input(){
    scanf("%d %d",&n,&m);
    for(int i=0;i<n;i++){
        scanf("%d %d",&a[i].x,&a[i].y);
    }
    for(int i=0;i<m;i++){
        scanf("%d %d",&b[i].x,&b[i].y);
    }
}

bool cmp(Node a1,Node a2){
    return a1.x < a2.x;
}

int get_num(Node a1,Node a2){
    int cnt = 0;
    for(int i=0;i<m;i++){
        if(b[i].x>=a1.x&&b[i].x<a2.x){
            if(cross(b[i],a1,a2)>0)
                cnt++;
        }
    } 
    return cnt;
}

void print(){
    for(int i=0;i<n;i++){
        printf("%d : %d %d\n",i,a[i].x,a[i]);
    }
}

void solve(){
    double ans = -1;
    printf("Case #%d: ",t);
    sort(a,a+n,cmp);
    //print();
    memset(num,0,sizeof(num));
    for(int i=0;i<n-1;i++){
        for(int j=i+1;j<n;j++){
            num[i][j] = get_num(a[i],a[j]);
        }
    }
    for(int i=0;i<n-2;i++){
        for(int j=i+1;j<n-1;j++){
            for(int k=j+1;k<n;k++){
                double area = fabs(cross(a[i],a[j],a[k]));
                int cnt = abs(num[i][k] - num[i][j] - num[j][k]);
                if(cnt==0)continue;
                if(area/cnt<ans||ans==-1){
                    ans = area / cnt;
                }
            }
        }
    }
    if(ans==-1)puts("-1");
    else printf("%.6lf\n",ans/2);
}

void File(){
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
}

int main(void){
    //File();
    int T;
    t = 0;
    scanf("%d",&T);
    while(T--){
        t ++;
        input();
        solve();
    }
    return 0;
}