wu


私信TA

用户名:cncfvc

访问量:215684

签 名:

读研狗没有时间刷题了~~

等  级
排  名 2
经  验 36073
参赛次数 8
文章发表 265
年  龄 25
在职情况 学生
学  校 电子科技大学
专  业 通信工程

  自我简介:

写代码 真好玩 ~

有点难理解。。。

参考代码:

#include <iostream>  
#include <cstdio>  
using namespace std;  
const int N=2;  
int m;  
struct Matrix{  
    int m[N][N];  
}matrix;  //定义表示矩阵的结构体
Matrix I={1,0,  
          0,1};  //单位矩阵
Matrix p;  
inline Matrix  multiply(Matrix a,Matrix b){  //inline表示内联函数  实现矩阵乘法
      int i,j,k;  
      Matrix re;  
      for(i=0;i<N;i++)  
      for(j=0;j<N;j++){  
         re.m[i][j]=0;  
         for(k=0;k<N;k++)  
           re.m[i][j]+=a.m[i][k]*b.m[k][j];  
           re.m[i][j]%=m;  
      }  
      return re;  
}  
inline Matrix quick_pow(int n){  
    Matrix re=p,b=I;  
    while(n>0){  
        if(n&1){  //当是奇数的时候
            b=multiply(b,re);  
        }  
        n=n>>1;  //向右移一位
        re=multiply(re,re);  
    }  
    return b;  
}  
  
int main()  
{  
   int b;  
   while(cin>>b>>m){  
       for(int i=0;i<2;i++)  
         for(int j=0;j<2;j++)  
           cin>>p.m[i][j];  
             
        matrix=quick_pow(b);  
          
       for(int i=0;i<2;i++){  
        for(int j=0;j<2;j++)  
         cout<<matrix.m[i][j]%m<<" ";    //坑  如果m=1且输出单位矩阵就要mod成0   
         cout<<endl;  
     }   
                       
   }  
    return 0;  
}


 

0.0分

1 人评分

  评论区