题意:由于人类对自然的疯狂破坏,人们意识到在大约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);
}
}