首页 > 代码库 > [CTSC1999]【网络流24题】星际转移

[CTSC1999]【网络流24题】星际转移

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

正解:网络判定+最大流。

可以二分答案,不过不如枚举时间,直接动态加边,然后在残量网络中增广。从源点向第0天的地球连容量为k的边,从第i天的星球向第i+1的星球连容量为inf的边,如果第i天到第i+1天有飞船经过,则连一条飞船容量的边。这样我们就能解决这个问题了。

 

 1 //It is made by wfj_2048~
 2 #include <algorithm>
 3 #include <iostream>
 4 #include <complex>
 5 #include <cstring>
 6 #include <cstdlib>
 7 #include <cstdio>
 8 #include <vector>
 9 #include <cmath>
10 #include <queue>
11 #include <stack>
12 #include <map>
13 #include <set>
14 #define inf (1<<30)
15 #define il inline
16 #define RG register
17 #define ll long long
18 #define c(t,i) ( t*(n+2)+i )
19 #define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
20 
21 using namespace std;
22 
23 struct edge{ int nt,to,flow,cap; }g[1000010];
24 
25 int head[100010],d[100010],q[100010],p[100010],r[100010],s[1010][1010],n,m,k,S,T,flow,num=1;
26 
27 il int gi(){
28     RG int x=0,q=1; RG char ch=getchar(); while ((ch<0 || ch>9) && ch!=-) ch=getchar();
29     if (ch==-) q=-1,ch=getchar(); while (ch>=0 && ch<=9) x=x*10+ch-48,ch=getchar(); return q*x;
30 }
31 
32 il void insert(RG int from,RG int to,RG int cap){ g[++num]=(edge){head[from],to,0,cap},head[from]=num; return; }
33 
34 il int bfs(RG int S,RG int T){
35     memset(d,0,sizeof(d)); RG int h=0,t=1; q[t]=S,d[S]=1;
36     while (h<t){
37     RG int x=q[++h];
38     for (RG int i=head[x];i;i=g[i].nt){
39         RG int v=g[i].to;
40         if (!d[v] && g[i].cap>g[i].flow){
41         q[++t]=v,d[v]=d[x]+1;
42         if (v==T) return 1;
43         }
44     }
45     }
46     return 0;
47 }
48 
49 il int dfs(RG int x,RG int T,RG int a){
50     if (x==T || !a) return a; RG int f,flow=0;
51     for (RG int i=head[x];i;i=g[i].nt){
52     RG int v=g[i].to;
53     if (d[v]==d[x]+1 && g[i].cap>g[i].flow){
54         f=dfs(v,T,min(a,g[i].cap-g[i].flow));
55         g[i].flow+=f,g[i^1].flow-=f;
56         flow+=f,a-=f; if (!a) return flow;
57     }
58     }
59     if (!flow) d[x]=-1; return flow;
60 }
61 
62 il int maxflow(RG int S,RG int T){ while (bfs(S,T)) flow+=dfs(S,T,inf); return flow; }
63 
64 il void work(){
65     n=gi(),m=gi(),k=gi(); S=100001,T=100002;
66     for (RG int i=1;i<=m;++i){
67     p[i]=gi(),r[i]=gi();
68     for (RG int j=1;j<=r[i];++j){
69         s[i][j]=gi(); if (s[i][j]==0) s[i][j]=n+1;
70         if (s[i][j]==-1) s[i][j]=n+2;
71     }
72     }
73     insert(S,n+1,k),insert(n+1,S,0);
74     for (RG int ans=1;ans<=(k+2)*(n+2);++ans){
75     for (RG int i=1;i<=n+2;++i) insert(c((ans-1),i),c(ans,i),inf),insert(c(ans,i),c((ans-1),i),0);
76     for (RG int i=1;i<=m;++i){
77         insert(c((ans-1),s[i][(ans-1)%r[i]+1]),c(ans,s[i][ans%r[i]+1]),p[i]);
78         insert(c(ans,s[i][ans%r[i]+1]),c((ans-1),s[i][(ans-1)%r[i]+1]),0);
79     }
80     insert(c(ans,(n+2)),T,inf),insert(T,c(ans,(n+2)),0);
81     if (maxflow(S,T)==k){ printf("%d\n",ans); return; }
82     }
83     printf("0"); return;
84 }
85 
86 int main(){
87     File("star");
88     work();
89     return 0;
90 }

 

[CTSC1999]【网络流24题】星际转移