1. 程式人生 > >poj(1273)(最大流)

poj(1273)(最大流)

#include"string.h"
#include"stdio.h"
#include"queue"
#define inf 9999999
using namespace std;
int r[300][300];//殘留網路,初始化為原圖
int pre[300];//記錄流的上一個結點,即從pre[i]流向i
int visit[300];記錄是否訪問
int n,m;
int bfs(int s,int t)//尋找一條從s到t的增廣路,若找到返回1
{
 int p,i;
 queue<int>q;
 memset(visit,0,sizeof(visit));//初始化必須放在這裡,因為在EK中要多次用到bfs,
 memset(pre,0,sizeof(pre));
 pre[s]=s;
 visit[s]=1;
 q.push(s);
 while(!q.empty())
 {
  p=q.front();
  q.pop();
  for(i=1;i<=m;i++)
  {
   if(visit[i]==0&&r[p][i])//滿足條件則入隊。
   {
    pre[i]=p;
    visit[i]=1;
    if(i==t)
     return 1;
    q.push(i);
   }
  }
 }
 return 0;
}
int EK(int s,int t))//s是源點t是匯點也是節點的最大值,計算從s到t的最大流
{
 int  d,i,flow=0,u;
 while(bfs(s,t))
 {
  d=inf;
  for(i=t;i!=s;i=pre[i])//找出所選路徑的最小r殘留值
   d=d<r[pre[i]][i]?d:r[pre[i]][i];
  for(i=t;i!=s;i=pre[i])//正向減去d,反向增加d值,相當於回朔。。
  {
   r[pre[i]][i]-=d;
   r[i][pre[i]]+=d;
  }
  flow+=d;
 }
 return flow;
}
int main()
{
 int i,j,k,h;
 while(scanf("%d%d",&n,&m)!=EOF)
 {
  int a,b,c;
  memset(r,0,sizeof(r));
  for(i=0;i<n;i++)
  {
   scanf("%d%d%d",&a,&b,&c);
   r[a][b]+=c;
  }
  printf("%d\n",EK(1,m));
 }
 return 0;
}