hdu 3936 FIB Query 菲波拉切数列的性质+矩阵乘法

参考：http://blog.csdn.net/jxy859/article/details/6686305

　　1.f(0)+f(1)+f(2)+…+f(n)=f(n+2)-1。 

　　2.f(1)+f(3)+f(5)+…+f(2n-1)=f(2n)。 

　　3.f(2)+f(4)+f(6)+…+f(2n) =f(2n+1)-1。 

　　4.[f(0)]^2+[f(1)]^2+…+[f(n)]^2=f(n)·f(n+1)。 

　　5.f(0)-f(1)+f(2)-…+(-1)^n·f(n)=(-1)^n·[f(n+1)-f(n)]+1。 

　　6.f(m+n-1)=f(m-1)·f(n-1)+f(m)·f(n)。 

        利用这一点，可以用程序编出时间复杂度仅为O（log n）的程序。 

　　7.[f(n)]^2=(-1)^(n-1)+f(n-1)·f(n+1)。 

　　8.f(2n-1)=[f(n)]^2-[f(n-2)]^2。 

　　9.3f(n)=f(n+2)+f(n-2)。 

　　10.f(2n-2m-2)[f(2n)+f(2n+2)]=f(2m+2)+f(4n-2m) [ n〉m≥-1,且n≥1] 

    11.f(2n+1)=[f(n)]^2+[f(n+1)]^2.
#include<iostream>
#include <stdlib.h>
#include <algorithm>
#include <stdio.h>
#include<vector>
using namespace std;
typedef __int64 i64;
const int M=1000000007;
void dfs(i64 a[][2],i64 b[][2])
{
    int tmp[2][2];
    tmp[0][0]=(a[0][0]*b[0][0])%M+(a[0][1]*b[1][0])%M;
    tmp[0][1]=(a[0][0]*b[0][1])%M+(a[0][1]*b[1][1])%M;
    tmp[1][0]=(a[1][0]*b[0][0])%M+(a[1][1]*b[1][0])%M;
    tmp[1][1]=(a[1][0]*b[0][1])%M+(a[1][1]*b[1][1])%M;
    a[0][0]=tmp[0][0]%M;
    a[0][1]=tmp[0][1]%M;
    a[1][0]=tmp[1][0]%M;
    a[1][1]=tmp[1][1]%M;
}
void dfs2(i64 a[][2],i64 b[][2])
{
    int tmp[2][2];
    tmp[0][0]=(a[0][0]*b[0][0])%M+(a[0][1]*b[1][0])%M;
    tmp[0][1]=(a[0][0]*b[0][1])%M+(a[0][1]*b[1][1])%M;
    tmp[1][0]=(a[1][0]*b[0][0])%M+(a[1][1]*b[1][0])%M;
    tmp[1][1]=(a[1][0]*b[0][1])%M+(a[1][1]*b[1][1])%M;
    b[0][0]=tmp[0][0]%M;
    b[0][1]=tmp[0][1]%M;
    b[1][0]=tmp[1][0]%M;
    b[1][1]=tmp[1][1]%M;
}
i64 getF(i64 x)
{
    if(x==0)return 0;
    i64 init[2][2]={1,0,0,1};
    i64 t[2][2]={1,1,1,0};
    while(x>1)
    {
        if(x&1)dfs2(t,init);
        dfs(t,t);
        x/=2;
    }
    dfs(t,init);
    return t[0][1];
}
int main()
{
    int Case;
    i64 a,b;
    scanf("%d",&Case);
    while(Case--)
    {
        scanf("%I64d%I64d",&a,&b);
        a--;
        i64 ans2=getF(2*b)*getF(2*b+1);
        i64 ans1=getF(2*a)*getF(2*a+1);
        cout<<((ans2-ans1)%M+M)%M<<endl;
    }
    return 0;
}