题目大意:给一个R*C大小的网格,网格上放了一些目标。可以在网格外发射子弹,子弹会沿着垂直或者水平的方向飞,并且打掉飞行路径上的所有目标,计算出至少需要多少子弹,各从哪些位置发射,才能把所有目标打掉。输入为r,c,n(目标的个数),下面n行是目标的坐标,输出最小值。
把行看成X集合,把列看成Y集合,那么点就是连接X与Y集合的一条边,现在相当于用最少的点来覆盖所有的边,即最小覆盖,最小覆盖=最大匹配,用匈牙利算法解决。
这题还要求打印出解集,那么我们在计算出最大匹配后,可以从X集合中未匹配的点开始扩展匈牙利树,在Y中被扩展到的点肯定得选,因为只有选择Y中这些点,X中未匹配的点才能被覆盖,而在X中,我们要选那些在扩展中没有被访问过的点,因为这些点是没有与Y中的那些扩展到的点相匹配,所以得选上,最终,打印出Y中被访问的点以及X中未被访问的点。
//匈牙利算法解决二分图的最小覆盖
/*
* ThinkingLion.cpp
*
* Created on: 2014年11月14日
* Author: dell
*/
#include<stdio.h>
#include<iostream>
#include<string>
#include<string.h>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<time.h>
#include<queue>
#include<stack>
#include<iterator>
#include<math.h>
#include<stdlib.h>
#include<limits.h>
//#define ONLINE_JUDGE
#define eps 1e-8
#define INF 0x7fffffff
#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define sfs(a) scanf("%s",a)
#define sf(a) scanf("%d",&a)
#define sfI(a) scanf("%I64d",&a)
#define pf(a) printf("%d\n",a)
#define pfs(a) printf("%s\n",a)
#define pfI(a) printf("%I64d\n",a)
#define sfd(a,b) scanf("%d%d",&a,&b)
#define for1(i,a,b) for(int i=(a);i<b;i++)
#define for2(i,a,b) for(int i=(a);i<=b;i++)
#define clr1(a) memset(a,0,sizeof(a))
#define clr2(a) memset(a,-1,sizeof(a))
#define ll __int64
const double PI = acos(-1.0);
template<class T> T gcd(T a, T b) {
return b ? gcd(b, a % b) : a;
}
template<class T> T lcm(T a, T b) {
return a / gcd(a, b) * b;
}
template<class T> inline T Min(T a, T b) {
return a < b ? a : b;
}
template<class T> inline T Max(T a, T b) {
return a > b ? a : b;
}
using namespace std;
#define N 1100
int r, c, n;
struct BPM{ //匈牙利算法解决二分图最大匹配(二分图有左右两个集合,分别为X,Y)
int n,m; //n代表行数,m代表列数
vector<int> G[N]; //邻接表存储边
int left[N]; //left[i]存储在Y集合中某个元素i在X集合中匹配的元素
int right[N]; //right[i]存储在X集合中某个元素i在Y集合中匹配的元素
bool T[N]; //标记X集合中的某个元素是否被访问
bool S[N]; //标记Y集合中的某个元素是否被访问
void Init(int n,int m){
this->n=n;
this->m=m;
for(int i=0;i<N;i++) G[i].clear();
}
void addEdge(int u,int v){
G[u].push_back(v);
}
bool match(int u){ //二分图最大匹配
S[u]=true;
for(int i=0;i<(int)G[u].size();i++){
int v=G[u][i];
if(!T[v]){
T[v]=true;
if(left[v]==-1 || match(left[v])){
left[v]=u;
right[u]=v;
return true;
}
}
}
return false;
}
int solve(){ //返回最大匹配数
memset(left,-1,sizeof left);
memset(right,-1,sizeof right);
int ans=0;
for(int i=0;i<n;i++){
memset(S,0,sizeof S);
memset(T,0,sizeof T);
if(match(i)) ans++;
}
return ans;
}
int mincover(vector<int> & X,vector<int>& Y){
int ans=solve();
memset(S,0,sizeof S);
memset(T,0,sizeof T);
for(int i=0;i<n;i++) //从X集合中选择还未匹配的点进行扩展
if(right[i] == -1)
match(i);
for(int i=0;i<n;i++) //X记录在X集合中选择哪些元素()
if(!S[i]) X.push_back(i);
for(int v=0;v<m;v++) //Y记录在Y集合中选择哪些元素
if(T[v]) Y.push_back(v);
return ans;
}
};
BPM solver;
int main() {
#ifndef ONLINE_JUDGE
freopen("in.txt", "r", stdin);
// freopen("out.txt","w",stdout);
#endif
while (scanf("%d%d%d", &r, &c, &n) != EOF && r + c + n) {
int x, y;
solver.Init(r,c);
vector<int>v1;
vector<int> v2;
for (int i = 0; i < n; i++) {
scanf("%d%d", &x, &y);
x--;y--;
solver.addEdge(x,y);
}
int ans=solver.mincover(v1,v2);
printf("%d",ans);
for(int i=0;i<(int)v1.size();i++) printf(" r%d",v1[i]+1);
for(int i=0;i<(int)v2.size();i++) printf(" c%d",v2[i]+1);
printf("\n");
// printf("%d\n", ans);
}
return 0;
}