HDU 4709 Herding

暴力+精度，精度把我坑了，狠狠的教训啊！

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4709

就是求所有点构成的三角形面积中，最小的面积。

直接上代码，一般控制好精度就没什么问题了。

#include<stdio.h>
#include<iostream>
using namespace std;
#include<queue>
#include<math.h>
#include<algorithm>
#include<string.h>

#define repA(p,q,i)  for( int (i)=(p); (i)!=(q); ++(i) )
#define repAE(p,q,i)  for( int (i)=(p); (i)<=(q); ++(i) )
#define repD(p,q,i)  for( int (i)=(p); (i)!=(q); --(i) )
#define repDE(p,q,i)  for( int (i)=(p); (i)>=(q); --(i) )
#define range 110
#define MAXN 4010*4010

double x[range],y[range];

double area(int i,int j,int k);

int main()
{
    int test,n;  scanf("%d",&test);
    while(test--)
    {
        scanf("%d",&n);
        repA(0,n,i)
            scanf("%lf%lf",&x[i],&y[i]);
        double minx = MAXN,s;
        repA(0,n,i)
        {
            repA(i+1,n,j)
            {
                repA(j+1,n,k)
                {
                    s = area(i,j,k) ;
                    if( s < 1e-6 )  continue;
                    minx = min(s,minx) ;
                }
            }
        }
        if( fabs( minx - MAXN ) < 1e-3 )  printf("Impossible\n");
        else printf("%.2lf\n",minx);
    }
}

double area(int i,int j,int k)
{
    return  fabs( (x[j]-x[i])*(y[k]-y[i]) - (y[j]-y[i])*(x[k]-x[i]) ) / 2 ;
}