題目鏈接：點擊查看
題目大意：給出一張 n 個點 m 條邊的有向圖，每個點可以生產(chǎn)至多 b[ i ] 個物品，每個物品花費為 a[ i ] ，可以出售至多 d[ i ] 個物品，每個物品售價為 c[ i ]，所有的物品都是一樣的，而運輸一個物品一公里需要花費一個單位的金錢，現(xiàn)在問如何調(diào)配，才能使得收益最大
題目分析：最大費用可行流，注意這里是可行流而不是最大流，如果是最大流的話下面這一組樣例跑不過去：
5 4
1 2 1 2
2 1 2 1
2 3 4 5
5 4 3 2
100 1 1 1
1 2 1
2 3 1
3 4 1
1 5 1
ans = 8
所以在這個題目中，最大費用的優(yōu)先級要高于最大流，所以spfa中的判斷條件是必須要有正貢獻(xiàn)的增廣路才可以增大流量，如果出現(xiàn)負(fù)貢獻(xiàn)的增廣路直接退出即可
建圖也非常簡單：
源點 -> 每個點，流量為 b[ i ] ，花費為 -a[ i ] 原圖中的 u 和 v ，u -> v ，流量為 inf ，花費為 -w 每個點 -> 匯點，流量為 d[ i ] ，花費為 c[ i ]
然后跑模板就好了
最后說點題外話，因為最小費用可行流是一個單獨的算法，一直在想為什么這個題不能用最小費用可行流的算法去解決，與前一個題目稍微比較了一下?洛谷 - P4043?，得出的結(jié)論是，最小費用可行流，實際上是用無源匯的可行流將其賦值上費用了而已，轉(zhuǎn)換后的圖還是以最大流為優(yōu)先級，而當(dāng)所有的邊流量都達(dá)到下限后，因為所有的邊費用都是正的，所以選擇的邊肯定越少費用越小，所以在之前的那個題目用最小費用可行流是可行的，而在這個題目中，因為邊權(quán)有負(fù)邊，所以不能用最大流來限制流量和費用，所以只能修改spfa達(dá)到以最大費用為優(yōu)先級的目的?
代碼：
//#pragma GCC optimize(2)
//#pragma GCC optimize("Ofast","inline","-ffast-math")
//#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<iostream>
#include<cstdio>
#include<string>
#include<ctime>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
#include<climits>
#include<queue>
#include<map>
#include<set>
#include<sstream>
#include<cassert>
#include<bitset>
#include<unordered_map>
using namespace std;typedef long long LL;typedef unsigned long long ull;const int inf=0x3f3f3f3f;const int N=510;//點const int M=1e4+100;//邊struct Edge
{int to,w,cost,next;
}edge[M];int head[N],cnt;void addedge(int u,int v,int w,int cost)
{edge[cnt].to=v;edge[cnt].w=w;edge[cnt].cost=cost;edge[cnt].next=head[u];head[u]=cnt++;edge[cnt].to=u;edge[cnt].w=0;edge[cnt].cost=-cost;edge[cnt].next=head[v];head[v]=cnt++;
}int d[N],incf[N],pre[N];bool vis[N];bool spfa(int s,int t)
{memset(d,0xcf,sizeof(d));memset(vis,false,sizeof(vis));memset(pre,-1,sizeof(pre));queue<int>q;q.push(s);vis[s]=true;incf[s]=inf;d[s]=0;while(!q.empty()){int u=q.front();q.pop();vis[u]=false;for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].to;int w=edge[i].w;int cost=edge[i].cost;if(!w)continue;if(d[v]<d[u]+cost){d[v]=d[u]+cost;pre[v]=i;incf[v]=min(incf[u],w);if(!vis[v]){vis[v]=true;q.push(v);}}}}
// return pre[t]!=-1;return d[t]>0;
}int update(int s,int t)
{int x=t;while(x!=s){int i=pre[x];edge[i].w-=incf[t];edge[i^1].w+=incf[t];x=edge[i^1].to;}return d[t]*incf[t];
}void init()
{memset(head,-1,sizeof(head));cnt=0;
}int solve(int st,int ed)
{int ans=0;while(spfa(st,ed))ans+=update(st,ed);return ans;
}int main()
{
#ifndef ONLINE_JUDGE
// freopen("data.in.txt","r",stdin);
// freopen("data.out.txt","w",stdout);
#endif
// ios::sync_with_stdio(false);int n,m,st=N-1,ed=st-1;while(scanf("%d%d",&n,&m)!=EOF){init();for(int i=1;i<=n;i++){int a,b,c,d;scanf("%d%d%d%d",&a,&b,&c,&d);addedge(st,i,b,-a);addedge(i,ed,d,c);}while(m--){int u,v,w;scanf("%d%d%d",&u,&v,&w);if(u!=v){addedge(u,v,inf,-w);addedge(v,u,inf,-w);}}printf("%d\n",solve(st,ed));}return 0;
}
?
總結(jié)
以上是生活随笔 為你收集整理的HDU - 6118 度度熊的交易计划(最大费用可行流) 的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔 網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔 推薦給好友。
