Cod sursa(job #2218176)

Utilizator PetyAlexandru Peticaru Pety Data 3 iulie 2018 15:46:32
Problema Iepuri Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.45 kb
#include <bits/stdc++.h>

using namespace std;

ifstream fin ("iepuri.in");
ofstream fout ("iepuri.out");

int sol[102];
const int MOD = 666013;
long long mat[4][4], rez[4][4];
int x, y, z, a, b, c, n, t;

void inmultire (long long a[4][4], long long b[4][4]) {
  long long aux[4][4];
  for (int i = 1; i <= 3; i++)
    for (int j = 1; j <= 3; j++)
      aux[i][j] = 0;
  for (int i = 1; i <= 3; i++)
    for (int j = 1; j <= 3; j++)
      for (int k = 1; k <= 3; k++)
        aux[i][j] = 1LL * (aux[i][j] + a[i][k] * b[k][j] % MOD) % MOD;
  for (int i = 1; i <= 3; i++)
    for (int j = 1; j <= 3; j++)
      a[i][j] = aux[i][j];
}


void put (int k) {
  while (k) {
    if (k & 1)
      inmultire(rez, mat);
    inmultire(mat, mat);
    k >>= 1;
  }
}

int main()
{
  fin >> t;
  for (int T = 1; T <= t; T++) {
    fin >> x >> y >> z >> a >> b >> c >> n;
    sol[0] = x; sol[1] = y; sol[2] = z;
    for (int i = 3; i <= n; i++) {
      sol[i] = (sol[i - 1] * a % MOD + sol[i - 2] * b % MOD + sol[i - 3] * c % MOD) % MOD;
    }
    mat[1][1] = mat[1][2] = mat[2][2] =  mat[3][1] = 0;
    mat[2][1] = mat[3][2] = 1;
    mat[1][3] = c;
    mat[2][3] = b;
    mat[3][3] = a;
    for (int i = 1; i <= 3; i++)
      for (int j = 1; j <= 3; j++)
        rez[i][j] = mat[i][j];
    put(n - 3);
    fout << 1LL * (1LL * rez[1][3] * x + 1LL * rez[2][3] * y + 1LL * rez[3][3] * z) % MOD << " " << sol[n] << "\n";
  }
  return 0;
}