Annales de mathématiques pures et appliquées/Tome 16/Analise indéterminée, article 1
# Annales de mathématiques pures et appliquées/Tome 16/Analise indéterminée, article 1

## ANALISE INDÉTERMINÉE.

### Résolution générale de l’équation indéterminée du premier degré à deux inconnues ;

La résolution générale de l’équation indéterminée du premier degré, à deux inconnues, paraît avoir été trouvée pour la première fois par les Indiens. Les commentateurs de Bhasker ou Bhascara-Acharya, attribuent cette découverte à Arya-Bhatta, géomètre indien, que l’on croit presque contemporain de Diophante ; mais les ouvrages de cet auteur ayant été perdus, il est difficile de juger du degré de généralité qu’il avait donné à sa solution. Cependant, nous possédons le traité d’algèbre de Bramegupta (qui, d’après les calculs de M. Bentley, a été composé au commencement du septième siècle de l’ère chrétienne) où l’on trouve la résolution générale de l’équation du premier degré à deux inconnues. Cette résolution est exposée aussi dans les ouvrages de Bhascara-Acharya, et dans ceux de tous les analistes indiens plus modernes. La méthode dont ils ont fait usage est semblable à celle que Bachet de Méziriac publia en France en 1624. On sait qu’elle consiste à réduire l’équation proposée ax+b=cy, à l’équation ax₁+1=y₁, et à résoudre celle-ci, en cherchant le plus grand commun diviseur entre a et c.

Lagrange a résolu l’équation dont il s’agit à l’aide des fractions continues, et M. Gauss l’a réduite à sa théorie des congruences ; mais toutes ces méthodes, qui dans le fond sont identiques entre elles, n’ont pas toute la généralité qu’on pourrait désirer. En effet, il est clair que les racines d’une équation à plusieurs inconnues, de même que celles d’une équation à une seule inconnue, doivent être fonctions de ses coefficiens exprimés généralement ; et cependant, même pour résoudre l’équation du premier degré à deux inconnues, qui est la plus simple de toutes, il faut connaître les coefficiens en nombres, ce qui montre combien les méthodes connues sont imparfaites.

La note que nous publions ici a pour objet de donner l’expression générale des racines entières d’une équation du premier degré à deux inconnues, en fonction de ses coefficiens. Elle est extraite d’un mémoire sur la théorie des nombres, présenté à l’Académie royale des sciences de Paris, et qui doit paraître dans le recueil des Savans étrangers. Notre méthode s’applique à toutes les équations indéterminées ; elle sert aussi à résoudre directement, et avec simplicité, les équations desquelles dépend la division du cercle, et à traiter beaucoup d’autres questions ; mais ces recherches ne sont pas de nature à trouver place ici, et nous les réservons pour une autre circonstance.

Étant donné l’équation à une seule inconnue

(1)xᵐ−1=0,

si l’on représente par Pₙ, Pₙ₋ₘ, Pₙ₋₂ₘ, Pₙ₋₃ₘ, … les sommes des puissances n. ⁱᵉᵐᵉˢ, (n−m) ⁱᵉᵐᵉˢ, (n−2m) ⁱᵉᵐᵉˢ, (n−3m) ⁱᵉᵐᵉˢ, … de ses racines, on aura

Pₙ=Pₙ₋ₘ=Pₙ₋₂ₘ=Pₙ₋₃ₘ=…;

de sorte que, si n est un multiple de m, on obtient Pₙ=m; et, dans le cas contraire, on trouve Pₙ=0. En exprimant les racines de l’équation (1) en fonctions circulaires, on aura

Pₙ={(Cos.0.ϖ/m+√(−1)Sin.0.ϖ/m )ⁿ+(Cos.2ϖ/m+√(−1)Sin.2ϖ/m )ⁿ+(Cos.4ϖ/m+√(−1)Sin.4ϖ/m )ⁿ+………………………+(Cos.2uϖ/m+√(−1)Sin.2uϖ/m )ⁿ+………………………+(Cos.2(m−1)ϖ/m+√(−1)Sin.2(m−1)ϖ/m )ⁿ}.

Si l’on transforme le second membre, au moyen de la relation connue (Cos.z+√(−1)Sin.z)ⁿ=Cos.nz+√(−1)Sin.nz, et qu’on néglige les imaginaires qui, dans le cas actuel, doivent nécessairement se détruire, on obtiendra

