Cod sursa(job #2462319)

Utilizator elenaisaiaElena Isaia elenaisaia Data 27 septembrie 2019 09:15:09
Problema Iepuri Scor 0
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.7 kb
#include <iostream>
#include <fstream>
#define mod 666013

using namespace std;

int x,y,z,a,b,c,n,p,m[3][3],q[3][3],rez;

void ori(int m1[3][3],int m2[3][3])
{
    int v1,v2,v3,v4,v5,v6,v7,v8,v9;
    v1=m1[0][0]*m2[0][0]+m1[0][1]*m2[1][0]+m1[0][2]*m2[2][0];
    v2=m1[0][0]*m2[0][1]+m1[0][1]*m2[1][1]+m1[0][2]*m2[2][1];
    v3=m1[0][0]*m2[0][2]+m1[0][1]*m2[1][2]+m1[0][2]*m2[2][2];
    v4=m1[1][0]*m2[0][0]+m1[1][1]*m2[1][0]+m1[1][2]*m2[2][0];
    v5=m1[1][0]*m2[0][1]+m1[1][1]*m2[1][1]+m1[1][2]*m2[2][1];
    v6=m1[1][0]*m2[0][2]+m1[1][1]*m2[1][2]+m1[1][2]*m2[2][2];
    v7=m1[2][0]*m2[0][0]+m1[2][1]*m2[1][0]+m1[2][2]*m2[2][0];
    v8=m1[2][0]*m2[0][1]+m1[2][1]*m2[1][1]+m1[2][2]*m2[2][1];
    v9=m1[2][0]*m2[0][2]+m1[2][1]*m2[1][2]+m1[2][2]*m2[2][2];
    m1[0][0]=v1;
    m1[0][1]=v2;
    m1[0][2]=v3;
    m1[1][0]=v4;
    m1[1][1]=v5;
    m1[1][2]=v6;
    m1[2][0]=v7;
    m1[2][1]=v8;
    m1[2][2]=v9;
}

void putere()
{
    while(n)
    {
        if(n%2)
        {
            n--;
            ori(m,q);
        }
        else
        {
            ori(m,m);
            n/=2;
        }
    }
}

void citire()
{
    ifstream fin("iepuri.in");
    ofstream fout("iepuri.out");
    fin>>p;
    for(int i=0;i<p;i++)
    {
        m[0][0]=q[0][0]=0;
        m[0][1]=q[0][1]=1;
        m[0][2]=q[0][2]=0;
        m[1][0]=q[1][0]=0;
        m[1][1]=q[1][1]=0;
        m[1][2]=q[1][2]=1;
        fin>>x>>y>>z>>a>>b>>c>>n;
        m[2][0]=q[2][0]=c;
        m[2][1]=q[2][1]=b;
        m[2][2]=q[2][2]=a;
        n-=3;
        putere();
        rez=(z*m[2][2]+y*m[2][1]+x*m[2][0])%mod;
        fout<<rez<<"\n";
    }
}

int main()
{
    citire();
    return 0;
}