计算凸包的周长与面积

 

先求凸包的顶点,再叉积求面积,两点间的距离公式求周长

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;

const double eps = 1e-8,pi=3.14159265;
// 若 一点的坐标在前面两点的右边,则更新右端点 
int n;

int dcmp(double x)  // 三态函数 ,在下面用于是否在精度 
{
	if(fabs(x)0?1:-1;
} 
#define Vector Point
struct Point 
{
	double x,y;
	inline Point(double x=0,double y=0):x(x),y(y){}
}p[10000+5],ch[10000+5];

bool myCmp(Point A,Point B){
	if(A.x!=B.x)return A.x1&&dcmp(Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2]))<=0)m--;
		ch[m++]=p[i];
	}
	int k=m;
	for(int i=n-2;i>=0;i--){  // 因 pn-1 已经遍历过,故从n-2开始 
		// 求上凸包,最终形成合围,为完整凸包
		while(m>k&&dcmp(Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2]))<=0)m--;
		ch[m++]=p[i]; 
	}
	if(n>1)m--;  // 因要形成凸包,故起始点重复,减1
	return m;  //返回顶点个数
}
double Dis(Point A,Point B){
	return sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y));
}
int main(){
	int m;
	cin >> n;
	for(int i=0;i> a >> b;
		p[i]=Point(a,b);
	}
	m=ConvexHull();
	
	//下面是计算凸包周长的 
/*	double ans=0.0; 
	for(int i=0;i 

你可能感兴趣的:(计算凸包的周长与面积)