最小圆覆盖模ssss板

# include<cstdio>
# include<algorithm>
# include<cmath>
using namespace std;
const int maxn=510;struct point{  //存点double x,y;point(double _x=0,double _y=0){x=_x;y=_y;}friend point operator + (const point &a,const point &b){return point(a.x+b.x,a.y+b.y);}friend point operator - (const point &a,const point &b){return point(a.x-b.x,a.y-b.y);}friend double operator ^ (point a,point b){return a.x*b.y-a.y*b.x;}
}dots[maxn];double distance11(point a,point b){  // 求两点之间的距离 return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}struct circle{   double r;point centre;circle(point _centre=point(0,0),double _r=0){centre=_centre;r=_r;}
}ansc;
struct triangle{point p[3];
};double trianglearea(triangle t){  //求三角形面积return fabs((t.p[1]-t.p[0])^(t.p[2]-t.p[0]))/2.0;
}circle circumcircleoftriangle(const triangle &t){   //三点确定圆心和坐标circle tmp;double a,b,c,c1,c2;double xa,ya,xb,yb,xc,yc;a=distance11(t.p[0],t.p[1]);b=distance11(t.p[1],t.p[2]);c=distance11(t.p[2],t.p[0]);tmp.r=a*b*c/trianglearea(t)/4.0;xa=t.p[0].x;ya=t.p[0].y;xb=t.p[1].x;yb=t.p[1].y;xc=t.p[2].x;yc=t.p[2].y;c1=(xa*xa+ya*ya-xb*xb-yb*yb)/2.0;c2=(xa*xa+ya*ya-xc*xc-yc*yc)/2.0;tmp.centre.x=(c1*(ya-yc)-c2*(ya-yb))/((xa-xb)*(ya-yc)-(xa-xc)*(ya-yb));tmp.centre.y=(c1*(xa-xc)-c2*(xa-xb))/((ya-yb)*(xa-xc)-(ya-yc)*(xa-xb));return tmp;
}//最小圆覆盖算法
circle mincircle2(int u,triangle t){ if(u==0)return circle(point(0,0),-2.0);else if(u==1)return circle(t.p[0],0);else if(u==2)return circle(point((t.p[0].x+t.p[1].x)/2.0,(t.p[0].y+t.p[1].y)/2.0),distance11(t.p[0],t.p[1])/2.0);else if(u==3)return circumcircleoftriangle(t);
}
void mincircle(int n,int u,triangle t){ansc=mincircle2(u,t);  //更新最小圆if(u==3)return ;for(int i=1;i<=n;i++){if(distance11(dots[i],ansc.centre)>ansc.r){t.p[u]=dots[i];mincircle(i-1,u+1,t);  //更新最小圆包含的点point temp=dots[i];for(int j=i;j>=2;--j){//将点重新排列dots[j]=dots[j-1];}dots[1]=temp;}}
}int main(){int n;
while(~scanf("%d",&n)){if(n==0)break;for(int i=1;i<=n;i++){scanf("%lf%lf",&dots[i].x,&dots[i].y);}triangle t;mincircle(n,0,t);printf("%.2f %.2f %.2f\n",ansc.centre.x,ansc.centre.y,ansc.r);
}}

阿里云建站—为企业提供互联网“快”服务

2020年因为一场突如其来的疫情，不少企业受到了严重冲击，疫情的冲击力对传统“纯线

下”行业的危害最大，互联网女皇玛丽·米克（MaryMeeker）4月17日发布了著名的年度互

联网趋势报告，报告中指出：拥有强大的互联网线上线下融合能力的企业在疫情中的表现最好，

线上线下融合的趋势已经存在一段时间了，但是疫情让这种需求变得更加的迫切。

如果你的企业完全依附于传统的、纯线下的经验模式，那么在2020年你将“必死无疑”，

一场巨大的，前所未有的互联网革命已经到来！

阿里云建站为了助力各行各业复工复产，疫情期间“马不停蹄”为数以万计的企业快速完成

建站，为他们朝着“线上线下融合”或者“纯线上”的互联网经营模式迈进，打下了坚实的基础。

“云·速成美站”模板建站1天就能上线，不懂技术没关系，打字就能建网站。千套网站模

板免费提供，百元就能建官网，一价全包，无任何隐形消费。

“云·企业官网”定制建站1周就能上线，高端量身定制建站，千元建官网无需自己动手，

建站专家1对1网站策划及设计，专业省心之选。

疫情，是一场大浪淘沙，每一次危机背后都隐藏着机会，危机越大，机会也越大。大环境

已经改变，如果你不努力不去改变来迎合这个大环境，那你必将被淘汰。

阿里云助力企业建站，优惠多多，福利多多，详情请点击如下链接

红包+折扣，阿里云上云大礼包！ https://www.aliyun.com/minisite/goods?userCode=nvs9lx2m

最小圆覆盖模ssss板相关推荐

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 ...

最小圆覆盖模ssss板

最小圆覆盖模ssss板相关推荐

最新文章

热门文章