【数学专题】矩阵乘法

整理的算法模板合集: ACM模板


在这里插入图片描述

luogu P3390 【模板】矩阵快速幂

在这里插入图片描述

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

using namespace std;
typedef long long ll;
const int N = 107, mod = 1e9 + 7;

int n, m;
int a[N][N];
int ans[N][N];
ll k;

void mul(int c[][N], int a[][N], int b[][N])
{
    int tmp[N][N] = {0};
    for(int i = 1; i <= n; ++ i){
        for(int j = 1; j <= n; ++ j){
            for(int k = 1; k <= n; ++ k){
                tmp[i][j] = (tmp[i][j] + (ll)a[i][k] * b[k][j]) % mod;
            }
        }
    }
    memcpy(c, tmp, sizeof tmp);
}

int main(){
    scanf("%d%lld", &n, &k);
    for(int i = 1; i <= n; ++ i)
        for(int j = 1; j <= n; ++ j)
            scanf("%d", &a[i][j]);

    for(int i = 1; i <= n; ++ i)
        ans[i][i] = 1;

    while(k){
        if(k & 1) mul(ans, ans, a);
        mul(a, a, a);
        k >>= 1;
    }
    for(int i = 1; i <= n;puts(""), ++ i)
        for(int j = 1; j <= n; ++ j)
            printf("%d ", ans[i][j]);
    return 0;
}

AcWing 1303. 斐波那契前 n 项和

在这里插入图片描述

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

using namespace std;
typedef long long ll;
const int N = 3;

int n, m;
int a[N][N];
int f[N];

void mul(int c[], int a[], int b[][N])
{
    int tmp[N] = {0};
    for(int j = 0; j < N; ++ j){//2的行
        for(int k = 0; k < N; ++ k){//1的列2的行,1的列=2的行
            //只有一行,求这一行的每一列
            tmp[j] = (tmp[j] + (ll)a[k] * b[k][j]) % m;
        }
    }
    memcpy(c, tmp, sizeof tmp);
}

void mul(int c[][N], int a[][N], int b[][N])
{
    int tmp[N][N] = {0};
    for(int i = 0; i < N; ++ i){
        for(int j = 0; j < N; ++ j){
            for(int k = 0; k < N; ++ k){
                tmp[i][j] = (tmp[i][j] + (ll)a[i][k] * b[k][j]) % m;
            }
        }
    }
    memcpy(c, tmp, sizeof tmp);
}

int main()
{
    scanf("%d%d", &n, &m);
    
    int f1[N] = {1, 1, 1};
    int a[N][N] = {
        0, 1, 0,
        1, 1, 1,
        0, 0, 1
    };
    n -- ;//这里是从第一项开始的,所以要--;
    
    while(n){
        if(n & 1) mul(f1, f1, a);
        mul(a, a, a);
        n >>= 1;
    }
    
    printf("%d\n", f1[2]);
    return 0;
}

AcWing 1304. 佳佳的斐波那契

在这里插入图片描述

在这里插入图片描述
在这里插入图片描述
计算公式为: T n = n ∗ S n − P n T_n = n * S_n - P_n Tn=nSnPn

//#pragma GCC optimize(2)
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;
typedef long long ll;
const int N = 4;

int n, m;
int f[N][N];

void mul(int c[][N], int a[][N], int b[][N])
{
    static int tmp[N][N];
    memset(tmp, 0, sizeof tmp);
    for(int i = 0; i < N; ++ i){
        for(int j = 0; j < N; ++ j){
            for(int k = 0; k < N; ++ k){
                tmp[i][j] = (tmp[i][j] + (ll)a[i][k] * b[k][j]) % m;
            }
        }
    }
    memcpy(c, tmp, sizeof tmp);
}

int main()
{
    scanf("%d%d", &n, &m);
    // {fn, fn+1, sn, pn}
    // pn = n * sn - tn
    int f1[N][N] = {1, 1, 1, 0};
    int a[N][N] = {
        {0, 1, 0, 0},
        {1, 1, 1, 0},
        {0, 0, 1, 1},
        {0, 0, 0, 1},
    };

    int k = n - 1;

    while(k){
        if(k & 1)mul(f1, f1, a);
        mul(a, a, a);
        k >>= 1;
    }
    cout << (((ll)n * f1[0][2] - f1[0][3] % m + m) % m) << endl;
    return 0;
}

AcWing 1305. GT考试

已标记关键词 清除标记
课程简介: 历经半个多月的时间,Debug亲自撸的 “企业员工角色权限管理平台” 终于完成了。正如字面意思,本课程讲解的是一个真正意义上的、企业级的项目实战,主要介绍了企业级应用系统中后端应用权限的管理,其中主要涵盖了六大核心业务模块、十几张数据库表。 其中的核心业务模块主要包括用户模块、部门模块、岗位模块、角色模块、菜单模块和系统日志模块;与此同时,Debug还亲自撸了额外的附属模块,包括字典管理模块、商品分类模块以及考勤管理模块等等,主要是为了更好地巩固相应的技术栈以及企业应用系统业务模块的开发流程! 核心技术栈列表: 值得介绍的是,本课程在技术栈层面涵盖了前端和后端的大部分常用技术,包括Spring Boot、Spring MVC、Mybatis、Mybatis-Plus、Shiro(身份认证与资源授权跟会话等等)、Spring AOP、防止XSS攻击、防止SQL注入攻击、过滤器Filter、验证码Kaptcha、热部署插件Devtools、POI、Vue、LayUI、ElementUI、JQuery、HTML、Bootstrap、Freemarker、一键打包部署运行工具Wagon等等,如下图所示: 课程内容与收益: 总的来说,本课程是一门具有很强实践性质的“项目实战”课程,即“企业应用员工角色权限管理平台”,主要介绍了当前企业级应用系统中员工、部门、岗位、角色、权限、菜单以及其他实体模块的管理;其中,还重点讲解了如何基于Shiro的资源授权实现员工-角色-操作权限、员工-角色-数据权限的管理;在课程的最后,还介绍了如何实现一键打包上传部署运行项目等等。如下图所示为本权限管理平台的数据库设计图: 以下为项目整体的运行效果截图: 值得一提的是,在本课程中,Debug也向各位小伙伴介绍了如何在企业级应用系统业务模块的开发中,前端到后端再到数据库,最后再到服务器的上线部署运行等流程,如下图所示:
©️2020 CSDN 皮肤主题: 撸撸猫 设计师:设计师小姐姐 返回首页