フィボナッチ数列II
制限時間:1000MS
メモリ制限:65536K
説明
1 1 2 3 5 8 13 21 34 55 89 144 ...の形式のシーケンスは、フィボナッチ数列のn番目の項を見つけます。
入力
n(1 〈n 〈2 31)
出力
数値は、フィボナッチ数列のn番目の項mod10000です。
サンプル入力
123456789
サンプル出力
4514
問題解決のアイデア
この質問は、元のタイトルのP1962フィボナッチ数(高速行列乗算+累乗)
で、数%の変更のみです。
ACコード
#include<cstdio>
using namespace std;
long long n,k;
struct node
{
long long oo[5][5];
}a;
node operator*(node x,node y)//重新定义“*”,变为矩阵乘法
{
node z;
for(long long i=1;i<=2;i++)//都改为2
for(long long j=1;j<=2;j++)
z.oo[i][j]=0;
for(long long o=1;o<=2;o++)
for(long long i=1;i<=2;i++)
for(long long j=1;j<=2;j++)
z.oo[i][j]=(z.oo[i][j]+x.oo[i][o]*y.oo[o][j])%10000;
return z;
}
node ksm(node x,long long k)//快速幂
{
node y=x;
while(k)
{
if(k%2)y=x*y;
x=x*x;
k/=2;
}
return y;
}
int main()
{
scanf("%lld",&n);
a.oo[1][1]=0;//初始值
a.oo[1][2]=a.oo[2][1]=a.oo[2][2]=1;
a=ksm(a,n);
printf("%lld",a.oo[1][1]);
return 0;
}