#include <cstdio>
int n, m;
struct Matrix {
int a[3][3];
Matrix operator*(const Matrix &b) {
Matrix res{};
for (int i = 1; i <= 2; i ++)
for (int j = 1; j <= 2; j ++)
for (int k = 1; k <= 2; k ++)
res.a[i][j] = (res.a[i][j] + 1LL * a[i][k] * b.a[k][j] % m) % m;
return res;
}
} ans, base;
int main() {
scanf("%d %d", &n, &m);
base.a[1][1] = base.a[1][2] = base.a[2][1] = 1, base.a[2][2] = 0;
ans.a[1][1] = ans.a[1][2] = 1;
if (n == 1)
return puts("1"), 0;
else if (n == 2)
return puts("2"), 0;
while (n) {
if (n & 1)
ans = ans * base;
base = base * base;
n >>= 1;
}
printf("%d", ans.a[1][1] - 1);
return 0;
}
一本通1643:【例 3】Fibonacci 前 n 项和
最新推荐文章于 2024-09-14 22:25:23 发布
文章详细介绍了如何使用C语言实现矩阵乘法的快速幂算法,通过递归和位操作优化计算过程,适用于处理较大规模的矩阵乘法问题。
摘要由CSDN通过智能技术生成