网络流24题13. 星际转移问题

星际转移问题

Description

由于人类对自然资源的消耗，人们意识到大约在 2300 年之后，地球就不能再居住了。
于是在月球上建立了新的绿地，以便在需要时移民。令人意想不到的是，2177 年冬由于未知的原因，地球环境发生了连锁崩溃，人类必须在最短的时间内迁往月球。现有 n 个太空站位于地球与月球之间，且有 m 艘公共交通太空船在其间来回穿梭。每个太空站可容纳无限多的人，而每艘太空船 i 只可容纳 H[i]个人。每艘太空船将周期性地停靠一系列的太空站，例如：(1，3，4)表示该太空船将周期性地停靠太空站 134134134…。每一艘太空船从一个太空站驶往任一太空站耗时均为 1。人们只能在太空船停靠太空站(或月球、地球)时上、下船。
初始时所有人全在地球上，太空船全在初始站。试设计一个算法，找出所有人尽快地全部转移到月球上的运输方案。

Input

第 1 行有 3 个正整数 n（太空站个数），m（太空船个数）和 k（需要运送的地球上的人的个数）。其中 1<=m<=13, 1<=n<=20, 1<=k<=50。
接下来的 m 行给出太空船的信息。第 i+1 行说明太空船 pi。第 1 个数表示 pi 可容纳的人数 Hpi；第 2 个数表示 pi 一个周期停靠的太空站个数 r，1<=r<=n+2；随后 r 个数是停靠的太空站的编号(Si1,Si2,…,Sir)，地球用 0 表示，月球用-1 表示。时刻 0 时，所有太空船都在初始站，然后开始运行。在时刻 1，2，3…等正点时刻各艘太空船停靠相应的太空站。人只有在 0,1,2…等正点时刻才能上下太空船。

Output

输出全部人员安全转移所需的时间。
若问题无解，输出0。

Sample Input

2 2 1
1 3 0 1 2
1 3 1 2 -1

Sample Output

5

题解

这题的建图有点烦的，不是很好想且细节比较麻烦，直接贴题解了。
【问题分析】

分层图网络流问题，枚举答案，构造网络流判定。

【建模方法】

首先判断从地球到月球是否存在一条路线，如果不存在那么无解，否则把每个太空站按照每天拆分成d个点， <i,d> 表示第i个站第d天。建立附加源S汇T，顺序枚举答案Day。

1、对于第Day天，从S到 <0,Day> 连接一条容量为无穷大的有向边。
2、从 <−1,Day> 到T连接一条容量为无穷大的有向边。
3、对于第i个太空船，设第Day-1天在a处，第Day天在b处，从 <a,Day−1> 到 <b,Day> 连接一条容量为该太空船容量的有向边。
4、对于第i个太空站，从 <i,Day−1 >到 <i,Day> 连接一条容量为无穷大的有向边。
5、求当前网络最大流，如果最大流量大于等于地球上人数K，停止枚举，当前Day值就是答案。

【建模分析】

我们把网络优化问题转化为枚举答案+可行性判定问题。枚举天数，按天数把图分层，因为乘船每坐一站天数都要增加一，把太空船航线抽象成图中的一条边，跨图的两层。由于太空船容量有限，边上也要加上容量限制。除了坐船以外，人还可以在某个空间站等待下一班太空船的到来，所以每个点要与下一层同一点连接一条容量为无穷的边。这样在层限制的图上求出的网络最大流就是在当前天数以内能够从地面到月球的最多的人数，该人数随天数递增不递减，存在单调性。所以可以枚举答案或二分答案，用网络流判定。网络流判定问题更适合枚举答案，而不是二分，因为新增一些点和边只需要在原有的基础上增广，不必重新求网络流。

#include
#include
#include
using namespace std;

const int N = 30, M = 1000000 + 10, D = 50, inf = 0x3f3f3f3f;
struct Edge{
    int fr, to, cap, flow;
}edg[M];
int hd[N * D], nxt[M];
int d[N * D], vis[N * D], q[N * D], dfn;
int s, t;
int n, m, k, ans, tot;
int a[N][D], h[N];
void insert(int u, int v, int w){
    edg[tot].fr = u, edg[tot].to = v, edg[tot].cap = w;
    nxt[tot] = hd[u], hd[u] = tot;
    tot++;
    edg[tot].fr = v, edg[tot].to = u;
    nxt[tot] = hd[v], hd[v] = tot;
    tot++;
}
bool bfs(){
    int head = 1, tail = 1;
    q[1] = s; vis[s] = ++dfn; d[s] = 0;
    while(head <= tail){
        int u = q[head++];
        for(int i = hd[u]; i >= 0; i = nxt[i]){
            Edge &e = edg[i];
            if(vis[e.to] == dfn || e.cap <= e.flow) continue;
            vis[e.to] = dfn;
            d[e.to] = d[u] + 1;
            q[++tail] = e.to;
        }
    }
    return vis[t] == dfn;
}
int dfs(int x, int a){
    if(x == t || a == 0) return a;
    int flow = 0, f;
    for(int i = hd[x]; i >= 0; i = nxt[i]){
        Edge &e = edg[i];
        if(d[e.to] == d[x] + 1 && (f = dfs(e.to, min(a, e.cap - e.flow))) > 0){
            flow += f;
            e.flow += f;
            edg[i^1].flow -= f;
            a -= f;
            if(a == 0) break;
        }
    }
    return flow;
}
void init(){
    memset(hd, -1, sizeof(hd));
    scanf("%d%d%d", &n, &m, &k);
    for(int i = 1; i <= m; i++){
        scanf("%d%d", &h[i], &a[i][0]);
        for(int j = 1; j <= a[i][0]; j++)
            scanf("%d", &a[i][j]);
    }
}
void work(){
    s = 0, t = 1;
    insert(s, 3, inf);
    insert(2, t, inf);
    for(int TIME = 1; TIME <= 50; TIME++){
        insert(s, TIME*(n+2)+3, inf);
        insert(TIME*(n+2)+2, t, inf);
        for(int i = 1; i <= m; i++){
            int yd = (TIME - 1) % a[i][0] + 1;
            int td = (yd != a[i][0]) ? (yd + 1) : (1);
            insert((TIME-1)*(n+2)+a[i][yd]+3, TIME*(n+2)+a[i][td]+3, h[i]);
        }
        for(int i = 2; i <= n + 3; i++)
            insert((TIME-1)*(n+2)+i, TIME*(n+2)+i, inf);
        while(bfs()) ans += dfs(s, inf);
        if(ans >= k){
            printf("%d\n", TIME);
            return;
        }
    }
    puts("0");
}
int main(){
    freopen("prog813.in", "r", stdin);
    freopen("prog813.out", "w", stdout);
    init();
    work();
    return 0;
}