HDU 4353

利用分式的性质可以很容易证明要求的是个三角形，这很简单。对于求三角形内的雷的个数，只需求出每条边上方有多少个雷，作一点运算即可。如

A,B,C（B是X轴坐标在中间的点），则AC（其上方的雷的个数）-AB-BC即可。

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cmath>
 using namespace std;

 struct node{
     int x,y;
 };
 node point[],mine[];
 int cpoint[][];

 void exchange(int &i,int &j){
     int tmp=i;
     i=j;
     j=tmp;
 }

 int count_point(int i,int j,int cm){
     int cp=;
     if(point[i].x>point[j].x)
     exchange(i,j);
     for(int k=;k<cm;k++){
         if(mine[k].x>=point[i].x&&mine[k].x<point[j].x){\\必须是开区间，否则会有重复计算，对于中间的一点。。
             if((point[i].x-mine[k].x)*(point[j].y-mine[k].y)-(point[j].x-mine[k].x)*(point[i].y-mine[k].y)>)
             cp++;
         }
     }
     return cp;
 }

 int main(){
     int T,n,m,kase=;
     scanf("%d",&T);
     while(T--){
         double ans=1e15;
         scanf("%d%d",&n,&m);
         for(int i=;i<n;i++){
             scanf("%d%d",&point[i].x,&point[i].y);
         }
         for(int i=;i<m;i++){
             scanf("%d%d",&mine[i].x,&mine[i].y);
         }
         memset(cpoint,,sizeof(cpoint));
         for(int i=;i<n;i++){
             for(int j=i+;j<n;j++)
             cpoint[i][j]=count_point(i,j,m);
         }
         int i,j,k;
         for(int ii=;ii<n;ii++){
             for(int jj=ii+;jj<n;jj++){
                 for(int kk=jj+;kk<n;kk++){
                     i=ii; j=jj;k=kk;
                     if(point[i].x>=min(point[j].x,point[k].x)&&point[i].x<=max(point[j].x,point[k].x))
                     exchange(i,i);
                     else if(point[j].x>=min(point[i].x,point[k].x)&&point[j].x<=max(point[i].x,point[k].x))
                     exchange(i,j);
                     else if(point[k].x>=min(point[i].x,point[j].x)&&point[k].x<=max(point[i].x,point[j].x))
                     exchange(i,k);
                     double area=fabs(((point[j].x-point[i].x)*(point[k].y-point[i].y)-(point[k].x-point[i].x)*(point[j].y-point[i].y))/2.0);
                     double cp=abs((cpoint[i][j]+cpoint[j][i])+(cpoint[i][k]+cpoint[k][i])-(cpoint[j][k]+cpoint[k][j]));
                     ans=min(ans,area/cp);
                 }
             }
         }
         printf("Case #%d: ",++kase);
         if(ans==1e15)
         printf("-1\n");
         else printf("%.6lf\n",ans);
     }
     return ;
 }