Cod sursa(job #2502099)

Utilizator adiaioanaAdia R. adiaioana Data 30 noiembrie 2019 13:04:08
Problema Iepuri Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.19 kb
#include <fstream>
#define ll long long
#define MOD 666013
using namespace std;
ifstream cin("iepuri.in");
ofstream cout("iepuri.out");
ll n,T,x,y,z,a,b,c,m[3][3],ans[3][3], rez[3][3];
void inm(ll mat1[3][3], ll mat2[3][3]);
void inm_mat_exp(ll exp);
int main()
{
    cin>>T;
    for(int i=1; i<=T; ++i)
    {
        cin>>x>>y>>z>>a>>b>>c>>n;

        for(int i=0; i<3; ++i)
            for(int j=0; j<3; ++j)
                rez[i][j]=m[i][j]=0;

        rez[0][0]=rez[1][1]=rez[2][2]=1;

        m[1][0]=m[2][1]=1;
        m[0][2]=c; m[1][2]=b; m[2][2]=a;

        inm_mat_exp(n-2);
        cout<<(x*rez[0][2]+y*rez[1][2]+z*rez[2][2])%MOD<<'\n';
    }
    return 0;
}

void inm_mat_exp(ll exp)
{
    while(exp)
    {
        if(exp%2)
            inm(rez,m);
        inm(m,m);
        exp/=2;
    }
}

void inm(ll mat1[3][3], ll mat2[3][3])
{
    for(int i=0; i<3; ++i)
        for(int j=0; j<3; ++j)
        {
            ans[i][j]=0;
            for(int l=0; l<3; ++l)
                ans[i][j]+=((mat1[i][l]*mat2[l][j])%MOD);
            ans[i][j]%=MOD;
        }

    for(int i=0; i<3; ++i)
        for(int j=0; j<3; ++j)
            mat1[i][j]=ans[i][j];
}