2018.10.10 loj#6002. 「网络流 24 题」最小路径覆盖（最大流）

传送门
一道求最小路径覆盖的板子题。

---

可以转化为求二分图的最大匹配。
如何转化呢？
在没有边的时候显然最小路径覆盖数等于点数。
假设当前加入一条边$(u,v)$。
相当于匹配了两条路径，于是假设现在匹配了$k$次，那么只需要$n-k$条路径。
代码：

#include<bits/stdc++.h>
#define N 505
#define M 20005
using namespace std;
inline int read(){
	int ans=0;
	char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
	return ans;
}
struct edge{int v,c,next;};
int n,m,tot=0,pred[N];
bool vis[N];
struct Dinic{
	int first[N],cur[N],cnt,d[N],s,t;
	edge e[M<<1];
	inline void init(){s=0,memset(first,-1,sizeof(first)),cnt=-1;}
	inline void addedge(int u,int v,int c){e[++cnt].v=v,e[cnt].c=c,e[cnt].next=first[u],first[u]=cnt;}
	inline void add(int u,int v,int c){addedge(u,v,c),addedge(v,u,0);}
	inline bool bfs(){
		queue<int>q;
		memset(d,-1,sizeof(d)),d[s]=0,q.push(s);
		while(!q.empty()){
			int x=q.front();
			q.pop();
			for(int i=first[x];~i;i=e[i].next){
				int v=e[i].v;
				if(~d[v]||!e[i].c)continue;
				d[v]=d[x]+1,q.push(v);
			}
		}
		return ~d[t];
	}
	inline int dfs(int p,int f){
		if(p==t||!f)return f;
		int flow=f;
		for(int&i=cur[p];~i;i=e[i].next){
			int v=e[i].v;
			if(!flow)break;
			if(e[i].c&&d[v]==d[p]+1){
				int tmp=dfs(v,min(flow,e[i].c));
				if(!tmp)d[v]=-1;
				e[i].c-=tmp,e[i^1].c+=tmp,flow-=tmp;
			}
		}
		return f-flow;
	}
	inline int solve(){
		int ret=0;
		while(bfs())memcpy(cur,first,sizeof(first)),ret+=dfs(s,0x3f3f3f3f);
		return ret;
	}
}dinic;
int main(){
	dinic.init(),n=read(),m=read(),dinic.t=n*2+1;
	for(int i=1;i<=n;++i)dinic.add(dinic.s,i,1),dinic.add(i+n,dinic.t,1);
	for(int u,v,i=1;i<=m;++i)u=read(),v=read(),dinic.add(u,v+n,1);
	int ans=n-dinic.solve();
	for(int i=0;i<=dinic.cnt;i+=2){
		if(dinic.e[i].c)continue;
		pred[dinic.e[i^1].v]=dinic.e[i].v-n,vis[dinic.e[i^1].v-n]=1;
	}
	for(int i=1;i<=n;++i){
		if(vis[i])continue;
		int k=i;
		while(k)printf("%d ",k),k=pred[k];
		puts("");
	}
	printf("%d",ans);
	return 0;
}