【网络流24题】家园

Luogu 2756 CodeVs 1034

题意：

地球和月球之间有n个无限容量的空间站，给出m个太空船的容量和行驶路线。问将K个人从地球运送到月球的最小时间。

题解：

这题很符合网络流最初的设定，人在地球月球空间站太空船之间流动，流动的过程中受到流量限制(太空船容量)

对于每个时间点给每个太空船和空间站建一个点，然后按照题意连边即可。

代码：

#include<cstdio>
#include<cstring>
#include<vector>
#define MAXN 20000+30
#define INF 0x3fffffff

using namespace std;


int m,n,K,Du[MAXN],dis[MAXN],dl[23333];
vector <int> f[MAXN],FB[MAXN],Flow[MAXN];
int cur[MAXN];

int a[30][30];
int Load[30],Num[30];

void ADD(int x,int y,int z)
{
	Du[x]++;Du[y]++;
	f[x].push_back(y);Flow[x].push_back(z);FB[x].push_back(Du[y]-1);
	f[y].push_back(x);Flow[y].push_back(0);FB[y].push_back(Du[x]-1);
}
bool BFS(int Begin,int End)
{
	int t=0,w=1,x,X;
	memset(dis,0xff,sizeof(dis));
	dis[Begin]=0;dl[1]=Begin;
	do
	{
		x=dl[++t];
		for(int i=0;i<Du[x];i++)
		{
			X=f[x][i];
			if(Flow[x][i]<=0||dis[X]>=0) continue;//不连通或已访问 
			dis[X]=dis[x]+1;dl[++w]=X;
		}
	}while(t<w);
	if(dis[End]>0) return true;
	else return false;
}
int Find(int x,int MFLOW,int y)
{
	if(x==y) return MFLOW;
	int X,h; 
	for(int i=cur[x];i<Du[x];i++)
	{
		cur[x]=i;
		X=f[x][i];
		if(Flow[x][i]>0&&dis[x]+1==dis[X]&&(h=Find(X,min(MFLOW,Flow[x][i]),y)))
		{
			Flow[x][i]-=h;
			Flow[X][FB[x][i]]+=h;
			return h;
		}
	}
	return 0;
}
int Solve(int x,int y)
{
	int X,Ans=0;
	while(1)
	{
		if(!BFS(x,y)) break;
		memset(cur,0,sizeof(cur));
		while((X=Find(x,0x7fffffff,y)))
			Ans+=X;
	}
		
	return Ans;
}

int main()
{
	scanf("%d%d%d",&n,&m,&K);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&Load[i],&Num[i]);
		for(int j=0;j<Num[i];j++)
			scanf("%d",&a[i][j]);
	}
	int x,S,T;
	for(int I=0;I<=20000/(n+m);I++)
	{
		S=0;T=(I+1)*(m+n)+1;
		
		for(int i=S;i<=T;i++)
			Du[i]=0,f[i].clear(),FB[i].clear(),Flow[i].clear();
		
		for(int i=0;i<=I;i++)
		{
			for(int j=1;j<=m;j++)
			{
				x=a[j][i%Num[j]];
				if(x==0) ADD(S,i*(n+m)+j,INF);
				else if(x==-1) ADD(i*(n+m)+j,T,INF);
				else ADD(i*(n+m)+j,i*(n+m)+m+x,INF),ADD(i*(n+m)+m+x,i*(n+m)+j,INF);
				
				if(i!=I)
					ADD(i*(n+m)+j,(i+1)*(n+m)+j,Load[j]);
			}
			for(int j=1;j<=n;j++)
				if(i!=I)
					ADD(i*(n+m)+m+j,(i+1)*(n+m)+m+j,INF);
		}
		if(Solve(S,T)>=K)
		{
			printf("%d",I);
			return 0;
		}
	}
	printf("%d",0);
	return 0;
}