省选专练USACO TRAIN5.1 圈奶牛

凸包入门题。

介绍一下我的算法：

Graham扫描法

利用栈的思想。

与Jarris步进法不同的是，这个算法不适用角度。

他利用了向量的叉乘。

根据有向的平行四边形来判断当前点是否更优。

OK看代码吧。

(*^▽^*)。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdio>
using namespace std;
const int N=10100;
const int INF=1e9+10;
int n;
double sqr(double sum){
	return sum*sum;
}
struct Points{
	double x,y;
	void read(){
		scanf("%lf%lf",&x,&y);
	}
	bool operator <(const Points &a)const{
		return x<a.x||(x==a.x&&y<=a.y);
	}
	Points operator -(const Points &a)const{
		return (Points){x-a.x,y-a.y};
	}
	double operator *(const Points &a)const{
		return x*a.y-y*a.x;
	}
	friend double dis(const Points &a,const Points &b){
		return sqrt(sqr(b.x-a.x)+sqr(b.y-a.y));
	}
}p[N],now;
int top=0;
Points q[N];
void push(Points now){
	while((now-q[top-1])*(q[top]-q[top-1])>0)top--;
	top++;
	q[top]=now;
}
void make_convex_hull(){
	q[0]=q[top=1]=p[1];
	for(int i=2;i<=n;i++){
		push(p[i]);
	}
	for(int i=n-1;i>=1;i--){
		push(p[i]);
	}
}
double len(){
	double ans=0;
	for(int i=1;i<top;i++){
		ans+=dis(q[i],q[i+1]);
	}
	return ans;
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		p[i].read();
	}
	sort(p+1,p+1+n);
	make_convex_hull();
	printf("%.2lf",len());
}