Pagini recente » Cod sursa (job #2685695) | Cod sursa (job #36834) | Cod sursa (job #1039347) | Cod sursa (job #392929) | Cod sursa (job #2593632)
#include <fstream>
using namespace std;
ifstream in("kfib.in");
ofstream out("kfib.out");
int mod=666013;
void produs(int a[2][2],int b[2][2]){
int aux[2][2];
for(int i=0; i<2; i++)
for(int j=0; j<2; j++)
{
aux[i][j]=0;
for(int k=0; k<=1; k++)
aux[i][j]=((aux[i][j]+1LL*a[i][k]*b[k][j])%mod);
}
for(int i=0; i<2; i++)
for(int j=0; j<2; j++)
a[i][j]=aux[i][j];
}
void putere(int a[2][2],int k){
int r[2][2]= {{1,0},{0,1}};
while(k!=0){
if(k%2==1)
produs(r,a);
produs(a,a);
k/=2;
}
for(int i=0;i<2;i++)
for(int j=0;j<2;j++)
a[i][j]=r[i][j];
}
int main(){
int k,a[2][2]={{0,1},{1,1}};
in>>k;
putere(a,k);
out<<a[1][0]%mod;
return 0;
}