矩阵乘

最新推荐文章于 2023-04-17 09:09:03 发布

czysjr

最新推荐文章于 2023-04-17 09:09:03 发布

阅读量389

点赞数

本文链接：https://blog.csdn.net/czysjr/article/details/39503341

版权

这篇博客回顾了矩阵乘法，强调了快速乘法在处理大整数乘法时的重要性。作者通过一道题目的例子说明如何使用快速乘法解决矩阵问题，并提到LC讲解中的一个简化技巧。同时，推荐了一篇关于构造矩阵的好文章，并提供了练习题目链接，包括对矩阵快速幂的运用。文章还讨论了如何处理矩阵运算的二分解法，特别是解决涉及大指数的矩阵求和问题。

摘要由CSDN通过智能技术生成

一年前这个时候学的矩阵乘法，但是当时好像就没怎么懂，前天碰到了这种题果然还是不会做，，，多做几道练习吧。

求K位的相邻两位之差不大于2的题。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#define mod 1000000007
using namespace std;
long long k;
inline void cheng(int a[11][11], int b[11][11]){
    int c[11][11] = {0};
    for(int i = 0; i < 10; i ++)
        for(int j = 0; j < 10; j ++)
            for(int k = 0; k < 10; k ++)
                c[i][j] = ((long long)a[i][k] * b[k][j] + c[i][j]) % mod;
    for(int i = 0; i < 10; i ++)
        for(int j = 0; j < 10; j ++)
            a[i][j] = c[i][j];  
}
int a[11][11], b[11][11], ans;
int main()
{
    freopen("sam.in", "r", stdin);
    freopen("sam.out", "w", stdout);
    scanf("%I64d", &k);
    long long n = k;
    if(n == 1){puts("10"); return 0;}
    for(int i = 0; i < 10; i ++)a[i][i] = 1;
    for(int i = 0; i < 10; i ++){
        for(int j = 0; j < 10; j ++)
            if(abs(j - i) <= 2)
                b[i][j] = 1;    
    }
    n--;
    while(n){
        if(n & 1)cheng(a, b);
        n = n >> 1;
        cheng(b, b);    
    }
    for(int j = 1; j < 10; j ++)
        for(int k = 0; k < 10; k ++)
            ans= (ans + a[k][j]) % mod;
    cout<<ans<<endl;
   // system("pause");
    return 0;
}

然后做了noi2012随机数生成器

这个只是一个简单的两位的矩阵但是两个longlong乘起来会爆，所以还需要用快速乘，以前没写过快速乘，不过还是很好写的，就是和快速幂一样的思想，从乘2到乘4再到乘8这样注意地把一个乘数分解就好了。但是我做这道题的时候还是不能自己构建出矩阵。。继续联系吧

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
using namespace std;
long long m, a, c, x0, n, q;
long long cheng1(long long x, long long y){
    long long ans = 0;
    for(; y; y>>=1, x = (x * 2) % m)
        if(y & 1) ans = (ans + x) % m;
    return ans;
}
inline void che