POJ 1113 Wall （凸包）

最新推荐文章于 2019-04-05 09:33:58 发布

原创最新推荐文章于 2019-04-05 09:33:58 发布 · 212 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#struct

ACM-计算几何专栏收录该内容

5 篇文章

订阅专栏

本文介绍了一个计算任意多边形（包括凹多边形）最小包围周长的方法，通过构造凸包并考虑固定距离的偏移，最终得出总周长。代码实现了关键步骤，包括点排序、凸包生成及周长计算。

给一个多边形（可能是凹的），然后给出一个最小周长围住这个多边形，围的墙至少要离这个多边形一定距离。

凸包，最后加上求周长即可。

#include<iostream>
#include<algorithm>
#include<cmath>
#define PI 3.141592653
using namespace std;

typedef struct{
	int x, y;
} POINT;

POINT ch[1000];

double dist(POINT a, POINT b){
	return sqrt(pow(a.x - b.x, 2.0) + pow(a.y - b.y, 2.0));
}

int crossProduct(POINT o, POINT a, POINT b){
	return (a.x - o.x) * (b.y - o.y) - (b.x - o.x) * (a.y - o.y);
}

bool cmp_position(POINT a, POINT b){
	return a.y < b.y || (a.y == b.y && a.x < b.x);
}

bool cmp_angle(POINT a, POINT b){
	int temp = crossProduct(ch[0], a, b);
	if(temp == 0)
		return (dist(a, ch[0]) < dist(b, ch[0]));
	else
		return temp > 0;
}

void print(POINT a[], int N){
	for(int i = 0; i < N; i++)
		cout << a[i].x << " " << a[i].y << endl;
}

int main(){
	int N, L;
	cin >> N >> L;
	for(int i = 0; i < N; i++)
		cin >> ch[i].x >> ch[i].y;
	swap(ch[0], *min_element(ch, ch + N, cmp_position));
	sort(ch + 1, ch + N, cmp_angle);
	ch[N] = ch[0];
	POINT s[1100];
	int top = 2;
	s[0] = ch[0];
	s[1] = ch[1];
	for(int i = 2; i <= N; i++){
		while(top >= 2 && crossProduct(s[top - 2], s[top - 1], ch[i]) <= 0){
			top--;
		}
		s[top++] = ch[i];
	}
	//print(s, top);
	double result = 0;
	for(int i = 1; i < top; i++)
		result += dist(s[i - 1], s[i]);
	result += PI * L * 2;
	cout << (int)(result + 0.5) << endl;
}