[Luogu P4843] [BZOJ 2502] 清理雪道

洛谷传送门

BZOJ传送门

题目描述

滑雪场坐落在FJ省西北部的若干座山上。

从空中鸟瞰，滑雪场可以看作一个有向无环图，每条弧代表一个斜坡（即雪道），弧的方向代表斜坡下降的方向。

你的团队负责每周定时清理雪道。你们拥有一架直升飞机，每次飞行可以从总部带一个人降落到滑雪场的某个地点，然后再飞回总部。从降落的地点出发，这个人可以顺着斜坡向下滑行，并清理他所经过的雪道。

由于每次飞行的耗费是固定的，为了最小化耗费，你想知道如何用最少的飞行次数才能完成清理雪道的任务。

输入输出格式

输入格式：

输入文件的第一行包含一个整数$n(2\le n\le 100)$ – 代表滑雪场的地点的数量。

接下来的$n$行，描述$1\sim n$号地点出发的斜坡，第$i$行的第一个数为$m_i(0\le m_i<n)$ ，后面共有$m_i$个整数，由空格隔开，每个整数$a_{i,j}$互不相同，代表从地点$i$下降到地点$a_{i,j}$的斜坡。每个地点至少有一个斜坡与之相连。

输出格式：

输出文件的第一行是一个整数$k$ – 直升飞机的最少飞行次数。

输入输出样例

输入样例#1：

8
1 3
1 7
2 4 5
1 8
1 8
0
2 6 5
0

输出样例#1：

4

解题分析

要强制每条边都有至少$1$的流量且为最小流， 上下界网络流打板即可。

代码如下：

#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <algorithm>
#include <queue>
#define R register
#define IN inline
#define W while
#define gc getchar()
#define INF 1e8
#define MX 105
template <class C>
IN void in(C &x)
{
	x = 0; R char c = gc;
	for (; !isdigit(c); c = gc);
	for (;  isdigit(c); c = gc)
	x = (x << 1) + (x << 3) + c - 48;
}
template <class C> IN C min(C a, C b) {return a < b ? a : b;}
int n, cnt, SS, TT, S, T, sum;
int head[MX], layer[MX], deg[MX];
struct Edge {int to, fl, nex;} edge[100500];
IN void add(R int from, R int to, R int fl)
{
	edge[++cnt] = {to, fl, head[from]}, head[from] = cnt;
	edge[++cnt] = {from, 0, head[to]}, head[to] = cnt;
}
namespace Dinic
{
	std::queue <int> q;
	IN bool BFS(R int st, R int ed)
	{
		std::memset(layer, 0, sizeof(layer));
		layer[st] = 1; R int now; q.push(st);
		W (!q.empty())
		{
			now = q.front(); q.pop();
			for (R int i = head[now]; ~i; i = edge[i].nex)
			{
				if (edge[i].fl && (!layer[edge[i].to]))
				layer[edge[i].to] = layer[now] + 1, q.push(edge[i].to);
			}
		} 
		return layer[ed];
	}
	int DFS(R int now, R int avai, R int ed)
	{
		if (now == ed) return avai;
		R int lef = avai, buf;
		for (R int i = head[now]; ~i; i = edge[i].nex)
		{
			if (edge[i].fl && layer[edge[i].to] == layer[now] + 1)
			{
				buf = DFS(edge[i].to, min(lef, edge[i].fl), ed);
				if (!buf) continue;
				lef -= buf, edge[i].fl -= buf, edge[i ^ 1].fl += buf;
				if (!lef) return avai;
			}
		}
		return avai - lef;
	}
	IN int solve(R int st, R int ed)
	{
		int ret = 0;
		W (BFS(st, ed)) ret += DFS(st, INF, ed);
		return ret;
	}
	IN void init()
	{
		solve(SS, TT);
		int res = edge[cnt].fl;
		edge[cnt].fl = edge[cnt - 1].fl = 0;
		for (R int i = head[SS]; i; i = edge[i].nex) edge[i].fl = edge[i ^ 1].fl = 0;
		for (R int i = head[TT]; i; i = edge[i].nex) edge[i].fl = edge[i ^ 1].fl = 0;
		printf("%d", res - solve(T, S));
	}
}
int main(void)
{
	int num, foo;
	std::memset(head, cnt = -1, sizeof(head));
	in(n); S = n + 1, T = n + 2, SS = n + 3, TT = n + 4;
	for (R int i = 1; i <= n; ++i)
	{
		in(num);
		add(S, i, INF), add(i, T, INF);
		for (R int j = 1; j <= num; ++j)
		in(foo), deg[i]--, deg[foo]++, add(i, foo, INF);
	}
	for (R int i = 1; i <= T; ++i)
	{
		if (deg[i] > 0) add(SS, i, deg[i]), sum += deg[i];
		else if (deg[i] < 0) add(i, TT, -deg[i]);
	}
	add(T, S, 1e9);
	Dinic::init();
}