deba@522: /* -*- mode: C++; indent-tabs-mode: nil; -*- deba@522: * deba@522: * This file is a part of LEMON, a generic C++ optimization library. deba@522: * alpar@956: * Copyright (C) 2003-2010 deba@522: * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport deba@522: * (Egervary Research Group on Combinatorial Optimization, EGRES). deba@522: * deba@522: * Permission to use, modify and distribute this software is granted deba@522: * provided that this copyright notice appears in all copies. For deba@522: * precise terms see the accompanying LICENSE file. deba@522: * deba@522: * This software is provided "AS IS" with no warranty of any kind, deba@522: * express or implied, and with no claim as to its suitability for any deba@522: * purpose. deba@522: * deba@522: */ deba@522: deba@522: #include deba@522: #include deba@522: #include deba@522: #include deba@522: deba@522: #include deba@522: #include deba@522: #include kpeter@672: #include deba@522: deba@522: #include "test_tools.h" deba@522: deba@522: using namespace lemon; deba@522: using namespace std; deba@522: deba@522: const char test_lgf[] = deba@522: "@nodes\n" deba@522: "label\n" deba@522: "0\n" deba@522: "1\n" deba@522: "2\n" deba@522: "3\n" deba@522: "4\n" deba@522: "5\n" deba@522: "6\n" deba@522: "7\n" deba@522: "8\n" deba@522: "9\n" deba@522: "@arcs\n" deba@522: " label cost\n" deba@522: "1 8 0 107\n" deba@522: "0 3 1 70\n" deba@522: "2 1 2 46\n" deba@522: "4 1 3 28\n" deba@522: "4 4 4 91\n" deba@522: "3 9 5 76\n" deba@522: "9 8 6 61\n" deba@522: "8 1 7 39\n" deba@522: "9 8 8 74\n" deba@522: "8 0 9 39\n" deba@522: "4 3 10 45\n" deba@522: "2 2 11 34\n" deba@522: "0 1 12 100\n" deba@522: "6 3 13 95\n" deba@522: "4 1 14 22\n" deba@522: "1 1 15 31\n" deba@522: "7 2 16 51\n" deba@522: "2 6 17 29\n" deba@522: "8 3 18 115\n" deba@522: "6 9 19 32\n" deba@522: "1 1 20 60\n" deba@522: "0 3 21 40\n" deba@522: "@attributes\n" deba@522: "source 0\n"; deba@522: kpeter@672: kpeter@672: void checkMinCostArborescenceCompile() kpeter@672: { kpeter@672: typedef double VType; kpeter@672: typedef concepts::Digraph Digraph; kpeter@672: typedef concepts::ReadMap CostMap; kpeter@672: typedef Digraph::Node Node; kpeter@672: typedef Digraph::Arc Arc; kpeter@672: typedef concepts::WriteMap ArbMap; kpeter@672: typedef concepts::ReadWriteMap PredMap; kpeter@672: kpeter@672: typedef MinCostArborescence:: kpeter@672: SetArborescenceMap:: kpeter@672: SetPredMap::Create MinCostArbType; kpeter@672: kpeter@672: Digraph g; kpeter@672: Node s, n; kpeter@672: Arc e; kpeter@672: VType c; kpeter@672: bool b; alpar@1257: ::lemon::ignore_unused_variable_warning(c,b); kpeter@672: int i; kpeter@672: CostMap cost; kpeter@672: ArbMap arb; kpeter@672: PredMap pred; kpeter@672: kpeter@672: MinCostArbType mcarb_test(g, cost); kpeter@672: const MinCostArbType& const_mcarb_test = mcarb_test; kpeter@672: kpeter@672: mcarb_test kpeter@672: .arborescenceMap(arb) kpeter@672: .predMap(pred) kpeter@672: .run(s); kpeter@672: kpeter@672: mcarb_test.init(); kpeter@672: mcarb_test.addSource(s); kpeter@672: mcarb_test.start(); kpeter@672: n = mcarb_test.processNextNode(); kpeter@672: b = const_mcarb_test.emptyQueue(); kpeter@672: i = const_mcarb_test.queueSize(); alpar@956: kpeter@672: c = const_mcarb_test.arborescenceCost(); kpeter@672: b = const_mcarb_test.arborescence(e); kpeter@672: e = const_mcarb_test.pred(n); kpeter@672: const MinCostArbType::ArborescenceMap &am = kpeter@672: const_mcarb_test.arborescenceMap(); kpeter@672: const MinCostArbType::PredMap &pm = kpeter@672: const_mcarb_test.predMap(); kpeter@672: b = const_mcarb_test.reached(n); kpeter@672: b = const_mcarb_test.processed(n); alpar@956: kpeter@672: i = const_mcarb_test.dualNum(); kpeter@672: c = const_mcarb_test.dualValue(); kpeter@672: i = const_mcarb_test.dualSize(i); kpeter@672: c = const_mcarb_test.dualValue(i); alpar@956: alpar@1257: ::lemon::ignore_unused_variable_warning(am); alpar@1257: ::lemon::ignore_unused_variable_warning(pm); kpeter@672: } kpeter@672: deba@522: int main() { deba@522: typedef SmartDigraph Digraph; deba@522: DIGRAPH_TYPEDEFS(Digraph); deba@522: deba@522: typedef Digraph::ArcMap CostMap; deba@522: deba@522: Digraph digraph; deba@522: CostMap cost(digraph); deba@522: Node source; deba@522: deba@522: std::istringstream is(test_lgf); deba@522: digraphReader(digraph, is). deba@522: arcMap("cost", cost). deba@522: node("source", source).run(); deba@522: deba@522: MinCostArborescence mca(digraph, cost); deba@522: mca.run(source); deba@522: deba@522: vector > > dualSolution(mca.dualNum()); deba@522: deba@522: for (int i = 0; i < mca.dualNum(); ++i) { deba@522: dualSolution[i].first = mca.dualValue(i); deba@522: for (MinCostArborescence::DualIt it(mca, i); deba@522: it != INVALID; ++it) { deba@522: dualSolution[i].second.insert(it); deba@522: } deba@522: } deba@522: deba@522: for (ArcIt it(digraph); it != INVALID; ++it) { deba@522: if (mca.reached(digraph.source(it))) { deba@522: double sum = 0.0; deba@522: for (int i = 0; i < int(dualSolution.size()); ++i) { deba@522: if (dualSolution[i].second.find(digraph.target(it)) deba@522: != dualSolution[i].second.end() && deba@522: dualSolution[i].second.find(digraph.source(it)) deba@522: == dualSolution[i].second.end()) { deba@522: sum += dualSolution[i].first; deba@522: } deba@522: } deba@522: if (mca.arborescence(it)) { kpeter@672: check(sum == cost[it], "Invalid dual solution"); deba@522: } kpeter@672: check(sum <= cost[it], "Invalid dual solution"); deba@522: } deba@522: } deba@522: deba@522: kpeter@672: check(mca.dualValue() == mca.arborescenceCost(), "Invalid dual solution"); deba@522: kpeter@672: check(mca.reached(source), "Invalid arborescence"); deba@522: for (ArcIt a(digraph); a != INVALID; ++a) { deba@522: check(!mca.reached(digraph.source(a)) || kpeter@672: mca.reached(digraph.target(a)), "Invalid arborescence"); deba@522: } deba@522: deba@522: for (NodeIt n(digraph); n != INVALID; ++n) { deba@522: if (!mca.reached(n)) continue; deba@522: int cnt = 0; deba@522: for (InArcIt a(digraph, n); a != INVALID; ++a) { deba@522: if (mca.arborescence(a)) { kpeter@672: check(mca.pred(n) == a, "Invalid arborescence"); deba@522: ++cnt; deba@522: } deba@522: } kpeter@672: check((n == source ? cnt == 0 : cnt == 1), "Invalid arborescence"); deba@522: } deba@522: deba@522: Digraph::ArcMap arborescence(digraph); kpeter@672: check(mca.arborescenceCost() == deba@522: minCostArborescence(digraph, cost, source, arborescence), kpeter@672: "Wrong result of the function interface"); deba@522: deba@522: return 0; deba@522: }