Cod sursa(job #2492043)

Utilizator mirceaisherebina mircea mirceaishere Data 13 noiembrie 2019 20:59:46
Problema Iepuri Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.3 kb
#include <fstream>
#include <iostream>
using namespace std;
ifstream fin("iepuri.in");
ofstream fout("iepuri.out");

int n, p, i, j, t, x, y, z, a, b, c;
long long f[3][3], constanta[3][3];


void InmultireMatrice(long long a[][3], long long b[][3]){
    long long c[3][3];
    for(int i=0; i<3; i++){
        for(int j=0; j<3; j++){
            c[i][j]=0;
            for(int k=0; k<3; k++)
                c[i][j]+=a[i][k]*b[k][j]%666013;
        }
    }
    for(int i=0; i<3; i++){
        for(int j=0; j<3; j++){
            a[i][j]=c[i][j]%666013;
        }
    }
}


int main(){
    fin>>t;
    for(; t>0; t--){
        fin>>x>>y>>z>>a>>b>>c>>p;
        constanta[0][0]=a;
        constanta[0][1]=b;
        constanta[0][2]=c;
        constanta[1][0]=constanta[2][1]=1;
        constanta[1][1]=constanta[1][2]=constanta[2][0]=constanta[2][2]=0;
        f[0][0]=a;
        f[0][1]=b;
        f[0][2]=c;
        f[1][0]=f[2][1]=1;
        f[1][1]=f[1][2]=f[2][0]=f[2][2]=0;
        p-=1;
        while(p!=0){
            if(p%2==1){
                InmultireMatrice(f, constanta);
            }
            InmultireMatrice(constanta, constanta);
            p=p/2;
        }
        fout<<( (f[2][0]*z)%666013 + (f[2][1]*y)%666013 + (f[2][2]*x)%666013 )%666013<<"\n";
    }
}