Pagini recente » Cod sursa (job #35960) | Cod sursa (job #997664) | Cod sursa (job #3180167) | Cod sursa (job #2181337) | Cod sursa (job #471506)
Cod sursa(job #471506)
#include <fstream>
#include <string.h>
using namespace std;
const char InFile[]="iepuri.in";
const char OutFile[]="iepuri.out";
const long long MOD=666013;
ifstream fin(InFile);
ofstream fout(OutFile);
long long T,N,A,B,C,X,Y,Z;
long long mat[4][4],rez[4][4],mA[4][4],mB[4][4];
void m_mul(int type)
{
long long mC[4][4];
memset(mC,0,sizeof(mC));
for(register int i=1;i<=3;++i)
{
for(register int j=1;j<=3;++j)
{
for(register int k=1;k<=3;++k)
{
mC[i][j]+=((mA[i][k]%MOD)*(mB[k][j]%MOD))%MOD;
mC[i][j]%=MOD;
}
}
}
if(type==1)
memcpy(rez,mC,sizeof(rez));
else
memcpy(mat,mC,sizeof(mat));
}
int main()
{
fin>>T;
for(register int i=0;i<T;++i)
{
fin>>X>>Y>>Z>>A>>B>>C>>N;
mat[1][1]=0;
mat[1][2]=1;
mat[1][3]=0;
mat[2][1]=0;
mat[2][2]=0;
mat[2][3]=1;
mat[3][1]=C;
mat[3][2]=B;
mat[3][3]=A;
memcpy(rez,mat,sizeof(rez));
--N;
for(;N;N>>=1)
{
if(N&1)
{
memcpy(mA,rez,sizeof(mA));
memcpy(mB,mat,sizeof(mB));
m_mul(1);
}
memcpy(mA,mat,sizeof(mA));
memcpy(mB,mat,sizeof(mB));
m_mul(2);
}
fout<<(((X*rez[1][1])%MOD+(Y*rez[1][2])%MOD)%MOD+(Z*rez[1][3])%MOD)%MOD<<"\n";
}
fin.close();
fout.close();
return 0;
}