HDU6080(很水的计算几何+floyd)

题面
题意就是平面上有n个村庄和m个点,在m个点中选最少的点使得这些点构成的凸包能包含n个村庄。

对于看到计几就是弃的选手,去漫展时把其他题都和别人讨论了出来了,只有这题没想过。

分析:
在m个点中选取点A和点B,对于向量AB,若n个村庄都在其左侧,这条边才可能是最终凸包上的向量。
考虑所有可能是凸包上的向量,问题就变成了选取最少的向量,使得其构成一个环,就是floyd了。

具体是这样的,枚举向量AB,再枚举n个村庄,若满足n个村庄都在其左侧,A向B连一条权为1的边,然后floyd求出最小环。
(用叉积判断点在向量的哪侧,分不清左右也没关系,把向量写出来,按一个确定的顺序叉,若有一个叉积<0,则不满足)

ans=m-最小环

复杂度O(n^2*m)代码习惯不好的我被卡常了,加了些正常的优化才过了。

#include <iostream>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>

using namespace std;
#define mmst(a, b) memset(a, b, sizeof(a))
#define mmcp(a, b) memcpy(a, b, sizeof(b))
#define imin(a, b) (((a)<(b)) ? (a) : (b))

typedef long long LL;

const int N=505,oo=1e7;

int n,m;
int a[N],b[N];
int x[N],y[N];
int f[N][N];

int main()
{
    while(scanf("%d",&n)==1)
    {

    for(int i=0;i<N;i++)
    for(int j=0;j<N;j++)
    f[i][j]=oo;

    for(int i=1;i<=n;i++)
    scanf("%d%d",&x[i],&y[i]);

    scanf("%d",&m);
    for(int i=1;i<=m;i++)
    scanf("%d%d",&a[i],&b[i]);

    for(int i=1;i<=m;i++)
    for(int j=1;j<=m;j++)
    if(j!=i)
    {
        bool flag=1;
        for(int k=1;k<=n;k++)
        {
            int cj=(a[j]-a[i])*(y[k]-b[i])-(b[j]-b[i])*(x[k]-a[i]);
            if(cj>0)
            {
                flag=0;
                break;
            }
        }
        if(flag)
        f[i][j]=1;
    }

    for(int k=1;k<=m;k++)
    for(int i=1;i<=m;i++)
    if(f[i][k]!=oo)
    for(int j=1;j<=m;j++)
    f[i][j]=imin(f[i][j],f[i][k]+f[k][j]);

    int ans=m+1;
    for(int i=1;i<=m;i++)
    ans=imin(ans,f[i][i]);

    if(ans!=m+1)
    cout<<m-ans<<endl;
    else
    cout<<"ToT"<<endl;
    }

    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值