矩阵运算，矩阵快速幂，模板

以你之姓@ 2022-06-11 05:41 221阅读 0赞

#include<cstdio>
    #include<algorithm>
    #include<cstring>
    #include<iostream>
    using namespace std;
    
    #define MOD 1000000007
    struct Mat{
        int n,m;   //行列
        int mat[19][19];
    };
    Mat operator *(Mat a, Mat b) {  //矩阵相乘  
    	Mat c;
    	memset(c.mat, 0, sizeof(c.mat));
    	c.n = a.n, c.m = b.m;
    	for (int i = 1; i <= a.n; i++) {
    		for (int j = 1; j <= b.m; j++) {
    			for (int k = 1; k <= a.m; k++) {
    				c.mat[i][j] += (a.mat[i][k] * b.mat[k][j]) % MOD;
    				c.mat[i][j] = (c.mat[i][j] % MOD + MOD) % MOD;
    			}
    		}
    	}
    	return c;
    }
    Mat operator +(Mat a, Mat b) { //矩阵相加  
    	Mat c;
    	memset(c.mat, 0, sizeof(c.mat));
    	c.n = a.n, c.m = a.m;
    	for (int i = 1; i <= a.n; i++) {
    		for (int j = 1; j <= a.m; j++) {
    			c.mat[i][j] = (a.mat[i][j] + b.mat[i][j]) % MOD;
    		}
    	}
    	return c;
    }
    Mat operator ^(Mat a,int k){   //矩阵幂
        Mat c;
        memset(c.mat,0,sizeof(c.mat));
        c.n = a.n,c.m = a.n;
        for(int i=1;i<=a.n;i++)c.mat[i][i] = 1;
        while(k){
            if(k&1){
                c = c*a;
            }
            a = a*a;
            k>>=1;
        }
        return c;
    }
    void out(Mat a){      //输出矩阵
        for(int i=1;i<=a.n;i++){
            for(int j=1;j<=a.m;j++){
                printf(j==a.m? "%d\n":"%d ",a.mat[i][j]);
            }
        }
    }
    int main(){
        Mat d;
        d.n=3,d.m=3;
        for(int i=1;i<=d.n;i++){
            for(int j=1;j<=d.m;j++){
                scanf("%d",&d.mat[i][j]);
            }
        }
        d=d^2;
        out(d);
    }

发表评论取消回复

表情：

评论列表（有 0 条评论，221人围观）

还没有评论，来说两句吧...

相关阅读

相关矩阵快速幂的模板（C++）

顾名思义，快速幂就是快速计算某个数的n次幂。其时间复杂度为 O(log₂N)，与朴素的O(N)相比效率有了极大的提高。前面学习了整数的快速幂：比如我们要求12^

迷南。/ 2023年03月06日 03:50/ 0 赞/ 2 阅读

相关矩阵快速幂（快速幂）模板题目Decoding Genome

这类题目的难点在于，数据量大，需要用到高精度，也就是快速幂矩阵的乘法与快速幂模板一： const int mod = 1e9+7; struct Matr

今天药忘吃喽~/ 2022年10月14日 13:51/ 0 赞/ 169 阅读

相关矩阵快速幂

昨天晚上矩阵小王子给我们讲了一下矩阵快速幂，学习了一下，写了一个模板。 1：思想矩阵快速幂的思想就是跟数的快速幂一样，假如我们要求2^11，次方，我们可以把 11 写成

Love The Way You Lie/ 2022年08月12日 18:30/ 0 赞/ 216 阅读

相关矩阵快速幂模板

int mod; const int NUM;//定义矩阵能表示的最大维数 int N;//N表示矩阵的维数，以下的矩阵加法、乘法、快速幂都是按N维矩阵

曾经终败给现在/ 2022年07月03日 15:26/ 0 赞/ 202 阅读

相关矩阵快速幂

include<cstdio> include<iostream> include<cmath> include<vector> inc

朴灿烈づ我的快乐病毒、/ 2022年06月18日 12:29/ 0 赞/ 252 阅读

相关矩阵快速幂基础模板

Think: 1注意记录数组初始化 2思想->大数取模 include <bits/stdc++.h> using namespace st

绝地灬酷狼/ 2022年06月16日 00:54/ 0 赞/ 227 阅读

相关矩阵运算，矩阵快速幂，模板

include<cstdio> include<algorithm> include<cstring> include<iostream>

以你之姓@/ 2022年06月11日 05:41/ 0 赞/ 222 阅读

相关快速幂，矩阵快速幂（模板）

1，整数快速幂： C++ Code <table style="width:100%;font-family:Consolas, 'Courier New';font-

左手的ㄟ右手/ 2022年06月11日 03:57/ 0 赞/ 297 阅读

相关矩阵快速幂

通常我们使用的快速幂是以二为底的，这次就遇到了一道以10为底的快速幂题目；先说下快速幂 long long power(long long a,long long

桃扇骨/ 2021年11月10日 14:31/ 0 赞/ 352 阅读

相关矩阵快速幂模板

const int maxn = 2; struct Matrix{ ll a[maxn][maxn]; void

小鱼儿/ 2021年09月19日 05:44/ 0 赞/ 384 阅读