【codevs1993】草地排水 最大流
【codevs1993】草地排水
题目描述 Description
在农夫约翰的农场上，每逢下雨，Bessie最喜欢的三叶草地就积聚了一潭水。这意味着草地被水淹没了，并且小草要继续生长还要花相当长一段时间。因此，农夫约翰修建了一套排水系统来使贝茜的草地免除被大水淹没的烦恼（不用担心，雨水会流向附近的一条小溪）。作为一名一流的技师，农夫约翰已经在每条排水沟的一端安上了控制器，这样他可以控制流入排水沟的水流量。
农夫约翰知道每一条排水沟每分钟可以流过的水量，和排水系统的准确布局（起点为水潭而终点为小溪的一张网）。需要注意的是，有些时候从一处到另一处不只有一条排水沟。
根据这些信息，计算从水潭排水到小溪的最大流量。对于给出的每条排水沟，雨水只能沿着一个方向流动，注意可能会出现雨水环形流动的情形。
输入描述 Input Description
第1行: 两个用空格分开的整数N (0 <= N <= 200) 和 M (2 <= M <= 200)。N是农夫John已经挖好的排水沟的数量，M是排水沟交叉点的数量。交点1是水潭，交点M是小溪。
第二行到第N+1行: 每行有三个整数，Si, Ei, 和 Ci。Si 和 Ei (1 <= Si, Ei <= M) 指明排水沟两端的交点，雨水从Si 流向Ei。Ci (0 <= Ci <= 10,000,000)是这条排水沟的最大容量。
输出描述 Output Description
输出一个整数，即排水的最大流量。
样例输入 Sample Input
5 4
1 2 40
1 4 20
2 4 20
2 3 30
3 4 10
样例输出 Sample Output
50
代码
ek算法
#include <cstdio> #include <cmath> #include <cstring> #include <ctime> #include <iostream> #include <algorithm> #include <set> #include <vector> #include <queue> #include <typeinfo> #include <map> #include <stack> typedef long long ll; #define inf 0x7fffffff using namespace std; inline ll read() { ll x=0,f=1; char ch=getchar(); while(ch<'0'||ch>'9') { if(ch=='-')f=-1; ch=getchar(); } while(ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); } return x*f; } //************************************************************************************** int pre[201]; int mp[201][201],flow[201][201]; int n,m; int father[201]; int ek(int s,int y) { int q[201]; int u; int f=0; while(1) { int di=1; int front=0; memset(pre,0,sizeof(pre)); pre[s]=0x7fffffff; q[front]=s; while(front<di) { u=q[front]; front++; for(int v=1;v<=n;v++) { if(pre[v]==0&&mp[u][v]>flow[u][v]) { father[v]=u; q[di++]=v; pre[v]=min(pre[u],mp[u][v]-flow[u][v]); } } } if(pre[y]==0)return f; for(u=y;u!=s;u=father[u]) { flow[father[u]][u]+=pre[y]; flow[u][father[u]]-=pre[y]; } f+=pre[y]; } } int main() { scanf("%d%d",&m,&n); for(int i=1;i<=m;i++) { int u,v,w; scanf("%d%d%d",&u,&v,&w); mp[u][v]+=w; } printf("%d\n",ek(1,n)); /* for(int i=1;i<=n;i++) { for(int j=1;j<=n;j++) printf("%d ",flow[i][j]); printf("\n"); }*/ return 0; }
dinic矩阵版
#include <cstdio> #include <cmath> #include <cstring> #include <ctime> #include <iostream> #include <algorithm> #include <set> #include <vector> #include <queue> #include <typeinfo> #include <map> #include <stack> typedef long long ll; using namespace std; inline ll read() { ll x=0,f=1; char ch=getchar(); while(ch<'0'||ch>'9') { if(ch=='-')f=-1; ch=getchar(); } while(ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); } return x*f; } //************************************************************************************** using namespace std; const int MAX=0x5fffffff;// int tab[250][250];//邻接矩阵 int dis[250];//距源点距离,分层图 int q[2000],h,r;//BFS队列 ,首,尾 int N,M,ANS;//N:点数;M,边数 int BFS() { int i,j; memset(dis,0xff,sizeof(dis));//以-1填充 dis[1]=0; h=0; r=1; q[1]=1; while (h<r) { j=q[++h]; for (i=1; i<=N; i++) if (dis[i]<0 && tab[j][i]>0) { dis[i]=dis[j]+1; q[++r]=i; } } if (dis[N]>0) return 1; else return 0;//汇点的DIS小于零,表明BFS不到汇点 } //Find代表一次增广,函数返回本次增广的流量,返回0表示无法增广 int find(int x,int low)//Low是源点到现在最窄的(剩余流量最小)的边的剩余流量 { int i,a=0; if (x==N)return low;//是汇点 for (i=1; i<=N; i++) if (tab[x][i] >0 //联通 && dis[i]==dis[x]+1 //是分层图的下一层 &&(a=find(i,min(low,tab[x][i]))))//能到汇点(a <> 0) { tab[x][i]-=a; tab[i][x]+=a; return a; } return 0; } int main() { int i,j,f,t,flow,tans; while (scanf("%d%d",&M,&N)!=EOF) { memset(tab,0,sizeof(tab)); for (i=1; i<=M; i++) { scanf("%d%d%d",&f,&t,&flow); tab[f][t]+=flow; } // ANS=0; while (BFS())//要不停地建立分层图,如果BFS不到汇点才结束 { while(tans=find(1,0x7fffffff))ANS+=tans;//一次BFS要不停地找增广路,直到找不到为止 } printf("%d\n",ANS); } }
dinic灵界表
#include <cstdio> #include <cmath> #include <cstring> #include <ctime> #include <iostream> #include <algorithm> #include <set> #include <vector> #include <queue> #include <typeinfo> #include <map> #include <stack> typedef long long ll; #define inf 0x7fffffff using namespace std; inline ll read() { ll x=0,f=1; char ch=getchar(); while(ch<'0'||ch>'9') { if(ch=='-')f=-1; ch=getchar(); } while(ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); } return x*f; } //************************************************************************************** struct ss { int to,v,next; } e[501]; int ans,q[501]; int head[501]; int vis[501]; int n,m; int t=2; void add(int u,int v,int w) { e[t].to=v; e[t].next=head[u]; e[t].v=w; head[u]=t++; } void insert(int u,int v,int w) { add(u,v,w); add(v,u,0); } bool bfs() { memset(vis,-1,sizeof(vis)); int di=1,front=0,now; q[front]=1; vis[1]=0; while(front<di) { now=q[front]; front++; if(front==501)front=0; for(int i=head[now]; i; i=e[i].next) { if(e[i].v&&vis[e[i].to]==-1) { vis[e[i].to]=vis[now]+1; q[di++]=e[i].to; if(di==501)di=0; } } } if(vis[n]==-1)return 0; else return 1; } int dfs(int x,int f) { if(x==n) return f; int w,used=0; for(int i=head[x]; i; i=e[i].next) { if(e[i].v&&vis[e[i].to]==vis[x]+1) { if(w=dfs(e[i].to,min(f,e[i].v))) { e[i].v-=w; e[i+1].v+=w; return w; } } } return 0; } void dinic() { while(bfs())ans+=dfs(1,inf); } int main() { scanf("%d%d",&m,&n); for(int i=1; i<=m; i++) { int a,b,c; scanf("%d%d%d",&a,&b,&c); insert(a,b,c); } dinic(); printf("%d\n",ans); return 0; }
