刷题集--家园

题意：由于人类对自然的疯狂破坏，人们意识到在大约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），目标是让所有的人尽快地全部转移到月球上。

暴力枚举答案（天数），然后更新图（算出当前每艘船从哪到哪，连边/每一点新建当前这天的点并接受上一天的流量代表人停留/原点向当天地球连边代表一直在地球上的人），当到达汇点（即月球）的人大于等于要求的，那么这一天就是最小天数。

代码：

#include<bits/stdc++.h>
using namespace std;
const int M=2000010,N=110;
int s,t,tot=-1;
int hd[M],nxt[M],las[M],to[M],dep[M],FL=0;
int fa[N],sz[N],mp[N][N],h[N];
void read(int &x)
{
    char c=getchar();x=0;int r=0;
    while(!isdigit(c))
    {
        if(c=='-')r=1;
        c=getchar();
    }
    while(isdigit(c))x=x*10+c-48,c=getchar();
    if(r)x=-x;
}
int find(int x)
{return x==fa[x]?x:fa[x]=find(fa[x]);}
void mg(int x,int y)
{
    x=find(x),y=find(y);
    if(x!=y)fa[x]=y;
    return;
}
void add(int u,int v,int w)
{
    nxt[++tot]=hd[u],to[tot]=v,las[tot]=w,hd[u]=tot;
    nxt[++tot]=hd[v],to[tot]=u,las[tot]=0,hd[v]=tot;
}
bool bfs()
{
    queue<int>q;int nw;
    memset(dep,0,sizeof dep);
    q.push(s),dep[s]=1;
    while(!q.empty())
    {
        nw=q.front(),q.pop();
        for(int i=hd[nw];i!=-1;i=nxt[i])
        {
            if(!dep[to[i]]&&las[i])
            {
                dep[to[i]]=dep[nw]+1;
                q.push(to[i]);
            }
        }
    }
    return dep[t];
}
int dfs(int pos,int flow)
{
    if(pos==t)return FL+=flow,flow;
    int tp;
    for(int i=hd[pos];i!=-1;i=nxt[i])
    {
        if(las[i]&&dep[to[i]]==dep[pos]+1)
        {
            tp=dfs(to[i],min(flow,las[i]));
            if(tp)
            {
                las[i]-=tp,las[i^1]+=tp;
                return tp;
            }
        }
    }
    return 0;
}
int dinic()
{
    int tp;
    while(bfs())
    {
        //for(int i=0;i<=n+5;i++)cur[i]=hd[i];
        while(dfs(s,2e9));
    }
}
int n,m,k;
int main()
{
    memset(hd,-1,sizeof hd);
    read(n),read(m),read(k);
    for(int i=0;i<=n+1;i++)fa[i]=i;
    for(int i=1;i<=m;i++)
    {
        read(h[i]),read(sz[i]);
        for(int j=0;j<sz[i];j++)
        {
            read(mp[i][j]);
            if(mp[i][j]==-1)mp[i][j]+=n+2;
            if(j)mg(mp[i][j-1],mp[i][j]);
        }
    }
    if(find(0)!=find(n+1))return puts("0"),0;
    s=n+2,t=n+1;
    for(int ans=1;;++ans)
    {
        add(s,(ans-1)*(n+5),2e9);
        for(int i=1;i<=m;i++)
        {
            int nw,nx;
            nw=mp[i][(ans-1)%sz[i]];
            nx=mp[i][ans%sz[i]];
            if(nw!=t)nw+=(ans-1)*(n+5);
            if(nx!=t)nx+=ans*(n+5);
            add(nw,nx,h[i]);
        }
        dinic();
        if(FL>=k){return printf("%d\n",ans),0;}
        for(int i=0;i<=n;i++)add(i+(ans-1)*(n+5),i+ans*(n+5),2e9);
    }
}