矩阵乘法，留待熟悉一下类再学习。

思路与二进制快速幂相类似，就是矩阵的乘法用一个函数表示，其它与快速幂几乎一模一样。

#include<stdio.h>

struct mat
{

       int matrix[2][2];

};
typedef struct mat Matrix;
Matrix mul(Matrix c, Matrix b)
{
Matrix result;
result.matrix[0][0] = (c.matrix[0][0]*b.matrix[0][0] + c.matrix[0][1]*b.matrix[1][0]) % 10000;
result.matrix[0][1] = (c.matrix[0][0]*b.matrix[0][1] + c.matrix[0][1]*b.matrix[1][1]) % 10000;
result.matrix[1][0] = (c.matrix[1][0]*b.matrix[0][0] + c.matrix[1][1]*b.matrix[1][0]) % 10000;
result.matrix[1][1] = (c.matrix[1][0]*b.matrix[0][1] + c.matrix[1][1]*b.matrix[1][1]) % 10000;
return result;
}
void qpow(int N)
{

        Matrix  c, d;

        c.matrix[0][0] = 1;c.matrix[0][1]=0;c.matrix[1][0]=0;c.matrix[1][1]=1;    //存a^b

        d.matrix[0][0] = 1;d.matrix[0][1]=1;d.matrix[1][0]=1;d.matrix[1][1]=0;    //存a的倍幂

        while (N > 0)

       {

               if (N & 1)     //或 if (b % 2 == 1)

                     c = mul(c,d);

              N = N >> 1;   //或  b = b / 2

              d = mul(d,d);

       }

       printf("%d\n",c.matrix[0][1]);

}
int main()
{

       int N;

       while(EOF!=scanf("%d",&N) && N!=-1)

       {

               qpow(N);

       }

       return 0;

}