HDU 2202 最大三角形 【凸包+枚举】

最大三角形

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

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

  
  

   
   3
3 4
2 6
3 7
6
2 6
3 9
2 0
8 0
6 6
7 7
  
  

Sample Output

  
  

   
   1.50
27.00
  
  

Author

Eddy

Recommend

lcy

数据比较多，全部枚举会超时。最大三角形三个顶点一定是在凸包的顶点上。

下面是代码。

#include <iostream>
#include <cstdio>
#include <math.h>
#include <algorithm>
#include <string.h>
#define maxn 50005
using namespace std;
struct point{
    int x,y;
}p[maxn];
point temp,pp;
int stack[maxn]={0};
int cross(point a,point b,point c)
{
    return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}
double d(point a,point b)
{
    return sqrt((double)(b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y));
}
double sss(double a,double b,double c)
{
    double t=(a+b+c)/2;
    double s=sqrt(t*(t-a)*(t-b)*(t-c));
    return s;
}
bool cmp (point a,point b)
{
    int jj=cross(p[0],a,b);
    if(jj>0) return true;
    else if(jj==0&&(d(p[0],a)<d(p[0],b)))
    return true;
    else return false;
}
int main()
{
    int total,i,k,j,top;
    double sum=0;
    while(scanf("%d",&total)!=EOF)
    {
        memset(p,0,sizeof(p));
        memset(stack,0,sizeof(stack));
        k=0;
        scanf("%d%d",&p[0].x,&p[0].y);
        pp.x=p[0].x;
        pp.y=p[0].y;
        for(i=1;i<total;i++)
        {
            scanf("%d%d",&p[i].x,&p[i].y);
            if(p[i].y<pp.y||(p[i].y==pp.y&&p[i].x<pp.x))
               {
                   pp.x=p[i].x;
                   pp.y=p[i].y;
                   k=i;
               }
        }
        p[k]=p[0];
        p[0]=pp;
        sort(p+1,p+total,cmp);
        stack[0]=0;
        stack[1]=1;
        top=1;
        for(i=2;i<total;i++)
        {
        while(top>0&&cross(p[stack[top-1]],p[stack[top]],p[i])<=0)
        top--;
        top++;
        stack[top]=i;
        }//类似栈的思想
        for(int i=0;i<=top-1;i++)
        {
            for(int j=0;j<=top;j++)
            {
                for(int k=0;k<=top;k++)
                if(j==i||j==k||i==k)
                continue;//不能构成三角形，跳过
            else if(sss(d(p[stack[i]],p[stack[j]]),d(p[stack[k]],p[stack[j]]),d(p[stack[k]],p[stack[i]]))>sum)
            sum=sss(d(p[stack[i]],p[stack[j]]),d(p[stack[k]],p[stack[j]]),d(p[stack[i]],p[stack[k]]));
            }
        }
        printf("%.2lf\n",sum);
        sum=0;

    }
}