[网络流24题] 洛谷P2764 最小路径覆盖问题 最大流

给出$n\le 150$个点，$m\le 6e3$条边的有向图。求有向图的最小路径覆盖。最小路径覆盖的定义在这里。
网络流的经典应用。如果不连边的话就是点数$n$，然后把所有的点拆成起点和终点：对有向图里面所有的边，都从起点连到终点。求最大匹配，多匹配一条边就对于路径数少一个。答案就是总点数$-$最大匹配数。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
namespace Dinic {
	const int N=1e3+3,M=1e4+3; 
	int front[N],head[N],dep[N],tot=-1,n,S,T;
	int to[N];
	bool tag[N];
	struct edge{ int v,nxt;ll w; }e[2*M];
	void init() {
		tot=-1;
		for(int i=1;i<=n;i++) front[i]=-1;
	}
	void add(int u,int v,ll w) {
		e[++tot].v=v,e[tot].w=w,e[tot].nxt=front[u],front[u]=tot;
		e[++tot].v=u,e[tot].w=0,e[tot].nxt=front[v],front[v]=tot;
	}
	bool bfs() {
		for(int i=1;i<=n;i++) dep[i]=0;
		queue<int> q;
		q.push(S);
		dep[S]=1;
		while(!q.empty()) {
			int u=q.front();
			q.pop();
			for(int i=front[u];i!=-1;i=e[i].nxt) {
				int v=e[i].v;
				if(dep[v]==0&&e[i].w) {
					dep[v]=dep[u]+1;
					q.push(v); 
				}
			}
		}
		return dep[T];
	}
	ll dfs(int cur,ll dis=1e18) {
		if(cur==T) return dis;
		ll flow=0,sum=0;
		for(int &i=head[cur];i!=-1;i=e[i].nxt) {
			int v=e[i].v;
			if(dep[v]==dep[cur]+1&&e[i].w) {
				flow=dfs(v,min(dis,e[i].w));
				dis-=flow,sum+=flow;
				e[i].w-=flow,e[i^1].w+=flow;
				if(!dis) break;
			}
		}
		return sum;
	}
	void solve() {
		ll ans=0;
		while(bfs()) {
			for(int i=1;i<=n;i++) head[i]=front[i];
			ans+=dfs(S);
		}
		n=(n-2)/2;
		for(int i=1;i<=n;i++) {
			for(int j=front[i];j!=-1;j=e[j].nxt) {
				if(e[j].v==S) continue;
				else if(e[j].w==0) {
					to[i]=e[j].v-n;
					break;
				}
			}
		}
		for(int i=1;i<=n;i++) {
			int cur=i;
			if(tag[cur]) continue;
			while(cur) { 
				printf("%d ",cur);
				cur=to[cur]; 
				tag[cur]=1;
			}
			if(to[i]) printf("\n"); 
		}
		printf("%lld\n",n-ans); 
	}
}
int main() {
	int n,m;
	scanf("%d%d",&n,&m); 
	int S=2*n+1;
	int T=2*n+2;
	Dinic::S=S;
	Dinic::T=T;
	Dinic::n=2*n+2;
	Dinic::init(); 
	for(int i=1;i<=n;i++) 
		Dinic::add(S,i,1);
	for(int i=1;i<=n;i++)
	 	Dinic::add(n+i,T,1);
	for(int i=1;i<=m;i++) {
		int u,v;
		scanf("%d%d",&u,&v);
		Dinic::add(u,v+n,1);
	}
	Dinic::solve();
	return 0;
}