IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

Calcul scientifique Python Discussion :

Sur la division euclidienne [Python 3.X]


Sujet :

Calcul scientifique Python

  1. #1
    Membre à l'essai
    Homme Profil pro
    statisticien (retraité)
    Inscrit en
    Mai 2017
    Messages
    14
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Val de Marne (Île de France)

    Informations professionnelles :
    Activité : statisticien (retraité)

    Informations forums :
    Inscription : Mai 2017
    Messages : 14
    Points : 10
    Points
    10
    Par défaut Sur la division euclidienne
    La division euclidienne de A par B (entiers relatifs), qui s’écrit A = BxQ + R, impose que le reste R soit positif ou nul (et bien sûr plus petit que la valeur absolue du diviseur B). Ceci implique que pour diviser A par B, on doit chercher le multiple de B qui se rapproche le plus de A, tout en restant inférieur à A, de manière à pouvoir rajouter un reste R positif (éventuellement nul).

    Ceci va donner, par exemple :
    On divise 103 par 13 -> 103 = 13 x 7 + 12 -> Q = 7 et R = 12
    On divise 103 par -13 -> 103 = -13 x (-7) + 12 -> Q = -7 et R = 12
    On divise -103 par -13 -> -103 = -13 x 8 + 1 -> Q = 8 et R = 1
    On divise -103 par 13 -> -103 = 13 x (-8) + 1 -> Q = -8 et R = 1

    La fonction Python « divmod » ne réalise pas la division euclidienne, puisque divmod(103, -13) = (-8,-1), avec un mauvais quotient et un reste négatif, ce qui est contraire à la division euclidienne.
    On peut bien sûr former le bon couple (-7,12) en constituant le couple (-divmod(103,13)[0], divmod(103,13)[1])… et se débrouiller dans chaque cas.

    Mais je n’arrive pas à savoir s’il existe une fonction Python qui donne directement le bon résultat, suivant les bonnes règles de la division euclidienne.

    Une idée ?

  2. #2
    Expert éminent
    Avatar de tyrtamos
    Homme Profil pro
    Retraité
    Inscrit en
    Décembre 2007
    Messages
    4 461
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Var (Provence Alpes Côte d'Azur)

    Informations professionnelles :
    Activité : Retraité

    Informations forums :
    Inscription : Décembre 2007
    Messages : 4 461
    Points : 9 248
    Points
    9 248
    Billets dans le blog
    6
    Par défaut
    Bonjour,

    Est-ce qu'une correction de ce genre conviendrait:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    def divmod2(x, y):
        q, r = divmod(x, y)
        if y<0:
            return q+1, r+abs(y)
        return q, r
    Un expert est une personne qui a fait toutes les erreurs qui peuvent être faites, dans un domaine étroit... (Niels Bohr)
    Mes recettes python: http://www.jpvweb.com

  3. #3
    Membre à l'essai
    Homme Profil pro
    statisticien (retraité)
    Inscrit en
    Mai 2017
    Messages
    14
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Val de Marne (Île de France)

    Informations professionnelles :
    Activité : statisticien (retraité)

    Informations forums :
    Inscription : Mai 2017
    Messages : 14
    Points : 10
    Points
    10
    Par défaut
    Oui, c'est tout à fait ça. Donc, il faut bien en passer par une modification de divmod pour aboutir à une véritable division euclidienne ; bizarre qu'une fonction directe n'ait pas été implémentée ; mais au fond, c'est sans doute un problème mineur...

    Merci pour cette réponse.

  4. #4
    Expert éminent
    Avatar de tyrtamos
    Homme Profil pro
    Retraité
    Inscrit en
    Décembre 2007
    Messages
    4 461
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Var (Provence Alpes Côte d'Azur)

    Informations professionnelles :
    Activité : Retraité

    Informations forums :
    Inscription : Décembre 2007
    Messages : 4 461
    Points : 9 248
    Points
    9 248
    Billets dans le blog
    6
    Par défaut
    En cas de besoin, on peut aussi faire ça en une seule ligne:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    divmod2 = lambda x, y : divmod(x, y) if y>0 else (x//y+1, x%y+abs(y))
    Un expert est une personne qui a fait toutes les erreurs qui peuvent être faites, dans un domaine étroit... (Niels Bohr)
    Mes recettes python: http://www.jpvweb.com

  5. #5
    Membre émérite

    Homme Profil pro
    Ingénieur calcul scientifique
    Inscrit en
    Mars 2013
    Messages
    1 229
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Alpes Maritimes (Provence Alpes Côte d'Azur)

    Informations professionnelles :
    Activité : Ingénieur calcul scientifique

    Informations forums :
    Inscription : Mars 2013
    Messages : 1 229
    Points : 2 328
    Points
    2 328
    Par défaut
    Petite citation de la doc :

    divmod(a, b)
    Prend deux nombres (non complexes) et donne leur quotient et reste de leur division entière sous forme d’une paire de nombres. Avec des opérandes de types différents, les règles des opérateurs binaires s’appliquent. Pour deux entiers le résultat est le même que (a // b, a % b). Pour des nombres à virgule flottante le résultat est (q, a % b), où q est généralement math.floor(a / b) mais peut valoir un de moins. Dans tous les cas q * b + a % b est très proche de a. Si a % b est différent de zéro, il a le même signe que b, et 0 <= abs(a % b) < abs(b).
    Et a//b c'est l'opérateur floor division. Donc le quotient entier de a par b est nécéssairement l'arrondi à la baisse du rapport a/b. Donc c'est normal.

    Et oui dans votre cas il faut en passer par un petit intermède.

  6. #6
    Expert éminent
    Avatar de fred1599
    Homme Profil pro
    Lead Dev Python
    Inscrit en
    Juillet 2006
    Messages
    3 823
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Meurthe et Moselle (Lorraine)

    Informations professionnelles :
    Activité : Lead Dev Python
    Secteur : Arts - Culture

    Informations forums :
    Inscription : Juillet 2006
    Messages : 3 823
    Points : 7 119
    Points
    7 119
    Par défaut
    On peut utiliser le module decimal.

    Il y a une méthode divmod qui fait le taf.

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    >>> import decimal
    >>> context = decimal.Context()
    >>> context.divmod(103, -13)
    (Decimal('-7'), Decimal('12'))
    Celui qui trouve sans chercher est celui qui a longtemps cherché sans trouver.(Bachelard)
    La connaissance s'acquiert par l'expérience, tout le reste n'est que de l'information.(Einstein)

+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. dépacement capacité (division euclidienne)
    Par Darksnakes dans le forum Débuter
    Réponses: 5
    Dernier message: 01/12/2008, 09h20
  2. Division euclidienne
    Par cixiii dans le forum R
    Réponses: 1
    Dernier message: 24/04/2008, 18h02
  3. Division euclidienne polynomiale (avec modulos)
    Par Batoche dans le forum Mathématiques
    Réponses: 5
    Dernier message: 18/12/2006, 14h03
  4. division euclidienne de polynôme
    Par gronaze dans le forum Mathématiques
    Réponses: 1
    Dernier message: 29/06/2006, 20h53

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo