费用流&网络流模版
费用流模版：
1 #include<cstdio> 2 #include<cstring> 3 #include<queue> 4 using namespace std; 5 6 const int Maxm=100000;//最大边数 7 const int Maxn=1000;//最大点数 8 struct Edge{ 9 Edge(){}; 10 Edge(int a,int b,int c,int d,int e){ 11 u=a; 12 v=b; 13 f=c; 14 w=d; 15 nxt=e; 16 } 17 int u,v,f,w,nxt;//U当前点 V来自点 F最大流量 W费用 NXT下一个点 18 }; 19 int cnt=1;//边计数 20 int inf=2147483647;//无限大 21 int g[Maxn+10];//点的边集的开始序号 22 Edge e[Maxm+10];//边集 23 int dist[Maxn+10];//费用 24 int src,sink;//源点与汇点 25 queue<int> que;//宽搜队列 26 bool inque[Maxn+10];//宽搜判断标志 27 int from[Maxn+10];//来源->用于计算费用 28 int ans=0;//存储最小费用 29 30 inline int remin(int a,int b){ 31 return a<b?a:b; 32 } 33 34 inline void insert(int u,int v,int f,int w){ 35 cnt++; 36 e[cnt]=Edge(u,v,f,w,g[u]); 37 g[u]=cnt;//增加一个边 38 } 39 40 inline void addEdge(int u,int v,int f,int w){ 41 insert(u,v,f,w);//插入正边 42 insert(v,u,0,-w);//插入反边 43 } 44 45 inline bool spfa(){ 46 while (!que.empty()) que.pop();//清空队列 47 for (int i=0;i<=sink;i++) dist[i]=inf;//清最大值 48 que.push(src); 49 inque[src]=true; 50 dist[src]=0;//加入源点 51 //标准SPFA计算最短路 流量作为通行标准，费用作为路径长度 52 while(!que.empty()){ 53 int now=que.front(); 54 que.pop(); 55 for (int i=g[now];i;i=e[i].nxt){ 56 if (e[i].f!=0 && dist[e[i].v]>dist[now]+e[i].w){ 57 dist[e[i].v]=dist[now]+e[i].w; 58 from[e[i].v]=i; 59 if (inque[e[i].v]==false){ 60 inque[e[i].v]=true; 61 que.push(e[i].v); 62 } 63 } 64 } 65 inque[now]=false; 66 } 67 if (dist[sink]==inf) return false;//无法在增广 68 return true; 69 } 70 71 inline void calcAns(){ 72 int minflow=inf; 73 for (int i=from[sink];i;i=from[e[i].u]) minflow=remin(minflow,e[i].f);//寻找整条路经的流量 74 for (int i=from[sink];i;i=from[e[i].u]) { 75 e[i].f-=minflow;//正边减流量 76 e[i^1].f+=minflow;//反边加流量 77 ans+=e[i].w*minflow;//计算费用 78 } 79 } 80 81 inline void minCostFlow(){ 82 while(spfa())calcAns(); 83 } 84 85 int main(){ 86 minCostFlow(); 87 return 0; 88 }
网络流模版：
1 #include<cstdio> 2 #include<cstring> 3 #include<queue> 4 using namespace std; 5 6 const int Maxm=100000;//最大边数 7 const int Maxn=1000;//最大点数 8 struct Edge{ 9 Edge(){}; 10 Edge(int a,int b,int c,int d,int e){ 11 u=a; 12 v=b; 13 f=c; 14 w=d; 15 nxt=e; 16 } 17 int u,v,f,w,nxt;//U当前点 V来自点 F最大流量 W费用 NXT下一个点 18 }; 19 int cnt=1;//边计数 20 int inf=2147483647;//无限大 21 int g[Maxn+10];//点的边集的开始序号 22 Edge e[Maxm+10];//边集 23 int dist[Maxn+10];//费用 24 int src,sink;//源点与汇点 25 queue<int> que;//宽搜队列 26 bool inque[Maxn+10];//宽搜判断标志 27 int from[Maxn+10];//来源->用于计算费用 28 int ans=0;//存储最小费用 29 30 inline int remin(int a,int b){ 31 return a<b?a:b; 32 } 33 34 inline void insert(int u,int v,int f,int w){ 35 cnt++; 36 e[cnt]=Edge(u,v,f,w,g[u]); 37 g[u]=cnt;//增加一个边 38 } 39 40 inline void addEdge(int u,int v,int f,int w){ 41 insert(u,v,f,w);//插入正边 42 insert(v,u,0,-w);//插入反边 43 } 44 45 inline bool spfa(){ 46 while (!que.empty()) que.pop();//清空队列 47 for (int i=0;i<=sink;i++) dist[i]=inf;//清最大值 48 que.push(src); 49 inque[src]=true; 50 dist[src]=0;//加入源点 51 //标准SPFA计算最短路 流量作为通行标准，费用作为路径长度 52 while(!que.empty()){ 53 int now=que.front(); 54 que.pop(); 55 for (int i=g[now];i;i=e[i].nxt){ 56 if (e[i].f!=0 && dist[e[i].v]>dist[now]+e[i].w){ 57 dist[e[i].v]=dist[now]+e[i].w; 58 from[e[i].v]=i; 59 if (inque[e[i].v]==false){ 60 inque[e[i].v]=true; 61 que.push(e[i].v); 62 } 63 } 64 } 65 inque[now]=false; 66 } 67 if (dist[sink]==inf) return false;//无法在增广 68 return true; 69 } 70 71 inline void calcAns(){ 72 int minflow=inf; 73 for (int i=from[sink];i;i=from[e[i].u]) minflow=remin(minflow,e[i].f);//寻找整条路经的流量 74 for (int i=from[sink];i;i=from[e[i].u]) { 75 e[i].f-=minflow;//正边减流量 76 e[i^1].f+=minflow;//反边加流量 77 ans+=e[i].w*minflow;//计算费用 78 } 79 } 80 81 inline void minCostFlow(){ 82 while(spfa())calcAns(); 83 } 84 85 int main(){ 86 minCostFlow(); 87 return 0; 88 }
