矩阵乘法的经典应用

入坑信奥的L同学

已于 2023-02-27 23:28:32 修改

阅读量603

点赞数

分类专栏：笔记算法进阶指南文章标签：矩阵乘法

于 2022-10-10 21:48:42 首次发布

本文链接：https://blog.csdn.net/m0_55341679/article/details/126695948

版权

~~矩乘是真的难qwq……~~
先简单介绍一下递推：递推算法是一种用若干步可重复运算来描述复杂问题的方法。递推是序列计算中的一种常用算法。通常是通过计算前面的一些项来得出序列中的指定项的值。比如斐波那契数列以及汉诺塔问题。

下面是正文部分：

例题1：AcWing 205.斐波那契
矩阵乘法大致要了解其运算规则，比如满足结合律、分配律之类的，此处可见蓝书P156。矩阵乘法用在这题上可以加速递推。由于矩阵乘法满足结合律，所以矩阵乘法可以用快速幂加速。具体算法流程可见蓝书P157。代码如下:

#include<iostream>
#include<cstring>
using namespace std;
const int mod=1e4;
void mul(int f[2],int a[2][2]){
   
    int c[2]={
   0,0};
    for(int j=0;j<2;j++){
   
        for(int k=0;k<2;k++){
   
            c[j]=(c[j]+(long long)f[k]*a[k][j])%mod;
        }
    }
    memcpy(f,c,sizeof c);
}
void mulself(int a[2][2]){
   
    int c[2][2]={
   0,0};
    for(int i=0;i<2;i++){
   
        for(int j=0;j<2;j++){
   
            for(int k=0;k<2;k++){
   
                c[i][j]=(c[i][j]+(long long)a[i][k]*a[k][j])%mod;
            }
        }
    }
    memcpy(a,c,sizeof c);
}
int main(){
   
    int n;
    while(cin>>n && n!=-1){
   
        int f[2]={
   0,1};
        int a[2][2]={
   {
   0,1},{
   1,1}};
        while(n){
   
            if(n&1) mul(f,a);
            n>>=1;
            mulself(a);
        }
        cout<<f[0]<<endl;
    }
}

例题2：luoguP7318「PMOI-4」人赢の梦
这题通过找规律可以发现最后的结果可以整理成一个仅与a，b有关的单项式。此外还可以发现a，b的次数与斐波那契数列存在这一些联系（容易找到），这样这个题目就转化为上面的题目了。
例题3：AcWing 225.矩阵幂求和
这题考察矩阵加法与乘法运算以及分治法。由于矩阵不满足可除性，所以不能用等比数列求和法，只能考虑采用分治法。我们由此联想到一个例题：
那么同样地，矩阵的“等比”数列也可以分治求解，但是略有不同。考虑将奇数减去1后转化为偶数的情况：

#include<iostream>
#include<cstring>
using namespace std;
const int N=35;
int a[N][N],self[N][N];
int n,m,mod;
void add(int b[][N],int c[][N]){
   
    for(int i=1;i<=n;i++){
   
        for(int j