1 files changed, 31 insertions, 0 deletions
diff --git a/test/math/binomial3.cpp b/test/math/binomial3.cpp
new file mode 100644
index 0000000..4a99689
--- /dev/null
+++ b/test/math/binomial3.cpp
@@ -0,0 +1,31 @@
+#include "../util.h"
+#include <math/extendedEuclid.cpp>
+#include <math/multInv.cpp>
+#include <math/binomial3.cpp>
+
+
+constexpr ll mod = 503;
+
+void stress_test() {
+	vector<ll> last = {1};
+	ll queries = 0;
+	for (ll i = 0; i < mod; i++) {
+		for (ll j = 0; j <= i; j++) {
+			ll got = calc_binom(i, j, mod);
+			ll expected = last[j];
+			if (got != expected) cerr << "calc_binom(" << i << ", " << j << "), got: " << got << ", expected: " << expected << FAIL;
+		}
+		queries += sz(last);
+
+		last.push_back(1);
+		for (ll j = i; j > 0; j--) {
+			last[j] = (last[j] + last[j - 1]) % mod;
+		}
+	}
+	cerr << "tested queries: " << queries << endl;
+}
+
+int main() {
+	stress_test();
+}
+