让我再讲一个故事吧。
又有一些小精灵要准备从银月城(S)迁徙到Nibel山(T)。
这两个地方之间的道路构成了一个网络。
每个道路都有它自己的容量,这决定了每天有多少小精灵可以同时从这儿通过。
和上一篇不同的是,由于上次迁徙的规模很大,
吸引了其它一些种族的注意,
这次每条道路都会有一些人/兽人/哥布林/...向精灵们征收过路费,
现在精灵们想知道,在花费最小的情况下,它们迁徙的速度最大是多少只每天。
费用流=最小费用最大流
在要求流最大的情况下要求费用最小,好像原来的isap已经派不上用场了呢!
让我们回到最朴实的EK算法上。
EK算法每一次只寻找一条增广路,
这带给它解决这一个方面的问题的得天独厚的优势。
这是原来的EK伪算法:
int BFS()
{
/*找到一条增广路*/
}
int ek()
{
/*对找到的增广路进行一系列处理*/
}
我们用BFS找增广路。
想象一下,
既然要求费用最小,
我们就把费用作为路径长度,
之后每一次跑一遍最短路,
那么就可以保证花费最小了!
所以,我们只要把原来的
BFS()
改成spfa()
或者dijkstra()
就好啦
ps.一般dijkstra只能跑不带负权边的图,
但是有一种特殊的技巧可以把边权魔改成正的。
以下是拿辣鸡spfa跑的费用流
#include<bits/stdc++.h>
using namespace std;
typedef pair<int,int> pii;
#define mp(a,b) make_pair(a,b)
#define ll first
#define rr second
inline int gotcha()
{
register int a=0,b=1,c=getchar();
while(!isdigit(c))b^=a=='-',c=getchar();
while(isdigit(c))a=a*10+c-48,c=getchar();
return b?a:-a;
}
const int _ = 5002 , __ = 50002<<1 , inf = 0x3f3f3f3f;
int to[__],ne[__],v[__],co[__],he[__]={0},ecnt=1;
int n,m,dis[_],pe[_],pv[_],S,T;
bool ed[_];
void adde(int a,int b,int c,int d){to[++ecnt]=b,v[ecnt]=c,co[ecnt]=d,ne[ecnt]=he[a],he[a]=ecnt;}
queue<int> q;
int spfa()
{
memset(dis,63,sizeof(dis)),memset(ed,0,sizeof(ed));
while(!q.empty())q.pop();
register int i,a;
q.push(S),ed[S]=1,dis[S]=0;
while(!q.empty())
{
a=q.front(),q.pop();ed[a]=0;
for(i=he[a];i;i=ne[i])
if(v[i]>0 && dis[to[i]]>dis[a]+co[i])
{
dis[to[i]]=dis[a]+co[i];
pe[to[i]]=i,pv[to[i]]=a;
if(!ed[to[i]])ed[to[i]]=1,q.push(to[i]);
}
}
return dis[T]<inf;
}
pii mfmc()
{
register int i,sco=0,sfl=0,flw;
while(spfa())
{
flw=inf;
for(i=T;i!=S;i=pv[i])flw=min(flw,v[pe[i]]);
for(i=T;i!=S;i=pv[i])v[pe[i]]-=flw,v[pe[i]^1]+=flw;
sco+=flw*dis[T],sfl+=flw;
}
return mp(sfl,sco);
}
int main()
{
register int i,j,k,a,b;
register pii tmp;
n=gotcha(),m=gotcha(),S=gotcha(),T=gotcha();
for(i=1;i<=m;i++)
{
j=gotcha(),k=gotcha(),a=gotcha(),b=gotcha();
adde(j,k,a,b),adde(k,j,0,-b);
}
tmp=mfmc();
printf("%d %d",tmp.ll,tmp.rr);
return 0;
}
这就不写伪代码了吧!?
以后补
原文地址:https://www.cnblogs.com/finder-iot/p/8409010.html
时间: 2024-10-03 10:24:13