HDU4720（最小圆覆盖问题）

题意：给出四个点的坐标，现在判断第四个点是否在最小覆盖圆里面，如果在则输出Safe（注意包括边界）,否则输出Danger.

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
using namespace std;
const int maxx=1005;
const int inf=0x3f3f3f3f;
struct node{double x,y;
};
node b[maxx];
node o;
double R;
double sqr(double x){return x*x;
}
double dis(node a,node b){return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
}
double incircle(node a){if(dis(o,a)<=R){return true;}return false;
}
node solve(double a,double b,double c,double d,double e,double f){double y=(f*a-c*d)/(b*d-e*a);double x=(f*b-c*e)/(a*e-b*d);return ((node){x,y});
}
int main(){int n;int t;scanf("%d",&t);int count=1;while(t--){for(int i=1;i<=4;i++){scanf("%lf %lf",&b[i].x,&b[i].y);}random_shuffle(b+1,b+n+1);R=0;for(int i=1;i<4;i++){if(!incircle(b[i])){o.x=b[i].x;o.y=b[i].y;R=0;for(int j=1;j<i;j++){if(!incircle(b[j])){o.x=(b[i].x+b[j].x)/2;o.y=(b[i].y+b[j].y)/2;R=dis(o,b[i]);for(int k=1;k<j;k++){if(!incircle(b[k])){o=solve(b[i].x-b[j].x,b[i].y-b[j].y,(sqr(b[j].x)-sqr(b[i].x)+sqr(b[j].y)-sqr(b[i].y))/2,b[i].x-b[k].x,b[i].y-b[k].y,(sqr(b[k].x)-sqr(b[i].x)+sqr(b[k].y)-sqr(b[i].y))/2);R=dis(b[i],o);}}}}}}double x=(b[4].x-o.x);double y=(b[4].y-o.y);double dist=sqrt(x*x+y*y);cout<<"Case #"<<count++<<": ";if(dist<=R){cout<<"Danger"<<endl;}else{cout<<"Safe"<<endl;}} return 0;
}

方法二：

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<iostream>
using namespace std;
const double sep=1e-3;
const double base=0.98;
const int maxx=1e3+5;
struct node{double x,y;
}t[maxx];
int n;
double x,y;
double dis(node a,node b){double x=(a.x-b.x);double y=(a.y-b.y);return sqrt(x*x+y*y);
}
void solve(){double ans=10000.0,temp=100.0,maxlen;node mid;mid.x=mid.y=0;int s=0;while( temp>sep ){maxlen=0;    for( int i=1 ; i<4 ; i++ ){if( maxlen<dis(mid,t[i]) ){maxlen=dis(mid,t[i]);s=i;}}ans=min(ans,maxlen);mid.x+=(t[s].x-mid.x)/maxlen*temp;mid.y+=(t[s].y-mid.y)/maxlen*temp;temp*=base;}double xx=mid.x-t[4].x;double yy=mid.y-t[4].y;double dist=sqrt(xx*xx+yy*yy);if(dist>ans){cout<<"Safe"<<endl;}else{cout<<"Danger"<<endl;}
}int main(){int count=1;int tt;scanf("%d",&tt);while(tt--){for( int i=1 ; i<=4 ; i++ ){scanf("%lf %lf",&t[i].x,&t[i].y);}cout<<"Case #"<<count++<<": ";solve();}return 0;
}

HDU4720（最小圆覆盖问题）相关推荐

HDU2215（最小圆覆盖问题）
题意:就是求最小圆覆盖问题方法一: #include<iostream> #include<algorithm> #include<cmath> #include ...
【hdoj】3007 Buried memory 【计算几何--最小圆覆盖】
传送门:Buried memory 苍天饶过谁,第三次在hdoj上交计算几何的题了,没一次是AC的. ┭┮﹏┭┮都是模板题啊,我都是抄板子的啊,为什么会这样,我怎么这么菜. 题意: 求最小圆覆盖的 ...
hdu 2215(最小圆覆盖)
解题思路:最小圆覆盖,注意由于树也有半径,所以要加上树的半径0.5 #include<iostream> #include<cmath> #include<cstdio& ...
POJ 2069最小球覆盖 HDU3007最小圆覆盖【模拟淬火算法】
POJ 2069最小球覆盖 1.给定N个三维点,要求覆盖这些点的最小球半径: 2.采用模拟淬火算法,随机选取一个点作为初始解,然后不断向当前最远的点靠近: 3.这是一个不断调整的过程,对应模拟淬火算法 ...
ZOJ - 1450 Minimal Circle HDU - 3007 Buried memory 最小圆覆盖模板【随机函数】【增量法】
题意给N个点,求最小的圆将这N个点全部覆盖,输出圆心坐标和半径分析最小的圆肯定落在三个点上,因此暴力枚举圆上的三个点即可,点增量算法O(n ^ 3),加入随机化,平均复杂度可以降到O(n^2) ...
P1742 最小圆覆盖
P1742 最小圆覆盖题意: 给出N个点,让你画一个最小的包含所有点的圆. 题解: 先说结论: 最优解的圆一定是在以某两个点连线为直径的圆或者某三个点组成的三角形的外接圆初始化将某个圆心定为第 ...
ZOJ 1450 Minimal Circle 点集的最小圆覆盖
From: http://blog.csdn.net/zmx354/article/details/17076267 给定一个点集,求出能覆盖点集内所有点的半径最小的圆.包含点在圆上的情况.个人感觉算 ...
luoguP1742 最小圆覆盖
最小圆覆盖首先没错,我是个蒟蒻.luogu 流程圆 C; for(i=1 to n) {if(P[i] 不在 C 内) {C = {P[i], 0};for(j=1 to i-1) {if(P[ ...
随机增量法：bzoj 1336 bzoj 1337 最小圆覆盖
1337: 最小圆覆盖 Time Limit: 1 Sec Memory Limit: 64 MB Submit: 1170 Solved: 573 [Submit][Status][Discus ...
python画图：小圆覆盖大圆问题
小圆覆盖大圆 MATLAB实现:https://blog.csdn.net/qq_41845823/article/details/120324970 用至少多少个半径为r的小圆覆盖半径为R的大圆是 ...

HDU4720（最小圆覆盖问题）

HDU4720（最小圆覆盖问题）相关推荐

最新文章

热门文章