HDU 4353
利用分式的性质可以很容易证明要求的是个三角形，这很简单。对于求三角形内的雷的个数，只需求出每条边上方有多少个雷，作一点运算即可。如
A,B,C（B是X轴坐标在中间的点），则AC（其上方的雷的个数）-AB-BC即可。
1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 #include <cmath> 5 using namespace std; 6 7 struct node{ 8 int x,y; 9 }; 10 node point[250],mine[550]; 11 int cpoint[250][250]; 12 13 void exchange(int &i,int &j){ 14 int tmp=i; 15 i=j; 16 j=tmp; 17 } 18 19 int count_point(int i,int j,int cm){ 20 int cp=0; 21 if(point[i].x>point[j].x) 22 exchange(i,j); 23 for(int k=0;k<cm;k++){ 24 if(mine[k].x>=point[i].x&&mine[k].x<point[j].x){\\必须是开区间，否则会有重复计算，对于中间的一点。。 25 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)>0) 26 cp++; 27 } 28 } 29 return cp; 30 } 31 32 int main(){ 33 int T,n,m,kase=0; 34 scanf("%d",&T); 35 while(T--){ 36 double ans=1e15; 37 scanf("%d%d",&n,&m); 38 for(int i=0;i<n;i++){ 39 scanf("%d%d",&point[i].x,&point[i].y); 40 } 41 for(int i=0;i<m;i++){ 42 scanf("%d%d",&mine[i].x,&mine[i].y); 43 } 44 memset(cpoint,0,sizeof(cpoint)); 45 for(int i=0;i<n;i++){ 46 for(int j=i+1;j<n;j++) 47 cpoint[i][j]=count_point(i,j,m); 48 } 49 int i,j,k; 50 for(int ii=0;ii<n;ii++){ 51 for(int jj=ii+1;jj<n;jj++){ 52 for(int kk=jj+1;kk<n;kk++){ 53 i=ii; j=jj;k=kk; 54 if(point[i].x>=min(point[j].x,point[k].x)&&point[i].x<=max(point[j].x,point[k].x)) 55 exchange(i,i); 56 else if(point[j].x>=min(point[i].x,point[k].x)&&point[j].x<=max(point[i].x,point[k].x)) 57 exchange(i,j); 58 else if(point[k].x>=min(point[i].x,point[j].x)&&point[k].x<=max(point[i].x,point[j].x)) 59 exchange(i,k); 60 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); 61 double cp=abs((cpoint[i][j]+cpoint[j][i])+(cpoint[i][k]+cpoint[k][i])-(cpoint[j][k]+cpoint[k][j])); 62 ans=min(ans,area/cp); 63 } 64 } 65 } 66 printf("Case #%d: ",++kase); 67 if(ans==1e15) 68 printf("-1\n"); 69 else printf("%.6lf\n",ans); 70 } 71 return 0; 72 }
