Vous ne trouvez pas de réponse à votre problème ? Alors posez la question dans le forum. Souvenez-vous qu'il n'y a jamais de question bête, mais rester dans l'ignorance parce que l'on n'ose pas poser une question, ça c'est une erreur !

RECHERCHE DES COEFFICIENTS DE BÉZOUT


Information sur la source

Catégorie :Maths & Algorithmes Classé sous : bezout, spécialité, pgcd, coefficients, euclide Niveau : Initié Date de création : 15/06/2008 Date de mise à jour : 19/06/2008 21:01:08 Vu / téléchargé: 2 559 / 55

Note :
Aucune note

Commentaire sur cette source (3)
Ajouter un commentaire et/ou une note

Description

Cliquez pour voir la capture en taille normale
C'est un tout petit programme avec une fonction récursive qui permet de rechercher les coefficients de Bézout.
Le théorème de Bézout dit que, si d=pgcd(a;b), alors il existe deux entiers relatifs u et v tels que a*u+b*v=d.
Par exemple, pgcd(45;8)=1, donc il existe deux entiers relatifs u et v tels que 45*u+8*v=1. Le programme en trouve deux :
u=-3 et v=17.
Il en existe d'autres, mais le programme détermine ceux qu'on obtient en remontant l'algorithme d'euclide.
 

Source

  • #include <iostream>
  • void Bezout(int rkm2,int rkm1,int ukm1=0,int vkm1=1,int ukm2=1,int vkm2=0){
  • /*r(k-4)=r(k-3)*q(k-2)+r(k-2) <=> r(k-2)=a*u(k-2)+b*v(k-2)
  • *r(k-3)=r(k-2)*q(k-1)+r(k-1) <=> r(k-1)=a*u(k-1)+b*v(k-1)
  • *r(k-2)=r(k-1)*q(k)+r(k) <=> r(k)=a*u(k)+b*v(k) avec u(k)=u(k-2)-qk*u(k-1) et v(k)=v(k-2)-qk*v(k-1)*/
  • int rk=rkm2%rkm1;
  • int qk=(rkm2-rk)/rkm1;
  • if(rk==0)
  • std::cout<<rkm1<<"=a*"<<ukm1<<"+b*"<<vkm1<<std::endl;
  • else
  • Bezout(rkm1,rk,ukm2-qk*ukm1,vkm2-qk*vkm1,ukm1,vkm1);
  • return;
  • }
  • int main(int argc,char* argv[]){
  • int a,b;
  • std::cout<<"a=";std::cin>>a;
  • std::cout<<"b=";std::cin>>b;
  • if(a<=0 || b<=0)
  • std::cout<<"On suppose a et b strictement positifs."<<std::endl;
  • else
  • Bezout(a,b);
  • return 0;
  • }
#include <iostream>

void Bezout(int rkm2,int rkm1,int ukm1=0,int vkm1=1,int ukm2=1,int vkm2=0){
/*r(k-4)=r(k-3)*q(k-2)+r(k-2) <=> r(k-2)=a*u(k-2)+b*v(k-2)
 *r(k-3)=r(k-2)*q(k-1)+r(k-1) <=> r(k-1)=a*u(k-1)+b*v(k-1)
 *r(k-2)=r(k-1)*q(k)+r(k)     <=> r(k)=a*u(k)+b*v(k) avec u(k)=u(k-2)-qk*u(k-1) et v(k)=v(k-2)-qk*v(k-1)*/
    int rk=rkm2%rkm1;
    int qk=(rkm2-rk)/rkm1;
    if(rk==0) 
        std::cout<<rkm1<<"=a*"<<ukm1<<"+b*"<<vkm1<<std::endl;
    else 
        Bezout(rkm1,rk,ukm2-qk*ukm1,vkm2-qk*vkm1,ukm1,vkm1);
    return;
}

int main(int argc,char* argv[]){
    int a,b;
    std::cout<<"a=";std::cin>>a;
    std::cout<<"b=";std::cin>>b;
    if(a<=0 || b<=0)
        std::cout<<"On suppose a et b strictement positifs."<<std::endl;
    else
        Bezout(a,b);
    return 0;
}

Fichier Zip

