Cod sursa(job #2399843)

Utilizator Anakin1001George Giorgiu Gica Anakin1001 Data 8 aprilie 2019 09:22:50
Problema Iepuri Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.38 kb
#include <fstream>
#define LL long long
#define mod 666013
using namespace std;
ifstream f("iepuri.in");
ofstream g("iepuri.out");
int b,n,i,k,T,t,e[1001];
LL p[4][4],v[4],d[4][4],x[4][4],a[4],S;
void init(LL a[4][4]){
    for(int i=1;i<=2;i++)
        a[i+1][i]=1;
    for(int i=1;i<=3;i++)
        a[i][3]=v[i];
}
void prod(LL a[4][4],LL b[4][4],LL c[4][4]){
    for(int i=1;i<=3;i++)
        for(int j=1;j<=3;j++)
            for(int k=1;k<=3;k++)
                c[i][j]=(c[i][j]+a[i][k]*b[k][j])%mod;

}
void swop(LL a[4][4],LL b[4][4]){
    for(int i=1;i<=3;i++)
        for(int j=1;j<=3;j++){
            a[i][j]=b[i][j];
            b[i][j]=0;
        }
}
void reset(){
    for(int i=1;i<=3;i++)
    for(int j=1;j<=3;j++){
        p[i][j]=x[i][j]=d[i][j]=0;
    }
}
int main()
{   f>>T;
    for(t=1;t<=T;t++){
        f>>a[1]>>a[2]>>a[3]>>v[3]>>v[2]>>v[1];
        f>>n;
        reset();
        init(x);
        p[1][1]=p[2][2]=p[3][3]=1;
        b=n;k=-1;
        while(b!=0){
            e[++k]=b%2;
            b=b/2;
        }
        while(k>=0){
            prod(p,p,d);
            swop(p,d);
            if(e[k]==1){
                prod(p,x,d);
                swop(p,d);
            }
            k--;
        }
        S=0;
        for(i=1;i<=3;i++)
            S=(S+a[i]*p[i][1])%mod;
        g<<S<<'\n';
    }
    return 0;
}