Cod sursa(job #1643330)

Utilizator gbibBacotiu Gabi gbib Data 9 martie 2016 18:31:50
Problema Al k-lea termen Fibonacci Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.02 kb
#include <bits/stdc++.h>
#define mod 1999999973
using namespace std;
ifstream in("kfib.in");
ofstream out("kfib.out");
void inmult(long long a[2][2], long long b[2][2]) //a<-a*b
{
    long long c[2][2];
    c[0][0]=0;
    c[0][1]=0;
    c[1][0]=0;
    c[1][1]=0;
    int i,k,j;
    for(i=0;i<=1;i++)
    {
        for(j=0;j<=1;j++)
        {
            for(k=0;k<=1;k++)
            {
                c[i][j]=(c[i][j]+(a[i][k]*b[k][j])%mod)%mod;
            }
        }
    }
    a[0][0]=c[0][0];
    a[0][1]=c[0][1];
    a[1][0]=c[1][0];
    a[1][1]=c[1][1];
}

void put(long long a[2][2],long long sol[2][2], long long p)
{
    long long bit;
    for(bit=0;(1<<bit)<=p;bit++)
    {
        if((1<<bit)&p)
        {
            inmult(sol,a);
        }
        inmult(a,a);
    }
}


int main()
{long long p;
long long a[2][2];
a[1][0]=a[0][1]=a[0][0]=1;
a[1][1]=0;
long long sol[2][2];
sol[1][1]=sol[0][0]=1;
sol[1][0]=sol[0][1]=0;
in>>p;
put(a,sol,p);
out<<sol[0][0]<<" ";
    return 0;
}