Le couloir - Solution| s4 - 4hs3 - (a2 - b2)2 = 0 |
Une fois cette équation résolue : a² + b² = u² + v² + 2d² = (u + v)² - 2uv + 2d² = (u + v)² - 2h(u + v) + 2d² donne la valeur de d :
| d = √(a² + b² - s² + 2hs)/2 |
Le problème est soluble si 1/h = 1/u + 1/v > 1/a + 1/b, puisque bien évidemment u et v sont < a et b !
Application numérique :
Comme la résolution de l'équation du quatrième degré n'est pas quelque chose de sympathique,
on peut utiliser une méthode par approximations successives, comme la méthode de Newton :
En partant de s0 on calcule sn+1 = sn - f(sn)/f '(sn)
|
s0 = a + b
s = (3s4 - 8hs3 + (a2 - b2)2) / (4s3 - 12hs2) |
Avec a = 3, b = 2, h = 1, s0=5, s converge vers 4,311185197... et d=1,231186...
Un script pour effectuer ces calculs :
Choisissons alors arbitrairement deux triplets de Pythagore primitifs
x² + y² = z² et x'² + y'² = z'² et soit p le PGCD de x et x', en multipliant le premier triplet par x/p et le second par x'/p :
(xx'/p)² + (yx'/p)² = (zx'/p)² et (x'x/p)² + (y'x/p)² = (z'x/p)²
Posons xx'/p = d/m, yx'/p = u/m, zx'/p = a/m, y'x/p = v/m et z'x/p = b/m pour obtenir un système "d'équations du couloir" avec a, b, u, v, d entiers
h = uv/(u + v) = mxx'yy' / p(xy' + x'y). Soit alors q = PGCD(xx'yy'/p², (xy' + x'y)/p)
En choisissant m = (xy' + x'y)/pq, h est aussi un entier. Tout ceci donne donc des données et une solution entières pour un problème du couloir.
|
(r,s) entiers arbitraires, premiers entre eux et de parité opposée, de même pour (r',s')
x = r² - s², y = 2rs, z = r² + s² et de même pour x', y', z' p = PGCD(x,x'), q = PGCD(xx'yy'/p², (xy' + x'y)/p), m = (xy' + x'y)/pq a = mzx'/p, b = mz'x/p h = xx'yy'/qp² d = mxx'/p
|
Par exemple, r = 2, s = 1, r' = 3, s' = 2 donne x = 3, y = 4, z = 5, x' = 5, y' = 12, z' = 13
p = PGCD(3,5) = 1, q = PGCD(720,56) = 8, m = 7
a = 175, b = 273, h = 90 et d = 105 :
Deux échelles de 1,75m et 2,73m se croisent à 90cm du sol dans un couloir de 1,05m de large.
On peut aussi échanger x' et y' soit x = 3, y = 4, z = 5, x' = 12, y' = 5, z' = 13
p = PGCD(3,12) = 3, q = PGCD(80,21) = 1, m = 21
a = 420, b = 273, h = 80 et d = 252
La plus petite possibilité est : a = 70, b = 119, h = 30, d = 56 (u = 42, v = 105)
Obtenue à partir des triplets de Pythagore primitifs (x et y inversés) :
(4, 3, 5) r = 2, s = 1 et (8, 15, 17) r' = 4, s' = 1.
p = PGCD(8,4) = 4 donne d = 8m, b = 10m, u = 6m, a = 17m, v = 15m
donc h = uv/(u + v) = m×6×15/21 = 30m/7 entier si m = 7 (multiple de 7, mais on cherche la plus petite valeur).
Un script pour obtenir de telles valeurs rationnelles
Choisir les valeurs de r, s≠r et r', s'≠r'
Le programme en déduit les valeurs par les formules ci-dessus.