标签：费用 num EK vis int 最小 read edge dis

最小费用最大流

这里只给出代码实现：

class flow_map
{
	private:
		int u;
		int s,t;
		int n,m;
		inline bool spfa()
		{
			memset(vis,0,sizeof(vis)),memset(dis,0x3f,sizeof(dis));
			queue<int> q;
			while(!q.empty()) q.pop();
			dis[s]=0,vis[s]=1,q.push(s),incf[s]=0x3f3f3f3f;
			register int u;
			while(!q.empty())
			{
				u=q.front(),q.pop(),vis[u]=0;
				for(register int j(head[u]);j;j=nex[j])
				{
					if(!flow[j]) continue;
					if(dis[to[j]]>dis[u]+cost[j])
					{
						dis[to[j]]=dis[u]+cost[j],
						incf[to[j]]=min(incf[u],flow[j]),
						pre[to[j]]=j;
						if(!vis[to[j]]) vis[to[j]]=1,q.push(to[j]);
					}
				}
			}
			return (dis[t]!=1061109567);
		}
		inline int against_edge(int id)
		{
			if(id&1) return (id+1);
			else return (id-1);
		}
	public:
		int edge_num;
		int to[maxm],head[maxn],nex[maxm],flow[maxm],cost[maxm];
		int dis[maxn];
		bool vis[maxn];
		int incf[maxn];
		int pre[maxn];
		int max_flow;
		inline void set()
		{
			read_(n),
			read_(m),
			read_(s),
			read_(t);
		}
		inline void add_edge(void)
		{
			++edge_num,
			read_(u),read_(to[edge_num]),
			read_(flow[edge_num]),
			read_(cost[edge_num]),
			nex[edge_num]=head[u],
			head[u]=edge_num;
		}
		inline void add_edge(int u,int v,int f,int c)
		{
			to[++edge_num]=v,
			nex[edge_num]=head[u],
			head[u]=edge_num,
			cost[edge_num]=c,
			flow[edge_num]=f;
		}
		inline void build()
		{
			for(register int i(0),u,v,w,c;i<m;++i)
				read_(u),
				read_(v),
				read_(w),
				read_(c),
				add_edge(u,v,w,c),
				add_edge(v,u,0,-c);
		}
		inline int max_flow__min_cost()
		{
			register int mxf(0),mic(0),x;
			while(spfa())
			{
				mxf+=incf[x=t],
				mic+=incf[t]*dis[t];
				while(x!=s)
					flow[pre[x]]-=incf[t],
					flow[against_edge(pre[x])]+=incf[t],
					x=to[against_edge(pre[x])];
			}
			return max_flow=mxf,mic;
		}
};

标签：费用,num,EK,vis,int,最小,read,edge,dis
来源： https://www.cnblogs.com/efX-bk/p/crazy-zip-ek-max-flow-min-cost.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。