Pagini recente » Cod sursa (job #2202936) | Cod sursa (job #285476) | Cod sursa (job #2988388)
#include <iostream>
#include <fstream>
using namespace std;
ifstream fin("kfib.in");
ofstream fout("kfib.out");
int mod=666013;
int m[2][2];
int sol[2][2];
void multiply(int a[2][2],int b[2][2]){
int c[2][2];
for(int i=0;i<2;i++){
for(int j=0;j<2;j++){
c[i][j]=0;
for(int k=0;k<2;k++)
c[i][j]+=1LL*a[i][k]*b[k][j];
c[i][j]%=mod;
}
}
for(int i=0;i<2;i++){
for(int j=0;j<2;j++){
a[i][j]=c[i][j];
}
}
}
int main(){
m[0][0]=0;
m[0][1]=m[1][0]=m[1][1]=1;
sol[0][0]=sol[1][0]=sol[1][1]=0;
sol[0][1]=1;
int p;
fin>>p;
while(p){
if(p%2){
multiply(sol,m);
}
multiply(m,m);
p/=2;
}
fout<<sol[0][0];
return 0;
}