Algorithmique & multiplications

Version imprimable

Bonjour, mon post est un problème que je n'arrive pas à résoudre.

J'ai 2 séries de nombres, en les multipliants entre-eux (1 d'une série multiplié avec 1 autre de l'autre série) doit me donner un nombre prédifini en résultat (ou s'en approcher le plus possible).

Pour moi cela ressemble à de la combinaison, mais je ne trouve pas de solution à mon problème.

Exemple de problème concret :

SÉRIE 1 : 8 ; 14.5 ; 52 ; 100 ; 100 ; 100 ; 100 ; 4.6 ; 8 ; 29 ; 100 ; 100 ; 100 ; 5 ; 9 ; 32 ; 100 ; 100 ; 8.25 ; 14.75 ; 54 ; 100 ; 18.25 ; 33 ; 100 ; 50 ; 90 ; 100 ; 72

SÉRIE 2 : 11.5 ; 12.75 ; 26 ; 78 ; 100 ; 100 ; 100 ; 6.75 ; 6.50 ; 14 ; 42 ; 100 ; 100 ; 7.25 ; 7.50 ; 14 ; 45 ; 100 ; 11.5 ; 11.75 ; 24 ; 68 ; 24 ; 25 ; 50 ; 64 ; 66 ; 100 ; 27

Il faut que je trouve les couples pour ces séries qui font 188.5 ± 1.

Il faudrait donc un algorithme qui calcule : 1er nombre de la série 1 × 1er nombre de la série 2, si c'est pas égal à 188.5±1 alors on continue, 1er nombre de la série 1 × 2ème nombre de la série 2 et à nouveau vérifier la condition, et ainsi de suite ...

Désolé je n'arrive pas trop à expliquer ce que je veux, j'espère que c'est assez clair.

À la fin il me faudrait les couples solutions de valeurs (nb série 1; nb série 2) qui donnent 188.5 ± 1.

Voici un début de code qui ne fonctionne pas, en espérant que vous m'indiquiez mes erreurs et comment l'améliorer :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 F1 = [8,14.5,52,100,100,100,100,4.6,8,29,100,100,100,5,9,32,100,100,8.25,14.75,54,100,18.25,33,100,50,90,100,72] F2 = [11.5,12.75,26,78,100,100,100,6.75,6.50,14,42,100,100,7.25,7.50,14,45,100,11.5,11.75,24,68,24,25,50,64,66,100,27] p=118.5 # résultat voulu e=1 # incertitude prod = [] for i in range (0,len(F1)): for j in range (0,len(F2)): if abs(p - F1(i) * F2(j)) <e : prod.append ((F1(i), F2(j))) print (prod)
Merci d'avance pour votre aide.

15/06/2019, 19h15
wiztricks

Salut,

Citation:

Envoyé par astronas

Voici un début de code qui ne fonctionne pas, en espérant que vous m'indiquiez mes erreurs et comment l'améliorer :

Le code que vous avez posté doit planter ligne 11 avec un beau message d'erreur qui devrait déjà vous indiquer quoi corriger ou vous poser des questions sur comment accéder au i-ème élément d'une liste et ouvrir votre tuto. préféré aux chapitres qui vont bien pour revoir tout çà.

- W

Ceci ?

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
F1 = [8,14.5,52,100,100,100,100,4.6,8,29,100,100,100,5,9,32,100,100,8.25,14.75,54,100,18.25,33,100,50,90,100,72]
F2 = [11.5,12.75,26,78,100,100,100,6.75,6.50,14,42,100,100,7.25,7.50,14,45,100,11.5,11.75,24,68,24,25,50,64,66,100,27]
 
 
p=118.5   # résultat voulu 
e=1   # incertitude 
 
prod = [] 
for i in range (len(F1)):
    for j in range (len(F2)):
        res = F1[i]*F2[j]
        if abs(p-res) < e:
            prod.append((F1[i],F2[j]))
print (prod)

