11. 运输问题1
★★☆ 输入文件:
maxflowa.in
输出文件:maxflowa.out
简单对比
时间限制:1 s 内存限制:128 MB【问题描述】一个工厂每天生产若干商品,需运输到销售部门进行销售。从产地到销地要经过某些城镇,有不同的路线可以行走,每条两城镇间的公路都有一定的流量限制。请你计算,在不考虑其它车辆使用公路的前提下,如何充分利用所有的公路,使产地运输到销地的商品最多,最多能运输多少商品。【输入格式】输入文件有若干行
第一行,一个整数n,表示共有n个城市$(2 \leq n \leq 100)$,产地是1号城市,销地是n号城市。
下面有n行,每行有n个数字。第p行第q列的数字表示城镇p与城镇q之间有无公路连接。数字为0表示无,大于0表示有公路,且该数字表示该公路流量。【输出格式】输出文件有一行
第一行,1个整数max,表示最大流量为max。【输入输出样例】输入文件名: maxflowa.in6
0 4 8 0 0 0
0 0 4 4 1 0
0 0 0 2 2 0
0 0 0 0 0 7
0 0 0 6 0 9
0 0 0 0 0 0输出文件名:maxflowa.out8
好了对于这么水的题我表示只是来存代码的(逃
注意读入建图的时候的反向边标记然后大力Dinic就好了w
#include <queue>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm> const int MAXV=;
const int MAXE=;
const int INF=0x7FFFFFFF; struct Edge{
int from;
int to;
int flow;
Edge* rev;
Edge* next;
};
Edge E[MAXE];
Edge* head[MAXV];
Edge* top=E; int v;
int depth[MAXV];
int flow[MAXV][MAXV]; int Dinic(int,int);
void Initialize();
void Insert(int,int);
int DFS(int,int,int);
bool BFS(int,int); int main(){
freopen("maxflowa.in","r",stdin);
freopen("maxflowa.out","w",stdout);
Initialize();
printf("%d\n",Dinic(,v));
return ;
} int Dinic(int s,int t){
int ans=;
while(BFS(s,t)){
ans+=DFS(s,INF,t);
}
return ans;
} int DFS(int s,int flow,int t){
if(s==t||flow==)
return flow;
int tmp=flow;
int k;
for(Edge* i=head[s];i!=NULL;i=i->next){
if(i->flow!=&&tmp!=&&depth[i->to]==depth[s]+){
k=DFS(i->to,std::min(tmp,i->flow),t);
if(k==){
depth[i->to]=;
continue;
}
tmp-=k;
i->flow-=k;
i->rev->flow+=k;
if(tmp==)
break;
}
}
return flow-tmp;
} bool BFS(int s,int t){
memset(depth,,sizeof(depth));
std::queue<int> q;
depth[s]=;
q.push(s);
while(!q.empty()){
s=q.front();
q.pop();
for(Edge* i=head[s];i!=NULL;i=i->next){
if(depth[i->to]==&&i->flow!=){
q.push(i->to);
depth[i->to]=depth[s]+;
if(i->to==t)
return true;
}
}
}
return false;
} void Initialize(){
scanf("%d",&v);
for(int i=;i<=v;i++){
for(int j=;j<=v;j++){
scanf("%d",&flow[i][j]);
}
}
for(int i=;i<=v;i++){
for(int j=;j<i;j++){
Insert(i,j);
}
}
} inline void Insert(int a,int b){
top->from=a;
top->to=b;
top->flow=flow[a][b];
top->next=head[a];
head[a]=top;
top->rev=top+;
top++;
top->from=b;
top->to=a;
top->flow=flow[b][a];
top->next=head[b];
head[b]=top;
top->rev=top-;
top++;
}
Backup