#include <bits/stdc++.h>
#include <fstream>
using namespace std;
typedef long long int64;
string __fname = "iepuri"; ifstream in (__fname + ".in"); ofstream out (__fname + ".out"); 
#define cin in 
#define cout out
#define sz(x)(int) x.size()
#define pb(x) push_back(x)
#define pii pair <int, int> 
#define pii64 pair <int64, int64>
#define ss cout << " "
#define nn cout << "\n"
#define ct(x) cout << (x)
#define cts(x) cout << (x) << " "
#define ctn(x) cout << (x) << "\n"
#define db(x) cout << "> " << #x << ": " << (x) << "\n"
#define qr queries()
#define yn(x) if (x) {ctn("Yes");}else {ctn("No");}
#define rt return
void solve(int);
void queries(){int n;cin >> n;for (int i = 1; i <= n; i++) solve(i);}
int64 ceildiv(int64 a, int64 b) {return a / b + !!(a % b);}
template<class T>T gcd (T a, T b){return (b ? gcd(b, a % b): a);}
template<class T>T lcm (T a, T b){return a * b / gcd(a, b);}
// // // // // // // // // // // // // // // // // // // // // // 
/*                  TEMPLATE - VANILLA                         */
// // // // // // // // // // // // // // // // // // // // // //
const int maxn = 200200;
const int64 mod = 666013;
const double pi = 3.14159265359;
const int ddx[] = {-1, -1, 0, 1, 1, 1, 0, -1};
const int ddy[] = {0, 1, 1, 1, 0, -1, -1, -1};
const int dx[] = {-1, 0, 1, 0};
const int dy[] = {0, 1, 0, -1};

void mul (int64 a[3][3], int64 b[3][3]) {
    int64 c[3][3];
    for (int i = 0; i < 3; i++)
        for (int j = 0; j < 3; j++) c[i][j] = 0;

    for (int i = 0; i < 3; i++){
        for (int j = 0; j < 3; j++) {
            for (int k = 0; k < 3; k++){
                c[i][j] += (a[i][k] * b[k][j]) % mod;
                c[i][j]%=mod;
            }
        }
    }
    for (int i = 0; i < 3; i++){
        for (int j = 0; j < 3; j++) {
            a[i][j] = c[i][j];
        }
    }
}

void pw (int64 a[3][3], int64 mc[3][3], int k) {
    if (k <= 1) return;
    pw (a, mc, k / 2);
    mul(a, a);
    if (k % 2 == 1) mul(a, mc);
}
void afis (int64 a[3][3]) {
    nn;
    for (int i = 0; i < 3; i++){
        for (int j = 0; j < 3; j++){
            cts(a[i][j]);
        }
        nn;
    }
}

void solve(int id){
    int64 x,y,z,a,b,c,n;
    cin >> x >> y >> z >> a >> b >> c >> n;
    int64 ma[3][3] = {
        {x, y, z},
        {0, 0, 0},
        {0, 0, 0}
    };

    int64 mb[3][3] = {
        {0, 0, c},
        {1, 0, b},
        {0, 1, a}
    };

    int64 mc[3][3] = {
        {0, 0, c},
        {1, 0, b},
        {0, 1, a}
    };
    pw(mb, mc, n);
    mul(ma, mb);
    ctn(ma[0][0]);

    rt;
}


int main(){
    ios_base::sync_with_stdio(0);cin.tie(0); ct(fixed); ct(setprecision(10)); 
    qr;
    rt 0;
}