#include <bits/stdc++.h>

std :: ifstream in ("iepuri.in");

std :: ofstream out ("iepuri.out");

const int DIM = 3;

const int mod = 666013;

int t;

int n;

int x;

int y;

int z;

int a;

int b;

int c;

int start[DIM][DIM];

int p[DIM][DIM];

int A[DIM][DIM];

void init()
{
    for(int i = 0; i < 3; i ++)
    {
        for(int j = 0; j < 3; j ++)
        {
            start[i][j] = 0;

            p[i][j] = 0;

            A[i][j] = 0;
        }
    }

    start[0][0] = z;

    start[0][1] = y;

    start[0][2] = x;

    for(int i = 0; i < 3; i ++)
    {
        p[i][i] = 1;
    }

    A[0][0] = a;

    A[1][0] = b;

    A[2][0] = c;

    A[0][1] = 1;

    A[1][2] = 1;
}

void multiply(int a[3][3], int b[3][3])
{
    int 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] = (1ll * c[i][j] + a[i][k] * b[k][j]) % mod;
            }
        }
    }

    for(int i = 0; i < 3; i ++)
    {
        for(int j = 0; j < 3; j ++)
        {
            a[i][j] = c[i][j];
        }
    }
}

void pow(int n)
{
    while(n)
    {
        if(n % 2)
        {
            multiply(p, A);
        }

        multiply(A, A);

        n /= 2;
    }
}

int main()
{

    in >> t;

    while(t --)
    {
        in >> x >> y >> z >> a >> b >> c >> n;

        init();

        pow(n - 2);

        multiply(start, p);

        out << start[0][0] << '\n';
    }

    return 0;
}