[蓝桥杯][2014年第五届真题]斐波那契--题解（快速矩阵幂）C++代码

题目描述
斐波那契数列大家都非常熟悉。它的定义是：

f(x) = 1 … (x=1,2)
f(x) = f(x-1) + f(x-2) … (x> 2)

对于给定的整数 n 和 m，我们希望求出：
f(1) + f(2) + … + f(n) 的值。但这个值可能非常大，所以我们把它对 f(m) 取模。
公式如下

但这个数字依然很大，所以需要再对 p 求模。
输入
输入为一行用空格分开的整数 n m p (0 < n, m, p < 10^18)
输出
输出为1个整数，表示答案
样例输入
2 3 5
样例输出
0
解析：
斐波那契求和公式：Sn=f(n+2)-1
本题具体公式：Sn=(f(n+2)-1)%f(m)%p

参考代码

#include<iostream>
#include<cstring>
using namespace std;
typedef long long LL;

//矩阵结构体
typedef struct m{
   
    LL a[2][2];
}matrix;
LL n,m,p;
LL fib_m;
//大数相乘取模
LL mul_mod(LL a,LL b,LL p){
   
    LL ans=0;
    while(b){
   
        if(b&1){
   
            ans+=a;
            ans%=p;
        }
        a<<=1;
        a%=p;
        b>>=1;
    }
    return ans;
}
//矩阵乘法
matrix mul_mat_1(matrix first,matrix second){
   
    matrix result;
    memset(result.a,0,sizeof(result.a));
    for(int i=0;i<2;<