POJ 1895 分层图网络流+输出路径

时间:2022-12-14 07:26:04

题意:
题目描述:在公元3141年人类的足迹已经遍布银河系。为了穿越那巨大的距离,人类发明了一种名为超时空轨道的技术。超时空轨道是双向的,连接两个星系,穿越轨道需要一天的时间。然而这个轨道只能同时给一艘飞船使用,也就是说,每条轨道每天只能有一艘飞船穿越。现在IBM公司要把K(K≤50)台超级计算机从地球运到Eisiem星系去,由于这些超级计算机个头巨大,一台计算机就要用一艘飞船来运。现在人类能够到达N(N≤50)个星系,拥有M(M≤200)条超时空轨道,太阳系的编号为S,Eisiem星系的编号为T。你需要求出至少需要几天才能将这些超级计算机全部运到目的地。注意,IBM公司是非常NB的公司,所有的超时空轨道都会优先给IBM公司使用。
思路:
POJ 1895 分层图网络流+输出路径
POJ 1895 分层图网络流+输出路径
From lydrainbowcat

//By SiriusRen
#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 5555
#define M 444444
int n,m,k,from[N],to[N],jy,rec[N][55],rect[N][55],inf=0x3f3f3f3f;
int first[N],next[M],v[M],w[M],tot,T=5551,day,vis[N],ans;
void Add(int x,int y,int z){w[tot]=z,v[tot]=y,next[tot]=first[x],first[x]=tot++;}
void add(int x,int y,int z){Add(x,y,z),Add(y,x,0);}
bool tell(){
memset(vis,-1,sizeof(vis)),vis[0]=0;
queue<int>q;q.push(0);
while(!q.empty()){
int t=q.front();q.pop();
for(int i=first[t];~i;i=next[i])
if(w[i]&&vis[v[i]]==-1)
vis[v[i]]=vis[t]+1,q.push(v[i]);
}return vis[T]!=-1;
}
int zeng(int x,int y){
if(x==T)return y;
int r=0;
for(int i=first[x];~i&&y>r;i=next[i])
if(w[i]&&vis[v[i]]==vis[x]+1){
int t=zeng(v[i],min(y-r,w[i]));
w[i]-=t,w[i^1]+=t,r+=t;
}
if(!r)vis[x]=-1;
return r;
}
void dfs(int x,int Day){
for(int i=first[x];~i;i=next[i])
if(w[i^1]){
w[i]++,w[i^1]--;
if(v[i]==T)return;
if(v[i]==x+n)dfs(v[i],Day+1);
else{
rec[Day][0]++;
rec[Day][rec[Day][0]]=jy;
rect[Day][rec[Day][0]]=(v[i]-1)%n+1;
dfs(v[i],Day+1);
}
break;
}
}
int main(){
memset(first,-1,sizeof(first));
scanf("%d%d%d%d%d",&n,&m,&k,&from[0],&to[0]);
for(int i=1;i<=m;i++){
scanf("%d%d",&from[i],&to[i]);
}add(0,from[0],k);
day=0;
while(++day){
for(int i=1;i<=m;i++)add(from[i]+(day-1)*n,to[i]+day*n,1),add(to[i]+(day-1)*n,from[i]+day*n,1);
for(int i=1;i<=n;i++)add(i+(day-1)*n,i+day*n,inf);
add(to[0]+day*n,T,inf);
while(tell())while(jy=zeng(0,inf))ans+=jy;
if(ans==k){printf("%d\n",day);break;}
}
for(jy=k;jy;jy--)dfs(0,0);
for(int i=1;i<=day;i++){
printf("%d ",rec[i][0]);
for(int j=1;j<=rec[i][0];j++)
printf("%d %d ",rec[i][j],rect[i][j]);
putchar('\n');
}
}

POJ 1895 分层图网络流+输出路径