2018.10.11【POJ2594】Treasure Exploration（最小链覆盖）

传送门

---

解析：

其实最小链覆盖就是允许节点相交的最小路径覆盖。

那么这个问题与最小路径覆盖的匹配模型的区别就在于，每个点可以直接和自己能到达的节点的后继匹配。

所以我们直接做一个$Floyed$求传递闭包，然后按照最小路径覆盖的方式建跑网络流就行了。

---

代码：

#include<iostream>
#include<cstdio>
#include<cctype>
#include<cstring>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const

inline int getint(){
	re int num;
	re char c;
	while(!isdigit(c=gc()));num=c^48;
	while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
	return num;
}

cs int N=501;
bool g[N][N];
int match[N];
bool vis[N];
int n,m;

inline bool find(int u){
	for(int re v=1;v<=n;++v){
		if(g[u][v]&&!vis[v]){
			vis[v]=true;
			if(match[v]==-1||find(match[v])){
				match[v]=u;
				return true;
			}
		}
	}
	return false;
}

signed main(){
	while(~scanf("%d%d",&n,&m)){
		if(n==0&&m==0)break;
		memset(g,0,sizeof g);
		for(int re i=1;i<=m;++i){
			int u=getint(),v=getint();
			g[u][v]=true;
		}
		
		for(int re k=1;k<=n;++k)
		for(int re i=1;i<=n;++i)
		if(i!=k)for(int re j=1;j<=n;++j)
		g[i][j]=g[i][j]||(g[i][k]&&g[k][j]);
		
		int ans=0;
		memset(match,-1,sizeof match);
		for(int re i=1;i<=n;++i){
			memset(vis,0,sizeof vis);
			if(find(i))++ans;
		}
		printf("%d\n",n-ans);
	}
	return 0;
}