矩阵快速幂_求一个邻接矩阵中走m步的方法数

最新推荐文章于 2024-08-22 15:56:58 发布

Z_X_B_233

最新推荐文章于 2024-08-22 15:56:58 发布

阅读量1k

点赞数

分类专栏：训练题文章标签：算法 matrix

本文链接：https://blog.csdn.net/cbcbcbz/article/details/70147507

版权

训练题专栏收录该内容

4 篇文章 0 订阅

订阅专栏

matrix multi(matrix a, matrix b)
{
    matrix tmp;
    for(int i = 0; i < 55; ++i) //55是base的范围
    {
        for(int j = 0; j < 55; ++j)
        {
            for(int k = 0; k < 55; ++k)
                tmp.m[i][j] = (tmp.m[i][j] + a.m[i][k] * b.m[k][j]) % Mod;
        }
    }
    return tmp;
}

matrix fast_mod(matrix base, int m,int n)  // 求矩阵 base 的m次幂 n为base的范围
{
    matrix ans;
    for(int i = 1; i <= n; i++)
        ans.m[1][i] = 1;
    while(m)
    {
        if(m & 1)  //实现 ans *= t; 其中要先把 ans赋值给 tmp，然后用 ans = tmp * t
            ans = multi(ans, base);
        base = multi(base, base);
        m >>= 1;
    }
    return ans;
}

ans[1][n+1]即为答案