Cod sursa(job #631445)

Utilizator sebii_cSebastian Claici sebii_c Data 8 noiembrie 2011 02:20:27
Problema Iepuri Scor 90
Compilator c Status done
Runda Arhiva de probleme Marime 1.1 kb
#include <stdio.h>
#include <string.h>
#define MOD 666013

int T[3][3], res[3][3], aux[3][3];

void matrix_mult(int (*A)[3], int (*B)[3], int (*C)[3])
{
    int i, j, k;
    memset(aux, 0, sizeof(aux));

    for (i=0; i<3; ++i)
	for (j=0; j<3; ++j)
	    for (k=0; k<3; ++k)
		aux[i][j] += (1ll*A[i][k]*B[k][j])%MOD;

    memcpy(C, aux, sizeof(aux));
} 
    
void fastexp(int exp)
{
    int i;
    for (i=0; (1<<i)<=exp; ++i) {
	if ((1<<i) & exp)
	     matrix_mult(res, T, res);
	matrix_mult(T, T, T);
    }
}

int main()
{
    freopen("iepuri.in", "r", stdin);
    freopen("iepuri.out", "w", stdout);
    int t, i;
    scanf("%d", &t);
    for (i=0; i<t; ++i) {
	int a, b, c, x, y, z;
	long n;
	scanf("%d %d %d %d %d %d %ld", &x, &y, &z, &a, &b, &c, &n);
	memset(T, 0, sizeof(T));
	memset(res, 0, sizeof(res));
	T[0][1] = 1;
	T[1][2] = 1;
	T[2][0] = c, T[2][1] = b, T[2][2] = a;
	res[0][0] = res[1][1] = res[2][2] = 1;
	fastexp(n-2);
	long long result = (1ll*res[2][0]*x%MOD + 1ll*res[2][1]*y%MOD + 1ll*res[2][2]*z%MOD)%MOD;
	printf("%lld\n", result);
    }
    return 0;
}