Erreur dans mon examen de TP, fonction PGCD

Erreur dans mon examen de TP, fonction PGCD - C - Programmation

Marsh Posté le 15-01-2004 à 16:59:44    

Code :
  1. #include <stdio.h>
  2. #include <stdlib.h>
  3. #include <string.h>
  4. #include <time.h>
  6. #define MAX 100 /* taille maximale des grands entiers manipulA~(c)s */
  7. #define BASE 10000
  9. typedef unsigned int CHIFFRE;
  10. typedef CHIFFRE BIGNUM[MAX];
  13. void affiche(BIGNUM z)
  14. {
  15. /* affiche a l'ecran le grand entier z */
  16.   int k;
  18. /* on cherche le chiffre le plus significatif */
  19.   k=MAX-1;
  20.   while ((k >=0) && (z[k] == 0))
  21.     k--;
  22. /* si le nombre n'est pas nul */   
  23.   if (k >= 0) 
  24.   {
  25.     /* on affiche le chiffre significatif */
  26.      printf("%u",z[k]);
  27.      k--;
  28.      while (k  >= 0)
  29.         {
  30. /* le format %.4u ecrit sur 4 caracteres en completant a gauche par des 0 */
  31. /* si necessaire */
  32.           printf("%.4u",z[k]);
  33.           k--;
  34.         } 
  35.    }
  36.   else
  37.      printf("0" ); 
  38.   printf("\n" );
  39. }
  43. void zero(BIGNUM z)
  44. {
  45. /* initialise le grand entier z a zero */
  46.   memset(z,0,sizeof(BIGNUM));
  47. }
  49. void alea(BIGNUM z, unsigned int n)
  50. {
  51. /* remplit les n premiers chiffres du grand entier z de fac,on aleatoire */
  52.   unsigned int j;
  53.  
  54.   zero(z);
  55.   srand(time(NULL));
  56.   for (j = 0; j < n; j++)
  57.      z[j] = rand() % BASE;
  58. }
  62. void copier(BIGNUM x, BIGNUM y)
  63. {
  64. /* copie la valeur du grand entier x dans le grand entier y */
  65.   memcpy(y,x,sizeof(BIGNUM));
  66. }
  69. /***** PARTIE A COMPLETER *************/
  71. void inc(BIGNUM z)
  72. {
  73. int i=0;
  74. while ( ((z[i] + 1) % BASE == 0)  &&  (i<MAX)  )
  75.        {  z[i] = 0;
  76.       i++;
  77.    } 
  78. z[i] += 1;
  80. }
  82. unsigned char plusgrand(BIGNUM z, BIGNUM y)
  83. { int k;
  85.   k=MAX-1;
  86.   while ((k >=0) && (z[k] == 0) && (y[k] == 0))
  87.     k--;
  89.        if ( z[k] >= y[k] )
  90.       return 1;
  91. else
  92.   return 0; 
  93. }
  97. void sub(BIGNUM z, BIGNUM y)
  98. {  int k, c;
  99.    BIGNUM x;                             
  100.    copier(y, x);
  101.  
  102.   for( k=0 ; k<MAX ; k++)
  103.      { 
  104.      if( z[k] >= x[k] )
  105.         z[k] -= x[k];
  107. else {
  108.        c = z[k] - x[k];
  109.        z[k] = BASE + c;
  110.  
  111.    if(k<MAX-1)
  112.       x[k+1] += 1;
  113. }
  114. }
  115. }
  119. void divise(BIGNUM a, BIGNUM b, BIGNUM q, BIGNUM r)
  120. { 
  121.    BIGNUM x;
  122.    copier(a, x);
  123.    zero(q);
  124.  
  125.    while ( plusgrand(x, b) ){
  126.            sub(x, b);
  127.            inc(q);
  128.            }
  129.    copier(x, r);
  130. }
  133. int notnul(BIGNUM a)
  134. {
  135.     int k = MAX-1;
  137. while (k >= 0){
  139.        if (a[k] != 0)
  140.        return 1;
  141.  
  142.    k--;
  143.       }   
  144. return 0;
  145. }
  148. void pgcd(BIGNUM a, BIGNUM b, BIGNUM d)
  149. {
  150.    if ( notnul(b) == 0 )
  151.         copier(a, d); 
  152.  
  153. else{
  154.        if (  plusgrand(a, b) )
  155.         {    sub(a, b);
  156. pgcd(a, b, d);
  157. }
  160. else  pgcd(b, a, d);
  162. }
  163. }
  166. int main()
  167. {
  168.   BIGNUM a,b,q,r,d;
  169.  
  171.   zero(a);
  172.   zero(b);
  173.   a[15] = 998;
  174.   a[1] = 12;
  175.   b[14] = 8;
  176.   b[0] = 3;
  177.   printf("a = " );
  178.   affiche(a);
  179.   printf("b = " );
  180.   affiche(b);
  181.   divise(a,b,q,r);
  182.   printf("q = " );
  183.   affiche(q);
  184.   printf("r = " );
  185.   affiche(r);
  186.   zero(a);
  187.   zero(b);
  188.   a[4] = 768;
  189.   a[1] = 12;
  190.   a[0] = 727;
  191.   b[3] = 2179;
  192.   b[2] = 2321;
  193.   b[0] = 428; 
  194.   printf("a = " );
  195.   affiche(a);
  196.   printf("b = " );
  197.   affiche(b);
  198.  
  199.   pgcd(a, b, d);
  200.   printf("d = " );
  201.   affiche(d);
  202.   return 0;
  203. }


 
