最大流(Dinic模版)——HDOJ 1532

来源:互联网 发布:域名实名不通过 编辑:程序博客网 时间:2024/06/07 05:55

HDOJ 1532 Drainage Ditches

这道题目是用来做Dinic算法的模版题目,而上一篇博客 最大流——HDOJ 3549 是用来作为EK算法的模版题目。

关于Dinic算法的讲解,网上有很多,我这里也不多讲什么。推荐一篇讲得比较容易懂的 http://comzyh.tk/blog/archives/568/#Dinic-Code

/*HDOJ 1532最大流经典入门题目 下面是Dinic算法,矩阵存储图,易于理解*/#include <iostream>using namespace std;#define min(a,b) (a)<(b)?a:b#define MAX 0x5fffffffint map[250][250];//邻接矩阵 int dis[250];//距源点距离,分层图 int que[2000],front,rear;//BFS队列 ,首,尾 int N,M,ANS;//N:点数;M,边数int BFS(){int i,j;memset(dis,0xff,sizeof(dis));//以-1填充 dis[1]=0;front=rear=0;que[0]=1;while (front <= rear){j=que[front++];for (i=1;i<=N;i++)if (dis[i]<0 && map[j][i]>0){dis[i]=dis[j]+1; que[++rear]=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 (map[x][i] >0 //联通 && dis[i]==dis[x]+1 //是分层图的下一层 &&(a=find(i,min(low,map[x][i]))))//能到汇点(a <> 0) {map[x][i] -= a;map[i][x] += a;return a;}return 0;}int main(){    int i,j,f,t,flow,tans;    while (scanf("%d%d",&M,&N)!=EOF){memset(map,0,sizeof(map));for(i=1;i<=M;i++){scanf("%d%d%d",&f,&t,&flow);map[f][t] += flow;}ANS=0;while(BFS())//要不停地建立分层图,如果BFS不到汇点才结束 {while(tans=find(1,0x7fffffff))ANS+=tans;//一次BFS要不停地找增广路,直到找不到为止 }printf("%d\n",ANS);    }    return 0;}


原创粉丝点击