summaryrefslogtreecommitdiff
path: root/content/math/binomial0.cpp
diff options
context:
space:
mode:
authorGloria Mundi <gloria@gloria-mundi.eu>2024-11-16 01:24:14 +0100
committerGloria Mundi <gloria@gloria-mundi.eu>2024-11-16 01:24:14 +0100
commit98567ec798aa8ca2cfbcb85c774dd470f30e30d4 (patch)
tree5113d5cc24d1ad5f93810b6442ce584a36950dc8 /content/math/binomial0.cpp
parentad3856a6b766087df0036de0b556f4700a6498c9 (diff)
parent8d11c6c8213f46f0fa19826917c255edd5d43cb1 (diff)
mzuenni tests
Diffstat (limited to 'content/math/binomial0.cpp')
-rw-r--r--content/math/binomial0.cpp14
1 files changed, 14 insertions, 0 deletions
diff --git a/content/math/binomial0.cpp b/content/math/binomial0.cpp
new file mode 100644
index 0000000..f37aea5
--- /dev/null
+++ b/content/math/binomial0.cpp
@@ -0,0 +1,14 @@
+constexpr ll lim = 10'000'000;
+ll fac[lim], inv[lim];
+
+void precalc() {
+ fac[0] = inv[0] = 1;
+ for (int i = 1; i < lim; i++) fac[i] = fac[i-1] * i % mod;
+ inv[lim - 1] = multInv(fac[lim - 1], mod);
+ for (int i = lim - 1; i > 0; i--) inv[i-1] = inv[i] * i % mod;
+}
+
+ll calc_binom(ll n, ll k) {
+ if (n < 0 || n < k || k < 0) return 0;
+ return (fac[n] * inv[n-k] % mod) * inv[k] % mod;
+}