hdu4353 Finding Mine三角形内的点数

http://acm.hdu.edu.cn/showproblem.php?pid=4353

题意：

　　求多边形面积和这个多边形内的金矿数的比值的最小值。

　　当xi<xj<xk时：

　　三角形内的点数=|ik上方的点-(ij上方的点+jk上方的点)|

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <algorithm>
 4 using namespace std;
 5 const int Ni = 205;
 6 const int Mi = 505;
 7 int n,m;
 8 double mx;
 9 struct node{
10     int x,y;//注意数据的范围，是否超int
11     bool operator < (const node &a) const
12     {
13         return x<a.x;
14     }
15 }a[Ni],b[Mi];
16 int sum[Ni][Ni];
17 inline double fabs(double a)
18 {
19     if(a>1e-8) return a;
20     if(-a>1e-8) return -a;
21     return 0;
22 }
23 bool cross(int i,int j,int k)
24 {
25     if(0<(a[j].x-a[i].x)*(b[k].y-a[i].y)-(a[j].y-a[i].y)*(b[k].x-a[i].x))
26         return 1;
27     else return 0;
28 }
29 inline int area(int i,int j,int k)
30 {
31     return (a[j].x-a[i].x)*(a[k].y-a[i].y)-(a[j].y-a[i].y)*(a[k].x-a[i].x);
32 }
33 inline bool judge(int i,int j,int k)
34 {
35     if(a[i].x<b[k].x&&b[k].x<=a[j].x){
36         if(cross(i,j,k)) return 1;
37     }
38     return 0;
39 }
40 void finds()
41 {
42     int i,j,k;
43     for(i=0;i<n;i++)
44     for(j=i+1;j<n;j++)
45     {
46         sum[i][j]=0;
47         for(k=0;k<m;k++) if(judge(i,j,k))
48             sum[i][j]++;
49     }
50     mx=1<<27;
51     for(i=0;i<n;i++)
52     for(j=i+1;j<n;j++)
53     for(k=j+1;k<n;k++)
54     {
55         int ans=sum[i][k]-sum[i][j]-sum[j][k];
56         double s=area(i,j,k);
57         if(ans&&(s=fabs(s/ans))<mx)
58         {
59              mx=s;
60         }
61     }
62 }
63 int main()
64 {
65     int cs=1,tcs,i;
66     scanf("%d",&tcs);
67     while(tcs--)
68     {
69         scanf("%d%d",&n,&m);
70         for(i=0;i<n;i++) scanf("%d%d",&a[i].x,&a[i].y);
71         sort(a,a+n);
72         for(i=0;i<m;i++) scanf("%d%d",&b[i].x,&b[i].y);
73         finds();
74         if(mx!=1<<27)
75             printf("Case #%d: %lf\n",cs++,mx/2);
76         else
77             printf("Case #%d: -1\n",cs++);
78     }
79     return 0;
80 }