Pagini recente » Cod sursa (job #2452857) | Cod sursa (job #1140322) | Cod sursa (job #589958) | Cod sursa (job #2445098) | Cod sursa (job #2282885)
#include <fstream>
using namespace std;
ifstream fin("kfib.in");
ofstream fout("kfib.out");
const int MOD = 666013;
int K;
int Z[2][2],Sol[2][2];
void Multiply(int A[2][2], int B[2][2])
{
int C[2][2];
for(int i = 0; i <= 1; ++i)
for(int j = 0; j <= 1; ++j)
{
C[i][j] = 0;
for(int k = 0; k <= 1; ++k)
C[i][j] = (C[i][j] + (1LL * A[i][k]*B[k][j]))%MOD;
}
for(int i = 0; i <= 1; ++i)
for(int j = 0; j <= 1; ++j)
A[i][j] = C[i][j];
}
int main()
{
Z[0][1] = Z[1][0] = Z[1][1] = 1;
Sol[0][1] = 1;
fin >> K;
while(K)
{
if(K % 2 == 1)
Multiply(Sol,Z);
Multiply(Z,Z);
K = K / 2;
}
fout << Sol[0][0];
return 0;
}