uva11419

题目链接:http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=560&page=show_problem&problem=2414

题目大意:给一个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;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值