[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天有飞船经过，则连一条飞船容量的边。这样我们就能解决这个问题了。

 //It is made by wfj_2048~
 #include <algorithm>
 #include <iostream>
 #include <complex>
 #include <cstring>
 #include <cstdlib>
 #include <cstdio>
 #include <vector>
 #include <cmath>
 #include <queue>
 #include <stack>
 #include <map>
 #include <set>
 #define inf (1<<30)
 #define il inline
 #define RG register
 #define ll long long
 #define c(t,i) ( t*(n+2)+i )
 #define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)

 using namespace std;

 struct edge{ int nt,to,flow,cap; }g[];

 int head[],d[],q[],p[],r[],s[][],n,m,k,S,T,flow,num=;

 il int gi(){
     RG int x=,q=; RG char ch=getchar(); while ((ch<'' || ch>'') && ch!='-') ch=getchar();
     if (ch=='-') q=-,ch=getchar(); while (ch>='' && ch<='') x=x*+ch-,ch=getchar(); return q*x;
 }

 il void insert(RG int from,RG int to,RG int cap){ g[++num]=(edge){head[from],to,,cap},head[from]=num; return; }

 il int bfs(RG int S,RG int T){
     memset(d,,sizeof(d)); RG int h=,t=; q[t]=S,d[S]=;
     while (h<t){
     RG int x=q[++h];
     for (RG int i=head[x];i;i=g[i].nt){
         RG int v=g[i].to;
         if (!d[v] && g[i].cap>g[i].flow){
         q[++t]=v,d[v]=d[x]+;
         if (v==T) return ;
         }
     }
     }
     return ;
 }

 il int dfs(RG int x,RG int T,RG int a){
     if (x==T || !a) return a; RG int f,flow=;
     for (RG int i=head[x];i;i=g[i].nt){
     RG int v=g[i].to;
     if (d[v]==d[x]+ && g[i].cap>g[i].flow){
         f=dfs(v,T,min(a,g[i].cap-g[i].flow));
         g[i].flow+=f,g[i^].flow-=f;
         flow+=f,a-=f; if (!a) return flow;
     }
     }
     if (!flow) d[x]=-; return flow;
 }

 il int maxflow(RG int S,RG int T){ while (bfs(S,T)) flow+=dfs(S,T,inf); return flow; }

 il void work(){
     n=gi(),m=gi(),k=gi(); S=,T=;
     for (RG int i=;i<=m;++i){
     p[i]=gi(),r[i]=gi();
     for (RG int j=;j<=r[i];++j){
         s[i][j]=gi(); if (s[i][j]==) s[i][j]=n+;
         if (s[i][j]==-) s[i][j]=n+;
     }
     }
     insert(S,n+,k),insert(n+,S,);
     for (RG int ans=;ans<=(k+)*(n+);++ans){
     for (RG int i=;i<=n+;++i) insert(c((ans-),i),c(ans,i),inf),insert(c(ans,i),c((ans-),i),);
     for (RG int i=;i<=m;++i){
         insert(c((ans-),s[i][(ans-)%r[i]+]),c(ans,s[i][ans%r[i]+]),p[i]);
         insert(c(ans,s[i][ans%r[i]+]),c((ans-),s[i][(ans-)%r[i]+]),);
     }
     insert(c(ans,(n+)),T,inf),insert(T,c(ans,(n+)),);
     if (maxflow(S,T)==k){ printf("%d\n",ans); return; }
     }
     printf(""); return;
 }

 int main(){
     File("star");
     work();
     return ;
 }