| Time Limit: 1000MS | Memory Limit: 65536K | |
| Total Submissions: 8650 | Accepted: 4285 |
Description
Obviously, it is desirable to use as few antennas as possible, but still provide coverage for each place of interest. We model the problem as follows: Let A be a rectangular matrix describing the surface of Sweden, where an entry of A either is a point of interest, which must be covered by at least one antenna, or empty space. Antennas can only be positioned at an entry in A. When an antenna is placed at row r and column c, this entry is considered covered, but also one of the neighbouring entries (c+1,r),(c,r+1),(c-1,r), or (c,r-1), is covered depending on the type chosen for this particular antenna. What is the least number of antennas for which there exists a placement in A such that all points of interest are covered?
Input
Output
Sample Input
2 7 9 ooo**oooo **oo*ooo* o*oo**o** ooooooooo *******oo o*o*oo*oo *******oo 10 1 * * * o * * * * * *
Sample Output
17 5 ///到现在还不是很理解如何 建的图 为啥要这样 先用网上的代码吧......,慢慢理解题目大意：
一个矩形中，有N个城市’*’，现在这n个城市都要覆盖无线，若放置一个基站，那么它至多可以覆盖相邻的两个城市。 问至少放置多少个基站才能使得所有的城市都覆盖无线？
解题思路：
思前想后，依稀可以认为是一道求二分图的最小路径覆盖问题
（注意不是最小点覆盖）
那么接下来需要确认的是，
究竟是求 有向二分图的最小路覆盖，还是求 无向二分图的最小路覆盖
因为有向和无向是截然不同的计算方法。
要确认是构造有向图，还是构造无向图，那么就需要先根据题意，看看构造二分图时所使用的方式，更适合构造哪一种二分图。
然后就进入了本题难点：如何构造二分图
首先要明确的是，输入的一堆“圈圈星星”可以看做是一张大地图，地图上有所有城市的坐标，但是这里有一个误区：不能简单地把城市的两个x、y坐标作为准备构造的二分图的两个顶点集。
城市才是要构造的二分图的顶点！
构造方法如下：
例如输入：
*oo
***
O*o
时，可以抽象为一个数字地图：
100
234
050
数字就是根据输入的城市次序作为该城市的编号，0代表该位置没有城市。
然后根据题目的“范围”规则，从第一个城市开始，以自身作为中心城市，向四个方向的城市进行连线（覆盖）
因此就能够得到边集：
e12 e21 e32 e43 e53
e23 e34
e35
可以看到，这些边都是有向边，但是每一条边都有与其对应的一条相反边。
即任意两个城市（顶点）之间的边是成对出现的
那么我们就可以确定下来，应该 构造无向二分图（其实无向=双向）
因为若要构造有向的二分图时，需要判断已出现的边，是很麻烦的工作
为了把有向图G构造为无向二分图，这里需要引入一个新名词“拆点”
其实就是把原有向图G的每一个顶点都”拆分（我认为复制更准确）”为2个点，分别属于所要构造的二分图的两个顶点集
例如在刚才的例子中抽出一条有向边e12举例说明：
复制顶点1和顶点2，使得1，2∈V1； 1’，2’∈V2 ，不难发现|V1|=|V2|
根据边e12和e21，得到无向二分图：
那么同理就可以得到刚才的例子的 无向二分图为：
再继而通过无向二分图，以V1的元素作为row，V2的元素作为col，构造 可达矩阵 存储到计算机
1’ 2’ 3’ 4’ 5’
1 F T F F F
2 T F T F F
3 F T F T T
4 F F T F F
5 F F T F F
接下来就是要求这个 无向二分图的最小路径覆盖 了
利用公式：
无向二分图的最小路径覆盖 = 顶点数 – 最大二分匹配数/2
顶点数：就是用于构造无向二分图的城市数，即进行“拆点”操作前的顶点数量
最大二分匹配书之所以要除以2，是因为进行了“拆点”擦奥做做使得匹配总数多了一倍，因此除以2得到原图的真正的匹配数
最后剩下的问题就是求最大二分匹配数了，用匈牙利算法，这就不多说了，参考POJ3041的做法，基本一摸一样。
从这道题得出了一个结论：
当二分图的两个顶点子集基数相等时，该二分图所有顶点的匹配数 等于 任意一个顶点子集匹配数的2倍
其实匈牙利算法解题是极为简单的，但是图论的难并不是难在解答，而是建图的过程，也难怪会有牛曰：用匈牙利算法，建图是痛苦的，最后是快乐的。
/
正式的定义，网上一大把，但他们的作用是让人看不懂……
二分图：把点分成两个集合X,Y，使得图的边的两个端点总是分别落在X和Y上，不会有X中的点连向X中的点，不会有Y中的点连向Y中的点
匹配：实质上是二分图中的一个边集，边集中出现的点不会重合，比如有a-b了，就不会有a-c了，要是有了a就重合了
最大匹配：这个边集的数目最大的那个匹配
匈牙利算法——
增广路：一条在X和Y之间交错的路径，【这条路上一条是匹配边，一条不是匹配边】，如此相交错，
其中第一条和最后一条边不是匹配边，（所以增广路的长度一定为奇数，不是匹配边的数目比是匹配边的数目多1），
按matrix67的神说法，当我们把这条路上不是匹配边的那些换成要匹配的，原来是匹配的换成不要匹配的，匹配数就+1
所以当有增广路存在时说明匹配数可以再增大
二分图中最小点覆盖等于最大匹配数
最小点覆盖：实质是个点集，点集里面的点能覆盖所有的边，最小点覆盖就是满足这个要求的点集中点数最小的那个
证明：所有的边分为匹配的（A）边和不是匹配的边（B），最小点覆盖的点集就是要每条匹配的边两端顶点中的一个，
比如现在有x1-y1属于A，x1-y2属于B，对于x1-y1这条匹配边取x1而不取y1，这样就能覆盖到x1-y2，即B也能覆盖到
二分图中最小边覆盖=顶点数-最小点覆盖（最大匹配）
最小边覆盖：实质是个边集，这个集合里的边能覆盖所有的点，最小边覆盖是满足这个要求的所有边集中边数最少的一个
这里顶点数等于总的顶点数，是二分图两边的顶点数，不是一边
证明：设最大匹配数为m，总顶点数为n。为了使边数最少，又因为一条边最多能干掉两个点，所以尽量用边干掉两个点
也就是取有匹配的那些边，当然这些边是越多越好，那就是最大匹配了，所以先用最大匹配数目的边干掉大多数点剩下的解决没有被匹配的点，就只能一条边干掉一个点了，设这些数目为a
显然，2m+a=n，而最小边覆盖=m+a，
所以最小边覆盖=（2m+a）-m=n-m
//自己照网上给的思路敲得 感觉好不规范.... #include<stdio.h> #include<string.h> int map[100][100]; int map1[100][100]; char str[100][100]; int book[100]; int match[100]; int dx[4]={1,0,-1,0}; int dy[4]={0,-1,0,1}; int v1,v2; int dfs(int x) { for(int y=1;y<=v2;y++) { if(map1[x][y]&&!book[y]) { book[y]=1; if(match[y]==0||dfs(match[y])) { match[y]=x; return 1; } } } return 0; } int main() { int t,m,n; int i,j; scanf("%d",&t); while(t--) { memset(match,0,sizeof(match)); memset(map,0,sizeof(map)); memset(map1,0,sizeof(map1)); scanf("%d%d",&n,&m); for(i=0;i<n;i++) scanf("%s",str[i]); int num=0; for(i=0;i<n;i++) { for(j=0;j<m;j++) { if(str[i][j]=='*') { map[i][j]=++num; } } } // printf("---%d\n",num); for(i=0;i<n;i++) for(j=0;j<m;j++) { if(map[i][j]) { for(int k=0;k<4;k++) { int xx=i+dx[k]; int yy=j+dy[k]; if(map[xx][yy]&&xx>=0&&yy>=0&&xx<n&&yy<m) { int a=map[i][j]; int b=map[xx][yy]; // printf("----%d %d %d %d %d %d\n",i,j,a,b); map1[a][b]=1; } } } } v1=v2=num; int sum=0; for(i=1;i<=v1;i++) { memset(book,0,sizeof(book)); if(dfs(i)) ++sum; } //for(i=1;i<=v1;i++) //{ // for(j=1;j<=v1-1;j++) // { // printf("%d ",map1[i][j]); // } // printf("%d\n",map1[i][j]); // } //printf("%d %d\n",v1,sum); printf("%d\n",v1-sum/2); } } //网上大神代码<strong> 優YoU <a target=_blank href="http://user.qzone.qq.com/289065406/blog/1299322779" target="_blank">http://user.qzone.qq.com/289065406/blog/1299322779</a></strong>//Memory Time //420K 16MS #include<iostream> using namespace std; int ipmap[41][11]; //标记存在城市'*'的位置，并依次记录城市的编号 int ip; //城市编号（最终是城市数量） int V1,V2; //二分图的两个顶点集 int M; //最大二分匹配 bool city[401][401]; //标记两个城市之间是否能连通 //通过“拆点”操作，把每一个城市拆分为2个，分别属于所构造的二分图的两个点集 bool vist[401]; int link[401]; int dire_r[4]={-1,1,0,0}; int dire_c[4]={0,0,-1,1}; //分别对应四个方位 上 下 左 右 /*Hungary Algorithm*/ bool dfs(int x) { for(int y=1;y<=V2;y++) if(city[x][y] && !vist[y]) { vist[y]=true; if(link[y]==0 || dfs(link[y])) { link[y]=x; return true; } } return false; } void search(void) { for(int x=1;x<=V1;x++) { memset(vist,false,sizeof(vist)); if(dfs(x)) M++; } return; } int main(void) { int test,h,w; cin>>test; while(test--) { /*Initial*/ memset(ipmap,0,sizeof(ipmap)); memset(city,false,sizeof(city)); memset(link,0,sizeof(link)); ip=0; M=0; /*Read in the maps*/ cin>>h>>w; int i,j; char temp; for(i=1;i<=h;i++) for(j=1;j<=w;j++) { cin>>temp; if(temp=='*') ipmap[i][j]=++ip; } /*Structure the Bipartite Graphs*/ for(i=1;i<=h;i++) for(j=1;j<=w;j++) if(ipmap[i][j]) for(int k=0;k<4;k++) { int x=i+dire_r[k]; int y=j+dire_c[k]; if(ipmap[x][y]) city[ ipmap[i][j] ][ ipmap[x][y] ]=true; //"拆点"操作是"顺便"被完成的 } //二分图构造完毕后，之后的问题就和POJ3041一样处理了 V1=V2=ip; /*增广轨搜索*/ search(); /*Output*/ cout<<ip-M/2<<endl; //无向二分图：最小路径覆盖数 = "拆点"前原图的顶点数 - 最大匹配数/2 } return 0; }
118
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
