bezout.c

2个数和多个数求最大公约数及 bezout 数组

#include <stdlib.h>
#include <stdio.h>

/* Ce fichier contient la procedure qui  calcule les parametre de bezout*/
int bezout(int a,int b, int *u, int *v)
{

  int c,q, r, f, i;
  int u0, v0;
    *u = 1;
    *v = 0; 	
    u0 = 0;
    v0 = 1;
/*on fait deux tests pour a==0 soit b==0*/
     
 if (a ==  0 )
  {
  *u = 1;
  *v = 0;
  return b;
  } 
  if (b ==  0 )
  {
  *u = 0;
  *v = 1;
  return a;
  } 
 
/* si  a inferieur  b,on change le position de a et b */
     
    if(a <b) 
    { 
        c = b;
        b = a; 
        a = c; 
    } 
 /* q est le quotient de a/b,r est le rest de a/b*/
   if ( a != 0 && b != 0)

   { 
    q = a/b;
    r = a%b;
  	
      
/* On fait une boucle pour calculer le q et r while  le reste  est different de zero. En suite on caulcule les parametres  de bezout. */   
   
    while (r!=0) 
      {    
	  f = *u - q * *v;
          i = u0 - q*v0;
         *u = *v;   
	 *v = f;
  	 u0 = v0 ;  
	 v0 = i;
          a = b; 
	  b = r;      
  	  r = a%b;
	  q = a/b ;        
     }

/*on fiche les parametres de bezout u et v*/ 
     printf("La  valeur du couple u et v est : (%3d , %3d) \n" ,f,i);	  
   
    }
  return b; /*b est le pgcd de a et b*/

}