Pagini recente » Cod sursa (job #2801149) | Cod sursa (job #516008) | Cod sursa (job #1672162) | Cod sursa (job #2901083) | Cod sursa (job #1819189)
#include <fstream>
using namespace std;
const int MOD = 666013;
ifstream f("iepuri.in");
ofstream g("iepuri.out");
int X, Y, Z, AA, BB, CC;
void mul(int A[3][3], int B[3][3]) // A <-- A*B
{
long long C[3][3];
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
{
C[i][j] = 0;
for(int k = 0; k < 3; k++)
C[i][j] += (long long)A[i][k] * B[k][j];
}
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
A[i][j] = C[i][j] % MOD;
}
void mul2(int A[3], int B[3][3]) // A <-- A*B
{
long long C[3];
for(int i = 0; i < 3; i++)
{
C[i] = 0;
for(int k = 0; k < 3; k++)
C[i] += (long long)A[k] * B[k][i];
}
for(int i = 0; i < 3; i++)
A[i] = C[i] % MOD;
}
int putere(int n)
{
int A[3][3] = {{0, 0, CC}, {1, 0, BB}, {0, 1, AA}};
int P[3] = {X, Y, Z};
while(n > 0)
{
while(n % 2 == 0)
{
mul(A, A);
n /= 2;
}
mul2(P, A);
n--;
}
return P[2];
}
int main()
{
int T, N;
f >> T;
while(T--)
{
f >> X >> Y >> Z >> AA >> BB >> CC;
f >> N;
g << putere(N - 2) << '\n';
}
f.close();
g.close();
return 0;
}