HDU2202

最新推荐文章于 2017-07-23 17:33:13 发布

克利斯提亚诺-梅西

最新推荐文章于 2017-07-23 17:33:13 发布

阅读量699

点赞数 1

分类专栏：邵大的学习

本文链接：https://blog.csdn.net/MyHeaven7/article/details/52193566

版权

邵大的学习专栏收录该内容

63 篇文章 0 订阅

订阅专栏

最大三角形

Time Limit: 5000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 3932 Accepted Submission(s): 1387

Problem Description

老师在计算几何这门课上给Eddy布置了一道题目，题目是这样的：给定二维的平面上n个不同的点，要求在这些点里寻找三个点，使他们构成的三角形拥有的面积最大。
Eddy对这道题目百思不得其解，想不通用什么方法来解决，因此他找到了聪明的你，请你帮他解决这个题目。

Input

输入数据包含多组测试用例，每个测试用例的第一行包含一个整数n，表示一共有n个互不相同的点，接下来的n行每行包含2个整数xi,yi，表示平面上第i个点的x与y坐标。你可以认为：3 <= n <= 50000 而且 -10000 <= xi, yi <= 10000.

Output

对于每一组测试数据，请输出构成的最大的三角形的面积，结果保留两位小数。
每组输出占一行。

Sample Input

Sample Output

题意：就是求所给点中，能组成三角形的最大的面积；

思路：利用三角形外接圆，即找到一个最小的圆能够包围所有的点（利用凸包），那么在这个圆上至少存在三个点以上，所以，此时就可以使用暴力求出最大三角形面积了；

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<iomanip>
using namespace std;
const int Max=50005;
struct Node
{
    int x;
    int y;
}a[Max],b[Max];
int cmp(Node m,Node n)
{
    if(m.x==n.x)
        return m.y<n.y;
    else
        return m.x<n.x;
}
int Count(Node i,Node j,Node k) //就算三角形的面积；
{
    int x1=i.x-j.x;
    int y1=i.y-j.y;
    int x2=k.x-j.x;
    int y2=k.y-j.y;
    return(x1*y2-x2*y1);
}
int convex(int n) //凸包计算；求出在圆上的点得数量m；和坐标。
{
    sort(a,a+n,cmp);
    int m=0,i,j,k;
    for(i=0;i<n;i++)
    {
        while(m>1&&Count(b[m-1],a[i],b[m-2])<=0)
            m--;
        b[m++]=a[i];
    }
         k=m;
    for(i=n-2;i>=0;i--)
    {
        while(m>k&&Count(b[m-1],a[i],b[m-2])<=0)
            m--;
        b[m++]=a[i];
    }
    if(n>1)
        m--;
    return m;
}
int main()
{
    int i,j,k,n,m;
    int  sum=0;
    while(cin>>n)
    {
        for(i=0;i<n;i++)
            cin>>a[i].x>>a[i].y;
            sum=0;
        m=convex(n);
        //cout<<m<<endl;
        for(i=0;i<m;i++)
            for(j=i+1;j<m;j++)
                for(k=j+1;k<m;k++)
        {
            sum=max(sum,Count(b[i],b[j],b[k]));
        }
        //sum=sum/2.0;
        printf("%.2lf\n",0.5*sum);
       // cout<<fixed<<setprecision(2)<<sum<<endl;
    }
    return 0;
}