Cod sursa(job #2150129)

Utilizator TherevengerkingSurani Adrian Therevengerking Data 3 martie 2018 11:52:09
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.92 kb
#include <bits/stdc++.h>

using namespace std;
ifstream fin("kfib.in");
ofstream fout("kfib.out");
#define modl 666013
void mcp(int a[4][4], int b[4][4])
{
    for(int i = 0; i < 2; ++i)
        for(int j = 0; j < 2; ++j)
            a[i][j] = b[i][j];
}
void prod(int a[4][4], int b[4][4], int c[4][4])
{
    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] = (c[i][j] % modl + (1LL * a[i][k] * b[k][j]) % modl) % modl;
        }
}
int a[4][4]
{
    {0, 1},
    {1, 1}
};
int ret[4][4]
{
    {1, 0},
    {0, 1}
};
int c[4][4];
int main()
{
    int n;
    fin >> n;
    n--;
    for(; n; n >>= 1)
    {
        if(n & 1)
        {
            prod(a, ret, c);
            mcp(ret, c);
        }
        prod(a, a, c);
        mcp(a, c);
    }
    fout << ret[1][1];
    return 0;
}