Cod sursa(job #2408882)

Utilizator Alex18maiAlex Enache Alex18mai Data 18 aprilie 2019 13:46:53
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.71 kb
//ALEX ENACHE

#include <vector>
#include <algorithm>
#include <queue>
#include <map>
#include <set>
#include <unordered_map>
#include <time.h>
#include <iomanip>
#include <deque>
#include <math.h>
#include <cmath>
#include <assert.h>
#include <stack>
#include <bitset>
#include <random>
#include <chrono>

using namespace std;

//#include <iostream>
#include <fstream>
ifstream cin ("kfib.in");ofstream cout ("kfib.out");

const long long MOD = 666013;

long long n[2][2];
long long ans[2][2];
long long aux[2][2];

void n_la_patrat(){
    for (int i=0; i<=1; i++){
        for (int j=0; j<=1; j++){
            for (int k=0; k<=1; k++){
                aux[i][j] += n[i][k] * n[k][j];
                aux[i][j] %= MOD;
            }
        }
    }
    for (int i=0; i<=1; i++){
        for (int j=0; j<=1; j++){
            n[i][j] = aux[i][j];
            aux[i][j] = 0;
        }
    }
}

void n_ori_ans(){
    for (int i=0; i<=1; i++){
        for (int j=0; j<=1; j++){
            for (int k=0; k<=1; k++){
                aux[i][j] += ans[i][k] * n[k][j];
                aux[i][j] %= MOD;
            }
        }
    }
    for (int i=0; i<=1; i++){
        for (int j=0; j<=1; j++){
            ans[i][j] = aux[i][j];
            aux[i][j] = 0;
        }
    }
}

void lgput(int p){
    while (p){
        if (p % 2){
            n_ori_ans();
        }
        n_la_patrat();
        p /= 2;
    }
}



int main() {

    //freopen("input", "r", stdin);freopen("output", "w", stdout);

    int p;
    cin>>p;

    n[0][0] = 0;
    n[0][1] = 1;
    n[1][0] = 1;
    n[1][1] = 1;

    for (int i=0; i<=1; i++){
        for (int j=0; j<=1; j++){
            ans[i][j] = n[i][j];
        }
    }

    lgput(p);

    cout<<ans[0][0];

    return 0;
}