summaryrefslogtreecommitdiff
path: root/test/graph/dinitzScaling.cpp
diff options
context:
space:
mode:
authorGloria Mundi <gloria@gloria-mundi.eu>2025-02-13 22:07:35 +0100
committerGloria Mundi <gloria@gloria-mundi.eu>2025-02-13 22:07:35 +0100
commit04ca8f7bd16c0c855f604188d617a1bf2e8eacfd (patch)
treec0b535dc4635bf7d99c714f030578c73ba023310 /test/graph/dinitzScaling.cpp
parent7dc90ef744cf16ac4b4cb4e7d22f4c4686ae7225 (diff)
rename Dinic to Dinitz
Diffstat (limited to 'test/graph/dinitzScaling.cpp')
-rw-r--r--test/graph/dinitzScaling.cpp61
1 files changed, 61 insertions, 0 deletions
diff --git a/test/graph/dinitzScaling.cpp b/test/graph/dinitzScaling.cpp
new file mode 100644
index 0000000..c06d486
--- /dev/null
+++ b/test/graph/dinitzScaling.cpp
@@ -0,0 +1,61 @@
+#include "../util.h"
+namespace dinitz {
+#include <graph/dinitzScaling.cpp>
+}
+
+namespace pushRelabel {
+#include <graph/pushRelabel.cpp>
+}
+
+void stress_test() {
+ ll queries = 0;
+ for (int tries = 0; tries < 20'000; tries++) {
+ int n = Random::integer<int>(2, 30);
+ int m = Random::integer<int>(n-1, max<int>(n, min<int>(500, n*(n-1) / 2 + 1)));
+
+ dinitz::adj.assign(n, {});
+ pushRelabel::adj.assign(n, {});
+
+ Graph<NoData, true> g(n);
+ g.erdosRenyi(m);
+ g.forEdges([](int a, int b){
+ ll w = Random::integer<ll>(1, 1'000'000'000'000ll);
+ dinitz::addEdge(a, b, w);
+ pushRelabel::addEdge(a, b, w);
+ });
+
+ ll got = dinitz::maxFlow(0, n - 1);
+ ll expected = pushRelabel::maxFlow(0, n - 1);
+
+ if (got != expected) cerr << "got: " << got << ", expected: " << expected << FAIL;
+ queries += n;
+ }
+ cerr << "tested random queries: " << queries << endl;
+}
+
+constexpr int N = 50000;
+constexpr int M = 200000;
+void performance_test() {
+ using namespace dinitz;
+ timer t;
+ Graph<NoData> g(N);
+ g.erdosRenyi(M);
+ adj.assign(N, {});
+ g.forEdges([](int a, int b){
+ ll w1 = Random::integer<ll>(1, 1'000'000'000'000ll);
+ ll w2 = Random::integer<ll>(1, 1'000'000'000'000ll);
+ addEdge(a, b, w1);
+ addEdge(b, a, w2);
+ });
+
+ t.start();
+ hash_t hash = maxFlow(0, N - 1);
+ t.stop();
+ if (t.time > 2000) cerr << "too slow: " << t.time << FAIL;
+ cerr << "tested performance: " << t.time << "ms (hash: " << hash << ")" << endl;
+}
+
+int main() {
+ stress_test();
+ performance_test();
+}