terça-feira, 21 de fevereiro de 2012

019-Equação Diofantina Linear

Em complemento aos artigos A combinação linear ax+byO algoritmo de Euclides, segue este sobre a resolução em [;\mathbb{Z};] da equação diofantina linear

[;ax+by=c;]

para [;a\neq 0;], [;b\neq 0;], [;c;], [; x;]e [;y;], todos [;\in \mathbb{Z};]
 Conforme já sabemos,

a) [;ax+by;] é múltiplo inteiro de [;d=mdc(a,b)>0;]; e
b) O algoritmo de Euclides nos permite achar [;d=mdc(a,b);] e ainda, [;s;] e [;t;] [;\in \mathbb{Z};], onde [;as+bt=d=mdc(a,b);] .

Vamos ver como podemos utilizar este conhecimento para resolver a equação em questão.

a) [;ax+by;] é múltiplo inteiro de [;mdc(a,b);].

Temos, então, um critério para saber se a equação [;ax+by=c;] tem solução em [;\mathbb{Z};]. A equação terá soluções inteiras, se e somente se, o [;mdc(a,b);] dividir [;c;].

Exemplo 1: Verificar se a equação [;91x+14y=70;] tem soluções em [;\mathbb{Z};]. Resolução: [;mdc(91,14)=7;] e [;7;] divide [;c=70;], logo, existem soluções inteiras para esta equação.

Exemplo 2: Verificar se a equação [;16x+24y=38;] tem soluções em [;\mathbb{Z};]. Resolução: [;mdc(16,24)=8;] e [;8;] não divide [;c=38;], logo, não existem soluções inteiras para esta equação.

Em consequência:

1- Se [;mdc(a,b)=1;], a equação [;ax+by=c;] sempre terá soluções em [;\mathbb{Z};]

Exemplo 3: [;6x+11y=100;], com [;mdc(6,11)=1;] 

2- A equação [;ax+by=1;] terá soluções em [;\mathbb{Z};] apenas se [;mdc(a,b)=1;]

Exemplo 4: [;7x+10y=1;] tem soluções porque [;mdc(7,10)=1;]

Exemplo 5: [;6x+15y=1;] não tem soluções porque [;mdc(6,15)=3\neq1;]

b) O algoritmo de Euclides nos permite achar [;mdc(a,b);]  e ainda, [;s;] e [;t;] [;\in \mathbb{Z};], onde [;as+bt=d=mdc(a,b);] 

Com isto, temos uma forma de achar uma solução inicial para [;ax+by=c;], com [;c\neq d;]. De fato, pois sabendo que [;c=\alpha d;] e [;\alpha=\fra{c}{d};], multiplicamos a primeira igualdade de b) por [;\alpha;] 

[;a(\alpha s)+b(\alpha t)=\alpha d=c;] 

e assim obtemos a solução inicial [;x=\alpha s;] e [;y=\alpha t;] para nossa equação.

Exemplo 6: Vamos achar uma solução inicial [;(x_0,y_0);] para a equação do exemplo 1:


[;91x+14y=70;] 

Com, [;a=91;] , [;b=14;] e [;c=70;]

Utilizamos o algoritmo de Euclides para calcular o [;mdc(91,14);] :

[;91=6.14+7;]  (*)
[;14=7.2+0;]

Logo, o último resto não-nulo é [;d=mdc(91,14)=7;]

Para [;91s+14t=d;] temos o equivalente [;91(1)+14(-6)=7;]obtido  por intermédio de (*).

Assim, [;s=1;] e [;t=-6;] com [;d=7;] e [;c=70;]

[;\alpha = \frac{c}{d}=\frac{70}{7}=10;]

[;x=\alpha s=10.1=10;] 
[;y=\alpha t=10.(-6)=-60;] 

Portanto, a equação [;91x+14y=70;] possui a solução inicial [;(x_0,y_0)=(10,-60);].

