[数学建模]网络最大流(自用,可能说的不明白)

FF算法(Ford-Fulkerson)
也称标号法

司书写Matlab可以直接用工具箱,所以只给了lingo程序

那么就只记记思路:

我们如果手算,步骤是这样的:
1.选定一条S到T的路,计算${\delta }_i$(${\delta }_i$是这条路中最小的边权)
2.更新这条路上每条边的权值(调整成参流)
比如:S$\to$A$\to$C$\to$T这条路,$\delta$=2,更新完之后SA为1,AC为0,CT为0
3.一二步重复执行,直到选的任一路的$\delta$都为0
4.这个网络的最大流就是$\Sigma$${\delta }_i$

那么用程序写出来呢?因为我们开始任意选的路会选错,所以不妨加条"后悔路"
也就是下图这样:

加入之后我们每次更新要注意把后悔路也更新了,更新方法是减$\delta$.
怎么理解呢?如果我们先选了S$\to$A$\to$C$\to$T这条路,CA就会被更新为2,TC被更新为2,这样相当于有"退回"的可能了.可以选S$\to$B$\to$C$\to$A$\to$T这种道路.

附C++代码:

#define readl(a) scanf("%lld",&a)
#define reads(a) scanf("%s",&a)
#define readc(a) scanf("%c",&a)
#define pb push_back
#define mem(a) memset(a,o,sizeof(a))
#define Buff ios::sync_with_stdio(false)
typedef long long ll;
using name space std;
const int INF = 2e9+7;
const int N = 1e5+7;
const int M = 1e6+7;
const int base = 100;

struct edge//边的结构体
{
	int to,nex,cap;//指向点,下一条边,容量
	edge(int to=0,int nex=0,int cap=0):to(to),nex(nex),cap(cap){}
};
class FF
{
public:
	int S,T,n,m,cnt;
	int head[N],vis[N];
	edge e[M];
	void addedge(int u,int v,int cap)
	{
		e[++cnt]=edge(v,head[u],cap);
		head[u]=cnt;
	}
	
	void buildGraph(int _n,int _m,int _S,int _T)//建图
	{
		n = _n;m = _m;
		S = _S;T = _T;
		mem(head);mem(vis);
		cnt = 1;
	
		for(int i=1;u,v,cap;i<=m;i++)
		{
			read(u);read(v);read(cap);
			addedge(u,v,cap);
			addedge(v,u,0);
		}
	}
	int dfs(int u.int flow)
	{
		if(u==T) return flow;
		vis[u] = 1;
		for(int i=head[u];i;i = e[i].nex)
		{
			int v = e[i].to;
			if(e[i].cap<=0||vis[v]) continue;
			int delta = dfs(v,min(flow,e[i].cap));
			if (delta<=0) continue;
			e[i].cap-=delta;
			e[i^1].cap+=delta;
			return maxFlow;
		}
	}
	int get_maxFlow()
	{
		int maxFlow = 0;
		while(delta = dfs(S,INF)) 
		{
			mem(vis);
			maxFlow += delta;
		}
		return maxFlow;
	}
};
FF ways;
signed main()
{
	int n,m, S, T;
	while(~scanf("%d %d %d %d",&n,&m,&S,&T))
	{
		ways.buildgraph(n,m,S,T);
		printf("%d\b",ways.get_maxFlow());
	}
	
}//源自b站某up