矩阵

最新推荐文章于 2020-08-12 12:14:34 发布

Albert_Bolt

最新推荐文章于 2020-08-12 12:14:34 发布

阅读量266

点赞数

分类专栏： ACM 文章标签： ACM 算法模板矩阵

本文链接：https://blog.csdn.net/Albert_Bolt/article/details/82425508

版权

ACM 同时被 3 个专栏收录

92 篇文章 1 订阅

订阅专栏

算法模板

13 篇文章 0 订阅

订阅专栏

矩阵

3 篇文章 0 订阅

订阅专栏

结构体

struct matrix{
    int a[maxn][maxn];
}m1,m2;

int n,k;//n*n矩阵的k次方

矩阵乘法

//矩阵乘法 
matrix mul(matrix m1,matrix m2)
{
    matrix result;
    memset(result.a,0,sizeof(result.a));
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            for(int k=1;k<=n;k++)
            {
                result.a[i][j]+=m1.a[i][k]*m2.a[k][j];
                result.a[i][j]%=9973;//模几根据题意 
            }
        }
    }
    return result;
}

矩阵加法

//矩阵加法 
matrix add(matrix m1,matrix m2)
{
    matrix result;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            result.a[i][j]=m1.a[i][j]+m2.a[i][j];
            result.a[i][j]%=9973;//模几根据题意
        }
    }
    return result;
}

矩阵快速幂

//矩阵快速幂 
matrix mpow(matrix m1,int k)
{
    if(k==1) return m1;//1次幂为原矩阵 
    matrix result=mpow(m1,k/2);//result=m1^(k/2)
    result=mul(result,result);//result=m1^k
    if(k%2!=0) result=mul(result,m1);//若k为奇数再乘一次m1 
    return result;
}