summaryrefslogtreecommitdiff
path: root/math/berlekampMassey.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 /math/berlekampMassey.cpp
parentad3856a6b766087df0036de0b556f4700a6498c9 (diff)
parent8d11c6c8213f46f0fa19826917c255edd5d43cb1 (diff)
mzuenni tests
Diffstat (limited to 'math/berlekampMassey.cpp')
-rw-r--r--math/berlekampMassey.cpp31
1 files changed, 0 insertions, 31 deletions
diff --git a/math/berlekampMassey.cpp b/math/berlekampMassey.cpp
deleted file mode 100644
index 29e084f..0000000
--- a/math/berlekampMassey.cpp
+++ /dev/null
@@ -1,31 +0,0 @@
-constexpr ll mod = 1'000'000'007;
-vector<ll> BerlekampMassey(const vector<ll>& s) {
- int n = sz(s), L = 0, m = 0;
- vector<ll> C(n), B(n), T;
- C[0] = B[0] = 1;
-
- ll b = 1;
- for (int i = 0; i < n; i++) {
- m++;
- ll d = s[i] % mod;
- for (int j = 1; j <= L; j++) {
- d = (d + C[j] * s[i - j]) % mod;
- }
- if (!d) continue;
- T = C;
- ll coef = d * powMod(b, mod-2, mod) % mod;
- for (int j = m; j < n; j++) {
- C[j] = (C[j] - coef * B[j - m]) % mod;
- }
- if (2 * L > i) continue;
- L = i + 1 - L;
- swap(B, T);
- b = d;
- m = 0;
- }
-
- C.resize(L + 1);
- C.erase(C.begin());
- for (auto& x : C) x = (mod - x) % mod;
- return C;
-}