思路与二进制快速幂相类似,就是矩阵的乘法用一个函数表示,其它与快速幂几乎一模一样。
#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;
}