Cod sursa(job #2046353)

Utilizator daniela12Sandu Daniela Teodora daniela12 Data 23 octombrie 2017 18:40:19
Problema Iepuri Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.19 kb
#include <iostream>
#include <fstream>
# define zi 666013
using namespace std;
ifstream f("iepuri.in");
ofstream g("iepuri.out");
long long int m[4][4], h[4][4];
void matrice (int a, int b, int c)
{
    for(int i=1;i<=3;++i)
        for(int j=1;j<=3;++j)
            m[i][j]=h[i][j]=0;;
    m[1][2]=m[2][3]=h[1][2]=h[2][3]=1;
    m[3][1]=h[3][1]=c;
    m[3][2]=h[3][2]=b;
    m[3][3]=h[3][3]=a;
}
void inm (long long int a[4][4], long long int b[4][4])
{
    int i, j, k, c[4][4]={0};
    for(i=1;i<=3;++i)
        for(j=1;j<=3;++j)
            for(k=1;k<=3;++k)
                c[i][j]=a[i][k]*b[k][j]+c[i][j];
    for(i=1;i<=3;++i)
        for(j=1;j<=3;++j)
            a[i][j]=c[i][j];
}
void putere (int n)
{   int p=1;
   while (p<n)
     if(p*2<=n)
           inm(m, m),p*=2;
        else
            inm(m,h), p++;
}
void putere2 (int n)
{
    while(n>1)
        inm(m, h), n--;
}
void citire ()
{
    int T,X,Y,Z,A,B,C,N;
    f>>T;
    while(T--)
    {
        f>>X>>Y>>Z>>A>>B>>C>>N;
        matrice(A,B,C);
        putere(N);
        g<<(m[1][1]*X+m[1][2]*Y+m[1][3]*Z)%zi<<"\n";
    }
    f.close();
    g.close();
}
int main()
{
    citire();
}