#include <iostream>
#include <cstdio>
using namespace std;
int const Mod = 666013;

struct Mat{
    int v [4][4];
    Mat (){
        for(int i = 1 ; i <= 3 ; ++ i)
            for(int j = 1 ; j <= 3 ; ++ j)
                v [i][j] = 0;
    }
    Mat operator * (Mat X) const{
        Mat r;
        for(int i = 1 ; i <= 3 ; ++ i)
            for(int j = 1 ; j <= 3 ; ++ j){
                int s = 0;
                for(int k = 1 ; k <= 3 ; ++ k)
                    s = (1LL * s + 1LL * v [i][k] * X.v [k][j]) % Mod;
                r.v [i][j] = s;
            }
        return r;
    }
    Mat operator << (int x) const{
        Mat c = *this , r;
        r.v [1][1] = r.v [2][2] = r.v [3][3] = 1;
        for(int i = 0 ; (1 << i) <= x ; ++ i){
            if ((1 << i) & x)
                r = r * c;
            c = c * c;
        }
        return r;
    }

}A , B;
int t , x , y , z , a , b , c , n;
int main()
{
    freopen ("iepuri.in" , "r" , stdin);
    freopen ("iepuri.out" , "w" , stdout);
    scanf ("%d" , &t);
    while (t --){
        scanf ("%d%d%d%d%d%d%d" , &x , &y , &z , &a , &b , &c , &n);
        Mat A , B;
        A.v [1][1] = x , A.v [1][2] = y , A.v [1][3] = z;
        B.v [1][1] = 0 , B.v [1][2] = 0 , B.v [1][3] = c;
        B.v [2][1] = 1 , B.v [2][2] = 0 , B.v [2][3] = b;
        B.v [3][1] = 0 , B.v [3][2] = 1 , B.v [3][3] = a;
        Mat r = (A * (B << (n - 2)));
        printf ("%d\n" , r.v [1][3]);
    }
    return 0;
}