九度OJ 1081 递推数列 -- 矩阵二分乘法

题目地址：http://ac.jobdu.com/problem.php?pid=1081

题目描述：

给定a0,a1,以及an=p*a(n-1) + q*a(n-2)中的p,q。这里n >= 2。 求第k个数对10000的模。

输入：

输入包括5个整数：a0、a1、p、q、k。

输出：

第k个数a(k)对10000的模。

样例输入：

20 1 1 14 5

样例输出：

8359

#include <stdio.h> #define MOD 10000 typedef struct matrix{    int a00, a01, a10, a11;}Matrix; void MatrixMul(Matrix * m, Matrix * n){    Matrix tmp;    tmp.a00 = (m->a00 * n->a00 + m->a01 * n->a10) % MOD;    tmp.a01 = (m->a00 * n->a01 + m->a01 * n->a11) % MOD;    tmp.a10 = (m->a10 * n->a00 + m->a11 * n->a10) % MOD;    tmp.a11 = (m->a10 * n->a01 + m->a11 * n->a11) % MOD;    *m = tmp;} int main(void){    int a0, a1, p, q, k;     while (scanf ("%d%d%d%d%d",            &a0, &a1, &p, &q, &k) != EOF){        if (k == 0) {printf ("%d\n", a0 % MOD); continue;}        if (k == 1) {printf ("%d\n", a1 % MOD); continue;}        Matrix pq = {p%MOD, q%MOD, 1, 0};        Matrix ans = {1, 0, 0, 1};        --k;        while (k >= 1){            if ((k & 1) == 1)                MatrixMul (&ans, &pq);            MatrixMul (&pq, &pq);            k = k >> 1;        }        printf ("%d\n", (ans.a00*a1+ans.a01*a0) % MOD);    }     return 0;}

参考资料：十个利用矩阵乘法解决的经典题目