目标：矩阵快速幂 掌握?
需要技能：矩阵乘法，c语言基础
要求能力：和我一样弱都能学会。
矩阵乘法大概就是：左行*右列（左边的矩阵第一行×右边第一列然后求和 就是 乘法的第一行第一个元素）
?
左边是行 右边是列 横是行 竖着是列
2行2列 *?2行2列 =?2行2列
?红框看为第一行和第一列 代表 答案的 第一行第一个元素
?绿框看为第一行和第二列 代表 答案的 第一行第二个元素
?蓝框看为第二行和第一列 代表 答案的 第二行第一个元素
?紫框看为第二行和第二列 代表 答案的 第二行第二个元素
#include<bits/stdc++.h> //万能头
using namespace std;//namespace
typedef long long ll;//定义ll 为long long 方便使用
#define N 111//最大数
const int M = 1e9+7;//Mod 需要模这个数
#define mod(x) ((x)%M )//定义 mod的函数
int n ;
struct mat{
ll m[N][N];
}unit;//矩阵结构体 以及单位矩阵unit (主对角为1其余为0)
mat operator * (mat a,mat b){ //重载operator 乘法
mat ret ;
ll val,i,j;
for(i = 0 ; i < n ; i++)
{
for(j= 0 ; j < n ; j++){
val = 0;
for(int k = 0 ; k < n;k++){
val+= mod((ll)a.m[i][k]*b.m[k][j]); //三个循环是求矩阵乘法 mod一下
}
ret.m[i][j] = mod(val);//也需要mod 一下 都是为了防止数太大
}
}
return ret;//返回 俩个矩阵乘法的结果
}
void init(){
for(int i = 0 ; i < N ;i++){
unit.m[i][i] = 1;
} //init 单位矩阵
}
mat quick_pow(mat a, ll n){//快速矩阵幂
mat res = unit ;
while(n){
if(n&1) res=res*a;
a = a*a;
n>>=1;
}
return res;
}
int main(){
cin >> n ;
init();
mat a,b;
cout<<"第一个"<<endl;
for(int i = 0 ; i < n ; i++){
for(int j = 0 ; j <n ;j++){
cin>>a.m[i][j];
}
}
//写入 矩阵
mat ans = quick_pow(a,2); //调用快速矩阵 2是幂
for(int i = 0 ; i < n ; i++){
for(int j = 0 ; j <n ;j++){
cout<<ans.m[i][j]<<" ";
}
cout<<endl;
}//输出结果
return 0;
}如果quick_pow看不懂，看看int如何快速幂，博客有介绍。