Cod sursa(job #1791572)

Utilizator RaduToporanRadu Toporan RaduToporan Data 29 octombrie 2016 14:51:42
Problema Iepuri Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.15 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,k;
    for (i=1; i<=3; i++)
        for (j=1; j<=3; j++)
        c[i][j]=0;
    for (i=1; i<=3; i++)
        for (j=1; j<=3; j++)
        for (k=1; k<=3; k++)
        c[i][j]=((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]+Y*b[2][1]+X*b[3][1]);
    }
    return 0;
}