Fibonacci[POJ3070]

最新推荐文章于 2021-04-25 19:21:33 发布

zsyz_lb2003

最新推荐文章于 2021-04-25 19:21:33 发布

阅读量130

点赞数 1

分类专栏：蓝书做题记录 caioj数论蓝书数论矩阵乘法

本文链接：https://blog.csdn.net/zyszlb2003/article/details/89514562

版权

蓝书做题记录同时被 3 个专栏收录

66 篇文章 2 订阅

订阅专栏

caioj数论

54 篇文章 0 订阅

订阅专栏

蓝书数论

29 篇文章 1 订阅

订阅专栏

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门

思路

通过矩阵乘法加速，
设
$F(0)=\begin{bmatrix}Fib_0=0,Fib_1=1\end{bmatrix},A=\begin{bmatrix}0&1\\1&1\end{bmatrix}$

$F(1)=F(0)*A=\begin{bmatrix}Fib_0*0+Fib_1*1=Fib_1,Fib_0*1+Fib_1*1=Fib_2\end{bmatrix}$

$. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .$

$F(n)=F(n-1)*A=\begin{bmatrix}Fib_n,Fib_{n+1}\end{bmatrix}$

$F(n)=F(n-2)*A^2$

$. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .$

$F(n)=F(0)*A^n$

$A^n$ 可用矩阵快速幂 $\log(n)$ 求解。

成功实现加速！

AC code

#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<cstring>
#define gc getchar()
using namespace std;
const int mod=10000;
inline void qr(int &x)
{
	x=0;int f=1;char c=gc;
	while(c<'0'||c>'9'){if(c=='-')f=-1;c=gc;}
	while(c>='0'&&c<='9'){x=x*10+(c^48);c=gc;}
	x*=f;
}
void mul(int f[2],int a[2][2])
{
	int c[2];memset(c,0,sizeof(c));
	for(int j=0;j<2;j++)
		for(int k=0;k<2;k++)
			c[j]=(c[j]+f[k]*a[k][j])%mod;
	memcpy(f,c,sizeof(c));
}
void muls(int a[2][2])
{
	int c[2][2];memset(c,0,sizeof(c));
	for(int i=0;i<2;i++)
		for(int j=0;j<2;j++)
			for(int k=0;k<2;k++)
				c[i][j]=(c[i][j]+a[i][k]*a[k][j])%mod;
	memcpy(a,c,sizeof(c));
}
int main()
{
	int n;
	while(qr(n),n!=-1)
	{
		int f[2]={0,1};
		int a[2][2]={{0,1},{1,1}};
		while(n)
		{
			if(n&1)mul(f,a);
			n>>=1;muls(a);
		}
		printf("%d\n",f[0]);
	}
	return 0;
}