LINUX.ORG.RU

История изменений

Исправление qulinxao, (текущая версия) :

#include <stdio.h>
 
void multiply(int F[2][2], int M[2][2]);
 
void power(int F[2][2], int n);
 
int fib(int n){
  if (n == 0)
    return 0;
  int F[2][2] = {{1,0},{0,1}};
  power(F, n);
  return F[0][0];
}

int E[2][2] = {{1,1},{1,0}};

void power(int F[2][2], int n){
  if( n == 0)
      return;
  
  if( n > 1){
      power(F, n/2);
      multiply(F, F);
  }
  
  if( n%2 )
     multiply(F, E);
}
 
void multiply(int F[2][2], int M[2][2]){
  int x =  F[0][0]*M[0][0] + F[0][1]*M[1][0];
  int y =  F[0][0]*M[0][1] + F[0][1]*M[1][1];
  int z =  F[1][0]*M[0][0] + F[1][1]*M[1][0];
  int w =  F[1][0]*M[0][1] + F[1][1]*M[1][1];
 
  F[0][0] = x;
  F[0][1] = y;
  F[1][0] = z;
  F[1][1] = w;
}
 
int main(){
  int n = 9;
  printf("%d", fib(9));
  getchar();
  return 0;
}

Исправление qulinxao, :

#include <stdio.h>
 
void multiply(int F[2][2], int M[2][2]);
 
void power(int F[2][2], int n);
 
int fib(int n){
  if (n == 0)
    return 0;
  int F[2][2] = {{1,0},{0,1}};
  power(F, n);
  return F[0][0];
}

int E[2][2] = {{1,1},{1,0}};

void power(int F[2][2], int n){
  if( n == 0)
      return;
  power(F, n/2);
  multiply(F, F);
  if( n%2)
     multiply(F, E);
}
 
void multiply(int F[2][2], int M[2][2]){
  int x =  F[0][0]*M[0][0] + F[0][1]*M[1][0];
  int y =  F[0][0]*M[0][1] + F[0][1]*M[1][1];
  int z =  F[1][0]*M[0][0] + F[1][1]*M[1][0];
  int w =  F[1][0]*M[0][1] + F[1][1]*M[1][1];
 
  F[0][0] = x;
  F[0][1] = y;
  F[1][0] = z;
  F[1][1] = w;
}
 
int main(){
  int n = 9;
  printf("%d", fib(9));
  getchar();
  return 0;
}

Исходная версия qulinxao, :

с просторов интернета:

#include <stdio.h>
 
void multiply(int F[2][2], int M[2][2]);
 
void power(int F[2][2], int n);
 
int fib(int n){
  if (n == 0)
    return 0;
  int F[2][2] = {{1,1},{1,0}};
  power(F, n-1);
  return F[0][0];
}

int E[2][2] = {{1,1},{1,0}};

void power(int F[2][2], int n){
  if( n == 0)
      return;
  power(F, n/2);
  multiply(F, F);
  if( n%2)
     multiply(F, E);
}
 
void multiply(int F[2][2], int M[2][2]){
  int x =  F[0][0]*M[0][0] + F[0][1]*M[1][0];
  int y =  F[0][0]*M[0][1] + F[0][1]*M[1][1];
  int z =  F[1][0]*M[0][0] + F[1][1]*M[1][0];
  int w =  F[1][0]*M[0][1] + F[1][1]*M[1][1];
 
  F[0][0] = x;
  F[0][1] = y;
  F[1][0] = z;
  F[1][1] = w;
}
 
int main(){
  int n = 9;
  printf("%d", fib(9));
  getchar();
  return 0;
}