Cod sursa(job #2326892)

Utilizator sebi110Ciobanu Sebastian sebi110 Data 24 ianuarie 2019 10:49:59
Problema Iepuri Scor 0
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.23 kb
#include <bits/stdc++.h>
#define MOD 666013
using namespace std;
typedef long long ll;
ll mat[4][4];
void mut(ll a[4][4],ll b[4][4])//a=b
{
    ll i,j;
    for(i=1;i<=3;i++)
        for(j=1;j<=3;j++)
            a[i][j]=b[i][j];
}
void inmult(ll a[4][4],ll b[4][4])
{
    ll rez[4][4],i,j,k;
    for(i=1;i<=3;i++)
    {
        for(j=1;j<=3;j++)
        {
            rez[i][j]=0;
            for(k=1;k<=3;k++)
            {
                rez[i][j]=(rez[i][j]+a[i][k]*b[k][j])%MOD;
            }
        }
    }
    mut(a,rez);
}
void ridic(ll k)
{
    if(k==1)
        return ;
    if(k%2==1)
    {
        ll aux[4][4];
        mut(aux,mat);
        ridic(k/2);
        inmult(mat,mat);
        inmult(mat,aux);
        return ;
    }
    ridic(k/2);
    inmult(mat,mat);
    return;
}
int main()
{
    freopen("iepuri.in","r",stdin);
    freopen("iepuri.out","r",stdout);
    ll t,x,y,z,a,b,c,n;
    cin>>t;
    while(t--)
    {
        cin>>x>>y>>z>>a>>b>>c>>n;
        mat[1][1]=mat[3][1]=mat[2][2]=mat[1][2]=0;
        mat[3][2]=mat[2][1]=1;
        mat[1][3]=c;mat[2][3]=b;mat[3][3]=a;
        ridic(n-2);
        cout<<(x*mat[1][3]+y*mat[2][3]+z*mat[3][3])%MOD<<'\n';
    }
    return 0;
}