Cod sursa(job #2055772)

Utilizator OlivianOlivian Dan Cretu Olivian Data 3 noiembrie 2017 18:35:52
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.97 kb
#include<cstdio>
#include<vector>
#define mod 666013
using namespace std;
int dim;
struct matrice{
    vector<vector<int> > v;
    matrice()
    {
        v.resize(dim+1,vector<int>(dim+1,0));
    }
    void unitate()
    {
        for(int i=1;i<=dim;i++) v[i][i]=1;
    }
    matrice operator *(const matrice &aux) const
    {
        matrice ret;
        for(int i=1;i<=dim;i++)
            for(int j=1;j<=dim;j++)
            for(int k=1;k<=dim;k++) ret.v[i][j]=(ret.v[i][j]+1LL*v[i][k]*aux.v[k][j])%mod;
            return ret;
    }
};
matrice rid_put(matrice m,int k)
{
    matrice ret;
    ret.unitate();
    for(;k;k>>=1)
    {
        if(k&1) ret=ret*m;
        m=m*m;
    }
    return ret;
}
int main()
{
    freopen("kfib.in","r",stdin);
    freopen("kfib.out","w",stdout);
    int k;
    scanf("%d",&k);
    dim=2;
    matrice m;
    m.v[1][2]=m.v[2][2]=m.v[2][1]=1;
    m.v[1][1]=0;
    m=rid_put(m,k-1);
    printf("%d",m.v[2][2]);

}