Cod sursa(job #3254495)

Utilizator Cezar1432Dogaru Cezar Cezar1432 Data 7 noiembrie 2024 17:23:50
Problema Iepuri Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.09 kb
#include <bits/stdc++.h>
#define MOD 666013
#define ll long long
using namespace std;
ifstream f("iepuri.in");
ofstream g("iepuri.out");
ll a[3][3],c[3][3],sol[3][3];
int t,A,B,C,X,Y,Z,N;
void inmultire(ll a[3][3], ll b[3][3]){

    ll c[3][3];
    int i,j,k;
    for(i=0;i<3;i++)
    for(j=0;j<3;j++){
        c[i][j]=0;
        for(k=0;k<3;k++)
            c[i][j]+=((a[i][k]%MOD)*(b[k][j]%MOD))%MOD;
        c[i][j]%=MOD;
    }
    for(i=0;i<3;i++)
        for(j=0;j<3;j++)
        a[i][j]=c[i][j];
}
void putere(ll a[3][3], int n){
    while(n){
        if(n%2==1)inmultire(sol,c),n--;
        else inmultire(c,c),n/=2;
    }
}
void rez(){
    for(int i=0;i<3;i++)
        for(int j=0;j<3;j++)
        c[i][j]=sol[i][j]=a[i][j]=0;
    f>>X>>Y>>Z>>A>>B>>C>>N;
    a[0][0]=X,a[0][1]=Y,a[0][2]=Z;

    sol[0][0]=sol[1][1]=sol[2][2]=1;
    c[1][0]=1;
    c[2][1]=1;
    c[0][2]=C;
    c[1][2]=B;
    c[2][2]=A;
    putere(sol,N);
    inmultire(a,sol);
    g<<a[0][0]<<'\n';


}
int main()
{
    f>>t;
    for(int i=1;i<=t;i++)
    rez();


    return 0;
}