hdu 3360 National Treasures 最小顶点覆盖(最大匹配)

来源:互联网 发布:js 迭代object 编辑:程序博客网 时间:2024/06/06 13:17
#include <cstdio>#include <cstring>#include <cmath>#include <vector>#include <string>#include <algorithm>#include <iostream>using namespace std;const int maxn=2505;const int maxc=55;vector<int>e[maxn];int vis[maxn],pre[maxn];int map[maxc][maxc];int dir[12][2]={{-1,-2},{-2,-1},{-2,1},{-1,2},{1,2},{2,1},{2,-1},{1,-2},{-1,0},{0,1},{1,0},{0,-1}};int find(int u)//判断增广路是否存在,匈牙利算法{    int i,j,v;    for(i=0;i<e[u].size();i++)    {        v=e[u][i];        if(!vis[v])        {            vis[v]=1;            if(pre[v]==-1||find(pre[v]))            {                pre[v]=u;                return 1;            }        }    }    return 0;}string a[maxn],b[maxn];int main(){    int n,m,tt=0;    while(scanf("%d%d",&n,&m)!=EOF)    {        if(n==0&&m==0)break;        int i,j,k;        for(i=1;i<=n*m;i++)            e[i].clear();        for(i=1;i<=n;i++)            for(j=1;j<=m;j++)            cin>>map[i][j];        int x,y,ans=0,t;        for(i=1;i<=n;i++)        {            for(j=1;j<=m;j++)            {                t=(i+j)%2;                if(map[i][j]==-1)continue;                for(k=0;k<12;k++)                {                    if(!(map[i][j]&(1<<k)))continue;                    x=i+dir[k][0];                    y=j+dir[k][1];                    if(x<=0||y<=0||x>n||y>m||map[x][y]==-1)continue;                    if(t)e[(i-1)*m+j].push_back((x-1)*m+y);                    else e[(x-1)*m+y].push_back((i-1)*m+j);                }            }        }        memset(pre,-1,sizeof(pre));        for(i=1;i<=n*m;i++)        {            memset(vis,0,sizeof(vis));            ans+=find(i);        }        printf("%d. %d\n",++tt,ans);    }    return 0;}/*    题意:在一个n*m的格子中,每个格子有一个数值,-1表示空,其余表示财宝。每个财宝的数值转换成二进制数,12个二进制位上数值,从右到左,第i个位是1表示图上相应第i序号位置需要有警卫。所有的要求位置有警卫财宝才安全。财宝可以被警卫替换。问至少需要替换多少财宝才能保证所有财宝的安全。    方法:需要警戒位置是财宝的讯号对财宝位置讯号建边。由于警戒位置与财宝位置的横纵坐标奇偶相反,可以建得二分图。对于所建图,根据题意就是找出最少的顶点使得剩余顶点覆盖所有的边,即最小顶点覆盖数为答案。    最大匹配:二分图G中,找出边数最大的子图M,使得M中各条边均无公共顶点,则M为最大匹配。可用匈牙利算法求得。    最小顶点覆盖:二分图G中,找出顶点数最少的子图M,使得M中所有的点可以覆盖G中所有的边(一个顶点可以覆盖与它相连的边)。    最小顶点覆盖=最大匹配*/