ACM1005：Number Sequence

#include<stdio.h>
#include<string.h>
#define SIZE 2
long long int ** multiple(long long int a[][SIZE], long long int b[][SIZE]);
int main()
{
	int A, B;
	long long int n;
	long long int factor[SIZE][SIZE];
	long long int ans[SIZE][SIZE];
 	while (scanf("%d%d%lld", &A, &B, &n) && (A||B||n))
	{
		memset(factor, 0, 2 * 2 * sizeof(int));
		memset(ans, 0, 2 * 2 * sizeof(int));
		factor[0][0] = A;
		factor[0][1] = B;
		factor[1][0] = 1;
		factor[1][1] = 0;
		//矩阵快速幂
		for (int i = 0; i < SIZE; i++)
			ans[i][i] = 1;
		n = n - 2;
		while (n > 0)
		{
			if (n & 1)
				memcpy(ans, multiple(ans, factor), 2 * 2 * sizeof(long long int));
			memcpy(factor, multiple(factor, factor), 2 * 2 * sizeof(long long int));
			n >>= 1;
		}
		printf("%lld\n", (ans[0][0] + ans[0][1]) % 7);
	}
	return 0;
}

long long int ** multiple(long long int a[][SIZE], long long int b[][SIZE])
{
	int i, j, k;
	long long int c[SIZE][SIZE];
	memset(c, 0, 2 * 2 * sizeof(long long int));
	for(i = 0; i < SIZE; i++)
		for(j = 0; j < SIZE; j++)
			for (k = 0; k < SIZE; k++)
			{
				c[i][j] += a[i][k] * b[k][j];
				c[i][j] = c[i][j] % 7;
			}
	return c;
}