Cod sursa(job #2612424)

Utilizator RaduToporanRadu Toporan RaduToporan Data 8 mai 2020 23:19:46
Problema Iepuri Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.17 kb
#include <cstdio>
 
int a[4][4],b[4][4],t,X,Y,Z,A,B,C,N;
 
void prod(int a[4][4], int b[4][4], int c[4][4])
{
    int i,j,k;
    for (i=1; i<=3; i++)
        for (j=1; j<=3; j++)
        {
            c[i][j]=0;
            for (k=1; k<=3; k++)
            c[i][j]=(1ll*((1ll*c[i][j])+(1ll*a[i][k]*b[k][j])))%666013;
        }
}
 
void atribuire(int a[4][4], int b[4][4])
{
    int i,j;
    for (i=1; i<=3; i++)
        for (j=1; j<=3; j++)
        a[i][j]=b[i][j];
}
 
void ridicare(int n, int b[4][4])
{
    int p[4][4],c[4][4];
    if (n==1) { atribuire(b,a); return; }
 
    ridicare(n/2,p);
    prod(p,p,b);
    if (n%2==1)
    {
        prod(a,b,c);
        atribuire(b,c);
    }
}
 
int main()
{
    freopen("iepuri.in","r",stdin);
    freopen("iepuri.out","w",stdout);
    scanf("%d",&t);
    for (int i=1; i<=t; i++)
    {
        scanf("%d%d%d%d%d%d%d",&X,&Y,&Z,&A,&B,&C,&N);
        a[1][1]=A; a[1][2]=1; a[1][3]=0;
        a[2][1]=B; a[2][2]=0; a[2][3]=1;
        a[3][1]=C; a[3][2]=0; a[3][3]=0;
        ridicare(N-2,b);
        printf("%d\n",((1ll*Z*b[1][1]%666013)+(1ll*Y*b[2][1])%666013+(1ll*X*b[3][1]))%666013);
    }
    return 0;
}