code ça marche mais se n'est pas une bonne execution

code ça marche mais se n'est pas une bonne execution - C++ - Programmation

Marsh Posté le 14-03-2012 à 15:34:29    

salut,
j'ai un grand besoin d'aide si quelqu'un a déja travailler sur le RSA et qu'il peut me guider je l'en remercie. [#ff0000][/#ff0000]
j'ai executer ce code mais ça ne génére pas les bonne clés.
POur info le fonctionnement et le suivant :
1. choisir p et q deux nombres premiers.
2. calculer n = p * q
3. calculer l'indicatrice d'Euler phi = (p - 1)(q - 1)
4. Choisir e (exposant public) tel qu'il soit premier avec phi
5. Trouver d (exposant prive) tel que : (d * e) = 1 (mod phi)
la cle publique est donc (e; n) et la cle privee est (d; n) ;
 
 
 

Code :
  1. #include "stdafx.h"
  2. #include <stdlib.h>
  3. #include <stdio.h>
  4. #include <time.h>
  5. #include <math.h>
  7. int myRandom();
  8. double isPremier(double);
  10. /**
  11. main()
  12. */
  13. int main(void) {
  14.    int _p, _q, _n, _indEuler;
  15.    int _e = 0;
  16.    int _d = 0;
  17.    int _flag = 0;
  18.  
  19.    do {
  20.                 _p = myRandom();
  21.                 _q = myRandom();
  22.    } while((isPremier(_p) == 0) || (isPremier(_q) == 0) );
  23.  
  24.    _n = _p * _q;
  25.    _indEuler = (_p - 1) * (_q - 1);
  26.  
  27.         while( !_flag ) {
  28.                 if( isPremier(_e) ) {
  29.                         if( (_e % _indEuler == 0) && (_e < _indEuler) ) {
  30.                                 _e++;
  31.                         } else {
  32.                                 // on a trouve l'exposant publique
  33.                                 // on sort du while
  34.                                 _flag = 1;
  35.                         }
  36.                 } else {
  37.                         _e++;
  38.                 }
  39.         }
  40.      
  41.         _flag = 0;
  42.         // on cherche l'exposant prive
  43.         while( !_flag ) {
  44.                 if( ( _e * _d - 1) % _n == 0 ) {
  45.                         // on a trouve
  46.                         _flag = 1;
  47.                 } else
  48.                         _d++;
  49.         }
  50.  
  51.         printf("p = %d\n", _p);
  52.         printf("q = %d\n", _q);
  53.         printf("n = %d\n", _n);
  54.         //printf("Indice Euler = %d\n", _indEuler);
  55.         printf("e = %d\n", _e);
  56.         printf("d = %d\n", _d);
  58.         return (0);
  59. }
  61. /**
  62.         myRandom()
  63. */
  64. int myRandom(void)
  65. {
  66.    static int first = 0;
  67.  
  68.    if (first == 0) {
  69. unsigned int t = (unsigned int) time(NULL);
  70. srand(t);
  71.      
  72.       first = 1;
  73.    }
  74.    return ( rand () % 10 );
  75. }
  77. /**
  78.         isPremier(int)
  79. */
  80. double isPremier(double _entier) {
  81.   static int _premiers[] = {2,3,5,7,11,13,17,19,23,29,31,37,41,
  82.                         43,47,53,59,61,67,71,73,79,83,89,97};
  83.   int i, n;
  84.   double d;
  87.   /* 1 n'est pas premier */
  88.   if (_entier == 1)
  89.   {
  90.     return 0;
  91.   }
  93.   n = sizeof (_premiers) / sizeof (*_premiers);
  95.   /* D'abord on regarde si n est divisible par les nombres premiers dans le tableau */
  96.   for (i = 0; i < n; i++) {
  97.         if (_entier == _premiers[i]) {
  98.       return 1;
  99.     }
  100.         if (((int)_entier) % _premiers[i] == 0) {
  101.       return 0;
  102.         }
  103.   }
  105.   /* Ensuite, on doit regarder a partir du dernier element du tableau+2 jusqu'a sqrt(_entier)... */
  106.   d = sqrt (_entier) + 0.5; /* Le 0.5 permet de tester si c'est un carre parfait... */
  107.   i = _premiers[i-1] + 2;
  109.   while (i < d) {
  110.     if (((int)_entier) % i == 0) {
  111.       return 0;
  112.     }
  113.     i += 2;
  114.   }
  116.   return 1;   
  117. }
  119. /**
  120.  
  121. */
  122. int pgcd(long _a, long _b)
  123. {
  124.     long _dividende = labs(_a); /* le _dividende contient la valeur absolue de a */
  125.     long _diviseur = labs(_b);  /* le _diviseur contient la valeur absolue de b  */
  126.     long _quotient;
  127.     long _reste;
  129.     int _fin = 0;
  131.     /*
  132.      * on ne calcule le pgcd de deux nombres que s'ils sont diffŽrents de zŽro
  133.      */
  134.     if(_a != 0 && _b != 0) {
  135.        while(!_fin) {
  136.            /* On applique la division euclidienne */
  137.            _reste = _dividende % _diviseur;
  138.            _quotient = (_dividende - _reste) / _diviseur;
  140.            /* Si le _reste est diffŽrent de 0, on continue l'algorithme */
  141.            if(_reste != 0) {
  142.                _dividende = _diviseur;
  143.                _diviseur = _reste;
  144.            }
  145.            else {
  146.                _fin = 1;
  147.            }
  148.        }
  149.     }
  150.     else {
  151.        /* Erreur ... */
  152.        _diviseur = 0;
  153.     }
  155.      
  156.     return _diviseur;
  157. }

[#ff0000][/#ff0000][#ff0e00][/#ff0e00]


Message édité par gilou le 14-03-2012 à 16:36:22
Reply

Marsh Posté le 14-03-2012 à 15:34:29   

Reply

Sujets relatifs:

Leave a Replay

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