网络流与费用流（下）费用流

最新推荐文章于 2023-01-07 16:14:47 发布

窝嘞割草

最新推荐文章于 2023-01-07 16:14:47 发布

阅读量471

点赞数

本文链接：https://blog.csdn.net/the301stdoub/article/details/51073297

版权

费用流

费用流算法

好像有好多高级算法的样子，但我现在只会一个最简单的spfa连续最短路增广，这个东西就是贪心啦，最短路使每次增广的流的费用最小，没啥好说的。

针对这个做法，要注意的地方是，无向图每条边要建4条边。。。还有如果费用有负值要加判负环，别的也没什么了，算法是比较脑残的，关键是建图。

const int maxn=1010;
const int maxm=50010;
const int INF=0x3f3f3f3f;

struct edge
{
	int v,w,c,next;
}e[maxm];
int h[maxn],ecnt;

inline void addedge(int u,int v,int w,int c)
{
	e[ecnt].v=v;
	e[ecnt].w=w;
	e[ecnt].c=c;
	e[ecnt].next=h[u];
	h[u]=ecnt++;
}

inline void Addedge(int u,int v,int w,int c)
{
	addedge(u,v,w,c);
	addedge(v,u,0,-c);
	
	addedge(v,u,w,c);//undirected
	addedge(u,v,0,-c);
}

int dis[maxn],pre[maxn],peg[maxn];
bool vis[maxn];
int flow,cost;

bool spfa(int s,int t)
{
	clr(vis,0);clr(pre,-1);
	clr(dis,0x3f);
	queue<int> q;
	q.push(s);
	vis[s]=1;
	dis[s]=0;
	while(!q.empty())
	{
		int u=q.front();q.pop();vis[u]=0;
		for(int i=h[u];~i;i=e[i].next)
		{
			int v=e[i].v;
			if(e[i].w && dis[u]+e[i].c<dis[v])
			{
				dis[v]=dis[u]+e[i].c;
				pre[v]=u;
				peg[v]=i;
				if(!vis[v])q.push(v);
				vis[v]=1;
			}
		}
	}
	if(dis[t]==INF)return 0;
	return 1;
}

void mincostflow(int s,int t)
{
	flow=cost=0;
	while(spfa(s,t))
	{
		int mf=INF;
		for(int u=t;~pre[u];u=pre[u])mf=min(mf,e[peg[u]].w);
		for(int u=t;~pre[u];u=pre[u])
		{
			e[peg[u]].w-=mf;
			e[peg[u]^1].w+=mf;
		}
		flow+=mf;
		cost+=mf*dis[t];
	}
}

其他方法留坑待填。。。