51Nod - 1242 斐波那契数列的第N项

本文介绍了一种高效计算斐波那契数列的方法——矩阵快速幂,并提供了完整的C++实现代码。该方法利用矩阵乘法和快速幂运算,有效解决了斐波那契数列的大数计算问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

斐波那契数列的定义如下:

F(0) = 0
F(1) = 1
F(n) = F(n - 1) + F(n - 2) (n >= 2)

(1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377, ...)
给出n,求F(n),由于结果很大,输出F(n) % 1000000009的结果即可。

Input 输入1个数n(1 <= n <= 10^18)。 Output 输出F(n) % 1000000009的结果。 Sample Input
11
Sample Output
89

矩阵快速幂模板。

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
typedef long long ll;
const int N = 2;
const int MOD = 1000000009;
struct Matrix
{
	ll m[N][N];
};
Matrix A = {1,1,1,0};
Matrix I = {1,0,0,1};
Matrix multi(Matrix a,Matrix b)
{
	Matrix c;
	
	for(int i = 0 ; i < N ; i++)
	{
		for(int j = 0 ; j < N ; j++)
		{
			c.m[i][j] = 0;
			for(int k = 0 ; k < N ; k++)
			{
				c.m[i][j] += a.m[i][k] * b.m[k][j] % MOD;
			}
			c.m[i][j] %= MOD;
		}
	} 
	return c;
}
Matrix power(Matrix A,ll k)
{
	Matrix ans = I ,p = A;
	while(k)
	{
		if(k & 1)
		{
			ans = multi(ans,p);
			k--;
		}
		k>>=1;
		p = multi(p,p);
	}
	return ans;
}
int main()
{
	ll n;
	while(~scanf("%lld",&n))
	{
		Matrix ans = power(A,n - 1);
		printf("%lld\n",ans.m[0][0]);
	}
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值