【ZOJ】3362 Beer Problem 最小费用流

2024-09-05 15:18
文章标签 最小 problem 费用 zoj beer 3362

本文主要是介绍【ZOJ】3362 Beer Problem 最小费用流,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

传送门:【ZOJ】3362 Beer Problem


题目分析:这道题本来应该很快就AC的,但是!因为我以前犯的一个致命错误导致我这题一天了到现在才调出来!唉。。失策。。貌似给的模板也有这个错误。。。马上就去改。。但是这个错误竟然还能过掉那么多的题。。害我还要一题一题的改回去。。


本题就是赤裸裸的最小费用流。

新建汇点t(源点即1),将所有的n-1个城市和汇点建边,容量为无穷大,费用为单位收益的相反数,然后所有的关系建立无向边,容量为题目中给出的,费用为给出的单位花费。跑一遍最小费用流以后,费用的相反数就是要求的最大收益。有没有发现我并没有说最小费用最大流?是的,最小费用的时候并不一定是最大流!由于每次找的增广路都是当前费用最小的(最短路),所以如果某一次得到的费用已经不小于0了,那么说明已经过了最优解了,可以直接break。只要在费用流算法中稍作修改即可。

今天各种卡题T U T不爱了。。。


代码如下:


#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;#define REP( i , a , b ) for ( int i = a ; i < b ; ++ i )
#define FOR( i , a , b ) for ( int i = a ; i <= b ; ++ i )
#define REV( i , a , b ) for ( int i = a ; i >= b ; -- i )
#define CLR( a , x ) memset ( a , x , sizeof a )
#define CPY( a , x ) memcpy ( a , x , sizeof a )typedef int type_c ;
typedef int type_f ;
typedef int type_w ;const int MAXN = 128 ;
const int MAXQ = 128 ;
const int MAXE = 10000 ;
const int INF = 0x3f3f3f3f ;struct Edge {int v , n ;type_c c ;type_w w ;Edge () {}Edge ( int v , type_c c , type_w w , int n ) : v ( v ) , c ( c ) , w ( w ) , n ( n ) {}
} ;struct Net {Edge E[MAXE] ;int H[MAXN] , cntE ;int d[MAXN] , cur[MAXN] , f[MAXN] ;int Q[MAXQ] , head , tail , inq[MAXN] ;int s , t ;type_f flow ;type_w cost ;int n , m ;void init () {cntE = 0 ;CLR ( H , -1 ) ;}void addedge ( int u , int v , type_c c , type_w w ) {E[cntE] = Edge ( v , c ,  w , H[u] ) ;H[u] = cntE ++ ;E[cntE] = Edge ( u , 0 , -w , H[v] ) ;H[v] = cntE ++ ;}int spfa () {head = tail = 0 ;CLR ( d , INF ) ;CLR ( inq , 0 ) ;Q[tail ++] = s ;f[s] = INF ;cur[s] = -1 ;d[s] = 0 ;while ( head != tail ) {int u = Q[head ++] ;if ( head == MAXQ )head = 0 ;inq[u] = 0 ;for ( int i = H[u] ; ~i ; i = E[i].n ) {int v = E[i].v ;if ( E[i].c && d[v] > d[u] + E[i].w ) {f[v] = min ( f[u] , E[i].c ) ;d[v] = d[u] + E[i].w ;cur[v] = i ;if ( !inq[v] ) {if ( d[v] < d[Q[head]] ) {if ( head == 0 )head = MAXQ ;Q[-- head] = v ;}else {Q[tail ++] = v ;if ( tail == MAXQ )tail = 0 ;}inq[v] = 1 ;}}}}if ( d[t] >= 0 )return 0 ;flow += f[t] ;cost += f[t] * d[t] ;for ( int i = cur[t] ; ~i ; i = cur[E[i ^ 1].v] ) {E[i].c -= f[t] ;E[i ^ 1].c += f[t] ;}return 1 ;}type_w mcmf () {flow = cost = 0 ;while ( spfa () ) ;return cost ;}void solve () {int u , v , c , w ;init () ;s = 1 , t = n + 1 ;FOR ( i , 2 , n ) {scanf ( "%d" , &w ) ;addedge ( i , t , INF , -w ) ;}while ( m -- ) {scanf ( "%d%d%d%d" , &u , &v , &c , &w ) ;addedge ( u , v , c , w ) ;addedge ( v , u , c , w ) ;}printf ( "%d\n" , -mcmf () ) ;}
} e ;int main () {while ( ~scanf ( "%d%d" , &e.n , &e.m ) )e.solve () ;return 0 ;
}


