Impl�mentation d'un algorithme foireuse

Bonjour,

Je fait un petit programme cryptant par courbe elliptique. Dans celui-ci je dois impl�menter l'exponentiation rapide. C'est un algo simple mais j'ai r�ussi � me planter et je ne vois pas o� :oops:...

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

template<typename T> std::pair<T, T> QuickExponentiation(T p, T a, T e, std::pair<T, T> point) { if( 1 < e) { if(p % 2 ==0) return AddPoints(p, a, QuickExponentiation(p, a, e/2, point), QuickExponentiation(p, a, e/2,point)); else return AddPoints(p, a, AddPoints(p, a, QuickExponentiation(p, a, (e-1)/2,point), QuickExponentiation(p, a, (e-1)/2, point)), point); } return point; }

---

La fonction AddPoints() fonctionne (v�rifi�), voici son prototye:
template<typename T>
std::pair< T, T> AddPoints(T p, T a, const std::pair<T, T>& firstPoint, const std::pair<T, T>& secondPoint);

Votre aide est la bienvenue.
Merci

Bonjour khazna,

Dans quelle mesure "�a ne marche pas", comment se comporte la fonction ?

Cdt.

Par exemple, lorsque je fais:

Code:

---

1 2	QuickExponentiation(29, 1, 2, std::pair<int, int>(28, 12))

---

�a me sort le point (7, 25) alors qu'en faisant l'addition avec AddPoints() �a me sort le bon r�sultat (-2, 22).
A ne rien y comprendre :aie:

Ton "e / 2" est certainement pass� en tant qu'entier � QuickExponentiation. C'est vraiment ce que tu veux ?

Oui, je crois :aie:

Euh.. J'ai peut �tre trouv� !!!

La fonction AddPoint prend comme param�tres des r�f�rences...

std::pair< T, T> AddPoints(T p, T a, const std::pair<T, T>& firstPoint, const std::pair<T, T>& secondPoint);

Or tu lui envoies un retour de fonction !!!

Utilise des variables tempons et rev�rifie...

C'est pas �a...
Raahhh foutue fonction :evilred:...

Sans savoir ce que sont sens�s faire a, p, e (celui ci, on peut encore deviner), AddPoint ou QuickExponentiation, difficile de dire ce qui ne fonctionne pas.

Point de perf : Pour un algorithme se disant quick, calculer deux fois la m�me valeur (QuickExponentiation(p, a, e/2, point)) syst�matiquement est �trange...

Que retourne exactement AddPoint ?

AddPoints retourne un point.
En fait p est un nombre premier.
Une courbe elliptique est d�finie par une �quation du type y�=x^3 + ax +b
Moi je me place dans le groupe Z/pZ(je vous passe un pdf expliquant �a).

Pour ce qui est du Quick, en fait j'essaie d�j� de le faire marcher, pour la vitesse je verrai apr�s :aie:...
J'ai fait un truc na�f qui fonctionne:

Code:

---

1 2 3 4 5 6 7 8 9 10 11

template<typename T> std::pair<T, T> QuickExponentiation(T p, T a, T e, std::pair<T, T> point) { std::pair<T, T> tmpPoint = point; for(T i = 0; i < e - 1; i++) { tmpPoint = AddPoints(p, a, tmpPoint, point); } return tmpPoint; }

---

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23

template<typename T> std::pair<T, T> QuickExponentiation(T p, T a, T e, std::pair<T, T> point) { std::pair<T, T> tmpPoint; std::pair<T, T> tmpPoint2;   if( 1 < e) { if(p % 2 ==0) { tmpPoint = QuickExponentiation(p, a, e/2, point); return AddPoints(p, a, tmpPoint , tmpPoint ); } else { tmpPoint = QuickExponentiation(p, a,(e-1)/2, point); tmpPoint2 = AddPoints(p, a, tmpPoint , tmpPoint ); return AddPoints(p, a, tmpPoint2 , point); } } return point; }

---

Et comme �a ?

Attention, indentation foireuse et trompeuse: Tu as oubli� les accolades (Edit: Corrig� depuis).

(m�me si dans le cas pr�sent, �a ne change pas grand-chose, puisque la section avant le else se termine par un return)

Toujours pas...
Je me demande bien pourquoi �a ne fonctionne pas: �a a l'air th�oriquement juste :D...

Citation:

---

Envoy� par M�dinoc

Attention, indentation foireuse et trompeuse: Tu as oubli� les accolades.

---

Mis � jour... (J'ai post� trop vite...) :mrgreen:

Apr�s cette apr�s midi de prog' endiabl� (bon pas toute l'aprem' non plus :aie:) j'ai r�ussi � finir mon crypteur par courbe elliptique (reste plus qu'a faire un prog' genre chat qui exploite ce qui a �t� fait et c'est bon) :yaisse2:.
N�anmmoins je vais chercher � r�soudre ce probl�me et si vous avez une piste je suis pr�t � l'explorer :) .
Merci pour votre aide, je l'aurais pas fini sans les membres de ce forum :king:.

Bon je viens de relire le code et j'ai trouv�: le test avec le modulo devait �tre e % 2 == 0
On a vraiment cherch� trop loin :aie: