Conversion float to integer

Version imprimable

08/10/2013, 13h44
Julien N

Conversion float to integer
Bonjour tout le monde,

Il y a un truc que je ne comprends pas avec la conversion d'un double (flottant) vers un entier avec la fonction int:
Code:

1 2 3 4 5 >>> n = 41063625 >>> print n**(1./3), int(n**(1./3)) 345.0 344 >>> print int(345.0) 345
8O

Python 2.7

Julien
08/10/2013, 14h00
Aegim

Bonjour,

A priori, ce serai du à une erreur de représentation du float, comme expliqué ici http://stackoverflow.com/questions/6...ersion-problem
08/10/2013, 14h16
VinsS
Salut,

En plus de la cause donnée par Aegim, il y a aussi une différence de comportement avec print
Code:

1 2 3 4 5 6 7 8 >>> d = 1./3 >>> d 0.3333333333333333 >>> 41063625 ** d 344.9999999999999 >>> print 41063625 ** d 345.0
08/10/2013, 15h27
Julien N
Hmm, c'est quelque chose que j'aurais du savoir. Merci à toi Aegim pour le lien.

Je vais jeter un oeil plus en détail au document suivant:
http://http://docs.oracle.com/cd/E19..._goldberg.html

@ VinsS
Oui c'est ça qui m'a perturbé, l'affichage de 345.0 (solution exacte).

Peut-être cela vous semblait-il evident au vu du code que j'ai fourni. Je cherche à verifier si le résultat d'une operation (en l'occurence la racine cubique d'un entier) est un entier. La methode que j'emploie en toujours la meme:
Code:

1 2 if n == int(n): print "OK c'est un entier"
Est-ce que par hasard vous connaissez une autre methode? Qui n'aurait pas le même problème.

Ciao

Ju

J'ai pu trouver ce code qui a priori donne le bon résultats sans passer par des floats :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
def cube_root(n): 
    "A modified Newton's Method solver for integral cube roots." 
    if int(n) != n : 
        raise ValueError("must provide an integer") 
    if n in (-1,0,1): return n 
    offset = (1,-1)[n > 0] 
    x = n/2 
    seen = set() 
    steps = 0 
    while 1: 
        y = x**3 
        if y == n: 
            return x 
        dydx = 3.0*x*x 
        x += int((n-y)/dydx)+offset 
        x = x or 1 
        if x in seen: 
            raise ValueError("%d is not a perfect cube" % n) 
        seen.add(x) 
        steps += 1

Le code vient de http://python.6.x6.nabble.com/cube-r....html#a1733529 et marche avec quelques essais. Ca pourrait donc être ce dont tu as besoin.

08/10/2013, 16h09
fred1599

racine cubique(8) -> 2 mais 2**3 -> 8

Pourquoi ne pas travailler avec les entiers plutôt que des flottants ?
08/10/2013, 16h39
Julien N

@Aegim
C'est une solution, mais assez... compliquée! :D

@fred1599
Bonne remarque. De temps en temps je resouds quelques excercices d'un jeu en ligne appelé "Poject Euler". ça me permet de me perfectionner et de découvrir des subtilités et methodes de calculs. Je crois que je ne suis pas le seul sur ce forum à le faire.

Le sujet que je traite est simple. On prend un entier n que l'on élève au cube. On cherche à trouver si parmi toutes les permutations possible du nombre resultant il existe d'autres racines cubiques entières et combien:

ex:

Citation:

345**3 --> 41063625
384**3 --> 56623104
405**3 --> 66430125
56623104 et 66430125 des permutations de 41063625

Bon le sujet de ce poste n'est pas la resolution du problème du jeu, hein, laissez le moi! :P

Ju

Pourquoi ne pas utiliser le module itertools, il est fait pour ça ?

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
>>> from itertools import permutations as perm
>>> def permut(n):
...     n = str(n)
...     size = len(n)
...     for i in perm(n, size):
...         yield ''.join(i)
... 
>>> print('\n'.join(permut(123)))
123
132
213
231
312
321

Bonne continuation...

08/10/2013, 16h56
Aegim

Je reconnais que c'est pas vraiment simple ^^

En plus simple (mais un peu brute), tu peux aussi faire dans un premier temps comme tu faisais initialement avec int. Et ensuite tu testes pour le résultat obtenu (344 dans ton exemple), le résultat +1 (ici 345)et le résultat -1 (343) si le cube est égale à ta valeur de départ. Comme ça en cas d'arrondi raté dans un sens ou dans l'autre tu t'y retrouves.
C'est loin d'être élégant comme solution, mais c'est plus simple ^^'
08/10/2013, 16h58
fred1599

Citation:

En plus simple (mais un peu brute), tu peux aussi faire dans un premier temps comme tu faisais initialement avec int. Et ensuite tu testes pour le résultat obtenu (345 dans ton exemple), le résultat +1 (ici 346)et le résultat -1 (343) si le cube est égale à ta valeur de départ. Comme ça en cas d'arrondi raté dans un sens ou dans l'autre tu t'y retrouves.
C'est loin d'être élégant comme solution, mais c'est plus simple ^^'

Non non c'est juste la bonne solution ;)

@Fred1599
J'utilise bien évidemment itertools ;)

@Aegim
J'ai continué à fouiller sur le net. On trouve pas mal de sujets abordant la racine carrée. J'ai retrouvé notamment un thread sur stackoverflow avec un algo similaire à celui que tu propose.

En voici un autre utilisant le module decimal qui après quelques tests semble aussi fonctionner:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 from decimal import Decimal,localcontext from math import exp,log def iscube(n): nd = Decimal(n) with localcontext() as ctx: ctx.prec = n.bit_length() i = int(exp(log(nd)*(1./3))) if i**3 != n: raise ValueError('input was not a perfect cube') return i
Il y a une astuce, visiblement le module decimal ne permet pas encore les puissances non entières, d'où le passage par le log...

Bref, je suis capable maintenant de verifier si une des permutations est également le resultat d'un entier au cube. Je vais garder en mémoire ce problème sur la precision des flottants. Je suis sûr que j'aurais un de ces quatre ce problème à nouveau.

Merci à tous les deux pour vos réponses rapides. A moins que vous ne souhaitiez rajouter quelque chose je passerais en résolu.

Ciao

Ju