【CTSC1999】巴蜀2367家园

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 　　文件第一行为三个正整数 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 　　文件只有一个数，若问题有解，输出完成全部人员安全转移的时刻，否则输出0。

把每个点按时间拆点，但是不用一下子全建出来，而可以枚举时间，每次加入当前时间对应的点，判断是否能满流。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int s=100005,t=100006,oo=0x3f3f3f3f;
int p[25][25],r[25],hp[25],
fir[100010],ne[2000010],to[2000010],w[2000010],
f[100010],que[100010],
fa[25],
n,m,k,tot;
void add(int u,int v,int x)
{
    tot++;
    ne[tot*2]=fir[u];
    fir[u]=tot*2;
    to[tot*2]=v;
    w[tot*2]=x;
    ne[tot*2+1]=fir[v];
    fir[v]=tot*2+1;
    to[tot*2+1]=u;
    w[tot*2+1]=0;
}
void init()
{
    int i,j;
    scanf("%d%d%d",&n,&m,&k);
    for (i=1;i<=m;i++)
    {
        scanf("%d%d",&hp[i],&r[i]);
        for (j=0;j<r[i];j++)
        {
            scanf("%d",&p[i][j]);
            if (p[i][j]==-1) p[i][j]=n+1;
        }
    }
    add(s,0,k);
    add(n+1,t,oo);
}
int find(int x)
{
    return x==fa[x]?x:fa[x]=find(fa[x]);
}
bool check()
{
    int i,j;
    for (i=0;i<=n+1;i++)
      fa[i]=i;
    for (i=1;i<=m;i++)
      for (j=1;j<r[i];j++)
        fa[find(p[i][j-1])]=find(p[i][j]);
    return find(0)!=find(n+1);
}
bool bfs()
{
    int i,u,v,hd,tl;
    memset(f,0,sizeof(f));
    f[s]=1;
    hd=tl=1;
    que[1]=s;
    while (hd<=tl)
    {
        u=que[hd++];
        for (i=fir[u];i;i=ne[i])
          if (w[i]&&!f[v=to[i]])
          {
            f[v]=f[u]+1;
            que[++tl]=v;
          }
    }
    return f[t];
}
int dfs(int u,int lim)
{
    int i,v,x,ret=0;
    if (u==t) return lim;
    for (i=fir[u];i&&ret<lim;i=ne[i])
      if (w[i]&&f[v=to[i]]==f[u]+1)
      {
        x=dfs(v,min(w[i],lim-ret));
        ret+=x;
        w[i]-=x;
        w[i^1]+=x;
      }
    if (!ret) f[u]=0;
    return ret;
}
bool ok(int x)
{
    int i,j;
    add((n+2)*x+n+1,t,oo);
    for (i=0;i<=n;i++)
      add((n+2)*(x-1)+i,(n+2)*x+i,oo);
    for (i=1;i<=m;i++)
      add((n+2)*(x-1)+p[i][(x-1)%r[i]],(n+2)*x+p[i][x%r[i]],hp[i]);/*
    printf("\n%d:\n",x);
    for (i=0;i<=t;i++)
      for (j=fir[i];j;j=ne[j])
        if (w[j]) printf("%d->%d:%d\n",i,to[j],w[j]);*/
    while (bfs()) while (dfs(s,oo));
    return w[fir[s]]==0;
}
int main()
{
    int ans;
    init();
    if (check())
    {
        printf("0\n");
        return 0;
    }
    for (ans=1;;ans++)
      if (ok(ans))
      {
        printf("%d\n",ans);
        return 0;
      }
}