【BZOJ1336】[Balkan2002]Alien最小圆覆盖 随机增量法

【BZOJ1336】[Balkan2002]Alien最小圆覆盖

Description

给出N个点,让你画一个最小的包含所有点的圆。

Input

先给出点的个数N,2<=N<=100000,再给出坐标Xi,Yi.(-10000.0<=xi,yi<=10000.0)

Output

输出圆的半径,及圆心的坐标

Sample Input

6
8.0 9.0
4.0 7.5
1.0 2.0
5.1 8.7
9.0 2.0
4.5 1.0

Sample Output

5.00
5.00 5.00

题解:特地学了一发随机增量法,期望复杂度据说是O(n),但不会证。

起初以1为圆心。先枚举一个点,如果这个点在圆外,则改为以哪个点为圆心;再枚举一个点,如果这个点再圆外,则圆心改为这两个点的中点;再枚举一个点,如果这个点再圆外,则将圆心改为三角形的内心。如何求内心?本人懒得推式子所以写的高斯消元。

由于是期望复杂度所以一开始要将原序列随机排序。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn=500010;
const double eps=1e-12;
struct point
{
	double x,y;
	point() {}
	point(double a,double b){x=a,y=b;}
}p[maxn],O;
int n;
double v[10][10],R;
point calc(point a,point b,point c)
{
	v[1][1]=2*b.x-2*a.x;
	v[1][2]=2*b.y-2*a.y;
	v[1][3]=b.x*b.x-a.x*a.x+b.y*b.y-a.y*a.y;
	v[2][1]=2*c.x-2*b.x;
	v[2][2]=2*c.y-2*b.y;
	v[2][3]=c.x*c.x-b.x*b.x+c.y*c.y-b.y*b.y;
	int i,j,k;
	for(i=1;i<=2;i++)
	{
		for(j=i+1;j<=2;j++)	if(fabs(v[j][i])>fabs(v[i][i]))	for(k=i;k<=3;k++)	swap(v[i][k],v[j][k]);
		double tmp=v[i][i];
		for(k=i;k<=3;k++)	v[i][k]/=tmp;
		for(j=1;j<=2;j++)	if(i!=j)	for(tmp=v[j][i],k=i;k<=3;k++)	v[j][k]-=tmp*v[i][k];
	}
	point ret(v[1][3],v[2][3]);
	return	ret;
}
double dis(point a,point b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int main()
{
	scanf("%d",&n);
	int i,j,k;
	for(i=1;i<=n;i++)	scanf("%lf%lf",&p[i].x,&p[i].y);
	random_shuffle(p+1,p+n+1);
	O=p[1],R=0;
	for(i=2;i<=n;i++)	if(dis(O,p[i])>R+eps)
	{
		O=p[i],R=0;
		for(j=1;j<i;j++)	if(dis(O,p[j])>R+eps)
		{
			O=point((p[i].x+p[j].x)/2,(p[i].y+p[j].y)/2),R=dis(O,p[i]);
			for(k=1;k<j;k++)	if(dis(O,p[k])>R+eps)
				O=calc(p[i],p[j],p[k]),R=dis(O,p[k]);
		}
	}
	printf("%.2lf\n%.2lf %.2lf",R,O.x,O.y);
	return 0;
}
时间: 2024-10-06 21:19:31

【BZOJ1336】[Balkan2002]Alien最小圆覆盖 随机增量法的相关文章

BZOJ 1336 Balkan2002 Alien最小圆覆盖 随机增量法

题目大意:求最小圆覆盖 随机增量法裸题 注意多输出几位小数不然过不去= = #include <cmath> #include <cstdio> #include <cstring> #include <iomanip> #include <iostream> #include <algorithm> #define M 100100 #define EPS 1e-7 using namespace std; struct Point

【bzoj1336/1337/2823】[Balkan2002]Alien最小圆覆盖 随机增量法

题目描述 给出N个点,让你画一个最小的包含所有点的圆. 输入 先给出点的个数N,2<=N<=100000,再给出坐标Xi,Yi.(-10000.0<=xi,yi<=10000.0) 输出 输出圆的半径,及圆心的坐标 样例输入 6 8.0 9.0 4.0 7.5 1.0 2.0 5.1 8.7 9.0 2.0 4.5 1.0 样例输出 5.00 5.00 5.00 题解 随机增量法求最小圆覆盖裸题 求法:设初始圆为某空圆,先枚举第一个点,如果不在当前圆内,则令当前圆为这一个点的最小圆

BZOJ 1337 最小圆覆盖 随机增量法

题意:链接 方法:随机增量法 解析: 首先把所有点打乱. 然后枚举第一个点,如果不在当前的圆内则把它设为圆心. 再枚举第二个点,如果不在当前的圆内则把圆设为其与第一个点的距离为直径的圆. 再枚举第三个点,如果不在当前的圆内则把圆设为这三个点构成三角形的外接圆. 然后最后就是答案了. 别问我这为什么是对的- -! 复杂度期望O(n),我是信了. 代码: #include <cmath> #include <cstdio> #include <iomanip> #inclu

[bzoj1336] [Balkan2002]Alien最小圆覆盖

最小圆覆盖..三个for是O(n)的QAQ..因为随机化后新的点不在当前圆内的几率不大.. 学习了下求中垂线的姿势... 1 #include<cstdio> 2 #include<cmath> 3 #include<iostream> 4 #include<cstdlib> 5 #include<algorithm> 6 #define d double 7 using namespace std; 8 const int maxn=1e5+2

hdu 3007【最小圆覆盖-随机增量法】

#include<iostream> #include<cstdio> #include<cmath> #include<algorithm> using namespace std; const int N=505; int n; double r; struct dian { double x,y; dian(double X=0,double Y=0) { x=X,y=Y; } dian operator + (const dian &a) c

BZOJ 1336 Balkan2002 Alien最小圆覆盖

题目大意:最小圆覆盖. 思路:再拍一份模板.做法见:http://blog.csdn.net/jiangyuze831/article/details/43950601 CODE: #define _CRT_SECURE_NO_WARNINGS #include <cmath> #include <cstdio> #include <cstring> #include <iomanip> #include <iostream> #include

bzoj2823: [AHOI2012]信号塔&amp;&amp;1336: [Balkan2002]Alien最小圆覆盖&amp;&amp;1337: 最小圆覆盖

首先我写了个凸包就溜了 这是最小圆覆盖问题,今晚学了一下 先随机化点,一个个加入 假设当前圆心为o,半径为r,加入的点为i 若i不在圆里面,令圆心为i,半径为0 再重新从1~i-1不停找j不在圆里面,令圆心为ij中点,直径为ij距离 再重新在1~j-1不停找k不在圆里面,三点可确定一圆,初中数学 复杂度看似O(n^3)实则O(n),好玄学 坑点:注意如果用点斜式表示方程有斜率为不存在的情况,需要特判 #include<cstdio> #include<iostream> #incl

【BZOJ-1336&amp;1337】Alie最小圆覆盖 最小圆覆盖(随机增量法)

1336: [Balkan2002]Alien最小圆覆盖 Time Limit: 1 Sec  Memory Limit: 162 MBSec  Special JudgeSubmit: 1573  Solved: 697[Submit][Status][Discuss] Description 给出N个点,让你画一个最小的包含所有点的圆. Input 先给出点的个数N,2<=N<=100000,再给出坐标Xi,Yi.(-10000.0<=xi,yi<=10000.0) Outpu

Bzoj 1336&amp;1337 Alien最小圆覆盖

1336: [Balkan2002]Alien最小圆覆盖 Time Limit: 1 Sec  Memory Limit: 162 MBSec  Special Judge Submit: 1473  Solved: 648 [Submit][Status][Discuss] Description Input 先给出点的个数N,2<=N<=100000,再给出坐标Xi,Yi.(-10000.0<=xi,yi<=10000.0) Output Sample Input 6 8.0