#二分图匹配，floyd，匈牙利算法#POJ 2594 Treasure Exploration

最新推荐文章于 2021-10-20 17:16:32 发布

lemondinosaur

最新推荐文章于 2021-10-20 17:16:32 发布

阅读量222

点赞数

分类专栏：二分图匹配最短路径最大流文章标签： POJ2594 Treasure Exploration

本文链接：https://blog.csdn.net/sugar_free_mint/article/details/103534155

版权

最短路径同时被 3 个专栏收录

58 篇文章 0 订阅

订阅专栏

最大流

20 篇文章 0 订阅

订阅专栏

二分图匹配

18 篇文章 0 订阅

订阅专栏

该博客介绍了如何利用Floyd算法找出无向图的传递闭包，将问题转化为最小不可相交路径覆盖，接着通过匈牙利算法求解最大匹配，最终得到最少需要的可重合路径数量。内容包括问题分析和解题代码。

摘要由CSDN通过智能技术生成

题目

给出一个无向图，问最少需要多少条可重合的路径覆盖所有点

分析

先跑一遍floyd传递闭包，把题目转换成最小不可相交路径覆盖，然后再跑二分图，最小不可相交路径覆盖=总点数-最大匹配即可

代码

#include <cstdio>
#include <cctype>
#include <cstring>
#define rr register
using namespace std;
bool a[501][501],cover[501];
int n,ans,lin[501];
inline signed iut(){
	rr int ans=0; rr char c=getchar();
	while (!isdigit(c)) c=getchar();
	while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
	return ans;
}
inline void print(int ans){
	if (ans>9) print(ans/10);
	putchar(ans%10+48);
}
inline bool dfs(int x){
	for (rr int i=1;i<=n;++i)
	if (a[x][i]&&!cover[i]){
		cover[i]=1;
		rr int q=lin[i];
		lin[i]=x;
		if (!q||dfs(q)) return 1;
	    lin[i]=q;
	}
	return 0;
}
signed main(){
	while (1){
		n=iut(),ans=0;
		rr int m=iut();
		if (!n&&!m) break;
		memset(a,0,sizeof(a));
		memset(lin,0,sizeof(lin));
		while (m--) a[iut()][iut()]=1;
		for (rr int k=1;k<=n;++k)
		for (rr int i=1;i<=n;++i) if (a[i][k])  
		for (rr int j=1;j<=n;++j) a[i][j]|=a[k][j];
	    for (rr int i=1;i<=n;++i){
	    	memset(cover,0,sizeof(cover));
	    	ans+=dfs(i);
	    }
	    print(n-ans),putchar(10);
	}
	return 0;
}