家园

题目描述 Description

由于人类对自然的疯狂破坏，人们意识到在大约2300年之后，地球不能再居住了，于是在月球上建立了新的绿地，以便在需要时移民。令人意想不到的是，2177年冬由于未知的原因，地球环境发生了连锁崩溃，人类必须在最短的时间内迁往月球。
现有n个太空站处于地球与月球之间（编号1..n），m艘公共交通太空船在其中来回穿梭，每个太空站Si可容纳无限的人，每艘太空船pi只可容纳Hpi人。对于每一艘太空船pi，将周期性地停靠一系列的太空站（Si1,Si2…Sir），如：（1，3，4）表示停靠太空站1 3 4 1 3 4 1 3 4 …。 任一艘太空船从任一个太空站驶往另一个任意的太空站耗时为1。人只能在太空船停靠太空站（或地球、月球）时上船或下船。初始时的人全在地球上，太空船全在初始站（太空船pi处于Si1），目标是让所有的人尽快地全部转移到月球上。

输入描述 Input Description

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

输出描述 Output Description

文件只有一个数，若问题有解，输出完成全部人员安全转移的时刻，否则输出0。

样例输入 Sample Input

2 2 1 

1 3 0 1 2

1 3 1 2 –1

样例输出 Sample Output

5

数据范围及提示 Data Size & Hint

 1<=m<=13, 1<=n<=20, 1<=k<=50。

此题又是最大流模型，但是区别是这题的图是变化的。

以时间为序，对于每一个时刻，我们都得到一张图，然后我们对此使用dinic算法求最大流。

可以发现一个站可以有多个状态，我们以(i,j)表示i站在时刻j的状态，很显然(i,j-1)状态可以流向(i,j)状态，且容量为无穷大。

假设飞船k,j-1时刻到j时刻从s站到了t站，并且飞船的承载力为flow，则我们在途中应添加边(s,j-1)-->(t,j)，容量为flow。

所以对于时间每增加1，原图会增加n+2个点，也就是每个站会产生一个新状态，对于新图我们用dinic求最大流直到没有增广路，也就是从上一个时刻过渡到该时刻，又有多少人到达月球。

这里假定测试数据最大时间是小于等于100的，当时间超过100就认为无法全部到达月球。

#include<cstdio>
#include<iostream>
#include<cstring>
#define Maxn 3100
using namespace std;

int flow[20],num[20],path[20][30];
const int inf=0x3f3f3f3f;
const int step=30;
struct line{
    int to,next,cap;
}p[Maxn*Maxn];
int head[Maxn];
int q[Maxn];
int d[Maxn];
int tot;
int src,t;
int n,m,k;
void addedge(int a,int b,int c){
    p[tot].to=b;
    p[tot].next=head[a];
    p[tot].cap=c;
    head[a]=tot++;
}
void insert(int a,int b,int c){
    addedge(a,b,c);
    addedge(b,a,0);
}
bool bfs(){
    memset(d,-1,sizeof d);
    int s=0,e=-1;
    q[++e]=src;
    d[src]=0;
    while(s<=e){
        int u=q[s++];
        for(int i=head[u];i!=-1;i=p[i].next){
            int v=p[i].to;
            if(d[v]==-1&&p[i].cap){
                d[v]=d[u]+1;
                q[++e]=v;
            }
        }
    }
    return d[t]!=-1;
}
int dfs(int u,int alpha){
    if(u==t) return alpha;
    int w,used=0;
    for(int i=head[u];i!=-1&&used<alpha;i=p[i].next){
        int v=p[i].to;
        if(p[i].cap&&d[v]==d[u]+1){
            w=dfs(v,min(alpha-used,p[i].cap));
            used+=w;
            p[i].cap-=w;
            p[i^1].cap+=w;
        }
    }
    if(!used) d[u]=-1;
    return used;
}
int dinic(int time){
    int ans=0;
    src=0,t=n+1+time*step;
    while(bfs())
        ans+=dfs(src,inf);
    return ans;
}
void add(int time){
    for(int i=0;i<n+2;i++)
        insert(i+(time-1)*step,i+time*step,inf);
    for(int i=0;i<m;i++)
        insert(path[i][(time-1)%num[i]]+(time-1)*step,path[i][time%num[i]]+time*step,flow[i]);
}
int main()
{
    int res;
    tot=0;
    memset(head,-1,sizeof head);
    cin>>n>>m>>k;
    for(int i=0;i<m;i++){
        cin>>flow[i]>>num[i];
        for(int j=0;j<num[i];j++){
            cin>>path[i][j];
            if(path[i][j]==-1)
                path[i][j]=n+1;
        }
    }
    for(res=1;res<=100;res++){
        add(res);
        k-=dinic(res);
        if(k<=0) break;
    }
    if(res<=100) cout<<res<<endl;
    else cout<<0<<endl;
	return 0;
}