旋转卡壳模版代码

本文提供了一种计算平面上给定点集所构成的最大三角形面积的算法实现。通过使用Graham扫描法找到凸包,然后遍历凸包上的点来确定最大的三角形面积。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

旋转卡壳模版代码
给定一个平面上的n个点,找出以这些为顶点的三角形,面积最大的一个,输出面积。

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

struct Point   
{
    double x,y;
};
double max(double a,double b)
{
	if(a>b)return a;
	return b;
}
bool mult(Point curr, Point end, Point begin)
{
    return (curr.x-begin.x)*(end.y-begin.y)
		>=(end.x-begin.x)*(curr.y-begin.y);
}
bool comp(const Point &a,const Point &b)
{
    return a.y<b.y||(a.y==b.y&&a.x<b.x);
}
int graham(Point pnt[], int n, Point res[])  //凸包模版代码
{
    int i,len,top=1;
    sort(pnt,pnt+n,comp);
    if(n==0)
    {
        return 0;
    }
    res[0]=pnt[0];
    if(n==1)
    {
        return 1;
    }
    res[1]=pnt[1];
    if(n==2)
    {
        return 2;
    }
    res[2]=pnt[2];
    for(i=2;i<n;i++)
    {
        while(top && mult(pnt[i],res[top],res[top-1]))
        {
            top--;
        }
        res[++top]=pnt[i];
    }
    len=top;
    res[++top]=pnt[n-2];
    for(i=n-3;i>=0;i--)
    {
        while(len!=top&&mult(pnt[i],res[top],res[top-1]))
        {
            top--;
        }         
        res[++top]=pnt[i];
    }
    return top;
}
double Area(Point a,Point b,Point c)
{
    return fabs((a.x-c.x)*(b.y-c.y)-(b.x-c.x)*(a.y-c.y))/2;
}

int main()
{
	Point pnt[400],res[400];
	int n;
	

	while(scanf("%d",&n)!=EOF)
	{
		for(int i=0;i<n;i++)
			scanf("%lf%lf",&pnt[i].x,&pnt[i].y);
		int top=graham(pnt,n,res);
		/*
		for(i=0;i<top;i++)
		{
		printf("%lf %lf\n",res[i].x,res[i].y);
		
		  }*/
		double ans=0.0;
		for(i=0;i<top;i++)
		{
			int j=(i+1)%top;
			int k=(j+1)%top;
			
			while(k!=i && Area(res[i],res[j],res[k])<Area(res[i],res[j],res[(k+1)%top]))
			{
				k=(k+1)%top;
			}
			if(k==i)
				continue;
			
			int kk=(k+1)%top;
			while(j!=kk && k!=i)
			{ 
				ans=max(ans,Area(res[i],res[j],res[k]));
				while(k!=i && Area(res[i],res[j],res[k])<Area(res[i],res[j],res[(k+1)%top]) )
					k=(k+1)%top;
					j=(j+1)%top;
			}
			
		}
		
		printf("%.2lf\n",ans);
		
	}
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值