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 averti
    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
    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 confirmé
    Avatar de tyrtamos
    Homme Profil pro
    Retraité
    Inscrit en
    Décembre 2007
    Messages
    4 486
    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 486
    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

  3. #3
    Membre averti
    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
    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 confirmé
    Avatar de tyrtamos
    Homme Profil pro
    Retraité
    Inscrit en
    Décembre 2007
    Messages
    4 486
    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 486
    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))

  5. #5
    Membre Expert

    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
    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 confirmé
    Avatar de fred1599
    Homme Profil pro
    Lead Dev Python
    Inscrit en
    Juillet 2006
    Messages
    4 049
    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 : 4 049
    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'))

+ 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