矩阵乘法模板代码

CI mod = 1e9 + 7;

struct matrix{ int a[maxm][maxm], n, m;};

matrix matrixMul(matrix p, matrix q){
    matrix res;
    res.n = p.n, res.m = q.m;
    f (i, 0, res.n - 1)
        f (j, 0, res.m - 1){
        	res.a[i][j] = 0;
            f (k, 0, p.m - 1)
                (res.a[i][j] += p.a[i][k] * q.a[k][j] % mod) %= mod;
		}
	return res;
}

matrix ksm(matrix a, int b){
    matrix ans = I;
    for (; b; b >>= 1, a = matrixMul(a, a))
    	if (b & 1)
    		ans = matrixMul(ans, a);
	return ans;
}

作者:DE_aemmprty

出处:https://www.cnblogs.com/aemmprty/p/17477287.html

版权:本作品采用「署名-非商业性使用-相同方式共享 4.0 国际」许可协议进行许可。

posted @   DE_aemmprty  阅读(27)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 微软正式发布.NET 10 Preview 1:开启下一代开发框架新篇章
· 没有源码,如何修改代码逻辑?
· DeepSeek R1 简明指南:架构、训练、本地部署及硬件要求
· NetPad:一个.NET开源、跨平台的C#编辑器
· PowerShell开发游戏 · 打蜜蜂
more_horiz
keyboard_arrow_up dark_mode palette
选择主题
点击右上角即可分享
微信分享提示