Cod sursa(job #1140321)

Utilizator span7aRazvan span7a Data 11 martie 2014 21:52:30
Problema Al k-lea termen Fibonacci Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 0.97 kb
#include<fstream>
#include<algorithm>
#include<cstring>
#define MOD 666013
using namespace std;
ifstream f("kfib.in");
ofstream g("kfib.out");
int sol[3][3],aux[3][3],k,auxiliar[3][3];
void inmulteste(int a[3][3],int b[3][3],int c[3][3])
{
    c[1][1]=a[1][1]*b[1][1]%MOD+a[1][2]*b[2][1]%MOD;
    c[1][2]=a[1][1]*b[1][2]%MOD+a[1][2]*b[2][2]%MOD;
    c[2][1]=a[2][1]*b[1][1]%MOD+a[2][2]*b[2][1]%MOD;
    c[2][2]=a[2][1]*b[1][2]%MOD+a[2][2]*b[2][2]%MOD;
}
void putere()
{
    sol[1][1]=sol[2][2]=1;
    aux[1][1]=aux[1][2]=aux[2][1]=1;
    while(k)
    {
        if(k%2==1)
        {
            memset(auxiliar,0,sizeof(auxiliar));
            inmulteste(sol,aux,auxiliar);
            memcpy(sol,auxiliar,sizeof(auxiliar));
        }
        memset(auxiliar,0,sizeof(auxiliar));
        inmulteste(aux,aux,auxiliar);k/=2;
        memcpy(aux,auxiliar,sizeof(auxiliar));
    }
        g<<sol[1][2]<<'\n';
}
int main()
{
    f>>k;
    putere();
    return 0;
}