初见安~这里是传送门:洛谷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——