关于网络流算法(2)

 1 /*
 2 不断寻找残量网络中汇点的可达路径，
 3 可达路径中的各条边每次迭代减少可达路径中所有边 capacity 最小值直到在残量网络中汇点不可达 
 4 */
 5 
 6 #include<iostream>
 7 #include<queue>
 8 using namespace std;
 9 const int arrsize = 10;
10 const int maxdata = INT_MAX;
11 
12 int cap[arrsize][arrsize];        //表示两个结点之间的最大可用流量 
13 int flow[arrsize],pre[arrsize];    //表示从原点到index结点的当前剩余可用流量 
14 queue<int> q;
15 
16 
17 int BFS(int src, int des)
18 //BFS函数，寻找可行通路 
19 {
20 //cout<<endl<<"进入BFS"<<endl; 
21     while(!q.empty())
22         q.pop();
23     
24     memset(pre,-1,sizeof(pre));
25     
26     pre[src] = 0;
27     flow[src] = maxdata;
28     
29     q.push(src);
30     while(!q.empty())
31     {
32         int index = q.front();
33         q.pop();
34         
35         if(index == des)
36             break;    
37         //层序寻找下一个可行通路
38         for(int i=0;i<arrsize;i++)
39         {
40             if(i!=src && cap[index][i]>0 && pre[i]==-1)
41             {
42                 pre[i]=index;                            //  保存当前结点的父结点 
43                 flow[i]=min(cap[index][i],flow[index]);    //★更新从源结点到当前结点的当前剩余可用流量 
44                 q.push(i);
45             }
46         }
47     }
48     if(pre[des]==-1)
49         return -1;
50     else
51         return flow[des];
52 }
53 
54 int maxFlow(int src, int des)
55 //s最大流函数，迭代残量网络 
56 {
57     int increasement,sumfolw = 0;
58     while((increasement = BFS(src,des)) != -1)
59     {        
60         int cur = des;
61         while(cur != src)
62         {
63             int last = pre[cur];
64             cap[last][cur] -= increasement;
65             cap[cur][last] += increasement;
66             cur=last; 
67         }
68         sumfolw += increasement;
69     }
70     return sumfolw;
71 }
72 
73 
74 
75 int main()
76 {
77     freopen("C:\\Users\\lenovo\\Desktop\\工作\\华为挑战赛\\数据_E-K.txt","r",stdin);
78     
79     for(short i=0;i<arrsize;i++)
80     for(short j=0;j<arrsize;j++)
81         cin>>cap[i][j];
82     
83     //打印邻接矩阵 
84     cout<<"邻接矩阵(有向图)为："<<endl;
85     for(short i=0;i<arrsize;i++)
86     for(short j=0;j<arrsize;j++)
87     {
88         cout<<cap[i][j]<<" ";
89         if(j==arrsize-1)
90             cout<<endl;
91     }
92     cout<<endl;
93     
94     memset(flow,0,sizeof(flow)); 
95     cout<<"最大流为："<<maxFlow(0,3)<<endl<<endl;
96 
97     fclose(stdin);
98 }