51nod-1358 浮波那契

思路：矩阵快速幂，对于3.4可乘以5变成17，同时n*=5, FB(n)=FB(n-5)+FB(n-17), 而n<=20时 FB(n)=1。 关键是对于初始化矩阵的构造，即e[m][m]*FB[n-1,n-2,...,n-17]=FB[n,n-1,...,n-16]

e[][]=就可以满足条件

Code ：

#include<iostream>
#include<cstring>
using namespace std;
typedef long long LL;

const int MOD=1e9+7;
const int MAX_S=20;
LL n,m;
LL e[MAX_S][MAX_S],sum[MAX_S][MAX_S];

LL Matrix_multiply(void);
void Mulrtiplt(LL a[][MAX_S],LL b[][MAX_S]);
int main()
{
	ios::sync_with_stdio(false);
	cin>>n;
	if(n<=4)	cout<<1<<endl;
	else	cout<<Matrix_multiply()<<endl;
	return 0;
}

LL Matrix_multiply(void)
{
	m=17;
	e[1][5]=e[1][m]=1;
	for(int i=1;i<=m;++i)
		e[i][i-1]=sum[i][i]=1;
	for(int i=1;i<=m;++i)
	{
		for(int j=1;j<=m;++j)
			cout<<e[i][j]<<" ";
		cout<<endl;
	}
	
	n=n*5-20;
	while(n){
		if(n%2)	 Mulrtiplt(sum,e);
		Mulrtiplt(e,e);
		n/=2;
	}
	LL ans=0;
	for(int i=1;i<=m;++i)
		ans=(ans+sum[1][i])%MOD;
	return ans;
}

void Mulrtiplt(LL a[][MAX_S],LL b[][MAX_S])
{
	LL d[MAX_S][MAX_S];
	memset(d,0,sizeof(d));
	for(int i=1;i<=m;++i)
		for(int k=1;k<=m;++k)
			for(int j=1;j<=m;++j)
				d[i][j]=(d[i][j]+a[i][k]*b[k][j])%MOD;
	for(int i=1;i<=m;++i)
		for(int j=1;j<=m;++j)
			a[i][j]=d[i][j];
}