矩阵快速幂: 因为结构体已经重载了 * 所以只需在原来普通的快速幂模板的基础上稍微修改一下。
其他操作主要是模拟。
#define _rep(i, a, b) for(int i = a; i <= b; i++)
#define _for(i, a, b) for(int i = a; i < b; i++)
#define M_sel *this
using namespace std;
typedef vector<int> vec_int;
typedef vector<vec_int> mat_int;
typedef long long ll;
const int Maxn = 1010;
const int Maxm = 1010;
const int Mod = 10000;
struct Matrix {
int n, m;
mat_int a;
Matrix(int n = 0, int m= 0) {
this->n = n; this->m = m;
vec_int tmp;
tmp.assign(m, 0);
_for(i, 0, n)
a.push_back(tmp);
}
Matrix operator + (const Matrix &b) const {
Matrix tmp(n, m);
_for(i, 0, n)
_for(j, 0, m)
tmp.a[i][j] = a[i][j] + b.a[i][j];
return tmp;
}
Matrix operator - (const Matrix &b) const {
Matrix tmp(n, m);
_for(i, 0, n)