739. [網(wǎng)絡流24題] 運輸問題

★★?? 輸入文件：tran.in?? 輸出文件：tran.out???簡單對比
時間限制：1 s?? 內存限制：128 MB

?問題描述：

?編程任務：
對于給定的m 個倉庫和n 個零售商店間運送貨物的費用，計算最優(yōu)運輸方案和最差運
輸方案。
?數(shù)據(jù)輸入：

?結果輸出:
程序運行結束時，將計算出的最少運輸費用和最多運輸費用輸出到文件tran.out中。
輸入文件示例 輸出文件示例
tran.in
2 3
220 280
170 120 210
77 39 105

150 186 122

tran.out

?

69140

?

對于所有數(shù)據(jù)：1<=N,M<=100

?

題解：
比較容易想到費用流。
建圖：
1>建立虛擬源S和匯T。
2>從S到m個倉庫建流量為a[i],費用為0的邊。
3>從n個商店向T建流量為b[i]，費用為0的邊。
4>從m個倉庫向n個商店建流量為無窮大，費用為c[i][j]的邊。
然后跑一邊最小費用，再跑一邊最大費用就好了。

AC代碼：

?

#include<cstdio> #include<cstring> #include<iostream> #define m(s,t) memset(s,t,sizeof s) #define R register #define inf 2139062143 using namespace std; int read(){R int x=0;bool f=1;R char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=0;ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}return f?x:-x; } const int N=310; struct node{int v,next,cap,cost; }e[N*N];int tot=1; int n,m,S,T,ans,head[N],pree[N],prev[N],flow[N],dis[N],q[N*10]; int a[N],b[N],c[N][N]; bool vis[N]; void add(int x,int y,int a,int b){e[++tot].v=y;e[tot].cap=a;e[tot].cost=b;e[tot].next=head[x];head[x]=tot; } void ins(int x,int y,int a,int b){add(x,y,a,b);add(y,x,0,-b); } void build(){S=0,T=m+n+1;for(int i=1;i<=m;i++) ins(S,i,a[i],0);for(int i=1;i<=n;i++) ins(i+m,T,b[i],0);for(int i=1;i<=m;i++){for(int j=1;j<=n;j++){ins(i,j+m,inf,c[i][j]);}} } void Cl(){tot=1;ans=0;m(head,0);m(pree,0);m(prev,0);m(flow,0); } bool spfa(int k){m(vis,0);k>0?m(dis,127):m(dis,128);int h=0,t=1;q[t]=S;dis[S]=0;vis[S]=1;flow[S]=inf;while(h!=t){int x=q[++h];vis[x]=0;for(int i=head[x];i;i=e[i].next){int v=e[i].v,cap=e[i].cap,cost=e[i].cost;if(cap>0&&((k>0&&dis[v]>dis[x]+cost)||(k<0&&dis[v]<dis[x]+cost))){dis[v]=dis[x]+cost;prev[v]=x;pree[v]=i;flow[v]=min(flow[x],cap);if(!vis[v]){vis[v]=1;q[++t]=v;}}}}if(k>0) return dis[T]<inf;else return dis[T]>0; } void work(){for(int i=T;i!=S;i=prev[i]){e[pree[i]].cap-=flow[T];e[pree[i]^1].cap+=flow[T];}ans+=flow[T]*dis[T]; } int main(){freopen("tran.in","r",stdin);freopen("tran.out","w",stdout);m=read();n=read();for(int i=1;i<=m;i++) a[i]=read();for(int i=1;i<=n;i++) b[i]=read();for(int i=1;i<=m;i++){for(int j=1;j<=n;j++){c[i][j]=read();}}build();while(spfa(1)) work();printf("%d\n",ans);Cl();build();while(spfa(-1)) work();printf("%d\n",ans);return 0; }

?

轉載于:https://www.cnblogs.com/shenben/p/6220688.html

總結

以上是生活随笔為你收集整理的739. [网络流24题] 运输问题的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內容還不錯，歡迎將生活随笔推薦給好友。