Cod sursa(job #471506)

Utilizator ChallengeMurtaza Alexandru Challenge Data 19 iulie 2010 11:50:08
Problema Iepuri Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.26 kb
#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;
}