HDU 2202(凸包內最大三角形面積)

先將凸包點求出來之後。

利用旋轉卡殼的方法,然後找不同兩點爲邊,在找其對踵點。三點求面積 取最大值就可以了。


#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
const int maxn=50001;
struct node
{
    double  x,y;
} st[maxn],num[maxn];
int n,k,e;
double dis(node p1,node p2)
{
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
int multi(node p1,node p2,node p3)//叉積
{
    return (p1.x-p3.x)*(p2.y-p3.y)-(p2.x-p3.x)*(p1.y-p3.y);
}
bool cmp(node p1,node p2)
{
    if(multi(p1,p2,num[0])>0) return true;
    if(multi(p1,p2,num[0])==0&&dis(p1,num[0])<dis(p2,num[0]))return true;
    return false;
}
void Graham()//凸包
{
    e=2;
    sort(num+1,num+n,cmp);
    st[0]=num[0];
    st[1]=num[1];
    st[2]=num[2];
    for(int i=3; i<n; i++)
    {
        while(e>1&&multi(num[i],st[e],st[e-1])>=0)
            e--;
        st[++e]=num[i];
    }
}
void rotating_calipers()
{
    e++;
    int q=1;
    int ans=0;
    st[e]=st[0];
    for(int i=0; i<e; i++)
    {
        q=1;
        for(int p=i+1; p<e; p++)
        {
            while(multi(st[i],st[q+1],st[p])>multi(st[i],st[q],st[p]))
                q=(q+1)%e;
            ans=max(ans,multi(st[i],st[q],st[p]));
        }
    }//求最大邊的方法,但是要考慮到不一定以p,p+1爲底的三角形,
    double sss=(double)ans/2;
    printf("%.2lf\n",sss);
}
int main()
{
    while(~scanf("%d",&n))
    {
        k=0;//K表示起始點
        for(int i=0; i<n; i++)
        {
            scanf("%lf %lf",&num[i].x,&num[i].y);
            if(num[i].x<num[k].x||(num[i].x==num[k].x&&num[i].y<num[k].y)) k=i;
        }
        swap(num[0],num[k]);
        //printf("%d\n",k);
        Graham();
        rotating_calipers();
    }
    return 0;
}


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章