Bonjour, voila je sors d'examen de TP de C
J'ai apparemennt tout fait mais lorsque j'execute la derniere fonction (pgcd) j'obtenais une segmentation fault (snif)
J'ai cherché tout le temps qui me restait mais j'ai pas trouvé d'ou venait l'erreur :/
Vous pouvez m'aider ? Pask du coup ca me turlupine, j'espère que ca vient pas d'autres fonctions au dessus ou quoi
Je précise que la main() et les fonctions affiche zero alea copier m'étaient données. voila
merci!


---------------
- mon feed-back
Reply

Marsh Posté le 15-01-2004 à 16:59:44   

Reply

Marsh Posté le 15-01-2004 à 19:55:38    

moi j'ai ça  
 

Code :
  1. int pgcd(int a, int b)
  2.     {
  3.     if (a%b == 0)
  4.         return b ;          // égal à 'b' si 'b' divise 'a'
  6.     return pgcd(b, a%b) ;   // sinon, égal un pgdc de 'b' et de 'a' mod 'b'
  7.     }

Reply

Marsh Posté le 15-01-2004 à 20:07:05    

jagstang a écrit :

moi j'ai ça  
 

Code :
  1. int pgcd(int a, int b)
  2.     {
  3.     if (a%b == 0)
  4.         return b ;          // égal à 'b' si 'b' divise 'a'
  6.     return pgcd(b, a%b) ;   // sinon, égal un pgdc de 'b' et de 'a' mod 'b'
  7.     }




 
 
Dans l'énnoncé ils demandaient de faire l'algo du PGCD avec les conditions suivantes :
 
pgcd(a, 0) = a
pgcd(a, b) = pgcd(b, a)  si a<b
pgcd(a, b) = pgcd(a-b, b) si a>=b
 
Vous la voyez ou l'erreur de segmentation dans le code :??:


---------------
- mon feed-back
Reply

Marsh Posté le 15-01-2004 à 23:11:55    

C'est bon j'ai trouvé... en fait c'était ma fonction plusgrand() qui foirait pff...
 
Voila le code correct :
 

Code :
  1. unsigned char plusgrand(BIGNUM z, BIGNUM y)
  2.   { int k;
  3.    
  4.       k=MAX-1;
  5.       while ((k >= 0) && (z[k] == 0) && (y[k] == 0))
  6.         k--;
  7.    
  8.       while((k >= 0) && (z[k] == y[k]))
  9.              k--;
  11.         if(k<0)
  12.         return 1;
  13.         else{
  14.                    if(z[k]>=y[k])
  15.                      return 1;
  16.                    else
  17.                      return 0;
  18.             }
  19.   }


Message édité par Zipo le 15-01-2004 à 23:39:02

---------------
- mon feed-back
Reply

Sujets relatifs:

Leave a Replay

Make sure you enter the(*)required information where indicate.HTML code is not allowed