poj3070--矩阵 快速幂

第一次写矩阵的题目，第一次写快速幂的题目，第一次进行移位操作，第一次进行按位与操作

//#define LOCAL
#include <stdio.h>
#include <string.h>

#define MAXN 20 + 10

int n;
int a11, a12, a21, a22;
int b11, b12, b21, b22;
int c11, c12, c21, c22;
int d11, d12, d21, d22;

int main()
{
#ifdef LOCAL
	freopen("input.txt", "r", stdin);
	//freopen("output.txt", "w", stdout);
#endif
	while(~scanf("%d", &n) && n != -1)
	{
		// 数据初始化
	 	a11 = 1;
	 	a12 = 1;
	 	a21 = 1;
	 	a22 = 0;
	 	b11 = 1;
	 	b12 = 0;
	 	b21 = 0;
	 	b22 = 1;
		// 算法主体
		if(n == 0)
		{
			printf("0\n");
			continue;
		}
			
		while(n)
		{
			if(n & 1)
			{
				c11 = (a11 * b11 + a12 * b21) % 10000;
				c12 = (a11 * b12 + a12 * b22) % 10000;
				c21 = (a21 * b11 + a22 * b21) % 10000;
				c22 = (a21 * b12 + a22 * b22) % 10000;
				b11 = c11;
				b12 = c12;
				b21 = c21;
				b22 = c22;
			}
			d11 = (a11 * a11 + a12 * a21) % 10000;
			d12 = (a11 * a12 + a12 * a22) % 10000;
			d21 = (a21 * a11 + a22 * a21) % 10000;
			d22 = (a21 * a12 + a22 * a22) % 10000;
			a11 = d11;
			a12 = d12;
			a21 = d21;
			a22 = d22;
			n >>= 1; 
		}
	 	printf("%d\n", b12); 
	}
	return 0;
}