daylight/math/matrix.hpp

• View this file on GitHub
• Last update: 2024-12-26 20:32:55+09:00
• Include: #include "daylight/math/matrix.hpp"

Depends on

• 基本テンプレート (daylight/base.hpp)

Verified with

• test/aoj/itp/1_7_D.test.cpp
• test/yukicoder/2441.test.cpp

Code

#include "daylight/base.hpp"
template<typename T>
vector<vector<T>> multiMatrix(const vector<vector<T>> &A,
							  const vector<vector<T>> &B) {
	int N = SZ(A);
	int M = SZ(B[0]);
	int L = SZ(A[0]);
	auto ret = make_v<T>(N, M);
	REP(i, N) {
		REP(j, M) {
			REP(k, L) {
				ret[i][j] += A[i][k] * B[k][j];
			}
		}
	}
	return ret;
}

template<typename T>
vector<vector<T>> powMatrix(const vector<vector<T>> &A,
							ll p) {
	auto ans = make_v<T>(SZ(A), SZ(A));
	REP(i, SZ(A)) ans[i][i] = 1;
	auto mul = A;
	for(; p > 0; p >>= 1) {
		if((p & 1) == 1) ans = multiMatrix(ans, mul);
		mul = multiMatrix(mul, mul);
	}
	return ans;
}

Traceback (most recent call last):
  File "/home/runner/.local/lib/python3.10/site-packages/competitive_verifier/oj_resolve/resolver.py", line 181, in resolve
    bundled_code = language.bundle(path, basedir=basedir)
  File "/home/runner/.local/lib/python3.10/site-packages/competitive_verifier/oj/verify/languages/cplusplus.py", line 252, in bundle
    bundler.update(path)
  File "/home/runner/.local/lib/python3.10/site-packages/competitive_verifier/oj/verify/languages/cplusplus_bundle.py", line 482, in update
    self.update(
  File "/home/runner/.local/lib/python3.10/site-packages/competitive_verifier/oj/verify/languages/cplusplus_bundle.py", line 477, in update
    raise BundleErrorAt(
competitive_verifier.oj.verify.languages.cplusplus_bundle.BundleErrorAt: daylight/base.hpp: line 103: unable to process #include in #if / #ifdef / #ifndef other than include guards

Back to top page