Cod sursa(job #774822)

Utilizator tester9x9Tester9x9 tester9x9 Data 6 august 2012 15:27:00
Problema Iepuri Scor 90
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.98 kb
#include <fstream>
#include <cstdio>
#define MOD 666013

using namespace std;


FILE* fin=fopen("iepuri.in","r");
FILE* fout=fopen("iepuri.out","w");
const int maxb=8192;
char buf[maxb];
int ptr=0;

inline int getint() {
    int nr=0;
    while(buf[ptr]<'0'||'9'<buf[ptr])
        if (++ptr>=maxb) fread(buf,maxb,1,fin),ptr=0;
    while ('0'<=buf[ptr]&&buf[ptr]<='9') {
        nr=nr*10+buf[ptr]-'0';
        if (++ptr>=maxb) fread(buf,maxb,1,fin),ptr=0;
    }
    return nr;
}


int A,B,C,X,Y,Z,N;
int M1[3][3],MX[3][3],AX[3][3],R[3][3];

void Multiply (int A[3][3], int B[3][3], int C[3][3])
{
    int i,j,k;
    for (i=0; i<3; i++)
        for (j=0; j<3; j++)
        {
            C[i][j]=0;
            for (k=0; k<3; k++)
            {
                C[i][j]+=1LL*A[i][k]*B[k][j]%MOD;
                if (C[i][j]>=MOD) C[i][j]-=MOD;
            }
        }
    for (i=0;i<3;i++)
        for (j=0;j<3;j++)
            A[i][j]=C[i][j];
}

void LgPow (int P)
{
    int i,j;
    for (i=0; i<3; i++)
        for (j=0; j<3; j++)
            R[i][j]=0;
    R[0][0]=R[1][1]=R[2][2]=1;
    for (i=1; i<=P; i<<=1)
    {
        if (i&P)
            Multiply(R,MX,AX);
        Multiply(MX,MX,AX);
    }
    for (i=0;i<3;i++)
        for (j=0;j<3;j++)
            MX[i][j]=R[i][j];
}

int Solve ()
{
    M1[0][0]=Z;
    M1[0][1]=Y;
    M1[0][2]=X;
    M1[1][0]=0;
    M1[1][1]=0;
    M1[1][2]=0;
    M1[2][0]=0;
    M1[2][1]=0;
    M1[2][2]=0;
    MX[0][0]=A;
    MX[0][1]=1;
    MX[0][2]=0;
    MX[1][0]=B;
    MX[1][1]=0;
    MX[1][2]=1;
    MX[2][0]=C;
    MX[2][1]=0;
    MX[2][2]=0;
    LgPow(N);
    Multiply(M1,MX,AX);
    return M1[0][0];
}


int T;

int main ()
{
    for (T=getint(); T; --T)
    {
        X=getint();Y=getint();Z=getint();
        A=getint();B=getint();C=getint();
        N=getint();
        N-=2;
        fprintf(fout,"%d\n",Solve());
    }
    fclose(fin);
    fclose(fout);
    return 0;
}