矩阵乘法和矩阵快速幂

矩阵乘法很久以前就学了

久到忘记了是什么时候

然而清楚地记得我不会

在集训的时候被废鼠姐姐一句话点醒(废鼠姐姐tql orz)

手动分割---------------------------------------------------------------------------------------------------------

设有矩阵A和B,A为p*n的矩阵,B为n*q的矩阵,则他们的乘积C是一个p*q的矩阵

C_{i,j}=\sum_{k=1}^{n}(A_{i,k}*B_{k,j})

如果还不懂,看图就好了

至于代码

先咕咕了(逃)

struct matrix{
	int r,c,m[N][N];
	matrix(){memset(m,0,sizeof m);}
	inline matrix operator * (const matrix& b) const {
		matrix ans;
		ans.r=r;
		ans.c=b.c;
		for(int i=1;i<=ans.r;++i)
			for(int j=1;j<=ans.c;++j)
				for(int k=1;k<=c;++k)
					ans.m[i][j]=(ans.m[i][j]+m[i][k]*b.m[k][j])%mod;
		return ans;
	}
};

猜你喜欢

转载自blog.csdn.net/Sakura_Chiyo/article/details/82049900