HDU4709 聚集（求三角形面积）

HDU4709 聚集（求三角形面积）

平面上给出N（1<=N<=100）个点，每个点给出的是浮点数的二位平面坐标（x，y）。现在要求的是这N个点中由最少3个点构成的图形中的最小面积。

输入：T（1<=T<=25）表输入实例个数。接下来每个实例第一行是N，接下来N行每行包括两个浮点数表坐标（X,Y）（-1000.0<=X,Y<=1000.0）.

输出：输出找到的最小面积值（保留小数点后两位），如果最小面积值是0.00，则输出“Impossible”

分析：本题只需枚举所有的3点，算他们构成的三角形面积即可。其中三角形面积用二维叉积计算。注意：一般浮点运算的题目都要设计eps来控制精度。

AC代码：

#include<cstdio>#include<cmath>#include<algorithm>using namespace std;const double eps=1e-8;struct point{    double x;    double y;    point(double _x=0,double _y=0){x=_x;y=_y;}    point operator -(const point &b)const    {        return point(x-b.x,y-b.y);    }    double operator ^(const point &b)const    {        return x*b.y-y*b.x;    }    void input()    {        scanf("%lf%lf",&x,&y);    }};int main(){    int t,n;    scanf("%d",&t);    point p[200];    while(t--)    {        scanf("%d",&n);        for(int i=0;i<n;i++)            p[i].input();        double min_area=1e9;        bool ok=false;        for(int i=0;i<n;i++)            for(int j=i+1;j<n;j++)                for(int k=j+1;k<n;k++)                {                    double area = fabs(((p[i]-p[j])^(p[i]-p[k]))/2);                    if(area<eps)continue;                    ok=true;                    min_area = min(min_area,area);                }        if(ok==false)printf("Impossible\n");        else printf("%.2lf\n",min_area);    }    return 0;}