Agora porque solução inicial? Ao contrário da equação [;axy+bx+cy+d=0;] que tem um número finito de soluções em [;\mathbb{Z};], a equação [;ax+by=c;] possui infinitas soluções inteiras. Isto deve-se ao fato do caráter aditivo desta equação. As equações diofantinas aditivas e multiplicativas mais simples são, respectivamente, [;x+y=a;] e [;xy=a;]. No primeiro caso, para qualquer valor [;x_i \in \mathbb{Z};] temos o correspondente [;y_i=a-x_i;] e, portanto, temos infinitas a soluções da forma [;(x_i, a-x_i);] . No segundo caso, as soluções de [;xy=a;] são, obviamente, os pares [;(d,d');] onde [;d;] e [;d';] [;\in \mathbb{Z};] e [;dd'=a;] e, como sabemos, existe um número finito de divisores de [;a;]. Em [;axy+bx+cy+d=0;], o monômio [;axy;] impõe um caráter multiplicativo a esta equação. Para a equação [;ax+by=c;], diferentemente de [;x+y=k;], uma inserção aleatória de variável inteira, por exemplo [;x=x_i;] [;\in \mathbb{Z};] , não garante um correspondente valor inteiro para [;y;] porque [;y=\frac{c-ax}{b};] e [;y;]será inteiro apenas para os [;c-ax;]múltiplos de [;b;]. Mas como temos infinitos múltiplos de [;b;], temos infinitas soluções para [;ax+by=c;]. Existe uma maneira de saber qual é o formato geral das soluções inteiras da equação diofantina linear. Para isto vamos usar o seguinte lema:

Sejam [;a,b,c \in \mathbb{Z};], de forma que [;mdc(a,b)=1;] . Se [;\frac{bc}{a};] for um inteiro, então [;a;] divide [;c;].

Demonstração: [;mdc(a,b)=1;] implica que existem [;x_0,y_0 \in \mathbb{Z};] onde [;ax_0+by_0=1;]. Decorre que [;(ac)x_0+(bc)y_0=c;] e como [;a;] divide [;bc;], por hipótese, vê-se facilmente que [;a;] divide [;c;].

SOLUÇÃO GERAL ( [;U=\mathbb{Z};] ) DA EQUAÇÃO [;ax+by=c;] .

Seja [;(x_0, y_0);] uma solução inicial e [;(x_1,y_1);] uma outra solução. Assim,

[;ax_0+by_0=ax_1+by_1=c \Rightarrow;]

[;a(x_0-x_1)=b(y_1-y_0);] 

Dividindo ambos os membros por [;mdc(a,b);], temos

[;a^'(x_0-x_1)=b^'(y_1-y_0);] (1) ,
 onde [;mdc(a^',b^')=1;] 

Pelas relações  [;y_1-y_0=\frac{a^'(x_0-x_1)}{b'};] , [;x_0-x_1=\frac{b^'(y_1-y_0)}{a'};] e pelo lema acima demonstrado, concluímos que [;b^';] divide [;x_0-x_1;] e [;a';] divide [; y_1-y_0;] . Logo, existem [;k;] e [;l;] [;\in Z;], tais que

[;x_0-x_1=lb^' \Rightarrow x_1=x_0-lb^';] (2)

[;y_1-y_0=ka^' \Rightarrow y_1=y_0+ka^';] (3)

Agora vamos provar que [;l=k;] . Para isso, basta substituir (2) e (3) em (1):

[;a^'(x_0-x_0+lb^')=b^'(y_0+ka^'-y_0) \Rightarrow;] 
[;a^'b^'l=a^'b^'k \Rightarrow;] 
[;l=k;]

Sendo [;d=mdc(a,b);] e lembrando que [;a^'=\frac{a}{d};] e [;b^'=\frac{b}{d};] concluímos que a solução geral em [;\mathbb{Z};] da equação diofantina linear [;ax+by=c;] é ( ver (2) e (3) ):

[;x=x_0-k\frac{b}{d};]

[;y=y_0+k\fra{a}{d};],

com [;k \in \mathbb{Z};]

ATALHOS PARA UMA SOLUÇÃO INICIAL [;(x_0,y_0);]

O algoritmo de Euclides nos permite, de modo geral, achar uma solução inicial para [;ax+by=c;].Mas, dependendo dos valores de [;a;] ,[;b;] e [;c;] podemos ter uma maneira rápida de conseguir [;(x_0,y_0);]. Temos três casos básicos:

1º Caso: Se [;c;] é múltiplo de [;a;] ou [;b;] então [;(x_0,y_0)=\left(\frac{c}{a},0\right);] ou [;(x_0,y_0)=\left(0,\frac{c}{b}\right);] , respectivamente.

Prova: 1)  [;c=an;], [;\rightarrow;]  [;ax+by=an \Rightarrow a.\frac{c}{a}+b.0=an \Rightarrow c+0=an;]
            2)  [;c=bn;] ,[;\rightarrow;] [;ax+by=bn \Rightarrow a.0+b.\frac{c}{b}=bn \Rightarrow 0+c=bn;]

Exemplo 7: Resolver em [;\mathbb{Z};] a equação [;3x+8y=15;]. Resolução: a equação tem solução pois  [;d=mdc(3,8)=1;] divide [;c=15;]. Como [;c;] é múltiplo de [;a=3;] , temos [;(x_0,y_0)=\left(\frac{c}{a},0\right);][;=\left(\frac{15}{3},0\right)=(5,0);].  Logo,

[;x=x_0-k\fra{b}{d}=5-8k;]

[;y=y_0+k\frac{a}{d}=3k;] 

[;k \in \mathbb{N};] 


2º Caso: Se [;c;] é múltiplo de [;a + b;] então [;(x_0,y_0)= \left(\frac{c}{a + b}, \frac{c}{a + b}\right);]

Prova:
 [;c=(a + b)n;] [;\rightarrow;] [;ax+by=(a+b)n \Rightarrow a.\frac{c}{a+b}+ b.\frac{c}{a+b};][;=(a+b).n \Rightarrow;][;\frac{c}{a + b}(a + b)=(a+b).n;][;\Rightarrow c=(a + b)n;]

