HUD 3007 Buried memory

最新推荐文章于 2018-02-24 20:01:20 发布

Cxiaokai

最新推荐文章于 2018-02-24 20:01:20 发布

阅读量998

点赞数

分类专栏：计算几何文章标签： ini struct c

本文链接：https://blog.csdn.net/Cxiaokai/article/details/6664595

版权

计算几何专栏收录该内容

20 篇文章 0 订阅

订阅专栏

本文介绍了如何使用最小的圆覆盖平面点集，首先通过求解点集的凸包，然后枚举凸包上的点找到最大距离。虽然旋转卡壳算法也可用于求最大距离，但在存在三点共线的情况下，可能需要额外计算以确保覆盖所有点。

摘要由CSDN通过智能技术生成

用最小的圆去覆盖平面点集。

思路：对平面的点集求凸包。枚举凸包上的点，找到最大距离。现在想想也可以用旋转卡壳来求最大距离。以最大距离的一半做为圆的半径。但是还有一种情况，就是有三个点的时候，很可能求出两点间的最大距离，以这个距离做为半径画出的圆，不能包围第三个点。这时还需要求一次距离。

求凸包+枚举凸包上的点

#include<stdio.h>
#include<math.h>
#include<iostream>
#include<algorithm>
using namespace std;
#define N 502
#define EPS 1e-8
struct point
{
    double x;
    double y;
}p[N];
int s[N],top;
int dblcmp(double d)
{
    if(fabs(d)<EPS)return 0;
    return d>0?1:-1;
}
double dis(point a,point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double det(double x1,double y1,double x2,double y2)
{
    return x1*y2-x2*y1;
}
double cross(point a,point b,point c)
{
    return det(a.x-c.x,a.y-c.y,b.x-c.x,b.y-c.y);
}
bool cmp1(point a,point b)
{
    if(a.y==b.y)return a.x<b.x;
    return a.y<b.y;
}
bool cmp2(point a,point b)
{
    double d=cross(a,b,p[0]);
    if(d>0)return true;
    else if(d<0)return false;
    {
        double dis1=dis(a,p[0]);
        double dis2=dis(b,p[0]);
        if(dis1<dis2)return true;
        else return false;
    }
}
void graham_scan(int n)
{
    sort(p,p+n,cmp1);
    sort(p+1,p+n,cmp2);
    top=-1;
    s[++top]=0;
    s[++top]=1;
    s[++top]=2;
    for(int i=3;i<n;i++)
    {
        while(top>=1&&dblcmp(cross(p[s[top]],p[i],p[s[top-1]]))<0)
        {
            top--;
        }
        s[++top]=i;
    }
}
int main()
{
    int n;
    while(scanf("%d",&n)&&n)
    {
        for(int i=0;i<n;i++)
        {
            scanf("%lf%lf",&p[i].x,&p[i].y);
        }
        graham_scan(n);
        double dd=-1;
        point center;
        for(int i=0;i<=top;i++)
        {
            for(int j=i+1;j<=top;j++)
            {
                if(dd<dis(p[s[i]],p[s[j]]))
                {
                    dd=dis(p[s[i]],p[s[j]]);
                    center.x=(p[s[i]].x+p[s[j]].x)/2;
                    center.y=(p[s[i]].y+p[s[j]].y)/2;
                }
            }
        }
        dd/=2;
        for(int i=0;i<=top;i++)
        {
            if(dd<dis(p[s[i]],center))
            {
                dd=dis(p[s[i]],center);
            }
        }
        printf("%.2lf %.2lf %.2lf\n",center.x,center.y,dd);
    }
    return 0;
}

直接枚举所有的点

#include<stdio.h>
#include<iostream>
#include<math.h>
using namespace std;
#define N 502
struct point
{
    double x;
    double y;
}p[N];
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()
{
    int n;
    while(scanf("%d",&n)&&n)
    {
        for(int i=0;i<n;i++)
        {
            scanf("%lf%lf",&p[i].x,&p[i].y);
        }
        double dd=-1;
        point center;
        for(int i=0;i<n;i++)
        {
            for(int j=i+1;j<n;j++)
            {
                if(dd<dis(p[i],p[j]))
                {
                    dd=dis(p[i],p[j]);
                    center.x=(p[i].x+p[j].x)/2;
                    center.y=(p[i].y+p[j].y)/2;
                }
            }
        }
        dd/=2;
        for(int i=0;i<n;i++)
        {
            if(dd<dis(center,p[i]))
            {
                dd=dis(center,p[i]);
            }
        }
        printf("%.2lf %.2lf %.2lf\n",center.x,center.y,dd);
    }
    return 0;
}