Pagini recente » Cod sursa (job #3258795) | Cod sursa (job #3210850) | Cod sursa (job #942205) | Cod sursa (job #1217992) | Cod sursa (job #2176007)
#include <cstdio>
using namespace std;
const int mod=666013,m=2;
struct matrice
{
int v[m][m];
matrice operator *(matrice a)
{
matrice b;
b.nula();
for(int i=0;i<m;i++)
for(int j=0;j<m;j++)
for(int k=0;k<m;k++)
b.v[i][j]=(b.v[i][j]+1LL*v[i][k]*a.v[k][j])%mod;
return b;
}
void nula()
{
for(int i=0;i<m;i++)
for(int j=0;j<m;j++) v[i][j]=0;
}
void unitate()
{
for(int i=0;i<m;i++)
for(int j=0;j<m;j++) v[i][j]=(i==j);
}
};
matrice rid_put(matrice a,int b)
{
matrice sol;
sol.unitate();
for(int i=1;i<=b;i<<=1)
{
if(b&i) sol=sol*a;
a=a*a;
}
return sol;
}
int main()
{
freopen("kfib.in","r",stdin);
freopen("kfib.out","w",stdout);
int k;
scanf("%d",&k);
if(k==0) {printf("0");return 0;}
matrice p;
p.v[0][0]=0;
p.v[0][1]=p.v[1][0]=p.v[1][1]=1;
matrice rez=rid_put(p,k-1);
printf("%d",rez.v[1][1]);
return 0;
}