矩阵模板
Posted passione-123456
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了矩阵模板相关的知识，希望对你有一定的参考价值。
struct mat { LL ma[sz][sz]; inline mat() { memset(ma, 0, sizeof ma); } inline mat operator+(const mat & T) const { mat res; for (int i = 0; i < sz; ++i) for (int j = 0; j < sz; ++j) res.ma[i][j] = (ma[i][j] + T.ma[i][j] + MOD) % MOD; return res; } inline mat operator-(const mat& T) const { mat res; for (int i = 0; i < sz; ++i) for (int j = 0; j < sz; ++j) res.ma[i][j] = (ma[i][j] - T.ma[i][j]) % MOD; return res; } inline mat operator*(const mat& T) const { mat res; for (int i = 0; i < sz; ++i) for (int j = 0; j < sz; ++j) for (int k = 0; k < sz; ++k) res.ma[i][j] = (res.ma[i][j] + ma[i][k] * T.ma[k][j]) % MOD; return res; } inline mat operator^(LL x) const { mat res, bas; for (int i = 0; i < sz; ++i) res.ma[i][i] = 1; for (int i = 0; i < sz; ++i) for (int j = 0; j < sz; ++j) bas.ma[i][j] = ma[i][j]; while (x) { if (x & 1) res = res * bas; bas = bas * bas; x >>= 1; } return res; } }
不错呀
以上是关于矩阵模板的主要内容，如果未能解决你的问题，请参考以下文章