想她一次就背十个单词,当我英语过六级后,我就去告诉她,我很在意她
一天一道数论题,当我可以秒杀数论题的时候,就开始做 DP
今日份快乐:洛谷 P1962(矩阵快速幂) 传送门
明天份快乐:洛谷 P1356(鸽舍原理) 传送门
题目:
请你求出 Fn mod 109+7 的值。
分析
根据数据范围,暴力走的话只能过 60%,如果想过 100% 的 263的话,必须有更快的方法去计算
斐波那契数列有这样的性质:
我们对上式进行递推,得
很明显,我们只需要计算出
我们有O(log2n)的算法去计算 x n-2 的结果,我们可以定义一个2*2的矩阵,对乘法进行重载,再跑一边普通的快速幂就可以
代码
#include <bits/stdc++.h>
using namespace std;
typedef struct name{
ll a[2][2];
name(){ // 单位矩阵
a[0][0] = 1;
a[0][1] = 0;
a[1][0] = 0;
a[1][1] = 1;
}
} node;
// 重载乘法
name operator * (const name &A, const name &B){
name c;
chushi(c.a, 0);
for(int i = 0; i < 2; i++){
for(int j = 0; j < 2; j++){
for(int k = 0; k < 2; k++){
c.a[i][j] += A.a[i][k] * B.a[k][j];
c.a[i][j] %= mod;
}
}
}
return c;
}
int main(){
ll n;
cin >> n;
if(n == 1) cout << 1 << endl; // 特殊处理 n = 1
else{
n--;
node a, b;
b.a[0][0] = b.a[0][1] = b.a[1][0] = 1;
b.a[1][1] =0;
while(n){
if(n & 1) a = a * b;
b = b * b;
n >>= 1;
}
b.a[0][0] = b.a[1][0] = 1;
b.a[0][1] = b.a[1][1] = 0;
a = a * b;
cout << (a.a[1][0]) % mod << endl;
}
return 0;
}
坚持的时候很狼狈,等成功以后,丑的还是丑的🤭