hdu 2855 Fibonacci Check-up

本题主要解决一个公式问题：

以下分析为转载：

={*()^k -()^k}

因为：

（a+1）^n=*a^k;

所以将上式化简成：

=（）^n+()^n

                  分子分母都乘*2    =()^n+()^n

                                         =()^2n+()^2n

                                         =F[ 2*N ]

再用矩阵快速幂解决就相当简单了。

下面是自己写的代码：

#include<iostream>
#include<cstdio>
#include<iostream>
using namespace std;
int m;
struct matrix
{
    int a[2][2];
};
matrix multi(matrix m3,matrix m4)
{
    matrix mm;
    for(int i=0;i<2;i++)
    {
        for(int j=0;j<2;j++)
        {
            mm.a[i][j]=0;
            for(int k1=0;k1<2;k1++)
            {
                mm.a[i][j]+=((m3.a[i][k1]%m*m4.a[k1][j]%m)%m);
            }
        }
    }
    return mm;
}
matrix cal(matrix m1,int k)
{
    matrix ma;
    ma.a[0][0]=1;
    ma.a[0][1]=0;
    ma.a[1][0]=0;
    ma.a[1][1]=1;
    while(k)
    {
        if(k%2)
            ma=multi(ma,m1);
        m1=multi(m1,m1);
        k/=2;
    }
    return ma;
}
int main()
{
    matrix m1;
    matrix ma;
    m1.a[0][0]=1;
    m1.a[0][1]=1;
    m1.a[1][0]=1;
    m1.a[1][1]=0;
    //(f(n),f(n-1))=((f(2),f(1))*m^n-2;

    int T,n;
    cin>>T;
    for(int i=0;i<T;i++)
    {
        cin>>n>>m;
        if(n==0)
            cout<<0<<endl;
        else if(n==1)
            cout<<1%m<<endl;
        else
        {
            ma=cal(m1,2*n-2);
            cout<<(ma.a[0][0]%m+ma.a[1][0]%m)%m<<endl;
        }

    }
    return 0;

}