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

来源:互联网 发布:ei数据库免费检索入口 编辑:程序博客网 时间:2024/06/01 12:57

题面
题意就是平面上有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;}