POJ - 3070 Fibonacci（快速幂）

题目大意：输入n求$F_{n}$的后四位，四位数中第一个非 0 数前面的 0 不需要输出
解题思路：题目给出了矩阵求法：
$\begin{equation} \left[\begin{matrix} F_{n+1}&F_{n}\\ F_{n}&F_{n-1} \end{matrix} \right] \end{equation}$ $= \begin{equation} \left[\begin{matrix} 1&1\\ 1&0 \end{matrix} \right] \end{equation}^{n}$
显然，要求$F_{n}$就要求得$\begin{equation} \left[\begin{matrix} 1&1\\ 1&0 \end{matrix} \right] \end{equation}^{n}$，然后取 [0][1] 位置上的数即可。
这里就用到快速幂，在过程中涉及到的矩阵乘法题目已经给出。
第一题快速幂，比较迷茫（躺）

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
const int mod = 10000;
struct node {
    int f[2][2];
};
node mul(node a, node b) {
    node tmp;
    for (int i = 0; i < 2; i++)
        for (int j = 0; j < 2; j++) {
            tmp.f[i][j] = 0;
            for (int k = 0; k < 2; k++)
                tmp.f[i][j] = (tmp.f[i][j] + a.f[i][k] * b.f[k][j]) % mod;
        }
    return tmp;
}
int pow(node &a, int n) {
    node tmp;
    tmp.f[0][1] = tmp.f[1][0] = 0;
    tmp.f[0][0] = tmp.f[1][1] = 1;
    while (n) {
        if (n&1) tmp = mul(tmp, a);
        a = mul(a, a);
        n >>= 1;
    }
    return tmp.f[0][1];
}
int main() {
    int n;
    while (scanf("%d", &n) && n != -1) {
        node F;
        F.f[0][0] = F.f[0][1] = F.f[1][0] = 1;
        F.f[1][1] = 0;
        printf("%d\n", pow(F, n));
    }
    return 0;
}