Cod sursa(job #1822366)

Utilizator GeoeyMexicanuBadita George GeoeyMexicanu Data 4 decembrie 2016 19:47:27
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.8 kb
#include <iostream>
#include <fstream>
#define MOD 666013

using namespace std;

ifstream f("kfib.in");
ofstream g("kfib.out");

int n;
long long a[4][4],s[4][4];
void inmultire(long long b[4][4],long long c[4][4])
{
    int i,j,k;
    long long d[4][4];
    for(i=1;i<=2;i++)
        for(j=1;j<=2;j++)
    {
        d[i][j]=0LL;
        for(k=1;k<=2;k++)
            d[i][j]=(0LL+d[i][j]+1LL*b[i][k]*c[k][j])%MOD;
    }
    for(i=1;i<=2;i++)
        for(j=1;j<=2;j++)
        b[i][j]=d[i][j];
}
int main()
{
    f>>n;
    a[1][2]=1;
    a[2][1]=1;
    a[2][2]=1;
    s[1][1]=0;
    s[1][2]=1;
    n--;
    while(n>0)
    {
        if(n%2!=0)
        {
            inmultire(s,a);
        }
        inmultire(a,a);
        n=n/2;
    }
    g<<s[1][2];
    return 0;
}