Pagini recente » Cod sursa (job #631229) | Cod sursa (job #1758628) | Cod sursa (job #200509) | Cod sursa (job #576545) | Cod sursa (job #2593438)
#include<bits/stdc++.h>
using namespace std;
const int MOD=666013;
long long mat[4][4],aux[4][4];
void solve(int a, int b, int c, int n){
n-=2;
for(int i=0;i<=3;i++)
for(int j=0;j<=3;j++)
aux[i][j]=0;
aux[0][0]=aux[1][1]=aux[2][2]=1;
mat[0][0]=mat[0][1]=0,mat[0][2]=c;
mat[1][0]=1,mat[1][1]=0,mat[1][2]=b;
mat[2][0]=0,mat[2][1]=1,mat[2][2]=a;
while(n){
if(n%2){
long long aux2[4][4];
for(int i=0;i<=3;i++)
for(int j=0;j<=3;j++)
aux2[i][j]=0;
for(int i=0;i<=2;i++)
for(int j=0;j<=2;j++)
for(int k=0;k<=2;k++)
aux2[i][j]=(aux2[i][j]+aux[i][k]*mat[k][j])%MOD;
for(int i=0;i<=2;i++)
for(int j=0;j<=2;j++)
aux[i][j]=aux2[i][j];
}
n/=2;
long long aux2[4][4];
for(int i=0;i<=3;i++)
for(int j=0;j<=3;j++)
aux2[i][j]=0;
for(int i=0;i<=2;i++)
for(int j=0;j<=2;j++)
for(int k=0;k<=2;k++)
aux2[i][j]=(aux2[i][j]+mat[i][k]*mat[k][j])%MOD;
for(int i=0;i<=2;i++)
for(int j=0;j<=2;j++)
mat[i][j]=aux2[i][j];
}
}
int main(){
freopen("iepuri.in","r",stdin);
freopen("iepuri.out","w",stdout);
int t;
scanf("%d", &t);
for(int i=1;i<=t;i++){
int x,y,z,a,b,c,n;
long long v[3];
scanf("%d%d%d%d%d%d%d", &x, &y, &z, &a, &b, &c, &n);
v[0]=x,v[1]=y,v[2]=z;
solve(a, b, c, n);
for(int j=0;j<=2;j++){
long long s=0;
for(int k=0;k<=2;k++)
s=(s+v[k]*aux[k][j])%MOD;
if(j==2)
printf("%d\n", s%MOD);
}
}
return 0;
}