Cod sursa(job #774811)

Utilizator tzipleatudTudor Tiplea tzipleatud Data 6 august 2012 15:21:38
Problema Iepuri Scor 90
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.55 kb
#include <fstream>
#include <cstring>
#define MOD 666013

using namespace std;

ifstream f("iepuri.in");
ofstream g("iepuri.out");


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 (f >> T; T; --T)
    {
        f >> X >> Y >> Z >> A >> B >> C >> N;
        N-=2;
        g << Solve() << '\n';
    }
    f.close();
    g.close();
    return 0;
}