Cod sursa(job #2612422)

Utilizator RaduToporanRadu Toporan RaduToporan Data 8 mai 2020 23:17:58
Problema Iepuri Scor 60
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.16 kb
#include <cstdio>

int a[4][4],b[4][4],i,j,k,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;
    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 (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",((Z*b[1][1]%666013)+(Y*b[2][1])%666013+(X*b[3][1]))%666013);
    }
    return 0;
}