Pₙ=Cos.0nϖ/m+Cos.2nϖ/m+Cos.4nϖ/m+…Cos.2unϖ/m+…Cos.2(m−1)nϖ/m,

c’est-à-dire,

Pₙ=∑ᵤ₌₀ᵘ⁼ᵐCos.2unϖ/m=(Sin⁡.2(n−n/2m )ϖ+Sin.nϖ/m)/2Sin.nϖ/m .

et la valeur de cette expression sera m ou zéro, suivant que le nombre n/m sera entier ou fractionnaire.

On sait qu’étant proposé de résoudre l’équation ax+b=cy, en nombres entiers, il suffit de trouver une valeur a de x, comprise entre zéro et c, car les autres s’en déduisent en ajoutant à celle-là un multiple quelconque de c; de sorte qu’on a, en général, x=α+cz, z étant un nombre entier quelconque.

Maintenant il faut observer que si, dans l’équation (2), on fait n=ax+b, m=c, et que l’on donne à x successivement toutes les valeurs 0, 1, 2, 3, …c−1, on obtiendra l’intégrale

∑ₓ₌₀ˣ⁼ᶜ{Cos.(0(ax+b)ϖ)/c+Cos.(2(ax+b)ϖ)/c+…+Cos.(2u(ax+b)ϖ)/c

+…+Cos.(2(c−1)(ax+b)ϖ)/c },

qui aura pour valeur c répété autant de fois que la quantité (ax+b)/c a de valeurs entières, lorsqu’on y fait x égal à un nombre entier moindre que c; d’où il suit que la formule

1/c ∑ₓ₌₀ˣ⁼ᶜ{Cos.(0(ax+b)ϖ)/c+Cos.(2(ax+b)ϖ)/c+…+Cos.(2u(ax+b)ϖ)/c

+…+Cos.(2(c−1)(ax+b)ϖ)/c };

exprimera le nombre des solutions de l’équation ax+b=cy, en supposant qu’on ne prenne pour x que des nombres entiers plus petits que c.

Si l’on considère le terme général de la série (3), on aura l’équation

1/c ∑ₓ₌₀ˣ⁼ᶜCos.(2u(ax+b)ϖ)/c=(Sin.2u/c (b+ac−1/2 a)ϖ−Sin.2u/c (b−1/2 a)ϖ)/2cSin.2uaϖ/c ,

dans le second membre de laquelle le numérateur est toujours zéro ; mais dont le dénominateur ne peut se réduire à zéro que lorsque a et c ont un diviseur commun, plus grand que l’unité, puisque u est toujours plus petit que c. Il résulte de là que, si a et c sont premiers entre eux, tous les termes de la série (3) s’évanouissent, excepté le premier, dont la valeur se réduit à

1/c ∑ₓ₌₀ˣ⁼ᶜCos.(0(ax+b)ϖ)/c=c/c=1.

Mais, si a et c ont un facteur commun g, on supposera a=mg, c=ng, et, en faisant u=n, on obtiendra

1/c ∑ₓ₌₀ˣ⁼ᶜCos.(2n(ax+b)ϖ)/c=(Sin.2n/c (b+ac−1/2 a)ϖ−Sin.2n/c (b−1/2 a)ϖ)/2cSin.2naϖ/c

=(Sin.(2(b+ang−1/2 a)ϖ)/g−Sin.2(b−1/2 a)ϖ/g)/2ngSin.aϖ/g .

Cette expression se réduit 0/0, en vertu de l’hypothèse a=mg. On devra donc différentier le numérateur et le dénominateur par rapport à a, pour en avoir la valeur déterminée, et l’on trouvera

(Sin⁡.2(b+ang−1/2 a) ϖ/g−Sin⁡.2(b−1/2 a) ϖ/g)/2ngSin.aϖ/g

=(2(ng−1/2 ) ϖ/g Cos⁡.2(b+ang−1/2 a) ϖ/g+ϖ/g Cos⁡.2(b−1/2 a) ϖ/g)/2ng.ϖ/g.Cos.aϖ/g

=2ngϖ.Cos⁡.2(b−1/2 a) ϖ/g /2ngϖ.Cos.aϖ/g =Cos⁡.2(b−1/2 a) ϖ/g /Cos.aϖ/g

=Cos⁡.2(b−mg/2 ) ϖ/g /Cos.mϖ=(Cos.2bϖ/g.Cos.mϖ+Sin.2bϖ/g.Sin.mϖ)/Cos.mϖ

=Cos.2bϖ/g.

Si, au lieu de prendre u=n, on fait, en général, u=en, e étant un nombre entier quelconque, on trouve

1/c ∑ₓ₌₀ˣ⁼ᶜCos⁡.2en.(ax+b)/c ϖ=Cos.2ebϖ/g ;

et, comme le nombre n est compris g−1 fois dans c−1, on pourra faire successivement e=0, 1, 2, 3…g−1; et la valeur de l’intégrale (3) sera exprimée (dans le cas actuel, où l’on suppose que a et c ont un commun diviseur g) par la série

1+Cos.2bϖ/g+Cos.4bϖ/g+…+Cos.2(g−1)ϖ/g,

dont la somme

(Sin⁡.2(b−1/2 b/g )ϖ+Sin.bϖ/g)/2Sin.bϖ/g ,

a pour valeur g; lorsque b/g est un nombre entier, et se réduit à zéro, dans le cas contraire.

De là résulte 1.ᵒ que l’équation ax+b=cy a toujours une solution entière, et plus petite que c, lorsque a et c n’ont d’autres diviseurs communs que l’unité ;

2.ᵒ Que si a et c ont un commun diviseur g, différent de l’unité, qui ne divise point b, cette équation n’admet aucune solution entière ;

3.ᵒ Qu’enfin, si b/g est un nombre entier, on trouvera pour x un nombre g de valeurs entières, plus petites que c, qui satisferont à l’équation proposée.

Puisque l’intégrale (3) représente le nombre des solutions entières de l’équation ax+b=cy, en prenant pour x des valeurs moindres que c, il est clair que la formule

1/c ∑ₓ₌₀ˣ⁼ᶜx{Cos.(0(ax+b)ϖ)/c+Cos.(2(ax+b)ϖ)/c+…

+Cos.(2u(ax+b)ϖ)/c+…+Cos.(2(c−1)(ax+b)ϖ)/c }

=1/c ∑ₓ₌₀ˣ⁼ᶜx{1+Cos⁡.2((ax+b)/c)ϖ+…+Cos⁡.2u((ax+b)/c)ϖ+…

+Cos⁡.2(c−1)((ax+b)/c)ϖ}

exprimera la somme des valeurs de x, entières et moindres que c, qui satisfont à l’équation ax+b=cy, lorsqu’elle est résoluble et que, lorsqu’elle ne l’est pas, cette intégrale se réduit à zéro.

Nous avons démontré que, si a et c ont un facteur commun, différent de l’unité, qui ne divise pas b, l’équation ax+b=cy n’admet aucune solution entière, et comme, si ce facteur commun divise aussi b, on peut toujours le supprimer, il sera permis dans en cas, de supposer que a et c sont premiers entre eux ; et alors on sera assuré qu’il existe toujours une valeur entière de x, comprise entre zéro et c y qui satisfait à l’équation dont il s’agit, et qu’il n’en existe qu’une seule.

Actuellement, pour trouver cette valeur de x, on considérera le terme général de l’intégrale (4), et on aura

1/c ∑ₓ₌₀ˣ⁼ᶜx.Cos⁡.2u (ax+b)/c ϖ

=((c−1)Sin⁡.2u(b+ca−1/2 a) ϖ/c−Sin⁡.2u(b+1/2 a) ϖ/c)/2cSin.uaϖ/c

+(Cos⁡.2u(ca+b−a) ϖ/c−Cos⁡.2u(b+a) ϖ/c)/c(2Sin.uaϖ/c )².

Il faudra faire successivement u=1, 2, 3, …(c−1), et ajouter au résultat le premier terme de la série (4) qui est

1/c ∑ₓ₌₀ˣ⁼ᶜx=c(c−1)/2c=c−1/c.

Puisque a et c sont premiers entre eux, et que u est plus petit que c, il s’ensuit que le dénominateur 2cSin.uaϖ/c du second membre de l’équation (5) ne pourra jamais s’évanouir ; on obtiendra, par conséquent, en faisant les réductions nécessaires,

((c−1)Sin⁡.2u(b+ca−1/2 a) ϖ/c−Sin⁡.2u(b+1/2 a) ϖ/c)/2c.Sin.uaϖ/c +(Cos⁡.2u(b+ca−a) ϖ/c−Cos⁡.2u(b+a) ϖ/c)/c(2Sin.uaϖ/c )² }=Sin⁡.2u(b−1/2 a) ϖ/c /2.Sin.uaϖ/c

