本题主要解决一个公式问题:
以下分析为转载:
={*()^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;
}