题目
分析
看似简单,但一看数据范围就知道不容易。
先看斐波那契数列的递推公式:
据此推出矩阵乘法:
进一步简化:
推出来后直接计算即可。
注意:矩阵乘法不满足交换律,不能搞反掉顺序!!!
代码
#include<iostream>
using namespace std;
typedef long long ll;
// 存矩阵
struct mat {
ll m[5][5];
};
ll n, m;
mat f, res;
// 矩阵乘法
void multiple(mat &x, mat &y, mat &z) {
memset(z.m, 0, sizeof z.m);
for (int i = 1; i <= 2; i++)
for (int j = 1; j <= 2; j++)
for (int k = 1; k <= 2; k++) {
z.m[i][k] += x.m[i][j] * y.m[j][k];
if (z.m[i][k] >= m) z.m[i][k] %= m;
}
}
void mpow(ll k) {
mat temp = f, t;
// 单位矩阵
res.m[1][1] = res.m[2][2] = 1;
while (k) {
if (k & 1) {
multiple(res, temp, t);
res = t;
}
multiple(temp, temp, t);
temp = t;
k >>= 1;
}
}
int main(){
cin >> n >> m;
f.m[1][1] = f.m[1][2] = f.m[2][1] = 1;
f.m[2][2] = 0;
if (n <= 2) cout << 1 << endl;
else {
mpow(n - 2);
ll ret = res.m[1][1] + res.m[1][2]; // 转移矩阵第一行都为1,直接加
if (ret >= m) ret %= m;
cout << ret << endl;
}
return 0;
}