USACO算法系列二十八——ditch

    题目：http://www.nocow.cn/index.php/Translate:USACO/ditch

    最近这段时间天气变化比较厉害，一不小心得了重感冒，每天头脑都昏昏沉沉的，然后新的一节是关于最大流的。以前还没有详细的看过这方面的资料，所以最近找了本算法导论，看了一下图论里面的最大流算法。觉得这本书讲的比网上的资料思路清晰一些，网上很多资料都是概述性的，知其然，不知所以然。所以我建议还是看看算法导论这本书，但是我发现这本书，还是有些错误的，第401页的图26-3 b)我觉得少了一条v3到v2的箭头，数值是5，然后画重点线的v3到v2的方向反了，数值是4。不知道是中文版+盗版的问题（穷啊），英文版就这样，另外405页上面的那段伪代码，也有一些需要商议的地方。

    总之，最大流算法的解法思想，网上随便搜索一下，应该很多资料。我就不多费唇舌了。

    这道题就是一个纯粹的最大流算法。

    我采用的是Ford-Fulkerson算法。不过里面的搜索增广路径使用的是简单的深度搜索而已，没有搜索完全部的路径，然后选择最大的那个路径。可惜的是在Test7的时候，超时了，还是有点晕啊，先休息一下。晚上再考虑。

    代码如下：

#include <iostream> #include <fstream> #define SIZE 201 #define LARGER 0X7FFFFFFF #define min(a,b) (a)<(b)?(a):(b) using namespace std; int n, m; int map[SIZE][SIZE] = {0}; int path[SIZE] = {0}; int result = 0; //读写文件流 ifstream fin("ditch.in"); ofstream fout("ditch.out"); bool isinpath(int pos,int length) { for (int i=0; i < length; i ++) { if (path[i] == pos) { return true; } } return false; } int ditchMap(int pos,int step, int ci) { if (pos == m) { for (int i=0; i < step-1; i ++) { map[path[i]][path[i+1]] -= ci; map[path[i+1]][path[i]] += ci; } //end for return ci; }//end if else { for (int i=1; i <= m; i ++) { if (map[pos][i] > 0) { if (isinpath(i,step)) { continue; } path[step] = i; int f= ditchMap(i, step+1, min(ci,map[pos][i])); if (f > 0) { return f; }//end if f>0 }//end if map[pos][i] }//end for }//end else return -1; } int main() { //输入数据 fin >> n >> m; for (int i=0; i < n; i ++) { int row =0; int col = 0; fin >> row >> col ; fin >> map[row][col]; } path[0] = 1; int temp = ditchMap(1, 1, LARGER); while (temp > 0) { result += temp; temp = ditchMap(1, 1, LARGER); } fout << result <<endl; return 0; }

    结果如下：

Executing... Test 1: TEST OK [0.011 secs, 3168 KB] Test 2: TEST OK [0.011 secs, 3168 KB] Test 3: TEST OK [0.011 secs, 3168 KB] Test 4: TEST OK [0.000 secs, 3168 KB] Test 5: TEST OK [0.000 secs, 3168 KB] Test 6: TEST OK [0.011 secs, 3168 KB] > Run 7: Execution error: Your program (`ditch') used more than the allotted runtime of 1 seconds (it ended or was stopped at 1.674 seconds) when presented with test case 7. It used 3168 KB of memory.