【网络流24题】No.19 负载平衡问题 (费用流)
【题意】
G 公司有 n 个沿铁路运输线环形排列的仓库, 每个仓库存储的货物数量不等。 如何用最
少搬运量可以使 n 个仓库的库存数量相同。搬运货物时,只能在相邻的仓库之间搬运。
输入文件示例
input.txt
5
17 9 14 16 4输出文件示例
output.txt
11
【分析】
其实我觉得这题可以贪心啊。。n^2贪心??、没细想。。
打的是费用流。。
大概这样建图:
懒得写了。。凌乱之美。。
求满流费用。。
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cstring> 4 #include<iostream> 5 #include<algorithm> 6 #include<queue> 7 #include<cmath> 8 using namespace std; 9 #define Maxn 1010 10 #define INF 0xfffffff 11 12 struct node 13 { 14 int x,y,f,o,c,next; 15 }t[Maxn*Maxn];int len; 16 int first[Maxn]; 17 18 int mymin(int x,int y) {return x<y?x:y;} 19 int mymax(int x,int y) {return x>y?x:y;} 20 21 void ins(int x,int y,int f,int c) 22 { 23 t[++len].x=x;t[len].y=y;t[len].f=f;t[len].c=c; 24 t[len].next=first[x];first[x]=len;t[len].o=len+1; 25 t[++len].x=y;t[len].y=x;t[len].f=0;t[len].c=-c; 26 t[len].next=first[y];first[y]=len;t[len].o=len-1; 27 } 28 29 int st,ed; 30 queue<int > q; 31 int dis[Maxn],pre[Maxn],flow[Maxn]; 32 bool inq[Maxn]; 33 bool bfs() 34 { 35 while(!q.empty()) q.pop(); 36 memset(dis,63,sizeof(dis)); 37 memset(inq,0,sizeof(inq)); 38 q.push(st);dis[st]=0;flow[st]=INF;inq[st]=1; 39 while(!q.empty()) 40 { 41 int x=q.front(); 42 for(int i=first[x];i;i=t[i].next) if(t[i].f>0) 43 { 44 int y=t[i].y; 45 if(dis[y]>dis[x]+t[i].c) 46 { 47 dis[y]=dis[x]+t[i].c; 48 pre[y]=i; 49 flow[y]=mymin(flow[x],t[i].f); 50 if(!inq[y]) 51 { 52 inq[y]=1; 53 q.push(y); 54 } 55 } 56 } 57 inq[x]=0;q.pop(); 58 } 59 if(dis[ed]>=INF-100000) return 0; 60 return 1; 61 } 62 63 void output() 64 { 65 for(int i=1;i<=len;i+=2) 66 printf("%d->%d %d %d\n",t[i].x,t[i].y,t[i].f,t[i].c); 67 printf("\n"); 68 } 69 70 int max_flow() 71 { 72 int ans=0,sum=0; 73 while(bfs()) 74 { 75 sum+=dis[ed]*flow[ed]; 76 ans+=flow[ed]; 77 int now=ed; 78 while(now!=st) 79 { 80 t[pre[now]].f-=flow[ed]; 81 t[t[pre[now]].o].f+=flow[ed]; 82 now=t[pre[now]].x; 83 } 84 } 85 return sum; 86 } 87 88 int n,a[Maxn]; 89 90 void init() 91 { 92 int h=0; 93 scanf("%d",&n); 94 st=n+1;ed=st+1; 95 len=0; 96 memset(first,0,sizeof(first)); 97 for(int i=1;i<=n;i++) 98 { 99 scanf("%d",&a[i]); 100 ins(st,i,a[i],0); 101 h+=a[i]; 102 } 103 h/=n; 104 for(int i=1;i<n;i++) ins(i,i+1,INF,1);ins(n,1,INF,1); 105 for(int i=2;i<=n;i++) ins(i,i-1,INF,1);ins(1,n,INF,1); 106 for(int i=1;i<=n;i++) ins(i,ed,h,0); 107 } 108 109 int main() 110 { 111 init(); 112 // output(); 113 // return 0; 114 int ans; 115 ans=max_flow(); 116 printf("%d\n",ans); 117 return 0; 118 }
2016-11-06 21:19:22