洛谷 P2936 [USACO09JAN]全流Total Flow
约翰总希望他的奶牛有足够的水喝,因此他找来了农场的水管地图,想算算牛棚得到的水的 总流量.农场里一共有N根水管.约翰发现水管网络混乱不堪,他试图对其进行简 化.他简化的方式是这样的:
两根水管串联,则可以用较小流量的那根水管代替总流量.
两根水管并联,则可以用流量为两根水管流量和的一根水管代替它们
当然,如果存在一根水管一端什么也没有连接,可以将它移除.
请写个程序算出从水井A到牛棚Z的总流量.数据保证所有输入的水管网络都可以用上述方法 简化.
输入输出格式
输入格式:
Line 1: A single integer: N
Lines 2..N + 1: Line i+1 describes pipe i with two letters and an integer, all space-separated: a_i, b_i, and F_i
输出格式:
Line 1: A single integer that the maximum flow from the well (‘A’) to the barn (‘Z’)
输入样例#1:
5
A B 3
B C 3
C D 5
D Z 4
B Z 6
输出样例#1:3
ISAP
所谓gap优化,就是dis=x的点都没有了以后,dis>x的点就都到不了了,
便可以直接结束。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
#if 0
Writers: Goes && G.S.M.
Have falled in love
#endif
const int M=888;
const int N=105;
int n,S=1,T=26;
struct ss{
int to,nex,va;
}edge[M*2];int head[N],ecnt=1;
void add(int x,int y,int va){
edge[++ecnt]=(ss){y,head[x],va};
edge[++ecnt]=(ss){x,head[y],0};
head[x]=ecnt-1;head[y]=ecnt;
}//网络流 式 建边
int dis[N],num[N],cur[N],fa[N];
void BFS(){//距离标记
for(int i=1;i<=80;i++) dis[i]=60;
queue<int> Q;Q.push(T);dis[T]=0;
while(!Q.empty()){
int pos=Q.front();Q.pop();
for(int i=head[pos];i;i=edge[i].nex)
if(edge[i^1].va&&dis[edge[i].to]>dis[pos]+1)
{
dis[edge[i].to]=dis[pos]+1;
Q.push(edge[i].to);
}
}return ;
}
int Out(){//增广(跟dinic一个道理)
int pos=T,f=0x3f3f3f3f;
while(pos!=S) pos=fa[pos],f=min(f,edge[cur[pos]].va);
pos=T;
while(pos!=S) pos=fa[pos],
edge[cur[pos]].va-=f,
edge[cur[pos]^1].va+=f;
return f;
}//就是这一波得到的能流过的量
int gswork(){
BFS();int pos=S,flow=0;//预处理距离标号
for(int i=1;i<=80;i++)//记距离的数量用于gap优化
num[dis[i]]++,cur[i]=head[i];
while(dis[S]<52){
if(pos==T) flow+=Out(),pos=S;
bool done =false;//标记是否前进
for(int i=cur[pos];i;i=edge[i].nex)
if(edge[i].va&&dis[pos]==dis[edge[i].to]+1)
{
fa[edge[i].to]=pos;cur[pos]=i;
pos=edge[i].to;done=true;i=0;//gsbreak;
}
if(!done)//如果往前走不动了
{
int m=51;
//更改逆向dis值
for(int i=head[pos];i;i=edge[i].nex)
if(edge[i].va) m=min(m,dis[edge[i].to]);
//gap优化
if(--num[dis[pos]]==0) break;
//更改相关值
num[dis[pos]=m+1]++;cur[pos]=head[pos];
if(pos!=S) pos=fa[pos];//回溯
}
}return flow;
}
inline void gsin(){
scanf("%d",&n);S=1,T=26;
for(int i=1;i<=n;i++){
char u,v;int x,y,va;
scanf("\n%c %c %d",&u,&v,&va);
x=u-64;y=v-64;add(x,y,va);
}return ;
}
int main()
{
gsin();
printf("%d",gswork());
return 0;
}