Pagini recente » Borderou de evaluare (job #1344170) | Cod sursa (job #861223) | Cod sursa (job #2191222) | Cod sursa (job #3248132) | Cod sursa (job #3231110)
#include <fstream>
using namespace std;
ifstream cin("kfib.in");
ofstream cout("kfib.out");
long long p,ans[1][2],m[2][2],MOD=666013;
void inmultire(long long a[2][2],long long b[2][2])
{
long long c[2][2];
for(int i=0;i<=1;i++)
for(int j=0;j<=1;j++)
{
c[i][j]=0;
for(int k=0;k<=1;k++)
c[i][j]=(c[i][j]+1LL*a[i][k]*b[k][j])%MOD;
}
for(int i=0;i<=1;i++)
for(int j=0;j<=1;j++)
a[i][j]=c[i][j];
}
void inmultire2(long long rez[1][2],long long a[2][2])
{
long long c[1][2];
int i=0;
for(int j=0;j<=1;j++)
{
c[i][j]=0;
for(int k=0;k<=1;k++)
c[i][j]=(c[i][j]+1LL*rez[i][k]*a[k][j])%MOD;
}
i=0;
for(int j=0;j<=1;j++)
rez[i][j]=c[i][j];
}
int main(void)
{
ans[0][1]=1;
m[0][1]=1;
m[1][0]=1;
m[1][1]=1;
cin>>p;
while(p)
{
if(p%2==1)
inmultire2(ans,m);
inmultire(m,m);
p=p/2;
}
cout<<ans[0][0];
return 0;
}