Fibonacci 第 n 项

这题是因为数据量很大，到了后面数组就没有办法保存斐波那契数了

考虑1×2的矩阵【 f [ n − 2 ] , f [ n − 1 ] 】。根据fibonacci数列的递推关系，我们希望通过乘以一个2 × 2 的矩阵，得到矩阵【 f [ n − 1 ] , f [ n ] 】 = 【 f [ n − 1 ] , f [ n − 1 ] + f [ n − 2 ] 】 很容易构造出这个2 × 2 矩阵A ，即：

所以，有【 f [ 1 ] , f [ 2 ] 】 × A ＝ 【 f [ 2 ] , f [ 3 ] 】
又因为矩阵乘法满足结合律，故有：
【 f [ 1 ] , f [ 2 ] 】 × A n − 1 = 【 f [ n ] , f [ n + 1 ] 】
这个矩阵的第一个元素即为所求。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll; 
ll n,mod;
struct nn
{
	ll x[3][3];
	nn()
	{
		memset(x,0,sizeof(x));
	}
}f;
nn solve(nn a,nn b)
{
	nn c;
	for(int i=1;i<=2;i++)
		for(int j=1;j<=2;j++)
			for(int k=1;k<=2;k++)
			{
				c.x[i][j]+=a.x[i][k]*b.x[k][j];
				c.x[i][j]%=mod;
			}
	return c;
}
nn solve1(nn a,nn b,int n,int m,int p)
{
	nn c;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=p;j++)
			for(int k=1;k<=m;k++)
			{
				c.x[i][j]+=a.x[i][k]*b.x[k][j];
				c.x[i][j]%=mod;
			}
	return c;
}
nn ac(int n)//快速幂 
{
	if(n==1)
		return f;
	nn t=ac(n/2);
	if(n%2==1)
		return solve(solve(t,t),f);
	else
		return solve(t,t);
}
int main()
{
	cin>>n>>mod;
	f.x[1][1]=1;
	f.x[1][2]=1;
	f.x[2][1]=1;
	nn r,p;
	p.x[1][1]=1;//斐波那契的第1、2项都是1
	p.x[2][1]=1;
	r=solve1(ac(n-2),p,2,2,1);//2*2大小的矩阵和2*1大小的矩阵
	cout<<r.x[1][1]%mod;
}