et partant

1/c ∑ₓ₌₀ˣ⁼ᶜx{1+Cos⁡.2((ax+b)/c)ϖ+Cos⁡.4((ax+b)/c)ϖ+…

+Cos⁡.2u((ax+b)/c)ϖ+…+Cos⁡.2(c−1)((ax+b)/c)ϖ}

=c−1/2+Sin⁡.2(b−1/2 a) ϖ/c /2Sin.aϖ/c +Sin⁡.4(b−1/2 a) ϖ/c /2Sin.2aϖ/c +…

+Sin⁡.2u(b−1/2 a) ϖ/c /2Sin.uaϖ/c +…+Sin⁡.2(c−1)(b−1/2 a) ϖ/c /2Sin.(c−1)aϖ/c

=c−1/2+1/2 ∑ᵤ₌₁ᵘ⁼ᶜ Sin⁡.2u(b−1/2 a) ϖ/c /Sin.uaϖ/c =α.

Cette formule très-simple donne pour α la plus petite valeur de x qui satisfasse à l’équation ax+b=cy, en nombres entiers ; et toutes les autres valeurs sont exprimées par l’équation x=α+cz, z étant un nombre entier quelconque.

Soit proposé, par exemple, de résoudre en nombres entiers, l’équation

3x+1=4y;

on aura, en comparant à l’équation générale ax+b=cy

a=3, b=1, c=4;

et par conséquent

α=4−1/2+1/2 ∑ᵤ₌₁ᵘ⁼⁴.Sin⁡.2u(1−3/2 ) ϖ/4 /Sin.3uϖ/4 =3/2−1/2 ∑ᵤ₌₁ᵘ⁼⁴.Sin.u ϖ/4 /Sin⁡.3u ϖ/4 ;

c’est-à-dire,

α=3/2−1/2 { Sin.ϖ/4 /Sin.3ϖ/4 +Sin.2ϖ/4 /Sin.6ϖ/4 +Sin.3ϖ/4 /Sin.9ϖ/4 }

=3/2−1/2 (1−1+1)=3/2−1/2=2/2=1;

et toutes les valeurs de x qui résolvent l’équation 3x+1=4y seront données par l’équation x=1+4z, comme on le sait d’ailleurs.

La valeur de α peut, en général, se calculer à l’aide des tables du sinus. Il est vrai que, par ce moyen, on n’obtiendra, le plus souvent, que des valeurs fractionnaires approchées ; mais comme, d’après ce qui précède, x ne peut avoir que des valeurs entières, on en trouvera la valeur exacte en substituant à cette valeur approchée le nombre entier le plus voisin.

On peut observer que, puisqu’on a

Sin.(2b−a) uϖ/c /Sin.auϖ/c =Sin.2buϖ/c.Cos.auϖ/c−Cos.2buϖ/c.Sin.auϖ/c /Sin.auϖ/c

=Sin.2buϖ/c.Cot.auϖ/c−Cos.2buϖ/c,

et que d’ailleurs

∑ᵤ₌₁ᵘ⁼ᶜCos.2buϖ/c=−1,

on pourra écrire

α=c−1/2+1/2 ∑ᵤ₌₁ᵘ⁼ᶜ.Sin.2buϖ/c.Cot.auϖ/c+1/2

=1/2 {c+∑ᵤ₌₁ᵘ⁼ᶜ.Sin.2buϖ/c.Cot.auϖ/c }.

On pourra faire usage de cette expression, aussi bien que de la précédente, pour résoudre l’équation proposée.

* ↑ Géomètre indien, de la ville de Bidder, sur la frontière septentrionale de l’Indostan, où il paraît qu’il enseignait les mathématiques, vers la fin du XII.ᵉ siècle. Voyez l’ouvrage de M. Hutton intitulé : Tracts on Mathematical, etc., 3 vol. in-8.ᵒ, Londres, 1812. J. D. G.
* ↑ Voy. entre autres, l’Introduction à l’analise infinitésimale d’Euler, tom. I, chap. XIV, n.ᵒ 260. J. D. G.
* ↑ Ces formules, étendues à un nombre quelconque d’équations du premier degré, entre un grand nombre d’inconnues, compléteraient la théorie exposée à la page 147 du III.ᵉ volume du présent recueil. J. D. G.
