ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

【计算几何】【凸包周长】【思维】城墙

2021-09-12 17:03:32  阅读:208  来源: 互联网

标签:return 周长 point int top Point stk 城墙 凸包


here

主要是看边界的情况,因为他说最边上的距离不能超过L这就意味着拐角处外面的城墙应该是属于弧状的。
这样最后的结果就会是凸包周长加上 2 Π L 2ΠL 2ΠL
套一下板子就好了。

ac code

#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];

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()
{
    double l;
	cin>>n>>l;
	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];

//    for(int i=0;i<=top;i++)
//    {
//    	cout<<stk[i].x<<" "<<stk[i].y<<endl;
//	}

    double fans=2.0*l*acos(-1);
    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("%.0f\n",fans/2);
		return;
	}
    printf("%.0f\n",fans);
}
signed main() {
	solve();
	return 0;
}

标签:return,周长,point,int,top,Point,stk,城墙,凸包
来源: https://blog.csdn.net/peacesir/article/details/120252689

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有