矩阵快速幂模板

最新推荐文章于 2024-10-25 17:33:50 发布

原创最新推荐文章于 2024-10-25 17:33:50 发布 · 129 阅读

0 ·

CC 4.0 BY-SA版权

ACM 专栏收录该内容

57 篇文章

订阅专栏

本文介绍了一种矩阵快速幂算法的实现方法，通过定义矩阵结构体并实现初始化、矩阵乘法和快速幂运算等功能，适用于解决特定类型的大规模计算问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

矩阵快速幂模板

const int maxarr=7;//矩阵大小

LL aa,bb;
struct node
{
    LL arr[maxarr][maxarr];
};
void init(node &a)//所有新建都要手动清空
{
    for(int i=0;i<maxarr;i++)
    {
        for(int j=0;j<maxarr;j++)
        {
            a.arr[i][j]=0;
        }
    }
}
node MatrixMultiplication(node &a,node &b)
{
    LL i,j,k;
    node c;
    init(c);
    for(i=0;i<maxarr;i++)
    {
        for(j=0;j<maxarr;j++)
        {
            for(k=0;k<maxarr;k++)
            {
                c.arr[i][j]=(c.arr[i][j]+a.arr[i][k]*b.arr[k][j])%mod;
            }
        }
    }
    return c;
}
node fast_power(node a,LL n)
{
    node b;
    init(b);
    for(int i=0;i<maxarr;i++)
    {
        b.arr[i][i]=1;
    }
    while(n>0)
    {
        if(n&1==1) b=MatrixMultiplication(a,b);
        a=MatrixMultiplication(a,a);
        n>>=1;
    }
    return b;
}

int main()
{
    node a,b;
    init(a);
    b=fast_power(a,n-2);
    return 0;
}