洛谷·[网络流24题] [CTSC1999]家园 / 星际转移问题

初见安~这里是传送门：洛谷P2754 [CTSC1999]家园 / 星际转移问题

题解

啊同样是题意比较友好的一个题呢。

可以看到每条航线都是一个有周期的环，也就是说每个时刻你面对的图都长得不一样。比较容易想到跑分层图最大流，只要最大流【也就是能运过去的最多人数】大于等于地球上人数那就一定可以运过去，再加一层费用作为时间求最小费用流。但是！如果我们要这么建图的话那层数必须是所有周期的公倍数。周期在15以内，其公倍数为：，这么多层的话空间是一定会炸的。

但其实我们可以发现，上述建图方式的话你走了多少层，那就花了多少时间。并且极端情况，13个点，每走一步都要等待12的时间，13*12，这个步数也是非常少的，换言之我们建这么多层图就可以跑出答案了。

这里用的是一层一层边走边建图的方式，这样就跑最大流就可以了。但是涉及到枚举答案所以要先保证答案存在，要用并查集检查一下图的连通性。代码可供参考。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define maxn 100005
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
int read() {
	int x = 0, f = 1, ch = getchar();
	while(!isdigit(ch)) {if(ch == '-') f = -1; ch = getchar();}
	while(isdigit(ch)) x = (x << 1) + (x << 3) + ch - '0', ch = getchar();
	return x * f;
}

struct edge {int to, w, nxt;} e[maxn];
int head[maxn], k = 0;
void add(int u, int v, int w) {
	e[k] = {v, w, head[u]}; head[u] = k++;
	e[k] = {u, 0, head[v]}; head[v] = k++;
}

int n, m, p, h[maxn], num[maxn], st[70][70], S, T;
int fa[maxn];
int get(int x) {return fa[x] == x? x : get(fa[x]);}

int dis[maxn];
queue<int> q;
bool bfs() {
	memset(dis, 0, sizeof dis); while(q.size()) q.pop();
	q.push(S); dis[S] = 1;
	while(q.size()) {
		register int u = q.front(), v; q.pop();
		for(int i = head[u]; ~i; i = e[i].nxt) {
			v = e[i].to; if(e[i].w && !dis[v]) {
				q.push(v), dis[v] = dis[u] + 1;
				if(v == T) return true;
			}
		}
	}
	return false;
}

int Dinic(int u, int flow) {
	if(u == T) return flow;
	register int rst = flow, k;
	for(int i = head[u], v; ~i; i = e[i].nxt) {
		v = e[i].to; if(e[i].w && dis[v] == dis[u] + 1) {
			k = Dinic(v, min(rst, e[i].w));
			if(!k) dis[v] = 0;
			e[i].w -= k, e[i ^ 1].w += k;
			rst -= k;
		}
	}
	return flow - rst;
}

signed main() {
	memset(head, -1, sizeof head);
	n = read(), m = read(), p = read(); S = 0, T = 100000;
	for(int i = 1; i <= n + 2; i++) fa[i] = i;
	for(int i = 1; i <= m; i++) {
		h[i] = read(); num[i] = read();
		for(int j = 0; j < num[i]; j++) {
			st[i][j] = read();
			if(st[i][j] < 0) st[i][j] = n + 2;
			if(!st[i][j]) st[i][j] = n + 1;
			if(j > 0) fa[get(st[i][j])] = get(st[i][j - 1]);
		}
	}
	
	if(get(n + 1) != get(n + 2)) {puts("0"); return 0;}
	register int mx = 0;
	for(int ans = 1; ; ans++) {
		add(S, ans * (n + 1), INF);
		for(int i = 1; i <= m; i++) {
			register int u = (ans - 1) % num[i], v = ans % num[i];
			if(st[i][u] == n + 2) u = T;
			else u = (ans - 1) * (n + 1) + st[i][u];
			if(st[i][v] == n + 2) v = T;
			else v = ans * (n + 1) + st[i][v];
			add(u, v, h[i]);
		}
		while(bfs()) mx += Dinic(S, INF);
		if(mx >= p) {printf("%d\n", ans);return 0;}
		for(int i = 1; i <= n + 1; i++) add((ans - 1) * (n + 1) + i, ans * (n + 1) + i, INF);
	}
	return 0;
}

迎评：）
——End——