首页 > 代码库 > POJ 1113
POJ 1113
求凸包,然后加上一个半径为L的圆周就好了。
因为每到一个凸包的顶点时,需要保持最小间距为L,则多个顶点围起来会有360度。
在输出时,不要把结果直接转化主INT,因为这个不是四舍五入。
#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <cmath>using namespace std;const int MAXN=1050;int n,l;int st[MAXN],stop,cnt;int ans[MAXN];struct point{ int x,y;}p[MAXN];bool cmp(point A,point B){ if(A.y<B.y) return true; else if(A.y==B.y){ if(A.x<B.x) return true; } return false;}bool multi(point a,point b,point c){ return (a.x-c.x)*(b.y-c.y)-(a.y-c.y)*(b.x-c.x)>0;}void slove(){ stop=0; cnt=0; st[stop++]=0; st[stop++]=1; for(int i=2;i<n;i++){ while(stop>1&&multi(p[st[stop-1]],p[i],p[st[stop-2]])) stop--; st[stop++]=i; } for(int i=0;i<stop;i++){ ans[cnt++]=st[i]; } stop=0; st[stop++]=n-1; st[stop++]=n-2; for(int i=n-3;i>=0;i--){ while(stop>1&&multi(p[st[stop-1]],p[i],p[st[stop-2]])) stop--; st[stop++]=i; } for(int i=0;i<stop;i++) ans[cnt++]=st[i];}double dist(point a,point b){ return sqrt((a.x-b.x)*(a.x-b.x)*1.0+(a.y-b.y)*(a.y-b.y)*1.0);}int main(){ while(scanf("%d%d",&n,&l)!=EOF){ for(int i=0;i<n;i++){ scanf("%d%d",&p[i].x,&p[i].y); } sort(p,p+n,cmp); slove(); double anl=4*acos(0.0)*l; for(int i=1;i<cnt;i++){ anl+=dist(p[ans[i-1]],p[ans[i]]); } printf("%.0lf\n",anl); } return 0;}
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。