Gym-101606B_计算几何

最新推荐文章于 2019-10-04 19:21:56 发布

ujn20161222

最新推荐文章于 2019-10-04 19:21:56 发布

阅读量220

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/ujn20161222/article/details/80149835

计算几何专栏收录该内容

2 篇文章

订阅专栏

本文通过C++实现了一个计算几何问题的解决方案，具体为：在平面上给定点集，找出一条线段使得所有其它点到该线段的距离的最大值与最小值之差最小。文章使用了枚举的方法来寻找最优线段，并利用向量运算和点到线的距离公式进行计算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

https://cn.vjudge.net/problem/Gym-101606B

枚举一个连线，算距离

#include<bits/stdc++.h>
using namespace std;

#define ll long long

ll rd(){
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

const int N=100+10;

struct P{
    int x,y;
    P operator-(P b){
        return {x-b.x,y-b.y};
    }
    double len(){return hypot(x,y);}
}p[N];

double cross(P a,P b){
    return a.x*b.y-a.y*b.x;
}

double dist(P p,P a,P b){
    return cross(p-a,b-a)/(b-a).len();
}

int main(){
    //freopen("in.txt","r",stdin);
    int n=rd();
    for(int i=0;i<n;i++)
    p[i].x=rd(),p[i].y=rd();
    double ans=1e100;
    for(int i=0;i<n;i++)
    for(int j=0;j<i;j++){
      //  cout<<"----"<<endl;
        double a=1e100,b=1e-100;
        for(int k=0;k<n;k++){
            a=min(a,dist(p[k],p[i],p[j]));
            b=max(b,dist(p[k],p[i],p[j]));
          //  cout<<a<<" "<<b<<" "<<k<<endl;
        }
        ans=min(ans,b-a);
        //cout<<i<<" "<<j<<endl;
        //cout<<a<<" "<<b<<" "<<ans<<endl;
    }
    printf("%.6f",ans);
    //cout<<ans<<endl;
}