[WIKIOI1034]家园(网络流)

题目描述 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。

别人的代码, 这个自己挊了好几次都挊不出来, 求助神犇去了.

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define inf 100000000
#define maxn 100000
#define maxe 1000000
struct node{
int from,to,flow;
node* next;
};
class dinic{
private:
node *mat[maxn],edge[maxe],*temp[maxn],*cur[maxn];
int numptr,n,flow,que[maxn],leval[maxn],path[maxn];
bool used[maxn];
bool bfs(int s,int t)
{
int head=0,tail=0,i,u;
node *ptr;
for(i=0;i<n;i++)leval[i]=-1;
que[tail++]=s;
leval[s]=1;
while(head!=tail)
{
u=que[head++];
for(ptr=mat[u];ptr;ptr=ptr->next)
{
if(leval[ptr->to]==-1&&ptr->flow>0)
{
leval[ptr->to]=leval[ptr->from]+1;
que[tail++]=ptr->to;
if(ptr->to==t)
{
return true;
}
}
}
}
return false;
}
public:
void init(int tempn)
{
flow=0;
int i;
n=tempn;
numptr=0;
for(i=0;i<n;i++)mat[i]=NULL;
}
void addedge(int from,int to,int tflow)
{
edge[numptr].from=from;
edge[numptr].to=to;
edge[numptr].flow=tflow;
edge[numptr].next=mat[from];
mat[from]=&edge[numptr++];
edge[numptr].from=to;
edge[numptr].to=from;
edge[numptr].flow=0;
edge[numptr].next=mat[to];
mat[to]=&edge[numptr++];
}
void solve(int s,int t)
{
int p,i,u,now,tempp;
node* ptr;
bool f;
while(bfs(s,t))
{
for(i=0;i<n;i++)temp[i]=mat[i],used[i]=true;
p=1;
path[p]=s;
while(p)
{
u=path[p];
if(u==t)
{
now=inf;
for(i=1;i<p;i++)
if(cur[path[i]]->flow<now)
now=cur[path[i]]->flow;
flow+=now;
for(i=1;i<p;i++)
{
cur[path[i]]->flow-=now;
ptr=edge+((cur[path[i]]-edge)^1);
ptr->flow+=now;
if(cur[path[i]]->flow==0)tempp=i;
}
p=tempp;
}
else{
f=false;
for(ptr=temp[u];ptr;ptr=ptr->next)
{
if(used[ptr->to]&&ptr->flow>0&&leval[ptr->from]+1==leval[ptr->to])
{
cur[u]=ptr;
temp[u]=ptr->next;
f=true;
path[++p]=ptr->to;
break;
}
}
if(!f)
{
p--;
used[u]=false;
}
}
}
}
}
int question()
{
return flow;
}
};
dinic M;
int a[50][50];
int p[50],num[50];
int main()
{
//freopen("in","r",stdin);
int i,j,k,n,m,temp1,temp2,nn;
while(scanf("%d%d%d",&n,&m,&k)!=EOF)
{
nn=n;
M.init(maxn);
for(i=0;i<m;i++)
{
scanf("%d%d",&p[i],&num[i]);
for(j=0;j<num[i];j++)
{
scanf("%d",&a[i][j]);
if(a[i][j]==-1)a[i][j]=maxn-1;
}
}
for(i=0;;i++)
{
for(j=0;j<m;j++)
{
temp1=a[j][i%num[j]];
if(temp1==maxn-1)continue;
temp2=a[j][(i+1)%num[j]];
if(temp2==0)continue;
if(temp1!=0)temp1+=i*nn;
if(temp2!=maxn-1)temp2+=(i+1)*nn;
M.addedge(temp1,temp2,p[j]);
}
if(i!=0)
{
for(j=1;j<=n;j++)
{
M.addedge(j+(i-1)*nn,j+i*nn,inf);
}
}
M.solve(0,maxn-1);
//cout<<i<<" "<<M.question()<<endl;
if(M.question()>=k)
{
cout<<i+1<<endl;
break;
}
if(i>50&&M.question()==0)
{
cout<<0<<endl;
break;
}
}
}
return 0;
}

相当纠结....