summaryrefslogtreecommitdiff
path: root/content/math/matrixPower.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/matrixPower.cpp
parentad3856a6b766087df0036de0b556f4700a6498c9 (diff)
parent8d11c6c8213f46f0fa19826917c255edd5d43cb1 (diff)
mzuenni tests
Diffstat (limited to 'content/math/matrixPower.cpp')
-rw-r--r--content/math/matrixPower.cpp14
1 files changed, 14 insertions, 0 deletions
diff --git a/content/math/matrixPower.cpp b/content/math/matrixPower.cpp
new file mode 100644
index 0000000..d981e6e
--- /dev/null
+++ b/content/math/matrixPower.cpp
@@ -0,0 +1,14 @@
+vector<mat> pows;
+
+void precalc(mat m) {
+ pows = {mat(sz(m.m), 1), m};
+ for (int i = 1; i < 60; i++) pows.push_back(pows[i] * pows[i]);
+}
+
+auto calc(ll b, vector<ll> v) {
+ for (ll i = 1; b > 0; i++) {
+ if (b & 1) v = pows[i] * v;
+ b /= 2;
+ }
+ return v;
+}