Cod sursa(job #1007230)

Utilizator Impaler_009Mihai Nitu Impaler_009 Data 8 octombrie 2013 16:40:30
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.93 kb
#include <fstream>
#include <cstring>

#define mod 666013
#define ll long long

using namespace std;

ifstream fin("kfib.in");
ofstream fout("kfib.out");

ll k;

ll u[2][2]={0,1,1,1};

void lgpow_matrix (ll m[][2], ll k)
{
    if (k==1)
    {
        memcpy (m,u,sizeof(u));
        return;
    }

    ll temp[2][2];
    memset (temp,0,sizeof(temp));

    lgpow_matrix (temp,k/2);

    for (int i=0; i<2; ++i)
      for (int j=0; j<2; ++j)
        for (int k=0; k<2; ++k)
        {
            m[i][j] = (m[i][j] + (temp[i][k]*temp[k][j])%mod ) %mod;
        }

    if (k%2)
    {
        ll temp = m[0][1];
        m[0][1] = (m[0][1] + m[0][0])%mod;
        m[0][0] = temp;

        temp = m[1][1];
        m[1][1] = (m[1][1] + m[1][0])%mod;
        m[1][0] = temp;
    }
}

int main()
{
    fin>>k;

    ll m[2][2];
    memset (m,0,sizeof(m));

    lgpow_matrix (m,k-1);

    fout<<m[1][1];
}