Cod sursa(job #1679873)

Utilizator Daria09Florea Daria Daria09 Data 8 aprilie 2016 12:16:55
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");
void product(long long a[][3],long long b[][3])
{
    long long c[3][3]; c[0][0]=0;c[0][1]=0;c[1][0]=0;c[1][1]=0;
    int i,j,k;
    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;
    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[][3],long long sol[][3], long long p)
{
    long long x;for(x=0;(1<<x)<=p;x++){if((1<<x)&p)product(sol,a);product(a,a);}
}
void solve()
{
    long long a[3][3],sol[3][3],n;
    a[1][0]=a[0][1]=a[0][0]=1;a[1][1]=0;
    sol[1][1]=sol[0][0]=1;sol[1][0]=sol[0][1]=0;
    f>>n;put(a,sol,n-1);
    g<<sol[0][0];
}
int main()
{
    solve();
    return 0;
}