Pagini recente » Cod sursa (job #1179538) | Cod sursa (job #1139622) | Cod sursa (job #1958415) | Cod sursa (job #1016852) | Cod sursa (job #1007230)
#include <fstream>
#include <cstring>
#define mod 666013
#define ll long long
using namespace std;
ifstream fin("kfib.in");
ofstream fout("kfib.out");
ll k;
ll u[2][2]={0,1,1,1};
void lgpow_matrix (ll m[][2], ll k)
{
if (k==1)
{
memcpy (m,u,sizeof(u));
return;
}
ll temp[2][2];
memset (temp,0,sizeof(temp));
lgpow_matrix (temp,k/2);
for (int i=0; i<2; ++i)
for (int j=0; j<2; ++j)
for (int k=0; k<2; ++k)
{
m[i][j] = (m[i][j] + (temp[i][k]*temp[k][j])%mod ) %mod;
}
if (k%2)
{
ll temp = m[0][1];
m[0][1] = (m[0][1] + m[0][0])%mod;
m[0][0] = temp;
temp = m[1][1];
m[1][1] = (m[1][1] + m[1][0])%mod;
m[1][0] = temp;
}
}
int main()
{
fin>>k;
ll m[2][2];
memset (m,0,sizeof(m));
lgpow_matrix (m,k-1);
fout<<m[1][1];
}