

此处借鉴大佬模板
struct Mat
{
LL m[101][101];
};//存储结构体
Mat a,e; //a是输入的矩阵,e是输出的矩阵
Mat Mul(Mat x,Mat y)
{
Mat c;
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j){
c.m[i][j] = 0;
}
}
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j){
for(int k=1;k<=n;++k){
c.m[i][j] = c.m[i][j]%mod + x.m[i][k]*y.m[k][j]%mod;
}
}
}
return c;
}
Mat pow(Mat x,LL y)//矩阵快速幂
{
Mat ans = e;
while(y){
if(y&1) ans = Mul(ans,x);
x = Mul(x,x);
y>>=1;
}
return ans;
}
可以参考原文
还有一个因为时间太久,忘记了出处,如有侵权,请及时联系,谢谢!
本文深入探讨了矩阵乘法和矩阵快速幂算法的实现细节,提供了具体的代码示例,适用于算法竞赛和高性能计算场景中矩阵操作的需求。
8091

被折叠的 条评论
为什么被折叠?



