最小圆覆盖模ssss板-优快云博客

本文链接：https://blog.youkuaiyun.com/weixin_44075041/article/details/99571859

# 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周就能上线，高端量身定制建站，千元建官网无需自己动手，