POJ 3932 Groundhog Build Home(最小圆覆盖)

转载请注明出处,谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526 by---cxlove

题目:用最小的圆覆盖所有的点

http://acm.hdu.edu.cn/showproblem.php?pid=3932

以下有两种方法。

首先是随机增量算法

------------------------------------------------------------------------------------

algorithm:

A、令Ci表示为前i个点的最小覆盖圆。当加入新点pi时如果pi不在Ci-1里那么pi必定在Ci的边界上。
B、再从新考虑这样一个问题,Ci为前i个点最小覆盖圆且p在Ci的的边界上!同理加入新点pi时如果p

i不在Ci-1里那么pi必定在Ci的边界上。这时我们就包含了两个点在这个最小圆的边界上。
C、再从新考虑这样一个问题,Ci为前i个点最小覆盖圆且有两个确定点再边界上!此时先让

O(N)的方法能够判定出最小圆。

------------------------------------------------------------------------------------

analysis:

现在来分析为什么是线性的。

C是线性的这是显然的。

B<-C的过程中。考虑pi 他在园内的概率为 (i-1)/i 。在圆外的概率为 1/i 所以加入pi的期望复杂度为:(1-i)/i*O(1) +(1/i)*O(i) {前者在园内那么不进入C,只用了O(1)。后者进入C用了O(i)的时间}这样分析出来,复杂度实际上仍旧

是线性的。

A<-B的过程中。考虑方法相同,这样A<-B仍旧是线性。于是难以置信的最小圆覆盖的复杂度变成了线性的。
-------------------------------------------------------------------------------------

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define LL long long
#define eps 1e-8
#define inf 999999.0
#define zero(a) fabs(a)0&&y+eps>0&&xR){
            central=p[i];R=0;
            for(int j=0;jR){
                    central.x=(p[i].x+p[j].x)/2;
                    central.y=(p[i].y+p[j].y)/2;
                    R=dist(central,p[j]);
                    for(int k=0;kR){
                            //3点确定圆
                            central=Circumcenter(p[i],p[j],p[k]);
                            R=dist(central,p[k]);
                        }
                }
        }
}
int main(){
    while(scanf("%lf%lf%d",&X,&Y,&n)!=EOF){
        for(int i=0;i

模拟退火算法。

随机选取若干个点,然后选定步长,从每个点随机走出去若干次,更新最优解。

将步长减小若干倍,直至达到一定的精度要求

好多地方随机,好多地方若干,这就和模拟退火的把握有关了。

在WA,TLE,AC之间徘徊,只刷到了800ms

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define LL long long
#define eps 1e-6
#define inf 999999.0
#define zero(a) fabs(a)0&&y+eps>0&&x0.001){
            for(int i=0;i





你可能感兴趣的:(POJ 3932 Groundhog Build Home(最小圆覆盖))