Pagini recente » Cod sursa (job #2830190) | Cod sursa (job #1394046) | Cod sursa (job #1568324) | Cod sursa (job #1711692) | Cod sursa (job #2855535)
#include <fstream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <queue>
#include <set>
#include <cstring>
#include <climits>
#include <iomanip>
#include <stack>
#include <cstdio>
#define MOD 666013
#define NMAX 4003
using namespace std;
int n;
int matriceInd[3][3],matConst[3][3];
int matriceInmult[3][3];
ifstream fin("expozitie.in");
ofstream fout("expozitie.out");
void inmultireMatrici(int a[][3],int b[][3],int rez[][3])
{
//initializarile
for(int i=0; i<=2; i++)
{
for(int j=0; j<=2; j++)
{
rez[i][j]=0;
}
}
for(int lin=1; lin<=2; lin++)
{
//linia din a
//inmultesc cu coloanele din b
for(int col=1; col<=2; col++)
{
for(int k=1; k<=2; k++)
{
rez[lin][col]+=1LL*(a[lin][k]*b[k][col])%MOD;
}
}
}
}
void putere(int k)
{
matriceInd[1][1]=0;
matriceInd[1][2]=1;
matriceInd[2][1]=1;
matriceInd[2][2]=1;
matConst[1][1]=0;
matConst[1][2]=1;
matConst[2][1]=1;
matConst[2][2]=1;
int aux[3][3];
while(k>0)
{
if(k%2==1)
{
inmultireMatrici(matriceInmult,matriceInd,aux);
memccpy(matriceInmult, aux, sizeof(aux));
}
inmultireMatrici(matriceInd,matConst,aux);
memccpy(matriceInd,aux,sizeof(aux));
k=k/2;
}
}
int main()
{
fin>>n;
putere(n-1);
fout<<matriceInmult[1][2]+matriceInmult[2][2];
return 0;
}