Pagini recente » Cod sursa (job #1099240) | Cod sursa (job #1750177) | Cod sursa (job #2986011) | Cod sursa (job #491051) | Cod sursa (job #2698665)
#include <fstream>
#define MOD 666013
#define ll long long
#include <string.h>
using namespace std;
ifstream f("iepuri.in");
ofstream g("iepuri.out");
ll a[3][3];
ll c[3][3];
int t,i;
ll x,y,z,A,B,C;
ll sol[3][3]= {{1,0,0},{0,1,0},{0,0,1}}; // matricea unitate
int n;
void mult(ll a[3][3],ll b[3][3])
{
ll c[3][3];
int i,j,k;
for(i=0; i<=2; i++)
{
for(j=0; j<=2; j++)
{
c[i][j]=0;
for(k=0; k<=2; k++)
c[i][j]+=a[i][k]*b[k][j];
c[i][j]=c[i][j] % MOD;
}
}
for(i=0; i<=2; i++)
for(j=0; j<=2; j++) a[i][j]=c[i][j];
}
void put(ll sol[3][3],int k)
{
while(k)
{
if(k%2==1)
{
mult(sol,c);
k--;
}
else
{
mult(c,c);
k/=2;
}
}
}
int main()
{
f>>t;
for(i=1; i<=t; i++)
{
f>>x>>y>>z>>A>>B>>C>>n;
memset(a,0,sizeof(a));
memset(c,0,sizeof(c));
memset(sol,0,sizeof(sol));
sol[0][0]=sol[1][1]=sol[2][2]=1;
a[0][0]=x; a[0][1]=y; a[0][2]=z;
c[0][2]=C; c[1][2]=B; c[2][2]=A;
c[1][0]=c[2][1]=1;
put(sol,n); // sol = C^n
mult(a,sol); // A= A*C^n
g<<a[0][0]<<'\n';
}
return 0;
}