斐波那契_矩阵乘法

Description

形如 1 1 2 3 5 8 13 21 34 55 89 144....的数列,求裴波拉契数列的第n项。

Input

 n (1〈 n 〈2^31）

Output

 一个数为裴波拉契数列的第n项mod 10;

题解

 a[1,1]:=0; a[1,2]:=1; a[2,1]:=1; a[2,2]:=1;

代码

type  arr=array[1..2,1..2] of longint;var  n:int64;  a,c:arr;procedure mi(a,b:arr);var  i,j,k:longint;begin  fillchar(c,sizeof(c),0);  for i:=1 to 2 do    for j:=1 to 2 do      for k:=1 to 2 do        c[i,k]:=(c[i,k]+a[i,j]*b[j,k]) mod 10;end;procedure main(n:int64);begin  if n<=1 then exit;  main(n div 2);  mi(c,c);  if odd(n) then mi(c,a);end;begin  readln(n); n:=n;a[1,1]:=0; a[1,2]:=1;  a[2,1]:=1; a[2,2]:=1;  c:=a;  main(n);  writeln((c[2,2]) mod 10);end.