Cod sursa(job #387188)

Utilizator Bogdan_CCebere Bogdan Bogdan_C Data 26 ianuarie 2010 23:11:10
Problema Al k-lea termen Fibonacci Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.8 kb
#include<cstdio>
#include<cstdlib>
#include<cstring>
#define FIN "kfib.in"
#define FOUT "kfib.out"
#define MOD 666013
using namespace std;
int id[2][2];
void multiply(int a[][2],int b[][2])
 {int rez[2][2];
  memset(rez,0,sizeof(rez));
    for(int i=0;i<2;i++)
      for(int j=0;j<2;j++)
        for(int k=0;k<2;k++)
          rez[i][j] =(rez[i][j]+(1LL * a[i][k]*b[k][j])%MOD)%MOD;
   memcpy(a,rez,sizeof(rez));  }
void exp(int a[][2],int k)
 {for(;k;k>>=1)
   {if(k&1) multiply(id,a);
     multiply(a,a);     
          }
printf("%d",id[1][0]);        
     }    
int main()
{freopen(FIN,"r",stdin);
freopen(FOUT,"w",stdout);
int k;
scanf("%d",&k);
int a[2][2];
a[0][0]=0;a[0][1]=1;a[1][0]=1;a[1][1]=1;
id[0][0]=id[1][1]=1;id[0][1] = id[1][0] = 0;
exp(a,k);  

    return 0;};