这篇关于【ZOJ】3362 Beer Problem 最小费用流的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1139322

相关文章

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n

poj 1287 Networking(prim or kruscal最小生成树)

题意给你点与点间距离,求最小生成树。 注意点是,两点之间可能有不同的路,输入的时候选择最小的,和之前有道最短路WA的题目类似。 prim代码: #include<stdio.h>const int MaxN = 51;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int P;int prim(){bool vis[MaxN];

poj 2349 Arctic Network uva 10369(prim or kruscal最小生成树)

题目很麻烦,因为不熟悉最小生成树的算法调试了好久。 感觉网上的题目解释都没说得很清楚,不适合新手。自己写一个。 题意:给你点的坐标,然后两点间可以有两种方式来通信:第一种是卫星通信,第二种是无线电通信。 卫星通信:任何两个有卫星频道的点间都可以直接建立连接,与点间的距离无关; 无线电通信:两个点之间的距离不能超过D,无线电收发器的功率越大,D越大,越昂贵。 计算无线电收发器D

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

poj 1734 (floyd求最小环并打印路径)

题意: 求图中的一个最小环,并打印路径。 解析: ans 保存最小环长度。 一直wa,最后终于找到原因,inf开太大爆掉了。。。 虽然0x3f3f3f3f用memset好用,但是还是有局限性。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#incl

hdu 1102 uva 10397(最小生成树prim)

hdu 1102: 题意: 给一个邻接矩阵,给一些村庄间已经修的路,问最小生成树。 解析: 把已经修的路的权值改为0,套个prim()。 注意prim 最外层循坏为n-1。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstri

poj 2175 最小费用最大流TLE

题意: 一条街上有n个大楼,坐标为xi,yi,bi个人在里面工作。 然后防空洞的坐标为pj,qj,可以容纳cj个人。 从大楼i中的人到防空洞j去避难所需的时间为 abs(xi - pi) + (yi - qi) + 1。 现在设计了一个避难计划,指定从大楼i到防空洞j避难的人数 eij。 判断如果按照原计划进行,所有人避难所用的时间总和是不是最小的。 若是,输出“OPETIMAL",若

poj 2135 有流量限制的最小费用最大流

题意: 农场里有n块地,其中约翰的家在1号地,二n号地有个很大的仓库。 农场有M条道路(双向),道路i连接着ai号地和bi号地,长度为ci。 约翰希望按照从家里出发,经过若干块地后到达仓库,然后再返回家中的顺序带朋友参观。 如果要求往返不能经过同一条路两次,求参观路线总长度的最小值。 解析: 如果只考虑去或者回的情况,问题只不过是无向图中两点之间的最短路问题。 但是现在要去要回

poj 3422 有流量限制的最小费用流 反用求最大 + 拆点

题意: 给一个n*n(50 * 50) 的数字迷宫,从左上点开始走,走到右下点。 每次只能往右移一格,或者往下移一格。 每个格子,第一次到达时可以获得格子对应的数字作为奖励,再次到达则没有奖励。 问走k次这个迷宫,最大能获得多少奖励。 解析: 拆点,拿样例来说明: 3 2 1 2 3 0 2 1 1 4 2 3*3的数字迷宫,走两次最大能获得多少奖励。 将每个点拆成两个