// 依赖的计算几何模板
inline int dcmp(double x){
constexpr double eps=1e-6;
if (x<-eps)return -1;
if (x>eps)return 1;
return 0;
}
struct Point {
double x,y;
Point()=default;
Point(double x,double y):x(x),y(y){}
Point operator+(const Point& p){
return Point(x+p.x,y+p.y);
}
Point operator-(const Point& p){
return Point(x-p.x,y-p.y);
}
// dot mul.
Point operator*(double k){
return Point(x*k,y*k);
}
// cross mul.
double operator*(const Point& p){
return x*p.y-y*p.x;
}
double dist(const Point& p){
const Point dp=*this-p;
return sqrt(dp.x*dp.x+dp.y*dp.y);
}
};
// Andrew算法
struct AndrewConvexHull {
// some define.
static const int N = 1e5+5;
int stk[N],tp;
bool used[N];
void cal(Point* p, int n) {
// init
tp=0;
for(int i=1;i<=n;i++)used[i]=0;
//
sort(p+1,p+1+n,[](const Point& a,const Point& b) {
if (a.x!=b.x)return a.x<b.x;
return a.y<b.y;
});
// 点1不标记为used
// 等到计算上凸壳的时候会被自动标记上
stk[++tp]=1;
for(int i=2;i<=n;i++){
// dcmp<0会保留共线的点
// 如果改成dcmp<=0则不会保留共线的点
while(tp>1&&dcmp((p[stk[tp]]-p[stk[tp-1]])*(p[i]-p[stk[tp]]))<=0) {
used[stk[tp]]=0;
tp--;
}
used[i]=1;
stk[++tp]=i;
}
int ltp=tp;
for(int i=n-1;i>=1;i--){
if(used[i])continue;
while(tp>ltp&&dcmp((p[stk[tp]]-p[stk[tp-1]])*(p[i]-p[stk[tp]]))<=0) {
used[stk[tp]]=0;
tp--;
}
used[i]=1;
stk[++tp]=i;
}
// stk[1]=stk[tp]=1,这里去掉一个重复的p[1]
if(tp>1)tp--;
}
}T;
模板题:
P2742 [USACO5.1] 圈奶牛Fencing the Cows /【模板】二维凸包
解法:
最小边长即凸包的边长。
计算出凸包之后,将凸包上相邻两点的距离算出并加入答案即可。
Code:
#include <bits/stdc++.h>
using namespace std;
#define X first
#define Y second
#define int long long
#define PI pair<int, int>
const int maxm = 2e5 + 5;
// 依赖的计算几何模板
inline int dcmp(double x){
constexpr double eps=1e-6;
if (x<-eps)return -1;
if (x>eps)return 1;
return 0;
}
struct Point {
double x,y;
Point()=default;
Point(double x,double y):x(x),y(y){}
Point operator+(const Point& p){
return Point(x+p.x,y+p.y);
}
Point operator-(const Point& p){
return Point(x-p.x,y-p.y);
}
// dot mul.
Point operator*(double k){
return Point(x*k,y*k);
}
// cross mul.
double operator*(const Point& p){
return x*p.y-y*p.x;
}
double dist(const Point& p){
const Point dp=*this-p;
return sqrt(dp.x*dp.x+dp.y*dp.y);
}
};
// Andrew算法
struct AndrewConvexHull {
// some define.
static const int N = 1e5+5;
int stk[N],tp;
bool used[N];
void cal(Point* p, int n) {
// init
tp=0;
for(int i=1;i<=n;i++)used[i]=0;
//
sort(p+1,p+1+n,[](const Point& a,const Point& b) {
if (a.x!=b.x)return a.x<b.x;
return a.y<b.y;
});
// 点1不标记为used
// 等到计算上凸壳的时候会被自动标记上
stk[++tp]=1;
for(int i=2;i<=n;i++){
// dcmp<0会保留共线的点
// 如果改成dcmp<=0则不会保留共线的点
while(tp>1&&dcmp((p[stk[tp]]-p[stk[tp-1]])*(p[i]-p[stk[tp]]))<=0) {
used[stk[tp]]=0;
tp--;
}
used[i]=1;
stk[++tp]=i;
}
int ltp=tp;
for(int i=n-1;i>=1;i--){
if(used[i])continue;
while(tp>ltp&&dcmp((p[stk[tp]]-p[stk[tp-1]])*(p[i]-p[stk[tp]]))<=0) {
used[stk[tp]]=0;
tp--;
}
used[i]=1;
stk[++tp]=i;
}
// stk[1]=stk[tp]=1,这里去掉一个重复的p[1]
if(tp>1)tp--;
}
}T;
int n;
Point a[maxm];
void solve(){
cin>>n;
for(int i=1;i<=n;i++){
cin>>a[i].x>>a[i].y;
}
T.cal(a,n);
double ans=0;
for(int i=1;i<=T.tp;i++){
int now=T.stk[i];
int pre=(i>1?T.stk[i-1]:T.stk[T.tp]);
ans+=a[now].dist(a[pre]);
}
printf("%.2f\n",ans);
}
signed main() {
// #define MULTI_CASE
ios::sync_with_stdio(0);
cin.tie(0);
#ifndef ONLINE_JUDGE
freopen("../in.txt", "r", stdin);
freopen("../out.txt", "w", stdout);
#endif
#ifdef MULTI_CASE
int T;
cin >> T;
while (T--)
#endif
solve();
return 0;
}