P2742【计算几何】【凸包周长double类型】[USACO5.1]圈奶牛Fencing the Cows /【模板】二维凸包

最新推荐文章于 2024-09-27 20:03:04 发布

while WA er

最新推荐文章于 2024-09-27 20:03:04 发布

阅读量119

点赞数

分类专栏：计算几何文章标签： c语言算法

本文链接：https://blog.csdn.net/peacesir/article/details/119480807

版权

计算几何专栏收录该内容

17 篇文章 0 订阅

订阅专栏

here

改了之前的模板，因为之前的模板我发现是int的读的时候。这里可以读入小数然后就出bug了

AC代码

#include <bits/stdc++.h>
#include <iostream>
#include<algorithm>
#include <cstdio>
#include <cmath>
#define re register int
#define int long long 
#define inf 0x3f
using namespace std;
const int N=1e5+5;
int n,top,ans;
struct Point{
	double x,y;
}point[N],stk[N];

inline double dis(Point a,Point b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

double cross(Point a,Point b,Point c){
  return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
}

void find_miny(){
   Point tmp=point[0];
   int flag=0;
   for(int i=1;i<n;i++){
     if(point[i].y<tmp.y||(point[i].y==tmp.y&&point[i].x<tmp.x)){
        tmp=point[i];
        flag=i;
     }
   }
   if(flag) swap(point[0],point[flag]);
}

bool cmp(Point a,Point b){
  double n=cross(point[0],a,b);
  if(n>0||(n==0&&dis(point[0],a)>dis(point[0],b))) return true;
  return false;
}

void Graham(){
  top=-1;
  stk[++top]=point[0]; stk[++top]=point[1];
  for(int i=2;i<n;i++){
    while(top&&cross(stk[top-1],stk[top],point[i])<0) top--;
    stk[++top]=point[i];
  }
}
void solve()
{
	cin>>n;
	top=0;
	memset(stk,0,sizeof(stk));
	memset(point,0,sizeof(point));
	for(int i=0;i<n;i++){
		cin>>point[i].x>>point[i].y;
    }
    
    find_miny();
    sort(point+1,point+n,cmp);
    Graham();
    stk[++top]=point[0];
    double fans=0;
    for(int i=0;i<top;i++)
    {
    	fans+=dis(stk[i],stk[i+1]);
	}
	if(n==1)
	{
		cout<<0<<endl;
		return;
	}
	if(top==2)
	{
		printf("%.2lf\n",fans/2);
		return;
	}
    printf("%.2lf\n",fans);
}
signed main() {
	solve();
	return 0;
}