题目描述
给定有向图 G=(V,E)G=(V,E) 。设 PP 是 GG 的一个简单路(顶点不相交)的集合。如果 VV 中每个定点恰好在PP的一条路上,则称 PP 是 GG 的一个路径覆盖。PP中路径可以从 VV 的任何一个定点开始,长度也是任意的,特别地,可以为 00 。GG 的最小路径覆盖是 GG 所含路径条数最少的路径覆盖。设计一个有效算法求一个 GAP (有向无环图) GG 的最小路径覆盖。
提示:设 V={1,2,...,n}V={1,2,...,n} ,构造网络 G_1={V_1,E_1}G1={V1,E1} 如下:
V_1={x_0,x_1,...,x_n}cup{y_0,y_1,...,y_n}V1={x0,x1,...,xn}∪{y0,y1,...,yn}
E_1={(x_0,x_i):iin V}cup{(y_i,y_0):iin V}cup{(x_i,y_j):(i,j)in E}E1={(x0,xi):i∈V}∪{(yi,y0):i∈V}∪{(xi,yj):(i,j)∈E}
每条边的容量均为 11 ,求网络 G_1G1 的 (x_0,y_0)(x0,y0) 最大流。
输入输出格式
输入格式:
第一行有 22 个正整数 nn 和 mm 。 nn 是给定 ext{GAP}GAP(有向无环图) GG 的顶点数, mm 是 GG 的边数。接下来的 mm行,每行有两个正整数 ii 和 jj 表示一条有向边 (i,j)(i,j)。
输出格式:
从第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 4 7 10 11 2 5 8 3 6 9 3
说明
1leq nleq 150,1leq mleq 60001≤n≤150,1≤m≤6000
由@FlierKing提供SPJ
这个题目有点难读,我读了挺久的,这个题目是让你找最少的路径使得所有的路径合在一起可以覆盖了整个图。
没有最少自然可以是n,就是所有点自成一条路径。
拆分题目有点难写,我开始根本就没有什么思路,后来看到二分图匹配有一个关于这个的定理,
DAG图的最小路径覆盖数 = 节点数 - 二分图的最大匹配
所以这个就可以求解,然后这个转化成二分图,就需要把一个点拆成两个点,然后再合并。
说多了也没用,自己看代码理解吧。
我开始看题解,以为只能用链式前向星建图,真的要哭了,不太会用,后来发现也可以用vector
#include <cstdio> #include <cstdlib> #include <cstring> #include <algorithm> #include <queue> #include <iostream> #include <vector> #define inf 0x3f3f3f3f using namespace std; const int maxn = 10050; int level[maxn], iter[maxn]; int head[maxn]; int n,m,s = 0, t = 0; struct node { int to, cap,flow; int nex; }exa[maxn<<1]; int cnt = 1; void add(int u,int v,int c) { exa[++cnt].to = v; exa[cnt].cap = c; exa[cnt].flow = 0; exa[cnt].nex = head[u]; head[u] = cnt; //printf("%d %d %d %d %d ",cnt,u, exa[cnt].to, exa[cnt].cap, exa[cnt].nex); exa[++cnt].to = u; exa[cnt].cap = 0; exa[cnt].flow = 0; exa[cnt].nex = head[v]; head[v] = cnt; //printf("%d %d %d %d %d ",cnt,v, exa[cnt].to, exa[cnt].cap, exa[cnt].nex); } void bfs(int s) { queue<int>que; que.push(s); memset(level, -1, sizeof(level)); level[s] = 0; while(!que.empty()) { int u = que.front(); que.pop(); //printf("%d %d %d ", u,head[u],exa[head[u]].nex); for(int i=head[u];i!=-1;i=exa[i].nex) { node &now = exa[i]; //printf("%d %d %d ", i, exa[i].to, exa[i].nex); if(level[now.to]<0&&now.cap>now.flow) { level[now.to] = level[u] + 1; que.push(now.to); } } } } int tag[maxn], to[maxn]; int dfs(int u,int v,int f) { if (u == v) return f; for(int &i=iter[u];i!=-1;i=exa[i].nex) { node &now = exa[i]; if(now.cap>now.flow&&level[now.to]>level[u]) { int d = dfs(now.to, v, min(f, now.cap - now.flow)); if(d>0) { to[u] = now.to; if (u != s) tag[now.to - n] = 1; now.flow += d; exa[i ^ 1].flow -= d; return d; } } } return 0; } void init() { memset(head, -1, sizeof(head)); } int dinic(int s,int t) { int flow = 0; while(1) { bfs(s); if (level[t] < 0) break; for (int i = s; i <= t; i++) iter[i] = head[i]; int f; while ((f = dfs(s, t, inf)) > 0) flow += f; } for(int i=1;i<=n;i++) { if(!tag[i]) { int now = i; printf("%d ", now); while(to[now]&&to[now]!=t) { printf("%d ", to[now] - n); now = to[now] - n; } printf(" "); } } return flow; } int main() { init(); cin >> n >> m; t = 2 * n + 1; for (int i = 1; i <= n; i++) add(s, i, 1); for (int i = 1; i <= n; i++) add(i+n, t, 1); for(int i=1;i<=m;i++) { int a, b; cin >> a >> b; add(a, b+n, 1); } int ans = dinic(s, t); printf("%d ",n - ans); return 0; }
#include <cstdio> #include <cstdlib> #include <cstring> #include <algorithm> #include <queue> #include <iostream> #include <vector> #define inf 0x3f3f3f3f using namespace std; const int maxn = 10050; int level[maxn], iter[maxn]; int n,m,s = 0, t = 0; struct node { int from,to, cap,flow; node(int from=0,int to=0,int cap=0,int flow=0):from(from),to(to),cap(cap),flow(flow){} }; vector<node>e; vector<int>G[maxn]; void add(int u,int v,int c) { e.push_back(node(u, v, c, 0)); e.push_back(node(v, u, 0, 0)); int len = e.size(); G[u].push_back(len - 2); G[v].push_back(len - 1); } void bfs(int s) { queue<int>que; que.push(s); memset(level, -1, sizeof(level)); level[s] = 0; while(!que.empty()) { int u = que.front(); que.pop(); for(int i=0;i<G[u].size();i++) { node &now = e[G[u][i]]; if(level[now.to]<0&&now.cap>now.flow) { level[now.to] = level[u] + 1; que.push(now.to); } } } } int tag[maxn], to[maxn]; int dfs(int u,int v,int f) { if (u == v) return f; for(int &i=iter[u];i<G[u].size();i++) { node &now = e[G[u][i]]; if(now.cap>now.flow&&level[now.to]>level[u]) { int d = dfs(now.to, v, min(f, now.cap - now.flow)); if(d>0) { to[u] = now.to; // printf("%d %d ", u, to[u]); if (u != s) tag[now.to - n] = 1; now.flow += d; e[G[u][i] ^ 1].flow -= d; return d; } } } return 0; } void init() { for (int i = 1; i <= n + 2; i++) G[i].clear(); e.clear(); } int dinic(int s,int t) { int flow = 0; while(1) { bfs(s); if (level[t] < 0) break; memset(iter, 0, sizeof(iter)); int f; while ((f = dfs(s, t, inf)) > 0) flow += f; } //cout << endl; for(int i=1;i<=n;i++) { if(!tag[i]) { int now = i; printf("%d ", now); while(to[now]&&to[now]!=t) { printf("%d ", to[now] - n); now = to[now] - n; } printf(" "); } } return flow; } int main() { init(); cin >> n >> m; t = 2 * n + 1; for (int i = 1; i <= n; i++) add(s, i, 1); for (int i = 1; i <= n; i++) add(i+n, t, 1); for(int i=1;i<=m;i++) { int a, b; cin >> a >> b; add(a, b+n, 1); } int ans = dinic(s, t); printf("%d ",n - ans); return 0; }