Pagini recente » Cod sursa (job #3334101) | Cod sursa (job #3322508) | Cod sursa (job #3315464) | Monitorul de evaluare | Cod sursa (job #3352233)
#include <fstream>
using namespace std;
long long m_inmultire[3][3];
void inmult_matrice(long long a[3][3], long long b[3][3])
{
long long rezultat[3][3] = {0};
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
for(int k = 0; k < 3; k++)
{
rezultat[i][j] += a[i][k] * b[k][j];
rezultat[i][j] = rezultat[i][j] % 666013;
}
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
b[i][j] = rezultat[i][j];
}
void exponentiere(long long m[3][3], int p)
{
long long rez[3][3];
for(int i = 0; i < 3; i ++)
for(int j = 0; j < 3; j ++)
{
if(i == j) rez[i][j] = 1;
else rez[i][j] = 0;
}
while(p > 0)
{
if(p % 2 == 1)
inmult_matrice(m, rez);
inmult_matrice(m, m);
p/=2;
}
for(int i = 0 ; i < 3; i ++)
for(int j = 0; j < 3; j ++)
m[i][j]=rez[i][j];
}
ifstream fin("iepuri.in");
ofstream fout("iepuri.out");
int main()
{
short int X, Y, Z, A, B, C, T;
int N;
fin >> T;
for(int i = 1; i <= T; i ++)
{
long long iepuri;
fin >> X >> Y >> Z >> A >> B >> C >> N;
long long m_inmultire[3][3] = { {A, B, C}, {1, 0, 0}, {0, 1, 0} };
exponentiere(m_inmultire, N-2);
iepuri = Z*m_inmultire[0][0] + Y*m_inmultire[0][1] + X*m_inmultire[0][2];
iepuri = iepuri % 666013;
fout << iepuri;
}
return 0;
}