week16 CSP-M4-B ZJM要抵御宇宙射线

ZJM要抵御宇宙射线

据传，2020年是宇宙射线集中爆发的一年，这和神秘的宇宙狗脱不了干系！但是瑞神和东东忙于正面对决宇宙狗，宇宙射线的抵御工作就落到了ZJM的身上。假设宇宙射线的发射点位于一个平面，ZJM已经通过特殊手段获取了所有宇宙射线的发射点，他们的坐标都是整数。而ZJM要构造一个保护罩，这个保护罩是一个圆形，中心位于一个宇宙射线的发射点上。同时，因为大部分经费都拨给了瑞神，所以ZJM要节省经费，做一个最小面积的保护罩。当ZJM决定好之后，东东来找ZJM一起对抗宇宙狗去了，所以ZJM把问题扔给了你~

Input

输入 第一行一个正整数N，表示宇宙射线发射点的个数
接下来N行，每行两个整数X,Y，表示宇宙射线发射点的位置

Output

输出包括两行
第一行输出保护罩的中心坐标x,y 用空格隔开
第二行输出保护罩半径的平方（所有输出保留两位小数，如有多解，输出x较小的点，如扔有多解，输入y较小的点）
无行末空格

Sample Input

5
0 0
0 1
1 0
0 -1
-1 0

Sample Output

0.00 0.00
1.00

我的思路：
这道题是要在给出的众多点中，选择一个点作为圆心，设置一定的半径，再以此作圆，要求这个圆要包含所有点。我的想法是可以遍历每一个点，并以这个点为圆心，再去找距离这个圆心最远的点，以此为半径；在找到以当前点为圆心的最大的半径后，再将此半径与之前记录的半径作比较，如果此时的半径更小则将结果更新为当前情况。遍历结束之后，也就得到最后的正确结果了。

我的总结：
注意是保护罩半径的“平方”！在比赛时，我直接给的半径，结果惨不忍睹。。。

我的代码：

#include<iostream>
#include<cmath>
#include <iomanip>
#include<algorithm>
using namespace std;

struct yuan{
long double x,y;
bool operator < (const yuan& aa) const
{
	return x==aa.x?(y<aa.y):(x<aa.x);
}
}yy[1005];

int n;
long double R=1e11,xa,ya,rr,rt;
int main()
{
	cin>>n;
	for(int i=0;i<n;i++) cin>>yy[i].x>>yy[i].y;
	sort(yy,yy+n);
	for(int i=0;i<n;i++)
	{
		rr=0;
		for(int j=0;j<n;j++)
		{
			if(j==i) continue;
			rt = pow(yy[i].x-yy[j].x,2)+pow(yy[i].y-yy[j].y,2);
			if(rt>rr) rr=rt;
		}
		if(rr<R)
		{
			R=rr;
			xa=yy[i].x;
			ya=yy[i].y;
		}
	}
	cout<<setiosflags(ios::fixed)<<setprecision(2)<<xa<<" "<<ya<<endl<<R;
	return 0;
}