网络流24题-太空飞行计划问题

太空飞行计划问题

题意 现有一个实验集合和一个仪器集合, 完成第i项实验都能得到$k_i$, 需要具仪器集合$c_i$, 现现给定若干实验, 若干仪器, 求最大收益, 并且要输出实验编号和仪器编号.

建图: 显然, 有仪器和实验两种集合, 考虑把最大收益转化成最大流, 但是实验是收益, 仪器是成本, 怎么转化呢, 每种实验有两种花费, 要么就直接买器材, 要么就直接不做这个实验. 如果给个前提 - 所有的实验我都做, 那么花费情况就变成, 我不做这个实验(放弃收益), 要么我买这个实验的所有仪器. 于是建立超级源点超级汇点, 二分图左半部分是实验, 边权是实验的收益, 右边半部分是仪器, 仪器连汇点, 边权是仪器的花费, 中间实验连仪器, 边权无穷大. 于是你比划比划就知道这个图的最小割就是在做所有实验下的最小花费, 于是所有的实验经费减去最小花费即为所求.

那么怎么输出方案呢, 这个题和这个题输出不同的地方就在于前者是输出左半部分, 然后输出右半部分, 而后者是以二元组的形式输出左半部分连着哪右半部分, 其实只用输出深度数组d[], 不为0的点就可以了, bfs都会给下一层的深度编号, dfs会把下一层不走的点变成0, 于是最后有深度的点都是一定会访问过的点

solution(输入是有毒)

#include <bits/stdc++.h>
using namespace std;
#define _rep(i, a, b) for (int i = (a); i <= (b); ++i)
#define _rev(i, a, b) for (int i = (a); i >= (b); --i)
#define _for(i, a, b) for (int i = (a); i < (b); ++i)
#define _rof(i, a, b) for (int i = (a); i > (b); --i)
#define oo 0x3f3f3f3f
#define ll long long
#define db double
#define eps 1e-8
#define bin(x) cout << bitset<10>(x) << endl;
#define what_is(x) cerr << #x << " is " << x << "s" << endl;
#define met(a, b) memset(a, b, sizeof(a))
#define all(x) x.begin(), x.end()
#define pii pair<int, int>
const int maxn = 120;
int n, m, pro[maxn], cnt = 1, d[maxn], head[maxn], take[maxn], s, t, sum;
char buf[10000];
struct node
{
	int nxt, to, cost;
} way[6000];
void addedge(int from, int to, int cost)
{
	way[++cnt].cost = cost;
	way[cnt].to = to;
	way[cnt].nxt = head[from];
	head[from] = cnt;
}
int nxt()
{
	int ret;
	scanf("%d", &ret);
	return ret;
}
bool bfs()
{
	met(d, 0);
	queue<int> q;
	q.push(s);
	d[s] = 1;
	while (!q.empty())
	{
		int cur = q.front();
		q.pop();
		for (int i = head[cur]; i; i = way[i].nxt)
		{
			int to = way[i].to, cost = way[i].cost;
			if (!d[to] && cost)
			{
				d[to] = d[cur] + 1;
				if (to == t)
					return 1;
				q.push(to);
			}
		}
	}
	return 0;
}
int dfs(int cur, int flow)
{
	if (cur == t)
		return flow;
	int rest = flow, k;
	for (int i = head[cur]; i && rest; i = way[i].nxt)
	{
		int to = way[i].to, cost = way[i].cost;
		if (d[to] == d[cur] + 1 && cost)
		{
			k = dfs(to, min(cost, rest));
			if (!k)
				d[to] = 0;
			way[i].cost -= k;
			way[i ^ 1].cost += k;
			rest -= k;
		}
	}
	return flow - rest;
}
void Dinic()
{
	int ans = 0;
	while (bfs())
		ans += dfs(s, oo);
	_rep(i, 1, m)if(d[i] > 0)cout << i << " ";
	cout << endl;
	_rep(i, m+1, m+n)if(d[i] > 0)cout << i - m <<" ";
	cout << endl;
	cout << sum - ans << endl;
}
int main()
{
	m = nxt(), n = nxt();
	t = n + m + 3;
	_rep(i, 1, m)
	{
		pro[i] = nxt();
		sum += pro[i];
		met(buf, 0);
		cin.getline(buf, 10000);
		int ulen = 0, tmp;
		while (sscanf(buf + ulen, "%d", &tmp) == 1)
		{
			addedge(i, tmp + m, oo), addedge(tmp + m, i, 0);
			if (tmp == 0)
				ulen++;
			else
			{
				while (tmp)
				{
					tmp /= 10;
					ulen++;
				}
			}
			ulen++;
		}
	}


	_rep(i, 1 + m, n + m)
		take[i] = nxt();

	_rep(i, 1, m) addedge(s, i, pro[i]), addedge(i, s, 0);
	_rep(i, 1 + m, m + n) addedge(i, t, take[i]), addedge(t, i, 0);
	Dinic();
}