【CTSC 1999】家园

题目描述
由于人类对自然资源的消耗，人们意识到大约在 2300 年之后，地球就不能再居住了。于是在月球上建立了新的绿地，以便在需要时移民。令人意想不到的是，2177 年冬由于未知的原因，地球环境发生了连锁崩溃，人类必须在最短的时间内迁往月球。现有 n 个太空站位于地球与月球之间，且有 m 艘公共交通太空船在其间来回穿梭。每个太空站可容纳无限多的人，而每艘太空船 i 只可容纳 H[i]个人。每艘太空船将周期性地停靠一系列的太空站，例如：(1，3，4)表示该太空船将周期性地停靠太空站 134134134…。每一艘太空船从一个太空站驶往任一太空站耗时均为 1。人们只能在太空船停靠太空站(或月球、地球)时上、下船。初始时所有人全在地球上，太空船全在初始站。试设计一个算法，找出让所有人尽快地全部转移到月球上的运输方案。
对于给定的太空船的信息，找到让所有人尽快地全部转移到月球上的运输方案。
输入格式：
第 1 行有 3 个正整数 n（太空站个数），m（太空船个数）和 k（需要运送的地球上的人的个数）。其中 n<=13 m<=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…等正点时刻才能上下太空船。
输出格式：
程序运行结束时，将全部人员安全转移所需的时间输出。如果问题
无解，则输出 0。

拆点！
建图很重要的套路。
一个点有多个不同信息一般要拆点。
这道题边就不是一成不变的，所以要拆点（对时间分层）。
把每个天空站拆成若干个，每个表示一个时刻。
然后二分答案，建图，跑最大流看看是否达到k。
也可以枚举答案，每次在上次的残余网络跑最大流，好像更快？
是否有解可以用并查集判，我为了方便如果枚举答案过大就直接break掉。

#include<cstdio>
#include<iostream>
#include<queue>
#include<cstring>
#define N 17
#define S 10001
#define T 10002
using namespace std;
struct node{
    int to,next,flow,from;
}ed[1000100];
int head[11000],deep[11000],size=-1,ind[110][110],a[110][110],curr[11000],max_flow,cnt[50],p[50],n,m,k;
void add(int from,int to,int flow)
{
    size++;
    ed[size].to=to;
    ed[size].from=from;
    ed[size].flow=flow;
    ed[size].next=head[from];
    head[from]=size;
}
bool bfs()
{
    queue <int> q;
    memset(deep,0,sizeof(deep));
    q.push(S);
    deep[S]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        if(u==T) return true;
        for(int i=head[u];~i;i=ed[i].next)
        {
            int v=ed[i].to;
            if(!deep[v]&&ed[i].flow)
            {
                q.push(v);
                deep[v]=deep[u]+1;
            }
        }
    }
    return false;
}
int dfs(int u,int cur)
{
    int rest=cur;
    if(u==T) return cur;
    for(int i=curr[u];~i;i=ed[i].next)
    {
        int v=ed[i].to;
        if(deep[v]==deep[u]+1&&ed[i].flow)
        {
            int new_flow=dfs(v,min(rest,ed[i].flow));
            ed[i].flow-=new_flow;
            ed[i^1].flow+=new_flow;
            rest-=new_flow;
            if(ed[i].flow) curr[u]=i;
            if(!rest) break;
        }
    }
    if(cur==rest) deep[u]=-1;
    return cur-rest;
}
int dinic()
{
    while(bfs())
    {
        for(int i=0;i<=10010;i++) curr[i]=head[i];
        max_flow+=dfs(S,1e9);
    }
    return max_flow;
}
void graph(int x)
{
    add(S,x*N,1e9);
    add(x*N,S,0);
    add(n+1+x*N,T,1e9);
    add(T,n+1+x*N,0);
    for(int i=1;i<=n;i++)
    {
        add(i+(x-1)*N,i+x*N,1e9);
        add(i+x*N,i+(x-1)*N,0);
    }
    for(int i=1;i<=m;i++)
    {
        int t=x%cnt[i];
        if(t!=0)
        {
            add(a[i][t-1]+(x-1)*N,a[i][t]+x*N,p[i]);
            add(a[i][t]+x*N,a[i][t-1]+(x-1)*N,0);
        }
        else
        {
            add(a[i][cnt[i]-1]+(x-1)*N,a[i][0]+x*N,p[i]);
            add(a[i][0]+x*N,a[i][cnt[i]-1]+(x-1)*N,0);
        }
    }
}
int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&p[i],&cnt[i]);
        for(int j=0;j<cnt[i];j++)
        {
            scanf("%d",&a[i][j]);
            if(a[i][j]==-1) a[i][j]=n+1;
        }
    }
    int ans=1;

    add(S,0,1e9);
    add(0,S,0);
    graph(1);
    while(dinic()<k&&ans<=500)
        graph(++ans);
    printf("%d",ans<500?ans:0);
    return 0;
}