Cod sursa(job #2739732)

Utilizator T_george_TGeorge Teodorescu T_george_T Data 9 aprilie 2021 20:12:35
Problema Iepuri Scor 90
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.14 kb
#include <bits/stdc++.h>

using namespace std;
ifstream in("iepuri.in");
ofstream out("iepuri.out");
const int MOD=666013;
long long p[4][4],r[4][4],rez[4][4];
void mult(long long a[][4],long long b[4][4]){
    int c[4][4];
    memset(c,0,sizeof(c));

    for(int i=1;i<=3;i++)
    for(int j=1;j<=3;j++){
    for(int k=1;k<=3;k++)
        c[i][j]=(c[i][j]+(a[i][k]*b[k][j])%MOD)%MOD;
    }
    for(int i=1;i<=3;i++)
    for(int j=1;j<=3;j++)
       a[i][j]=c[i][j];
}
void logp(long long n){

    for(long long i=0;(1<<i)<=n;i++){
         if((1<<i)&n)
            mult(r,p);
         mult(p,p);
    }
     for(int i=1;i<=3;i++)
        for(int j=1;j<=3;j++)
        p[i][j]=r[i][j];
}
int main()
{
    long long t,x,y,z,a,b,c,n;
    for(in>>t;t>0;t--){
        in>>x>>y>>z>>a>>b>>c>>n;
        memset(r,0,sizeof(r));
        for(int i=1;i<=3;i++)
            r[i][i]=1;
        p[1][1]=0,p[1][2]=0,p[1][3]=c;
        p[2][1]=1,p[2][2]=0,p[2][3]=b;
        p[3][1]=0,p[3][2]=1,p[3][3]=a;
        logp(n);
        rez[1][1]=x,rez[1][2]=y,rez[1][3]=z;
        mult(rez,r);
        out<<rez[1][1]<<'\n';
    }
    return 0;
}