Cod sursa(job #3197582)

Utilizator sildegilici silca sildeg Data 27 ianuarie 2024 10:23:31
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.26 kb
#include <fstream>
using namespace std;
ifstream fin("kfib.in");
ofstream fout("kfib.out");
const int MOD = 666013;
int Sol[2][2],M[2][2];

void Multiply2(int A[1][2],int B[2][2])
{
    long long C[1][2];
    for(int i = 0; i < 1; ++i)
        for(int j = 0; j < 2; ++j)
            {
                C[i][j] = 0;
                for(int k = 0; k < 2; ++k)
                    C[i][j] += 1LL*A[i][k]*B[k][j];
                C[i][j] %= MOD;
            }
    for(int i = 0; i < 1; ++i)
        for(int j = 0; j < 2; ++j)
            A[i][j] = C[i][j];
}

void Multiply(int A[2][2],int B[2][2])
{
    long long C[2][2];
    for(int i = 0; i < 2; ++i)
        for(int j = 0; j < 2; ++j)
            {
                C[i][j] = 0;
                for(int k = 0; k < 2; ++k)
                    C[i][j] += 1LL*A[i][k]*B[k][j];
                C[i][j] %= MOD;
            }
    for(int i = 0; i < 2; ++i)
        for(int j = 0; j < 2; ++j)
            A[i][j] = C[i][j];
}

int main()
{
    int K;
    fin >> K;

    Sol[0][1] = 1;

    M[0][1] = M[1][0] = M[1][1] = 1;

    while(K)
    {
        if(K%2 == 1)
            Multiply2(Sol,M);
        Multiply(M,M);
        K = K/2;
    }

    fout << Sol[0][0]<< "\n";

    return 0;
}