2019-02-24 二维凸包

二维凸包:

#include <bits/stdc++.h>
using namespace std;
const int maxn=100000+10;
int n,top=2;
double ans;
struct node{//存储直角坐标系的x和y值
	double x,y;
};
node p[maxn],z[maxn];
bool cmp(node u,node v){//比较大小
	return u.x==v.x?u.y<v.y:u.x<v.x;//还需要判断u.x与v.x是否相等，如果相等则以y排序。
}
double xl(node a,node b){
	if(abs(a.x-b.x)<0.0000000001)//double类型不能直接比较，但可以相减判断精度
		return 1000000000.00;//跟INF差不多
	else	
		return (a.y-b.y)/(a.x-b.x);//计算斜率
}
bool cmpxl(node a,node b,node c){
	return xl(a,b)<xl(b,c);//判断是否斜率减小，即能否更换
}
double dist(node a,node b){//计算两点之间的距离
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));//勾股定理计算距离
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%lf%lf",&p[i].x,&p[i].y);//输入每个数的坐标
	sort(p+1,p+n+1,cmp);//以x值排序（从左到右）
	z[1]=p[1],z[2]=p[2];//z为一个栈,先设其坐标为最左边的数
	for(int i=3;i<=n;i++){//从前往后
		while(top>1&&cmpxl(p[i],z[top],z[top-1]))//判断斜率是否一直上升
			top--;
		top++;
		z[top]=p[i];//入栈
	}
	for(int i=1;i<top;i++)//从最左边的点开始，计算下半部分的最大值
		ans+=dist(z[i],z[i+1]);
	z[1]=p[n],z[2]=p[n-1];//z为一个栈，设其坐标为最右边的数
	top=2;
	for(int i=n-2;i>=1;i--){//从后往前
		while(top>1&&cmpxl(p[i],z[top],z[top-1]))//判断斜率是否一直上升
			top--;
		top++;
		z[top]=p[i];//入栈
	}
	for(int i=1;i<top;i++)//从最右边的值开始，计算上半部分的最大值
		ans+=dist(z[i],z[i+1]);
	printf("%.2lf",ans);//输出
	return 0;
}

Wall
重要注释有注明，其他注释如上

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn=100000+10;
int n,top=2,k;
double ans;
struct node{
	double x,y;
};
node p[maxn],z[maxn];
bool cmp(node u,node v){
	return u.x==v.x?u.y<v.y:u.x<v.x;//还需要判断u.x与v.x是否相等，如果相等则以y排序。
}
double xl(node a,node b){
	if(abs(a.x-b.x)<0.0000000001)
		return 1000000000.00;
	else	
		return (a.y-b.y)/(a.x-b.x);
}
bool cmpxl(node a,node b,node c){
	return xl(a,b)<xl(b,c);
}
double dist(node a,node b){
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int main(){
	scanf("%d%d",&n,&k);
	for(int i=1;i<=n;i++)
		scanf("%lf%lf",&p[i].x,&p[i].y);
	sort(p+1,p+n+1,cmp);
	z[1]=p[1],z[2]=p[2];
	for(int i=3;i<=n;i++){
		while(top>1&&cmpxl(p[i],z[top],z[top-1]))
			top--;
		z[++top]=p[i];
	}
	for(int i=1;i<top;i++)
		ans+=dist(z[i],z[i+1]);
	z[1]=p[n],z[2]=p[n-1];
	top=2;
	for(int i=n-2;i>=1;i--){
		while(top>1&&cmpxl(p[i],z[top],z[top-1]))
			top--;
		z[++top]=p[i];
	}
	for(int i=1;i<top;i++)
		ans+=dist(z[i],z[i+1]);
	ans=ans+k*2*3.1415926535897932384626;//最终答案只要凸包加上以k为半径的圆的周长就行了。
	printf("%d",(int)(ans+0.5));//确保四舍五入
	return 0;
}