问题 D: LiMn2O4的数学之路（矩阵快速幂求斐波那契数列）

题目描述

长期以来，LiMn2O4常常会因为自己成为不了数学选手而苦恼，并曾经有过一次转型的尝试。那是一个月黑风高的十一长假，LiMn2O4拿起了厚厚的《具体数学》。路过的practer看见了，practer说你有没有办法求一下这个公式的值：

LiMn2O4看了眼公式，说这个值有无理数，怎么表示？可以做？于是practer给了LiMn2O4第另一个公式：

求第一个公式减第二个公式的值。LiMn2O4稍加思考后觉得这个太简单了，他还要看《具体数学》，没有时间。现在交给聪明的你来解决这个问题，请你求出新的公式的第N项。

输入

第一行有一个正整数K。
接下来有K行，每行有一个数N。
数据保证1≤K≤100 00,0≤N≤1 000 000 000。结果对1 000 000 007取模。

输出

输出k行，第i行代表第i个数据的的结果

样例输入
复制样例数据 2
1
1024

样例输出
1
754854590

利用矩阵乘法求斐波那契数列：

斐波那契数列的数列形式递推公式转变为矩阵形式的递推公式

f[i] = 1 * f[i-1]+1 * f[i-2] ；f[i-1] = 1 * f[i-1] + 0 * f[i-2];

注意：因为复杂度是（logn）数据范围可以很大

#include <iostream>
#include <string>
#include <cstdio>
#include<algorithm>
using namespace std;
const int mod=1000000007;

struct mat///矩阵mat
{
    long long x[3][3];
    mat()///构造函数，初始化x[][]
    {
        for(int i=0;i<3;i++) for(int j=0;j<3;j++)
            x[i][j]=0;
    }
};

mat muti(mat a,mat b)///矩阵乘法
{
    mat c;
    for(int i=1;i<=2;i++)
    {
        for(int j=1;j<=2;j++)
        {
            for(int k=1;k<=2;k++)///矩阵a的第i行的第1到k个元素和矩阵b的第j列的第1到k个元素相乘，再相加
            {
                c.x[i][j]=(c.x[i][j]+a.x[i][k]*b.x[j][k])%mod;
            }
        }
    }
    return c;
}

mat qmod(mat &a,int n)///矩阵快速幂
{
    mat c;
    c.x[1][1]=0,c.x[1][2]=1,c.x[2][1]=0,c.x[2][2]=1;///单位矩阵,类似于快速幂中将sum初始化为1
    while(n)
    {
        if(n%2==1)  c=muti(c,a);
        a=muti(a,a);
        n>>=1;
    }
    return c;
}

int main()
{
    int k;cin>>k;
    while(k--)
    {
        mat a;
        a.x[1][1]=1,a.x[1][2]=1,a.x[2][1]=1,a.x[2][2]=0;
        mat c,d;
        int n;
        cin>>n;
        c=qmod(a,n);///n次方，就是斐波那契数列的第n项
        cout<<c.x[1][1]<<endl;///输出c11即可
    }
}