Pagini recente » Cod sursa (job #3267933) | Cod sursa (job #3285851) | Cod sursa (job #3267921) | Cod sursa (job #3283586) | Cod sursa (job #3284403)
#include <fstream>
#define int long long
using namespace std;
ifstream fin("kfib.in");
ofstream fout("kfib.out");
int n, i, j, mod=666013, a[3][3], b[3][3];
void inmultire(int a[3][3], int b[3][3])
{
int c[3][3];
int i, j;
for(i=1; i<=2; i++)
for(j=1; j<=2; j++)
c[i][j]=0;
for(i=1; i<=2; i++)
for(j=1; j<=2; j++)
for(int k=1; k<=2; k++)
c[i][j]=(c[i][j]+a[i][k]*b[k][j])%mod;
for(i=1; i<=2; i++)
for(j=1; j<=2; j++)
a[i][j]=c[i][j];
}
void putere(int pow)
{
while(pow)
{
if(pow%2==1)
inmultire(a, b);
inmultire(b, b);
pow/=2;
}
}
signed main()
{
fin>>n;
a[1][1]=0, a[1][2]=1, a[2][1]=1, a[2][2]=1;
b[1][1]=0, b[1][2]=1, b[2][1]=1, b[2][2]=1;
putere(n);
b[1][1]=1, b[1][2]=1;
b[2][1]=0, b[2][2]=0;
inmultire(a, b);
fout<<a[1][2];
return 0;
}