/Users/jst/src/grafalgo/cpp/graphAlgorithms/maxFlo/checkMf.cpp

#include "stdinc.h"
#include "Mflograph.h"
#include "List.h"

using namespace grafalgo;

main() {
        vertex u,v; edge e; int sum;
        Mflograph fg; fg.read(cin);

        // verify that capacity constraints and min flow requirements 
        // are respected
        for (e = fg.first(); e != 0; e = fg.next(e)) {
                u = fg.tail(e); v = fg.head(e); string s;
                if (fg.f(u,e) < 0)
                        cout << "Negative flow on edge " 
                             << e << "=" << fg.edge2string(e,s) << endl;
                if (fg.f(u,e) > fg.cap(u,e))
                        cout << "Flow exceeds capacity on edge "
                             << e << "=" << fg.edge2string(e,s) << endl;
                if (fg.f(u,e) < fg.minFlo(e))
                        cout << "Flow less than min flow requirement on edge "
                             << e << "=" << fg.edge2string(e,s) << endl;
        }

        // verify that flow at each node is balanced
        for (u = 1; u <= fg.n(); u++) { 
                if (u == fg.src() || u == fg.snk()) continue;
                sum = 0;
                for (e = fg.firstAt(u); e != 0; e = fg.nextAt(u,e)) {
                        sum -= fg.f(u,e);
                }
                if (sum != 0) cout << "Vertex " << u << " is not balanced\n";
        }

        // Verify that the flow is maximum by computing hop-count distance
        // over all edges that are not saturated. If sink ends up with a
        // distance smaller than n, then there is a path to sink with
        // non-zero residual capacity.
        int *d = new int[fg.n()+1];
        for (u = 1; u <= fg.n(); u++) d[u] = fg.n();
        d[fg.src()] = 0;
        List q(fg.n()); q.addLast(fg.src());
        while (!q.empty()) {
                u = q.first(); q.removeFirst();
                for (e = fg.firstAt(u); e != 0; e = fg.nextAt(u,e)) {
                        v = fg.mate(u,e);
                        if (fg.res(u,e) > 0 && d[v] > d[u] + 1) {
                                d[v] = d[u] + 1; q.addLast(v);
                        }
                }
        }
        if (d[fg.snk()] < fg.n()) printf("Not a maximum flow\n");
}