Pagini recente » Cod sursa (job #1670235) | Cod sursa (job #2817156) | Cod sursa (job #2738788) | Cod sursa (job #507078) | Cod sursa (job #2462319)
#include <iostream>
#include <fstream>
#define mod 666013
using namespace std;
int x,y,z,a,b,c,n,p,m[3][3],q[3][3],rez;
void ori(int m1[3][3],int m2[3][3])
{
int v1,v2,v3,v4,v5,v6,v7,v8,v9;
v1=m1[0][0]*m2[0][0]+m1[0][1]*m2[1][0]+m1[0][2]*m2[2][0];
v2=m1[0][0]*m2[0][1]+m1[0][1]*m2[1][1]+m1[0][2]*m2[2][1];
v3=m1[0][0]*m2[0][2]+m1[0][1]*m2[1][2]+m1[0][2]*m2[2][2];
v4=m1[1][0]*m2[0][0]+m1[1][1]*m2[1][0]+m1[1][2]*m2[2][0];
v5=m1[1][0]*m2[0][1]+m1[1][1]*m2[1][1]+m1[1][2]*m2[2][1];
v6=m1[1][0]*m2[0][2]+m1[1][1]*m2[1][2]+m1[1][2]*m2[2][2];
v7=m1[2][0]*m2[0][0]+m1[2][1]*m2[1][0]+m1[2][2]*m2[2][0];
v8=m1[2][0]*m2[0][1]+m1[2][1]*m2[1][1]+m1[2][2]*m2[2][1];
v9=m1[2][0]*m2[0][2]+m1[2][1]*m2[1][2]+m1[2][2]*m2[2][2];
m1[0][0]=v1;
m1[0][1]=v2;
m1[0][2]=v3;
m1[1][0]=v4;
m1[1][1]=v5;
m1[1][2]=v6;
m1[2][0]=v7;
m1[2][1]=v8;
m1[2][2]=v9;
}
void putere()
{
while(n)
{
if(n%2)
{
n--;
ori(m,q);
}
else
{
ori(m,m);
n/=2;
}
}
}
void citire()
{
ifstream fin("iepuri.in");
ofstream fout("iepuri.out");
fin>>p;
for(int i=0;i<p;i++)
{
m[0][0]=q[0][0]=0;
m[0][1]=q[0][1]=1;
m[0][2]=q[0][2]=0;
m[1][0]=q[1][0]=0;
m[1][1]=q[1][1]=0;
m[1][2]=q[1][2]=1;
fin>>x>>y>>z>>a>>b>>c>>n;
m[2][0]=q[2][0]=c;
m[2][1]=q[2][1]=b;
m[2][2]=q[2][2]=a;
n-=3;
putere();
rez=(z*m[2][2]+y*m[2][1]+x*m[2][0])%mod;
fout<<rez<<"\n";
}
}
int main()
{
citire();
return 0;
}