【poj3070】 Fibonacci_4165 51.fibonacci(poj3070)-CSDN博客

本文介绍了一种使用矩阵快速幂的方法来高效地求解Fibonacci数列的第n项。通过定义特定的矩阵并利用快速幂运算，可以在对数时间内得到结果。

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

题意：用矩阵乘法求fibonacci数列的第n项。

Solution
　　矩乘入门题啊，题目把题解已经说的很清楚里= =。
　　矩乘其实很简单，通过自己YY或者是搜索对于一个递推公式求出它所对应的矩阵，然后套个快速幂就可以迅速求解第n项。

代码：

// poj3070
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define LL long long
#define inf 2147483640
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;

int n,a[2][2],b[2][2];

void mul(int a[2][2],int b[2][2],int ans[2][2]) {
    int t[2][2];
    for (int i=0;i<2;i++)
        for (int j=0;j<2;j++) {
            t[i][j]=0;
            for (int k=0;k<2;k++) t[i][j]=(t[i][j]+a[i][k]*b[k][j])%10000;
        }
    for (int i=0;i<2;i++)
        for (int j=0;j<2;j++) ans[i][j]=t[i][j];
}
void pow(int k) {
    while (k) {
        if (k&1) mul(a,b,b);
        k>>=1;
        mul(a,a,a);
    }
}
int main() {
    while (scanf("%d",&n)!=EOF) {
        if (n==-1) break;
        a[0][0]=a[0][1]=a[1][0]=1;a[1][1]=0;
        b[0][0]=b[1][1]=1;
        b[1][0]=b[0][1]=0;
        pow(n);
        printf("%d\n",b[1][0]);
    }
    return 0;
}