739. [网络流24题] 运输问题
生活随笔
收集整理的這篇文章主要介紹了
739. [网络流24题] 运输问题
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
739. [網絡流24題] 運輸問題
★★?? 輸入文件:tran.in?? 輸出文件:tran.out???簡單對比
時間限制:1 s?? 內存限制:128 MB
?問題描述:
?編程任務:
對于給定的m 個倉庫和n 個零售商店間運送貨物的費用,計算最優運輸方案和最差運
輸方案。
?數據輸入:
?結果輸出:
程序運行結束時,將計算出的最少運輸費用和最多運輸費用輸出到文件tran.out中。
輸入文件示例 輸出文件示例
tran.in
2 3
220 280
170 120 210
77 39 105
150 186 122
tran.out
?
69140
?
對于所有數據: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题] 运输问题的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 《HTML5+CSS3网页设计入门必读》
- 下一篇: java PreparedStatem