洛谷OJ：P2764 最小路径覆盖问题（网络流）

题目描述

给定有向图 G=(V,E) 。设 P 是 G 的一个简单路(顶点不相交)的集合。如果 V 中每个定点恰好在P的一条路上，则称 P 是 G 的一个路径覆盖。P中路径可以从 V 的任何一个定点开始，长度也是任意的，特别地，可以为 0 。G 的最小路径覆盖是 G 所含路径条数最少的路径覆盖。设计一个有效算法求一个 DAG (有向无环图) G 的最小路径覆盖。

提示：设 V={1,2,...,n} ，构造网络 G1​={V1​,E1​} 如下：

V1​={x0​,x1​,...,xn​}∪{y0​,y1​,...,yn​}

E1​={(x0​,xi​):i∈V}∪{(yi​,y0​):i∈V}∪{(xi​,yj​):(i,j)∈E}

每条边的容量均为 1 ，求网络 G1​ 的 (x0​,y0​) 最大流。

输入格式

第一行有 2 个正整数 n 和 m 。 nn 是给定GAP(有向无环图) G 的顶点数， m 是 G 的边数。接下来的 m行,每行有两个正整数 i 和 j 表示一条有向边 (i,j)。

输出格式

从第1 行开始，每行输出一条路径。文件的最后一行是最少路径数。

输入输出样例

输入 #1复制

11 12
1 2
1 3
1 4
2 5
3 6
4 7
5 8
6 9
7 10
8 11
9 11
10 11

输出 #1复制

1 4 7 10 11
2 5 8
3 6 9
3

说明/提示

1≤n≤150,1≤m≤6000

由@FlierKing提供SPJ

思路：利用网络流解决最少路径覆盖的问题一般建图方式是将每个点拆分成两个点【例如点的标号为x,则我们可以拆分成x和x+】，分别与源点和汇点相连，如下图所示【出处：巨巨的知乎】：

建图前：                               

 建图后： 

在该图中，我们令每一条从源点指向汇点的边的边长为1，跑一遍最大流，便能得到最大合并路径数，再用点数去减即得最小路径覆盖数。这几乎是显然的：从A点到B'点的每一条流，都代表着一次合并。而从源点只给每个点输送1单位流量，又保证了每个点只被经过一次。【同样摘自上述链接】，对于本题输出路径的操作，只需额外开一个nxt数组即可，而起点x一定是edges[x+n][t]=1的点。【t是汇点】。

#include<set>
#include<map>
#include<queue>
#include<vector>
#include<string>
#include<math.h>
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<functional>
using namespace std;
#define maxn 550
#define ll long long
#define inf 1e18
int n,m,s,t,lv[maxn],cur[maxn],nxt[maxn];
ll edges[maxn][maxn];
//lv是每个点的层数
bool bfs(){
	memset(lv,-1,sizeof(lv));
	lv[s]=0;
	queue<int>q;
	q.push(s);
	while(!q.empty()){
		int now=q.front();
		q.pop();
		for(int i=1;i<=2*n+1;i++){
			if(edges[now][i]>0 && lv[i]==-1){
				lv[i]=lv[now]+1;
				q.push(i);
			}
		}
	}
	return lv[t]!=-1;
}
ll dfs(int p=s,ll flow=inf){
	if(p==t)
		return flow;
	ll mn=flow;
	for(int i=1;i<=2*n+1;i++){
		if(edges[p][i]>0 && lv[i]==lv[p]+1){
			ll c=dfs(i,min(edges[p][i],mn));
			mn-=c;
			edges[p][i]-=c;
			edges[i][p]+=c;
			if(c && p) nxt[p]=i-n;
		}
	}
	return flow-mn;
}
ll dinic(){
	ll maxFlow=0;
	while(bfs())
		maxFlow+=dfs();
	return maxFlow;
}
int main(void){
	int x,y;
	scanf("%d%d",&n,&m);
	t=2*n+1;
	for(int i=1;i<=n;i++){
		edges[s][i]=1;
		edges[i+n][t]=1;
	}
	for(int i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		edges[x][y+n]=1;
	}
	ll pathNum=n-dinic();
	for(int i=1;i<=n;i++){
		if(edges[i+n][t]){
			int now=i;
			printf("%d",now);
			while(nxt[now]){
				now=nxt[now];
				printf(" %d",now);
			}
			printf("\n");
		}
	}
	printf("%lld\n",pathNum);
	return 0;
}