HDU 1045 Fire Net (二分图匹配)

题意:给出n*n的图,要在空白格内填小圆圈,圆圈两两之间不能在同一行、列,除非有黑格子挡住,求最多填多少个圆圈。

题解:二分图匹配
对每一行、列能填的连续白格子整体编号,分别用 l x [ i ] [ j ] 、 l y [ i ] [ j ] lx[i][j]、ly[i][j] lx[i][j]ly[i][j]表示。
对于点 ( i , j ) (i, j) (i,j),将 l x [ i ] [ j ] lx[i][j] lx[i][j] l y [ i ] [ j ] ly[i][j] ly[i][j]建边,表示这一行、一列,都只能有这一个点(圆圈),跑匈牙利算法求出的最大匹配就是能放置的最多圆圈个数。

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<fstream>
#include<set>
#include<map>
#include<sstream>
#include<iomanip>
#define ll long long
using namespace std;
/* ***********************************************************
//二分图匹配(匈牙利算法的 DFS 实现) (邻接矩阵形式)
//初始化: g[][] 两边顶点的划分情况
//建立 g[i][j] 表示 i->j 的有向边就可以了,是左边向右边的匹配
//g 没有边相连则初始化为 0
//uN 是匹配左边的顶点数, vN 是匹配右边的顶点数
//调用:res=hungary(); 输出最大匹配数
//优点:适用于稠密图, DFS 找增广路,实现简洁易于理解
//时间复杂度:O(VE)
//*************************************************************/
//顶点编号从 1 开始的
const int MAXN = 10;
int uN, vN;//u,v 的数目,使用前面必须赋值
int g[MAXN][MAXN];//邻接矩阵
int linker[MAXN];
bool used[MAXN];
bool dfs(int u) {
    for (int v = 1; v <= vN; v++)
        if (g[u][v] && !used[v]) {
            used[v] = true;
            if (linker[v] == -1 || dfs(linker[v])) {
                linker[v] = u;
                return true;
            }
        }
    return false;
}
int hungary() {
	int res = 0;
	memset(linker, -1, sizeof(linker));
	for (int u = 1; u <= uN; u++) {
		memset(used, false, sizeof(used));
		if (dfs(u)) res++;
	}
	return res;
}
int n, lx[11][11], ly[11][11];
char s[11][11];
int main() {
	while (~scanf("%d", &n) && n) {
		uN = vN = 0;
		memset(g, 0, sizeof(g));
		for (int i = 1; i <= n; i++) scanf("%s", s[i] + 1);
		for (int i = 1; i <= n; i++) {
			for (int j = 1; j <= n; j++) {
				if (s[i][j] == 'X') continue;
				if (j == 1 || s[i][j - 1] == 'X') lx[i][j] = ++uN;
				else lx[i][j] = lx[i][j - 1];
			}
		}
		for (int j = 1; j <= n; j++) {
			for (int i = 1; i <= n; i++) {
				if (s[i][j] == 'X') continue;
				if (i == 1 || s[i - 1][j] == 'X') ly[i][j] = ++vN;
				else ly[i][j] = ly[i - 1][j];
			}
		}
		for (int i = 1; i <= n; i++) {
			for (int j = 1; j <= n; j++) {
				if (s[i][j] == 'X') continue;;
				g[lx[i][j]][ly[i][j]] = 1;
			}
		}
		printf("%d\n", hungary());
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值