Cod sursa(job #1796975)

Utilizator ajeccAjechiloae Eugen ajecc Data 3 noiembrie 2016 21:55:02
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 2.33 kb
#include <bits/stdc++.h>
#define for0(i, n) for(int i = 0; i < n; i++)
#define for1(i, n) for(int i = 1; i <= n; i++)
#define pb push_back
#define mp make_pair
#define all(v) v.begin(), v.end()
#define V vector<int>
#define VP vector<pair<int, int> >
#define clr(A, x) memset(A, x, sizeof(A))
#define cpy(A, B) memcpy(A, B, sizeof(B))
#define g(s) getline(cin, s) ///ai grija la fin/cin ///
#define FASTIO ios_base::sync_with_stdio(0)
const long long INFLL = (1LL<<62);
const int INFINT = 2000000000;
typedef long long ll;
typedef unsigned long long ull;
using namespace std;
/*template <typename T>
string to_string(const T& n){
    ostringstream os;
    os << n;
    return os.str();
}*/

/*void invers_modular(int a, int b, int &d, int &x, int &y)
{
    if(!b)
    {
        d=a;
        x=1;
        y=0;
        return ;
    }
    int x0, y0;
    invers_modular(b, a%b, d, x0, y0);
    x=y0;
    y=x0-a/b*y0;
}*/ // daca x<0 se aduna cu mod pana e mai mare, x fiind rezultatul

/*ull putere(ull baza, ull exponent, ull MOD)
{
    if(exponent == 0) return 1;
    if(exponent % 2 == 0) return putere((baza * baza) % MOD, exponent / 2, MOD) % MOD;
    return ((baza % MOD) * (putere(baza, exponent - 1, MOD) % MOD) % MOD);
}*/
ifstream fin("kfib.in"); /// modifica cu numele corespunzator
ofstream fout("kfib.out"); /// modifica cu numele corespunzator

const int MOD = 666013;

ll fibo[4], aux[4];

void putere(ll fibo[], int n)
{
    if(n == 1)
        fibo[1] = fibo[2] = fibo[3] = 1;
    else
    {
        if(n % 2)
        {
            putere(fibo, n - 1);
            aux[0] = fibo[1];
            aux[1] = (fibo[0] + fibo[1]) % MOD;
            aux[2] = fibo[3];
            aux[3] = (fibo[2] + fibo[3]) % MOD;
            for0(i, 4)
                fibo[i] = aux[i];
        }
        else
        {
            putere(fibo, n / 2);
            aux[0] = (fibo[0] * fibo[0] + fibo[1] * fibo[2]) % MOD;
            aux[1] = (fibo[0] * fibo[1] + fibo[1] * fibo[3]) % MOD;
            aux[2] = (fibo[0] * fibo[2] + fibo[2] * fibo[3]) % MOD;
            aux[3] = (fibo[1] * fibo[2] + fibo[3] * fibo[3]) % MOD;
            for0(i, 4)
                fibo[i] = aux[i];
        }
    }
}

int main()
{
    int k;
    fin >> k;
    putere(fibo, k + 1);
    fout << fibo[0];

    return 0;
}