Pagini recente » Cod sursa (job #2720530) | Cod sursa (job #990751) | Cod sursa (job #317324) | Cod sursa (job #1268935) | Cod sursa (job #1385712)
#include <iostream>
#include <fstream>
using namespace std;
ifstream f("kfib.in");
ofstream g("kfib.out");
const int mod=666013;
int k,a[3][3],b[3][3],c[3][3];
void citire(){f>>k;}
void inm(int x[3][3],int y[3][3],int z[3][3])
{
for(int i=1;i<3;i++)
for(int j=1;j<3;j++)
for(int h=1;h<3;h++)
z[i][j]=(z[i][j]+x[i][h]*y[h][j])%mod;
}
void afisare(int x[3][3])
{
for(int i=1;i<3;i++)
{for(int j=1;j<3;j++)g<<x[i][j]<<" ";
g<<endl;}
}
void putere(int p)
{ int i;
for(i=0;(1<<i)<=p;i++)
{if((1<<i)&p)
{c[1][1]=0;c[1][2]=0;c[2][1]=0;c[2][2]=0;
inm(a,b,c);
a[1][1]=c[1][1];a[1][2]=c[1][2];a[2][1]=c[2][1];a[2][2]=c[2][2];
}
c[1][1]=0;c[1][2]=0;c[2][1]=0;c[2][2]=0;
inm(b,b,c);
b[1][1]=c[1][1];b[1][2]=c[1][2];b[2][1]=c[2][1];b[2][2]=c[2][2];
}
}
int main()
{
citire();
a[1][1]=1;
a[2][2]=1;
b[1][1]=0;
b[1][2]=1;
b[2][1]=1;
b[2][2]=1;
putere(k-1);
g<<a[2][2];
return 0;
}