Exemplo 8: Resolver em [;\mathbb{Z};] a equação [;2x+7y=99;]. Resolução: a equação tem solução pois [;d=mdc(2,7)=1;] divide [;c=99;]. Como [;c;] é múltiplo de [;a+b=2+7=9;], temos[;(x_0,y_0)=\left(\frac{c}{a+b},\frac{c}{a+b}\right)=\left(\frac{99}{9},\frac{99}{9}\right)=(11,11);]. Assim,

[;x=x_0-k\fra{b}{d}=11-7k;]

[;y=y_0+k\frac{a}{d}=11+2k;] 

 [;k \in \mathbb{N};] 


3º Caso: Se [;c;] é múltiplo de [;|a-b|;] então [;(x_0,y_0)= \left(-\frac{c}{|a - b|}, +\frac{c}{|a - b|}\right);], se [;a<b;] e [;(x_0,y_0)= \left(+\frac{c}{|a - b|}, -\frac{c}{|a - b|}\right);], se [;a>b;].

Prova: análoga ao do 2º caso para cada desigualdade.

Exemplo 9: Resolver em [;\mathbb{Z};] a equação [;6x+30y=24;]. Resolução: [;c=24;] é múltiplo de [;d=mdc(6,30)=6;] o que diz que a equação tem solução. Como [;c;] é múltiplo de [;|a-b|=|6-30|=|-24|=24;] e [;a=6 < b=30;], temos [;(x_0,y_0)= \left(-\frac{c}{|a - b|}, +\frac{c}{|a - b|}\right);][;=\left(-\frac{24}{24},+\frac{24}{24}\right)=(-1,1);]. Portanto,

[;x=x_0-k\fra{b}{d}=-1-5k;]

[;y=y_0+k\frac{a}{d}=1+k;]

 [;k \in \mathbb{N};] 


EXEMPLO 10


Valdir, com R$ [;100,00;] no bolso, deseja comprar ingressos de teatro e de cinema que estão em promoção, para fins de sorteio em uma festa. Sabendo que o ingresso do teatro custa R$ [;10,00;] e o do cinema custa  R$ [;15,00;], de quantas formas possíveis Valdir pode comprar ingressos de teatro e cinema com R$ [;100,00;]?

RESOLUÇÃO

A questão equivale a resolver a equação diofantina [;10x+15y=100;], sendo [; x;]a quantidade de ingressos de teatro e [;y;] a quantidade de ingressos de cinema. Temos, então, [;a=10;], [;b=15;] e [;c=100;].


Temos, ainda, [;d=mdc(10,15)=5;] que, sendo divisor de [;c=100;], nos garante que a equação tem solução em [;\mathbb{Z};]. No entanto, procuramos soluções positivas e verificaremos sua existência analisando as soluções gerais.

Como [;c=100;] é múltiplo de [;a=10;], uma primeira solução em [;\mathbb{Z};] seria [;(x_0,y_0)=\left(\frac{c}{a},0 \right)=(10,0);] ( conforme o 1º caso de atalho para obter uma solução inicial ).As outras soluções inteiras são representadas por

[;x=x_0-k \left(\frac{b}{d}\right) \Rightarrow x=10-3k;] 
[;y=y_0+k\left(\frac{a}{d}\right) \Rightarrow y=2k;] 

Para obter soluções positivas temos que ter [;10-3k>0;] e [;2k>0;] e isto equivale à [;k<3,33...;] e [;k>0;]. Mas como [;k;]é natural, temos [;0<k<4;]. Logo, para [;k=1;], [;k=2;] e [;k=3;], as condições do problema de Valdir são satisfeitas. Então ele possui [;3;] modos de comprar ingressos de teatro e cinema com R$ [;100,00;] naquelas condições de preço:

Para [;k=1;] :

[;x=10-3.1=7;] ingressos de teatro;
[;y=2.1=2;] ingressos de cinema.

Para [;k=2;]:

[;x=10-3.2=4;] ingressos de teatro;
[;y=2.2=4;] ingressos de cinema.

Para [;k=3;]:

[;x=10-3.3=1;] ingresso de teatro;
[;y=2.3=6;] ingressos de cinema.

E assim, Valdir fica satisfeito de poder escolher a alternativa que mais lhe convier.



REFERÊNCIA BIBLIOGRÁFICA

- Teoria dos Números, de  Salahoddin Shokranian, Marcus Soares, Hemar Godinho; Editora UNB, 1998;
- O que é a Matemática?, de Richard Courant e Herbert Robbins, Editora Moderna, 2000.

Imagem: http://www.overmundo.com.br/overblog/do-pio-de-ave-ao-piano-de-nazareth


2 comentários:

  1. Já conhecia a forma de resolver as equações diofantinas, o que me chamou a atenção foi os casos particulares, me passou despercebido. Com um post muito bem elaborado e escrito, restou pouco para outros blogueiros publicar sobre tal assunto. Parabéns e que continue nos supreendendo com posts interessantes.

    ResponderExcluir
  2. Obrigado, Prof! Na verdade eu sempre achava chato utilizar o algoritmo de Euclides (rs). Uma vez vi um post seu muito interessante sobre as equações do tipo [;x^2+y^2=a;]. Valeu!

    ResponderExcluir