Cod sursa(job #2502046)

Utilizator adiaioanaAdia R. adiaioana Data 30 noiembrie 2019 12:46:27
Problema Iepuri Scor 0
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.19 kb
#include <fstream>
#define MOD 666013
using namespace std;
ifstream cin("iepuri.in");
ofstream cout("iepuri.out");
int n,T,x,y,z,a,b,c,m[5][5],v[5][5],ans[5][5], rez[5][5];
void inm(int mat1[5][5], int mat2[5][5]);
void inm_mat_exp(int 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(int exp) //v*m
{
    while(exp)
    {
        if(exp%2==1)
            inm(rez,m);
        inm(m,m);
        exp/=2;
    }
}

void inm(int mat1[5][5], int mat2[5][5])
{
    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], ans[i][j]=0;
}