Cod sursa(job #2087140)

Utilizator PinkiePie1189Preoteasa Mircea-Costin PinkiePie1189 Data 12 decembrie 2017 23:20:12
Problema Iepuri Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.8 kb
#include<stdio.h>
#define MOD 666013
FILE*fin,*fout;
void inmultire(long long a[4][4],long long b[4][4]);
void lgput(int b);
void reset(long long a[4][4]);
long long X[4][4],A[4][4],rez[4][4],r[4][4];
int main()
{
    fin=fopen("iepuri.in","r");
    fout=fopen("iepuri.out","w");
    int T;
    fscanf(fin,"%d",&T);
    for(int i=1;i<=T;i++)
    {
        int x,y,z,a,b,c,n;
        fscanf(fin,"%d%d%d%d%d%d%d",&x,&y,&z,&a,&b,&c,&n);
        reset(X);
        for(int i=1;i<=2;i++)
        {
            X[i][i+1]=1;
        }
        X[3][1]=c;X[3][2]=b;X[3][3]=a;
        reset(A);
        A[1][1]=x,A[2][1]=y,A[3][1]=z;
        if(n==0)
        {
            fprintf(fout,"%d",x);
        }
        else
        {
            lgput(n);
            inmultire(X,A);
            fprintf(fout,"%lld\n",X[1][1]);
        }
    }
    fclose(fin);
    fclose(fout);
    return 0;
}
void inmultire(long long a[4][4],long long b[4][4])
{
    reset(rez);
    for(int i=1;i<=3;i++)
    {
        for(int j=1;j<=3;j++)
        {
            for(int k=1;k<=3;k++)
            {
                rez[i][j]=(rez[i][j]%MOD+((a[i][k]%MOD)*(b[k][j]%MOD))%MOD)%MOD;
            }
        }
    }
    for(int i=1;i<=3;i++)
    {
        for(int j=1;j<=3;j++)
        {
            a[i][j]=rez[i][j];
        }
    }
}
void lgput(int b)
{
    reset(r);
    for(int i=1;i<=3;i++)
    {
        r[i][i]=1;
    }
    while(b>1)
    {
        if(b%2==0)
        {
            inmultire(X,X);
        }
        else
        {
            inmultire(r,X);
            inmultire(X,X);
        }
        b/=2;
    }
    inmultire(X,r);
}
void reset(long long a[4][4])
{
    for(int i=1;i<=3;i++)
    {
        for(int j=1;j<=3;j++)
        {
            a[i][j]=0;
        }
    }
}