很变态,被自己板子的特殊性折磨了四个小时,警钟长鸣。
每条边必须要经过至少 1 1 1 次,可以想到是上下界网络流。
每次可以降落到滑雪场的任何一个地点,且具有多个滑雪终点,可以想到是有源汇。
求最少的飞行次数,可以想到是最小流。
即本题为有源汇上下界最小流。
建图:
可参考:网络流中有关上下界限制的学习笔记
超级源点 s 1 s1 s1 向所有地点连下界为 0 0 0,上界为 i n f inf inf 的边。
由所有地点向超级汇点 t 1 t1 t1 连下界为 0 0 0 ,上界为 i n f inf inf 的边。
原图中存在的所有边,下界为 1 1 1,上界为 i n f inf inf。
由超级汇点 t 1 t1 t1 向超级源点 s 1 s1 s1 连一条下界为 0 0 0,上界为 i n f inf inf 的边,让原图的源汇点流量平衡。
对原图中的每个点 x x x,若 i n x < 0 in_x<0 inx<0 ,向平衡汇点 t 2 t2 t2 连容量为 − i n x -in_x −inx 的边。
对原图中的每个点 x x x,若 i n x > 0 in_x>0 inx>0 ,由平衡源点 s 2 s2 s2 连容量为 i n x in_x inx 的边。
跑一遍平衡源点 s 2 s2 s2 到平衡汇点 t 2 t2 t2 的最大流,则可行流为从 t 2 t2 t2 到 s 2 s2 s2 边上经过的流量。
去掉原图汇点 t 1 t1 t1 到原图源点 s 1 s1 s1 间的 i n f inf inf 边,求从原图汇点 t 1 t1 t1 到原图源点 s 1 s1 s1 间的最大流。
则最小流为两个流量相减。
参考代码
#include <bits/stdc++.h>
#define itn int
#define int long long
#define endl "\n"
#define PII pair<int, int>
using namespace std;
const int N = 2e5 + 10;
const int M = 5e4 + 10;
const itn inf = 0x3f3f3f3f;
const int mod = 1e9 + 7;
// const int mod = 998244353;
struct Edge {
int from, to, cap, flow;
Edge(int f, int t, int c, int fl) {
from = f;
to = t;
cap = c;
flow = fl;
}
};
struct Dinic {
int n, m, s, t; //结点数,边数(包括反向弧),源点编号和汇点编号
vector<Edge> edges; //边表。edge[e]和edge[e^1]互为反向弧
vector<int> G[N]; //邻接表,G[i][j]表示节点i和第j条边在e数组中的序号
bool vis[N]; // BFS使用
int d[N]; //从起点到i的距离
int cur[N]; //当前弧下标
void clear_all(int n) {
for (int i = 0; i < n; i++)
G[i].clear();
edges.clear();
}
void clear_flow() {
int len = edges.size();
for (int i = 0; i < len; i++)
edges[i].flow = 0;
}
void add_edge(int from, int to, int cap) {
edges.push_back(Edge(from, to, cap, 0));
edges.push_back(Edge(to, from, 0, 0));
m = edges.size();
G[from].push_back(m - 2);
G[to].push_back(m - 1);
}
bool BFS() {
memset(vis, 0, sizeof(vis));
queue<int> q;
q.push(s);
d[s] = 0;
vis[s] = 1;
while (!q.empty()) {
int x = q.front();
q.pop();
int len = G[x].size();
for (int i = 0; i < len; i++) {
Edge& e = edges[G[x][i]];
if (!vis[e.to] && e.cap > e.flow) {
vis[e.to] = 1;
d[e.to] = d[x] + 1;
q.push(e.to);
}
}
}
return vis[t];
}
int DFS(int x, int a) {
if (x == t || a == 0)
return a;
int flow = 0, f, len = G[x].size();
for (int& i = cur[x]; i < len; i++) {
Edge& e = edges[G[x][i]];
if (d[x] + 1 == d[e.to] &&
(f = DFS(e.to, min(a, e.cap - e.flow))) > 0) {
e.flow += f;
edges[G[x][i] ^ 1].flow -= f;
flow += f;
a -= f;
if (a == 0)
break;
}
}
return flow;
}
int maxflow(int s, int t) {
this->s = s;
this->t = t;
int flow = 0;
while (BFS()) {
memset(cur, 0, sizeof(cur));
flow += DFS(s, inf);
}
return flow;
}
int mincut() { // call this after maxflow
int ans = 0;
int len = edges.size();
for (int i = 0; i < len; i++) {
Edge& e = edges[i];
if (vis[e.from] && !vis[e.to] && e.cap > 0)
ans++;
}
return ans;
}
void reduce() {
int len = edges.size();
for (int i = 0; i < len; i++)
edges[i].cap -= edges[i].flow;
}
} dinic;
int n, m;
int down[N], up[N], in[N];
int s1, s2, t1, t2;
void ins(int u, int v, int l, int r) {
dinic.add_edge(u, v, r - l);
in[u] -= l;
in[v] += l;
}
void solve() {
cin >> n;
s1 = n + 1, t1 = n + 2; //原图原汇
s2 = n + 3, t2 = n + 4; //平衡原汇
for (int i = 1; i <= n; i++) {
int x;
cin >> x;
dinic.add_edge(s1, i, inf);
dinic.add_edge(i, t1, inf);
for (int j = 1; j <= x; j++) {
int v;
cin >> v;
//上限inf 下限1
in[i]--;
in[v]++;
dinic.add_edge(i, v, inf);
}
}
for (int i = 1; i <= t1; i++) {
if (in[i] > 0)
dinic.add_edge(s2, i, in[i]);
else if (in[i] < 0)
dinic.add_edge(i, t2, -in[i]);
}
dinic.add_edge(t1, s1, inf);
int len = dinic.edges.size();
int ans1 = dinic.maxflow(s2, t2);
//以下为自用dinic板子特殊性 reduce将所有容量-流量 clear将所有流量清空
dinic.reduce();
ans1 = dinic.edges[len - 1].cap;
dinic.clear_flow();
dinic.edges[len - 2].cap = dinic.edges[len - 1].cap = 0;
int ans2 = dinic.maxflow(t1, s1);
// cout << ans1 << " " << ans2 << endl;
cout << ans1 - ans2 << endl;
}
signed main() {
ios_base::sync_with_stdio(false), cin.tie(0), cout.tie(0);
cout << fixed << setprecision(12);
// init();
int T = 1;
// cin >> T;
while (T--) {
solve();
}
return 0;
}