summaryrefslogtreecommitdiff
path: root/test/math/discreteLogarithm.cpp
diff options
context:
space:
mode:
authormzuenni <mzuenni@users.noreply.github.com>2024-07-28 22:54:40 +0200
committerGitHub <noreply@github.com>2024-07-28 22:54:40 +0200
commit8d11c6c8213f46f0fa19826917c255edd5d43cb1 (patch)
tree96d75baff33d5a04b5a60f1a41f514a26c716874 /test/math/discreteLogarithm.cpp
parent8c33b4e0d3030cfed17fc64b4fe41133339f6d87 (diff)
Test (#4)
* update * moved content in subdir * rename file * add test setup * add test setup * add github action * automaticly test all cpp files * timeout after 10s * setulimit and dont zero memory * test build pdf * install latexmk * update * update * ngerman * fonts * removed old code * add first test * added tests * test in sorted order * more tests * simplified test * more tests * fix suffix tree * fixes and improvements * done ust lst directly * fix swap * add links to pdf * fix constants * add primorial * add comment * various improvements * more tests * added missing stuf * more tests * fix tests * more tests * more tests * more tests * fix recursion? * test trie * more tests * only use python temporarily for listings * only use python temporarily for listings * more tests * fix longestCommonSubstring * more tests * more tests * made code more similiar * fix? * more tests * more tests * more tests * add ahoCorasick test + limit 4GB stack size * more tests * fix test * add additional test * more tests * more tests * fix? * better fix * fix virtual tree * more tests * more tests * recursive closest pair * more tests * decrease limit * new tests * more tests * fix name * more tests * add test * new test * more tests * more tests * more tests * more tests * new test and content * new code * new code * larger tests * fix and test * new test * new test * update pdf * remove comments * new test * more tests * more testcases * more tests * increased limit * more tests * more tests * more tests * new tests * more tests * shortened code * new test * add basic tests for bigint * more tests * removed old files * new test * ignore some files * more auto more ccw * fix test * more tests * fix * new tests * more tests * more tests * stronger test * actually verify delaunay... * more tests * fix header * more tests * run tests parallel? * test parralel? * add --missing * separate workflows * test * is the pdf checked? * separate workflows * fix workflow * more workflows --------- Co-authored-by: Yidi <noob999noob999@gmail.com>
Diffstat (limited to 'test/math/discreteLogarithm.cpp')
-rw-r--r--test/math/discreteLogarithm.cpp64
1 files changed, 64 insertions, 0 deletions
diff --git a/test/math/discreteLogarithm.cpp b/test/math/discreteLogarithm.cpp
new file mode 100644
index 0000000..0f9eecf
--- /dev/null
+++ b/test/math/discreteLogarithm.cpp
@@ -0,0 +1,64 @@
+#include "../util.h"
+#include <math/modPowIterativ.cpp>
+#include <math/legendre.cpp>
+
+ll overwrite = 0;
+ll getMemory(ll /**/) {return overwrite - 1;} //dlog code adds one...
+#define sqrtl getMemory
+#include <math/discreteLogarithm.cpp>
+#undef sqrtl
+
+template<typename F>
+void stress_test(F&& f) {
+ ll work = 0;
+ for (ll tries = 0; tries < 3'000; tries++) {
+ ll p = Random::prime<ll>(1'000);
+ overwrite = f(p);
+ ll a = Random::integer<ll>(1, p);
+ vector<bool> naive(p);
+ for (ll i = 0, j = 1; i < p; i++, j = (j * a) % p) {
+ naive[j] = true;
+ }
+ for (ll b = 0; b < p; b++) {
+ ll got = dlog(a, b, p);
+ if (got < -1 || got >= p) cerr << "error: out of range" << FAIL;
+ if ((got >= 0) != naive[b]) {
+ cerr << a << " " << b << " " << p << endl;
+ cerr << got << endl;
+ cerr << "error" << FAIL;
+ }
+ if (got >= 0 && powMod(a, got, p) != b) {
+ cerr << a << "^" << got << " = " << powMod(a, got, p) << " != " << b << " % " << p << endl;
+ cerr << "error: wrong" << FAIL;
+ }
+ work++;
+ }
+ }
+ cerr << "stress tested: " << work << endl;
+}
+
+constexpr int N = 25;
+constexpr ll mod = 1'394'633'899;
+void performance_test() {
+ timer t;
+ hash_t hash = 0;
+ overwrite = sqrt(mod);
+ for (int operations = 0; operations < N; operations++) {
+ ll a = Random::integer<ll>(1, mod);
+ ll b = Random::integer<ll>(0, mod);
+ t.start();
+ hash += dlog(a, b, mod);
+ t.stop();
+ }
+ if (t.time > 500) cerr << "too slow: " << t.time << FAIL;
+ cerr << "tested performance: " << t.time << "ms (hash: " << hash << ")" << endl;
+}
+
+
+int main() {
+ stress_test([](ll p){return sqrtl(p);});
+ stress_test([](ll p){return min<ll>(10, p - 1);});
+ stress_test([](ll p){return min<ll>(p - 1, sqrtl(p) + 100);});
+ performance_test();
+}
+