矩阵快速幂模板

最新推荐文章于 2022-08-27 21:03:02 发布

outer_form

最新推荐文章于 2022-08-27 21:03:02 发布

阅读量340

点赞数

分类专栏：矩阵文章标签： c++ 矩阵

本文链接：https://blog.csdn.net/outer_form/article/details/50557642

版权

矩阵专栏收录该内容

11 篇文章 0 订阅

订阅专栏

void mul(matrix a,matrix b,matrix c){
    matrix d={{0}};
    int i,j,l;
    for(i=0;i<=k;i++)
        for(j=0;j<=k;j++)
            for(l=0;l<=k;l++)
                d[i][j]=(d[i][j]+a[i][l]*b[l][j])%p;
    memcpy(c,d,sizeof d);
}
void quick_pow(matrix a,LL b,matrix c){
    matrix d={{0}},t;
    for(int i=0;i<=k;i++)
        d[i][i]=1;
    memcpy(t,a,sizeof t);
    while(b){
        if(b&1)
            mul(d,t,d);
        mul(t,t,t);
        b>>=1;
    }
    memcpy(c,d,sizeof d);
}