https://vjudge.net/problem/POJ-1459
题解转载自：優YoU http://user.qzone.qq.com/289065406/blog/1299339754
解题思路
多源多汇最大流问题
题目给出很多都是废话，特别是符号s(u)，d(u)，Con还有那条公式都别管，混淆视听
难点在于构图
电站p(u)均为源点，用户c(u)均为汇点，中转站当普通点处理
第一个误区是例图， 结点 和 边 都有x/y（流量和容量），这个很容易使人产生矛盾（因为学习最大流问题是，只有 边 才有流量和容量。
但是不难发现，题目所给的例图中有多个源点，多个汇点，多个普通点，只有源点和汇点才标有 x/y，普通点没有标x/y，而且所给出的所有边都有x/y。 这无疑在促使我们对图做一个变形： 建议一个超级源s，一个超级汇t，使s指向所有源点，并把源点的 容量y 分别作为这些边的 容量，使所有汇点指向t，并把汇点的容量y分别作为这些边的 容量，然后本来是源点和汇点的点，全部变为普通点。这样就把“多源多汇最大流”变形为“单源单汇最大流”问题。
第二个误区就是流量值。 学习最大流问题时，会发现边上的流量值是给定初始值的，但是这题的输入只有容量，没有流量，很多人马上感觉到无从入手。其实边上的流量初始值为多少都没有所谓，解最大流需要用到的只有容量。但是一般为了方便起见， 会把所有边的流量初始化为0。这样做有一个最大的好处，就是可以回避 反向弧 的存在，这个下面详细叙述。
本题中要注意的是：
1、 如果输入中，某一点上有环，就无视掉。环是否存在不影响最终结果。
2、 一般两点之间都是单边，一旦存在双边也没有问题，因为由定义知两个方向的容量一定相等（其实不相等也无妨，因为其中有一条为 反向弧，前面已经提到 反向弧 是可以直接回避、无视的，因此反向弧上的容量为多少就不重要了），而且在寻找增广路的标号过程中，搜索的是未标号的点，就是说(u,v)这条弧即使是双向的，但一旦从u到达v后，就不能回头了，因为两者都被标记了，即另外一条弧就不起任何作用了。
下面详细说说为什么能够回避反向弧。
首先需要明确，任意一个点j上记录的信息有：
1、 寻找增光路时，除超级源s外，增广路上任一点j都有一个唯一的前驱i（i被记录在j）
2、 min{从i到j的容流差,l(vi)}
3、 构图时，除超级汇t外，图上任一点j都会直接指向一些点(这些点作为后继点，同在记录在j)
从这个特点可以知道，从超级源开始寻找增广路时，万一遇到双向边，正向弧，反向弧自动被回避。万一遇到单向边，如果是非饱和正向弧，就会继续走下去；如果是反向弧，这条弧必然是 零弧（每条边初始化流量均为0），定义知如果增广路有反向弧，它必须要是非零弧，而且由于反向弧每次都不会经过，所以在改进增广路时反向弧上的流量也不会被改变，永远为0，也就与最终结果无关了
最后当无法寻找增广路时，最大流就是与超级源s直接关联的边上的 流量之和
#include<algorithm> #include<iostream> #include<cstring> #include <sstream> #include<vector> #include<cstdio> #include<cmath> #include<queue> #include<stack> #include<map> #include<set> #include <string> #define MP make_pair #define ll long long #define CLR(a, b) memset(a, b, sizeof(a)) const ll mod=2147493647; using namespace std; const double eps=1e-8; const int MAXN=110; const int INF=0x3f3f3f3f; int maze[MAXN][MAXN]; int gap[MAXN],dis[MAXN],pre[MAXN],cur[MAXN]; int sap(int st,int ed,int nodenum) { memset(cur,0,sizeof(cur)); memset(dis,0,sizeof(dis)); memset(gap,0,sizeof(gap)); int u=pre[st]=st,maxflow=0,aug=-1; gap[0]=nodenum; while(dis[st]<nodenum) { loop: for(int v=cur[u]; v<nodenum; v++) { if(maze[u][v]&& dis[u]==dis[v]+1) { if(aug==-1||aug>maze[u][v]) aug=maze[u][v]; pre[v]=u; u=cur[u]=v; if(v==ed) { maxflow+=aug; for(u=pre[u]; v!=st; v=u,u=pre[u]) { maze[u][v]-=aug; maze[v][u]+=aug; } aug=-1; } goto loop; } } int mindis=nodenum-1; for(int v=0; v<nodenum; v++) { if(maze[u][v]&&mindis>dis[v]) { cur[u]=v; mindis=dis[v]; } } if((--gap[dis[u]])==0) break; gap[dis[u]=mindis+1]++; u=pre[u]; } return maxflow; } int main() { int n,np,nc,m; int u,v,w; while(cin>>n>>np>>nc>>m) { int s=n,t=n+1; char ch; memset(maze,0,sizeof(maze)); for(int i=0; i<m; i++) { while(cin>>ch,ch!='('); cin>>u>>ch>>v>>ch>>w; maze[u][v]=w; } for(int i=0; i<np; i++) { while(cin>>ch,ch!='('); cin>>u>>ch>>w; maze[s][u]=w; } for(int i=0; i<nc; i++) { while(cin>>ch,ch!='('); cin>>u>>ch>>w; maze[u][t]=w; } cout<<sap(n,n+1,n+2)<<endl; } return 0; }