J'ai voulu voir jusqu'où on pouvais aller par curiosité (et c'est pas beau à voir :aie:)
Code:

1 2 3 4 5 6 >>> F1 = [8,14.5,52,100,100,100,100,4.6,8,29,100,100,100,5,9,32,100,100,8.25,14.75,54,100,18.25,33,100,50,90,100,72] >>> F2 = [11.5,12.75,26,78,100,100,100,6.75,6.50,14,42,100,100,7.25,7.50,14,45,100,11.5,11.75,24,68,24,25,50,64,66,100,27] >>> p = 118.5 >>> e = 1 >>> sum(([(a, b) for a in F1 if p-e <= a*b <= p+e] for b in F2), []) [(18.25, 6.5)]
ça m'intrigue et je me demande s'il n'y a pas moyen d'utiliser numpy/scipy

16/06/2019, 21h16
lg_53

si bien sûr que vous pouvez utiliser numpy. Votre code sera plus court et s'éxécutera plus rapidement.

Mais bon si P1 et P2 ne font qu'une vingtaine de valeurs, vous ne verez pas vraiment la différence. Sur des milliers de valeurs, oui.
17/06/2019, 00h05
flapili

Citation:

Envoyé par lg_53

si bien sûr que vous pouvez utiliser numpy. Votre code sera plus court et s'éxécutera plus rapidement.

Mais bon si P1 et P2 ne font qu'une vingtaine de valeurs, vous ne verez pas vraiment la différence. Sur des milliers de valeurs, oui.

je pense que tu m'a confondu avec l'auteur du topic :roll:
Sinon, difficile de faire plus court qu'une ligne :aie:

Bonjour,

Code:

1
2
3
4
5
6
7
F1 = [8,14.5,52,100,100,100,100,4.6,8,29,100,100,100,5,9,32,100,100,8.25,14.75,54,100,18.25,33,100,50,90,100,72]
F2 = [11.5,12.75,26,78,100,100,100,6.75,6.50,14,42,100,100,7.25,7.50,14,45,100,11.5,11.75,24,68,24,25,50,64,66,100,27]
 
for x in F1:
    for y in F2:
        if 188.5-1<=x*y<=188.5+1:
            print([x, y])

Ce qui donne:

Code:

1
2
[29, 6.5]
[14.75, 12.75]

Et en une ligne:

Code:

print([[x, y] for x in F1 for y in F2 if 188.5-1<=x*y<=188.5+1])

qui affiche:

Code:

[[29, 6.5], [14.75, 12.75]]

Citation:

Envoyé par tyrtamos

Code:

print([[x, y] for x in F1 for y in F2 if 188.5-1<=x*y<=188.5+1])

qui affiche:

Code:

[[29, 6.5], [14.75, 12.75]]

J'ai cru que je m'était planté mais non p vaut bien 118.5 et non 188.5 :P
edit: en faite c'est l'auteur qui s'est planté entre son explication et son code, on a qu'a dire qu'on a tous les deux raison :chin:

par curiosité j'ai voulu voir si les perfs étaient raisonnables:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 >>> import timeit >>> import numpy >>> a = numpy.random.rand(10000)*100 >>> b = numpy.random.rand(10000)*100 >>> p = 118.5 >>> e = 1 >>> def t1(): return sum(([(x, y) for x in a if p-e <= x*y <= p+e] for y in b), []) >>> def t2(): return [[x, y] for x in a for y in b if p-e <=x*y<= p+e] >>> timeit.timeit(t1, number=1) 54.252527317000045 >>> timeit.timeit(t2, number=1) 50.68921082600002
Je m'avoue vaincu :hola:

Maintenant je me demande vraiment s'il y a moyen d'utiliser numpy pour ça ? peu être avec un numpy.outer puis parcourir l'array 2D et si la valeurs est dans l'intervalle on affiche F1[i],F2[j] où i et j sont respectivement les indices de l'array 2D

edit 2: je doit être mauvais avec numpy
Code:

1 2 3 4 5 >>> def t3(): return [i for i, x in numpy.ndenumerate(numpy.outer(a, b)) if p-e <= x <= p+e] >>> timeit.timeit(t3, number=1) 82.1113105259999

Non en numpy ce sera plutot qqch comme ceci :
Code:

1 2 def t4(): return np.where( np.abs(np.outer(a,b)-p)<=e )
A noter que ta fonction t3, et ma fonction t4, ne retourne pas les couple de valeur à multiplier, mais les indices pour retrouver les couples dans les tableaux. Pour obtenir les couples on pourrait faire comme ceci :
Code:

1 2 3 4 5 F1 = np.array(a) F2 = np.array(b) def t5(): ### Bon ok, c'est en 2 lignes... Il faut bien factoriser le calcul de idxs ! idxs = np.where( np.abs(np.outer(F1,F2)-p)<=e ) return list(zip( F1[idxs[0]], F2[idxs[1]] ))
C'est 10 fois plus rapide chez moi.

20/06/2019, 14h14
astronas

Merci pour vos réponses tout le monde. Bonne fin de semaine à tous :zoubi:

P.S : je me suis bien planté dans le code, p est bien égal à 188.5 et non 118.5 :oops: