最小花费最大流模板

自己写的最小花费最大流模板

#include<iostream>#include<cstdio>#include<cmath>#include<cstring>#include<queue>using namespace std;const int maxn=10010;//点数const int maxx=100100;//边数，注意为实际上边数的两倍const int INF=1<<30;int to[maxx],nex[maxx],cap[maxx],flow[maxx],cost[maxx];int head[maxn];int pre[maxn],dis[maxn];int eid;void init(){    eid=0;    memset(head,-1,sizeof(head));}//初始化void addEdge(int u,int v,int ca,int co)//注意是u到v的有向边，ca是容量，co是花费{    to[eid]=v;    cap[eid]=ca;    cost[eid]=co;    flow[eid]=0;    nex[eid]=head[u];    head[u]=eid++;    to[eid]=u;    cap[eid]=0;    cost[eid]=-co;    flow[eid]=0;    nex[eid]=head[v];    head[v]=eid++;}bool inQ[maxn];bool spfa(int s,int e,int n)//s点到e点的是否存在最短路，n为点的数量，点可以从0开始，即0到n{    queue<int> que;    memset(inQ,false,sizeof(inQ));    for(int i=0;i<=n;i++)        dis[i]=INF;    dis[s]=0,inQ[s]=true;    pre[s]=-1;    que.push(s);    while(!que.empty())    {        int u=que.front();        que.pop();        inQ[u]=false;        for(int i=head[u];~i;i=nex[i])        {            int v=to[i];            if(cap[i]-flow[i]>0&&dis[v]>dis[u]+cost[i])            {                dis[v]=dis[u]+cost[i];                pre[v]=i;//这是为了修改边的容量                if(!inQ[v])                    que.push(v),inQ[v]=true;            }        }    }    return dis[e]!=INF;}int minCostMaxFlow(int s,int e,int &minCost,int n)//最小流，s点到e的最小花费最大流，n含义同上{    int ans=0;    while(spfa(s,e,n))    {        int mint=INF;        for(int i=pre[e];~i;i=pre[to[i^1]])        {            if(mint>cap[i]-flow[i])                mint=cap[i]-flow[i];        }        ans+=mint;        for(int i=pre[e];~i;i=pre[to[i^1]])        {            flow[i]+=mint;flow[i^1]-=mint;            minCost+=mint*cost[i];        }    }    return ans;}