summaryrefslogtreecommitdiff
path: root/math/berlekampMassey.cpp
diff options
context:
space:
mode:
authormzuenni <michi.zuendorf@gmail.com>2022-06-27 17:19:28 +0200
committermzuenni <michi.zuendorf@gmail.com>2022-06-27 17:19:28 +0200
commit5ab8a5088b729a9953b8dff1b2a985dc8fb2098b (patch)
treeed40d6936c0e9eee40ba62751cbf99ecddbaddc2 /math/berlekampMassey.cpp
parentadabbad9c51cf7cd3874bfde8eac1fbcf84fec10 (diff)
updated tcr
Diffstat (limited to 'math/berlekampMassey.cpp')
-rw-r--r--math/berlekampMassey.cpp30
1 files changed, 30 insertions, 0 deletions
diff --git a/math/berlekampMassey.cpp b/math/berlekampMassey.cpp
new file mode 100644
index 0000000..4999254
--- /dev/null
+++ b/math/berlekampMassey.cpp
@@ -0,0 +1,30 @@
+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;
+ B = T;
+ b = d;
+ m = 0;
+ }
+
+ C.resize(L + 1);
+ C.erase(C.begin());
+ for (auto& x : C) x = (mod - x) % mod;
+ return C;
+} \ No newline at end of file