凸包算法，并求凸包的周长

最新推荐文章于 2022-03-14 15:11:17 发布

自爄創煇熿

最新推荐文章于 2022-03-14 15:11:17 发布

阅读量754

点赞数

分类专栏：计算几何

本文链接：https://blog.csdn.net/chl_3205/article/details/8836094

版权

该程序实现了一个二维点集的Graham扫描算法，找到点集的凸包，并计算凸包的周长。输入包括点的数量和一个额外的参数，最后输出凸包周长（包括点集的周长和一个圆的周长）。程序通过叉积判断点的相对位置，按斜率排序，并使用动态规划方法构建凸包。

摘要由CSDN通过智能技术生成

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=1005;
const double PI=acos(-1);//是元周率的值。
struct point
{
	double x;
	double y;
};
point p[maxn];
point s[maxn];
double cross(point p0,point p1,point p2)  //叉积.
{
	return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
double dis(point a,point b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool cmp(point p1,point p2)  //按照斜率排序。
{
	if(cross(p1,p2,p[0])>0||(cross(p1,p2,p[0])==0&&(dis(p1,p[0])<dis(p2,p[0]))))
		return true;
    return false;
}
double max(double a,double b)
{
	if(a>b) return a;
    return b;
}
int graham(int n) //求凸包的函数，所求的点在s数组中。
{
	int top,i,tmp=0;
	for(i=1;i<n;i++)  //找基点。
	{
		if((p[i].y<p[tmp].y)||((p[i].y==p[tmp].y)&&(p[i].x<p[tmp].x)))
			tmp=i;
	}
	swap(p[0],p[tmp]);
	sort(p