Pour les "Membres Club", vous pouvez télécharger directement un fichier contenu dans le zip sans télécharger le zip en entier !
  •   bezout
    • coefficients de Bezout.ex_Télécharger ce fichier [Réservé aux membres club]38 400 octets
    • main.cppTélécharger ce fichier [Réservé aux membres club]Voir ce fichier938 octets

Télécharger le zip

Historique

19 juin 2008 21:01:08 :
J'ai utilisé le modulo à la place du calcul de la partie entiere... J'avais betement oublié l'usage du modulo % qui donne le reste de la division euclidienne !

Commentaires et avis

signaler à un administrateur
Commentaire de coucou747 le 16/06/2008 09:37:19

salut

http://www.codyx.org/snippet_coefs-bezout-equations-diophantiennes_607.aspx

c'est plus un snippet qu'une source...

Sinon, tu devrais renvoyer les coeficients, plutot que de les afficher, car la, si tu dois faire une autre fonction qui utilise ces coefs, bah tu ne peux pas.

signaler à un administrateur
Commentaire de Pamaury le 17/06/2008 21:54:25

Il y a déjà une source avec plein de commentaires sur comment le coder plus proprement:
http://www.cppfrance.com/codes/ALGORITHME-EUCLIDE-ETENDU_36400.aspx

signaler à un administrateur
Commentaire de linkid le 23/06/2008 11:14:10

Il est aussi en Python et en code casio 35+ pour ceux qui sont intéressés...
http://www.pythonfrance.com/codes/BEZOUT-ALGORITHME-EUCLIDE-ETENDU_45696.aspx

Ajouter un commentaire

Discussions en rapport avec ce code source dans le forum

Probleme sur un programme qui calcule un pgcd [ par El loco ] Voila j ai un probleme sur le programme suivant, il marche correctement avec une grande serie de nombres mais quand je tape 15 et 32 il me donne un pg Format de données [ par leskritiques ] Voici le programme que je veux modofier :#include "chaine.h"#include "entreeSortie.h"#include "chaine.cpp"#include "entreeSortie.cpp"int pgcd (const i PGCD? [ par bethoring ] Tout d'abord bonsoir, je voulais savoir si dans math.h yavé une fonction pour faire un pgcd ( jai regarder dedans et il me semble pas menfin on sai ja algorithme d'euclide [ par dOsSpr0uTosS ] salut à tous, je debute en C++(mais alors vraiment c'est à dire que j'ai commencer y'a une semaine et j'ai pas appris grand chose ...) et donc pour m' PROGRAMME en lige de commande C++ ( ou C... ) [ par tekila_bandita ] Bonjourje pense connaitre la reponse, mais existe-t-il un moyen de faire des programmes en lignes de commandes du style foncion?? Quand je dis fonctio calcule de coefficients de corrélation et corélation partielle [ par mohamedabid ] mohamedabidBonjours Je cherche une implémentation sur java ou c++ des fonctions de corrélations d'une série temporelleExiste-il une fonction prédéfini Programmation langage C [ par Ligorin ] Bonjour, je souhaiterai avoir votre aide,en langage C ,sur le problème de monnaie,en générale. Soit a+4b+6c=8. Objectif : Min a+b+c La solution e


Nos sponsors

Sondage...

CalendriCode

Juillet 2009
LMMJVSD
  12345
6789101112
13141516171819
20212223242526
2728293031  

Consulter la suite du CalendriCode

Téléchargements

Logiciels à télécharger sur le même thème :

Comparez les prix Nouvelle version

Photothèque Nouveau !



Développement réalisé par Nicolas SOREL (Nix) avec l'aide de : Cyril DURAND et Emmanuel (EBArtSoft), Merci à Vincent pour ses précieux conseils
CodeS-SourceS.com© Toute reproduction même partielle est interdite sauf accord écrit du Webmaster
CodeS-SourceS.com© est une marque déposée tous droits réservés
Temps d'éxécution de la page : 0,484 sec

Google Coop CodeS-SourceS Google Coop CodeS-SourceS


Certaines images présentes sur le site (notament certains avatars) sont issues des collections IconShock, donc si vous souhaitez utiliser ces icons vous devez les acheter, ne les copiez pas et ne utilisez pas dans vos sites et applications sans les avoir commandé.