第二场 hdu 6050 Funny Function

Function Fx,yFx,ysatisfies: 

For given integers N and M,calculate Fm,1Fm,1 modulo 1e9+7. 

InputThere is one integer T in the first line. 
The next T lines,each line includes two integers N and M . 
1<=T<=10000,1<=N,M<2^63. 
OutputFor each given N and M,print the answer in a single line. 
Sample Input

2
2 2
3 3

Sample Output

2
33


题目大意：按照题目所给的方式求F[m][1]?
解题思路：对于任意i>=1，当j>=3时，有  通过归纳法可以得到  进而推导出  通过矩阵快速幂求解。
这个思路能够推出矩阵A、B0、B1  

 然后能够用快速幂模板很容易求出F[m][1]


AC代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <bits/stdc++.h>
#define N 101000
#define ll long long
using namespace std;
const int SMod=1e9+7;
struct Matrix
{
    ll m[2][2];
};
Matrix Mul(Matrix a,Matrix b)
{
    Matrix c;
    memset(c.m,0,sizeof(c.m));
    for(int i=0; i<2; i++)
        for(int j=0; j<2; j++)
            for(int k=0; k<2; k++)
                c.m[i][j]+=((a.m[i][k]*b.m[k][j])%SMod+SMod)%SMod;
    return c;
}
Matrix fastm(Matrix a,ll n)
{
    Matrix res;
    memset(res.m,0,sizeof(res.m));
    res.m[0][0]=res.m[1][1]=1;
    while(n)
    {
        if(n&1)
            res=Mul(res,a);
        n>>=1;
        a=Mul(a,a);
    }
    return res;
}
int main()
{
    Matrix A,B0,B1;
    ll t,n,m;
    //freopen("1006.in","r",stdin);
    //freopen("1006.out","w",stdout);
    while(~scanf("%lld",&t))
    {
        while(t--)
        {
            A.m[0][0]=0,A.m[0][1]=1;
            A.m[1][0]=2,A.m[1][1]=1;
            B0.m[0][0]=1,B0.m[0][1]=0;
            B0.m[1][0]=0,B0.m[1][1]=1;
            B1.m[0][0]=-1,B1.m[0][1]=1;
            B1.m[1][0]=2,B1.m[1][1]=0;
            scanf("%lld%lld",&n,&m);
            A=fastm(A,n);
//            for(int i=0;i<2;i++)
//            {
//                for(int j=0;j<2;j++)
//                printf("%lld ",A.m[i][j]);
//                printf("
");
//            }
            if(n%2==0)
            {
                for(int i=0; i<2; i++)
                    for(int j=0; j<2; j++)
                        A.m[i][j]-=B0.m[i][j];
            }
            else
            {
                for(int i=0; i<2; i++)
                    for(int j=0; j<2; j++)
                        A.m[i][j]-=B1.m[i][j];
            }
            A=fastm(A,m-1);
            printf("%lld
",(A.m[0][0]+A.m[0][1]+SMod)%SMod);
        }
    }
    return 0;
}