首页 > 代码库 > 矩阵乘法
矩阵乘法
#include<bits/stdc++.h>
using namespace std;
const int N=2;
const int MOD=10000;
struct MAT{int a[N][N];};
MAT operator*(MAT x, MAT y){
MAT z;
for(int i=0;i<N;++i)
for(int j=0;j<N;++j){
z.a[i][j]=0;
for(int k=0;k<N;++k)
z.a[i][j]=(z.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;
}
return z;
}
MAT pow(int n){
MAT ret,s;
ret.a[0][0]=ret.a[1][1]=1;
ret.a[1][0]=ret.a[0][1]=0;
s.a[0][0]=s.a[0][1]=s.a[1][0]=1;
s.a[1][1]=0;
for(;n;n>>=1,s=s*s)if(n&1)ret=ret*s;
return ret;
}
int main(){
int n;
while(cin>>n,n>=0){
MAT ans=pow(n);
cout<<ans.a[0][1]<<endl;
}
}
矩阵乘法
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。