[模板]矩阵乘法（斐波那契数列）

洛谷P1962

题目背景

大家都知道，斐波那契数列是满足如下性质的一个数列：

• f(1) = 1

• f(2) = 1

• f(n) = f(n-1) + f(n-2) (n ≥ 2 且 n 为整数)

题目描述

请你求出 f(n) mod 1000000007 的值。

输入输出格式

输入格式：

·第 1 行：一个整数 n

输出格式：

第 1 行： f(n) mod 1000000007 的值

输入输出样例

输入样例#1： 

5

输出样例#1： 

5

输入样例#2： 

10

输出样例#2： 

55

说明

对于 60% 的数据： n ≤ 92

对于 100% 的数据： n在long long(INT64)范围内。

矩阵的模板并不难，主要就是如何画出有用的矩阵

另外，还需特别记住一点：矩阵乘法不满足交换律，但满足结合律。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int mod=1e9+7;
struct matrix
{    
    ll m[11][11];
    }a,b,ans;
ll n;

inline matrix mul(matrix a,matrix b)
{
    matrix ans;
    for(int i=1;i<=2;i++)
        for(int j=1;j<=2;j++)
        {
            ans.m[i][j]=0;
            for(int k=1;k<=2;k++)
                ans.m[i][j]=(ans.m[i][j]+a.m[i][k]*b.m[k][j])%mod;
            }
    return ans;
    }
inline matrix matpow(matrix a,ll k)
{
    matrix base=a,ans=a;k--;
    while(k)
    {
        if(k&1) ans=mul(ans,base);
        base=mul(base,base);
        k>>=1;
        }
    return ans;
    }
int main()
{
    scanf("%lld",&n);
    if(n<=2){printf("1
");return 0;}
    a.m[1][1]=1;a.m[1][2]=1;a.m[2][1]=1;a.m[2][2]=0;
    b.m[1][1]=1;b.m[2][1]=1;//前面的是f(2) 后面的是f(1);
    a=matpow(a,n-1);
    ans=mul(a,b);
    cout<<ans.m[2][1]<<endl;//输出的是下面的f(n)  上面的是f(n+1)
    return 0;
    }