Pagini recente » Cod sursa (job #2852975) | Cod sursa (job #1321009) | Cod sursa (job #43795) | Cod sursa (job #817707) | Cod sursa (job #1165320)
#include <iostream>
#include <fstream>
#include <cstring>
using namespace std;
ifstream f("kfib.in");
ofstream g("kfib.out");
int rez[2][2];
void mul(int ac[2][2], int b[2][2])
{
int i,j,k;
for(i=0;i<=1;i++)
for(j=0;j<=1;j++){
rez[i][j]=0;
for(k=0;k<=1;k++)
{
rez[i][j]+=(1LL*ac[i][k]*b[k][j])%666013;
rez[i][j]%=666013;
}
}
memcpy(ac, rez, sizeof(rez));
}
void pow(int aa[2][2], int n)
{
int p[2][2];
p[0][0]=p[0][1]=p[1][0]=1;
p[1][1]=0;
while(n)
{
if(n&1)
mul(p,aa);
mul(aa,aa);
n/=2;
}
memcpy(aa, p, sizeof(aa));
}
int main()
{
int a[2][2],n;
a[0][0]=a[0][1]=a[1][0]=1;
a[1][1]=0;
f>>n;
pow(a,n-2);
g<<a[0][0]<<"\n";
return 0;
}