poj2594 Treasure Exploration(有向图可相交最小路径覆盖/floyd传递闭包)

题目

DAG，n(n<=500)个点，m(m<=5000)条边，

求最小路径覆盖的条数，路径允许相交

思路来源

https://www.cnblogs.com/hujunzheng/p/3919231.html

题解

图片来自思路来源，如果允许路径相交，就把间接相连的路径也连边

这样如果要走形如1-2-5，2-3-4这样相交的路径，就认为是走1-2-5，3-4不相交的路径即可

floyd求传递闭包，即更新一下任两点的连通性，连通就连边，再跑最大匹配

代码

#include<iostream>
#include<cstdio> 
#include<cstring>
using namespace std;
const int N=505;
const int M=5005;
bool vis[N*2];
bool dp[N][N];
int head[N*2],cnt,ans;
int link[N*2]; 
int t,n,m,u,v;
struct edge{int to,next,w;}e[N*N*2];
void init()
{
	memset(head,0,sizeof head);
	memset(dp,0,sizeof dp);
	cnt=0;
}
void add(int u,int v)
{
	e[++cnt].to=v;
	e[cnt].next=head[u];
	head[u]=cnt;
}
bool dfs(int u)
{
	for(int i=head[u];i;i=e[i].next)
	{
		int v=e[i].to;
		if(!vis[v])
		{
			vis[v]=1;
			if(link[v]==-1||dfs(link[v]))
			{
				link[v]=u;
				return 1;
			}
		}
	}
	return 0;
}
int hungary(int n)
{
	int res=0;
	memset(vis,0,sizeof vis);
	memset(link,-1,sizeof link);
	for(int u=1;u<=n;++u)
	{
		memset(vis,0,sizeof vis);
		res+=dfs(u);
	}
	return res;
}
int main()
{
	while(~scanf("%d%d",&n,&m))
	{
		if(!n&&!m)break;
		init();
		for(int i=1;i<=m;++i)
		{
			scanf("%d%d",&u,&v);
			dp[u][v]=1;
		}
		for(int k=1;k<=n;++k)
		{
			for(int i=1;i<=n;++i)
			{
				for(int j=1;j<=n;++j)
				{
					dp[i][j]|=dp[i][k]&dp[k][j];
				}
			}
		}
		for(int i=1;i<=n;++i)
		{
			for(int j=1;j<=n;++j)
			{
				if(i==j)continue;
				if(dp[i][j])add(i,j);
			}
		}
		ans=hungary(n);
		printf("%d\n",n-ans);
	}
	return 0;
}