【每日一题】洛谷 P1962 (矩阵快速幂与斐波那契数列)

想她一次就背十个单词，当我英语过六级后，我就去告诉她，我很在意她
一天一道数论题，当我可以秒杀数论题的时候，就开始做 DP

---

今日份快乐：洛谷 P1962(矩阵快速幂) 传送门

明天份快乐：洛谷 P1356(鸽舍原理) 传送门

---

题目：

请你求出 F_n mod 10⁹+7 的值。

分析

根据数据范围，暴力走的话只能过 60%，如果想过 100% 的 2⁶³的话，必须有更快的方法去计算
斐波那契数列有这样的性质：

我们对上式进行递推，得

很明显，我们只需要计算出

我们有O(log₂n)的算法去计算 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;
}
 

坚持的时候很狼狈，等成功以后，丑的还是丑的🤭