Pagini recente » Cod sursa (job #851704) | Cod sursa (job #2496666) | Cod sursa (job #254699) | Cod sursa (job #400492) | Cod sursa (job #2809382)
#include <bits/stdc++.h>
#define MOD 666013
using namespace std;
ifstream fin("kfib.in");
ofstream fout("kfib.out");
int n,i,j,mat[5][5],c[5][5],aux[5][5];
void inmultire(int a[][5],int b[][5],int c[][5])
{
int i,j,k;
for(i=0;i<2;i++)
for(j=0;j<2;j++)
for(k=0;k<2;k++)
c[i][j]=(c[i][j]+a[i][k]*b[k][j]%MOD)%MOD;
}
int main()
{
fin>>n;
mat[0][0]=mat[1][1]=1;
c[0][1]=c[1][0]=c[1][1]=1;
while(n)
{
if(n%2==1)
{
memset(aux,0,sizeof(aux));
inmultire(mat,c,aux);
memcpy(mat,aux,sizeof(aux));
}
memset(aux,0,sizeof(aux));
inmultire(c,c,aux);
memcpy(c,aux,sizeof(aux));
n/=2;
}
fout<<mat[0][1];
return 0;
}