Pagini recente » Cod sursa (job #987439) | Cod sursa (job #3136036) | Cod sursa (job #2848225) | Cod sursa (job #1361019) | Cod sursa (job #2711453)
#include <bits/stdc++.h>
using namespace std;
ifstream f("kfib.in");
ofstream g("kfib.out");
const int MOD=666013;
int n,a[3][3],sol[3][3];
void initialization()
{
a[0][0]=0;
a[0][1]=a[1][0]=a[1][1]=1;
sol[1][0]=sol[1][1]=0;
sol[0][0]=sol[0][1]=1;
}
void multiplication(int a[3][3],int b[3][3])
{
int ans[3][3];
memset(ans,0,sizeof(ans));
for(int i=0; i<2; i++)
for(int j=0; j<2; j++)
for(int k=0; k<2; k++)
ans[i][j]=(ans[i][j]+1LL*a[i][k]*b[k][j])%MOD;
memcpy(a,ans,sizeof(ans));
}
void solve(int p)
{
while(p)
{
if(p%2)
multiplication(sol,a);
p/=2;
multiplication(a,a);
}
}
int main()
{
f>>n;
initialization();
solve(n-1);
g<<sol[0][0];
return 0;
}