Pagini recente » Cod sursa (job #1289998) | Cod sursa (job #1441933) | Cod sursa (job #3184795) | Cod sursa (job #2875435) | Cod sursa (job #1140321)
#include<fstream>
#include<algorithm>
#include<cstring>
#define MOD 666013
using namespace std;
ifstream f("kfib.in");
ofstream g("kfib.out");
int sol[3][3],aux[3][3],k,auxiliar[3][3];
void inmulteste(int a[3][3],int b[3][3],int c[3][3])
{
c[1][1]=a[1][1]*b[1][1]%MOD+a[1][2]*b[2][1]%MOD;
c[1][2]=a[1][1]*b[1][2]%MOD+a[1][2]*b[2][2]%MOD;
c[2][1]=a[2][1]*b[1][1]%MOD+a[2][2]*b[2][1]%MOD;
c[2][2]=a[2][1]*b[1][2]%MOD+a[2][2]*b[2][2]%MOD;
}
void putere()
{
sol[1][1]=sol[2][2]=1;
aux[1][1]=aux[1][2]=aux[2][1]=1;
while(k)
{
if(k%2==1)
{
memset(auxiliar,0,sizeof(auxiliar));
inmulteste(sol,aux,auxiliar);
memcpy(sol,auxiliar,sizeof(auxiliar));
}
memset(auxiliar,0,sizeof(auxiliar));
inmulteste(aux,aux,auxiliar);k/=2;
memcpy(aux,auxiliar,sizeof(auxiliar));
}
g<<sol[1][2]<<'\n';
}
int main()
{
f>>k;
putere();
return 0;
}