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中所有的边(一个顶点可以覆盖与它相连的边)。 最小顶点覆盖=最大匹配*/