Equação de Pell

Na matemática, mais especificamente na Teoria dos Números, a equação de Pell (também chamada de equação de Pell-Fermat) é a equação:

x 2 d y 2 = 1 {\displaystyle x^{2}-dy^{2}=1}

Onde x {\displaystyle x} e y {\displaystyle y} são números inteiros e d {\displaystyle d} um número natural.

Esta equação foi nomeada em homenagem ao matemático inglês John Pell, foi estudada por Brahmagupta no século VII e por Fermat no século XVII.[1]

Introdução

As equações de Pell-Fermat são estudadas há milênios na Índia e na Grécia. Eles tinham uma grande interesse particularmente no caso de d {\displaystyle d} = 2 uma vez que sua solução forneciam uma boas aproximações racionais de 2 x / y . {\displaystyle {\sqrt {2}}\approx x/y.} Baudhayana encontrou os pares (17,12) e (577, 408) forneciam muito boas a aproximações 2 577 / 408. {\displaystyle {\sqrt {2}}\approx 577/408.} Já Arquimedes usou a equação no caso de n = 3 e obteve a aproximação 3 1351 / 780 {\displaystyle {\sqrt {3}}\approx 1351/780} . com Brahmagupta , que desenvolveu o método chakravala para resolver a equação de Pell e outras equações indeterminadas quadráticas em sua Brahma Sphuta Siddhanta em 628, cerca de mil anos antes da época de Pell. O nome de Pell, nestas equações, ocorre devido a um erro de Euler atribuindo ao matemático inglês John Pell (1610-1685) o estudo da mesma. Aparentemente foi Lord Brouncker (1620-1684) o primeiro matemático europeu moderno a estudar as equações de Pell-Fermat.[2]


Soluções

Note que

 
  
    
      
        
          x
          
            2
          
        
        
        d
        
          y
          
            2
          
        
        =
        1
        
        (
        x
        
        
          
            d
          
        
        y
        )
        (
        x
        +
        
          
            d
          
        
        y
        )
        =
        1
      
    
    {\textstyle x^{2}-dy^{2}=1\Longleftrightarrow (x-{\sqrt {d}}y)(x+{\sqrt {d}}y)=1}
  

Se d N {\displaystyle {\sqrt {d}}\in \mathbb {N} } , isto é, se d {\displaystyle d} é um quadrado perfeito, então

 
  
    
      
        (
        x
        
        
          
            d
          
        
        y
        )
        (
        x
        +
        
          
            d
          
        
        y
        )
        =
        1
        
        
          
            {
            
              
                
                  x
                  
                  
                    
                      d
                    
                  
                  y
                  =
                  1
                
              
              
                
                  x
                  +
                  
                    
                      d
                    
                  
                  y
                  =
                  1
                
              
            
            
          
        
        
        
          ou
        
        
        
          
            {
            
              
                
                  x
                  
                  
                    
                      d
                    
                  
                  y
                  =
                  
                  1
                
              
              
                
                  x
                  +
                  
                    
                      d
                    
                  
                  y
                  =
                  
                  1
                
              
            
            
          
        
      
    
    {\displaystyle (x-{\sqrt {d}}y)(x+{\sqrt {d}}y)=1\Longrightarrow {\begin{cases}x-{\sqrt {d}}y=1\\x+{\sqrt {d}}y=1\end{cases}}\quad {\text{ou}}\quad {\begin{cases}x-{\sqrt {d}}y=-1\\x+{\sqrt {d}}y=-1\end{cases}}}
  

Como d N {\displaystyle {\sqrt {d}}\in \mathbb {N} } , então existe n {\displaystyle n} natural tal que n = d {\displaystyle n={\sqrt {d}}} . Assim, no primeiro caso acima, temos que

 
  
    
      
        
          
            {
            
              
                
                  x
                  
                  n
                  y
                  =
                  1
                
              
              
                
                  x
                  +
                  n
                  y
                  =
                  1
                
              
            
            
          
        
        
        (
        x
        
        n
        y
        )
        +
        (
        x
        +
        n
        y
        )
        =
        1
        +
        1
      
    
    {\displaystyle {\begin{cases}x-ny=1\\x+ny=1\end{cases}}\Longrightarrow (x-ny)+(x+ny)=1+1}
  

Assim,

 
  
    
      
        (
        x
        +
        x
        )
        +
        (
        n
        y
        
        n
        y
        )
        =
        2
        
        x
        =
        1.
      
    
    {\displaystyle (x+x)+(ny-ny)=2\Longrightarrow x=1.}
  

Substituindo x = 1 {\displaystyle x=1} em uma das equações do sistema, teremos que y = 0 {\displaystyle y=0} .

Resolvendo o caso

 
  
    
      
        
          
            {
            
              
                
                  x
                  
                  n
                  y
                  =
                  
                  1
                
              
              
                
                  x
                  +
                  n
                  y
                  =
                  
                  1
                
              
            
            
          
        
      
    
    {\displaystyle {\begin{cases}x-ny=-1\\x+ny=-1\end{cases}}}
  

Teremos x = 1 {\displaystyle x=-1} e y = 0 {\displaystyle y=0} . Assim, os pares ( x , y ) = ( 1 , 0 ) {\displaystyle (x,y)=(1,0)} e ( x , y ) = ( 1 , 0 ) {\displaystyle (x,y)=(-1,0)} são ditos soluções triviais.

Joseph Louis Lagrange provou que, contanto que d {\displaystyle d} não seja um quadrado perfeito, a equação de Pell tem infinitas soluções inteiras distintas.

Os inversíveis em Z [ d ] {\displaystyle Z[{\sqrt {d}}]}

Como x 2 d y 2 = ( x + d y ) ( x d y ) {\displaystyle x^{2}-dy^{2}=(x+{\sqrt {d}}y)(x-{\sqrt {d}}y)} , então iremos usar os números da forma a + d b {\displaystyle a+{\sqrt {d}}b} , números em Z [ d ] {\displaystyle \mathbb {Z} [{\sqrt {d}}]} , para resolver a equação.

Definição

Definiremos a norma N : R R {\displaystyle N:\mathbb {R} \rightarrow \mathbb {R} } da seguinte forma

 
  
    
      
        N
        (
        u
        )
        =
        u
        
          
            
              u
              ¯
            
          
        
        .
      
    
    {\displaystyle N(u)=u{\bar {u}}.}
  

Onde u = x + d y {\displaystyle u=x+{\sqrt {d}}y} e u ¯ = x d y {\displaystyle {\bar {u}}=x-{\sqrt {d}}y} .


Proposição

Seja u = x + d y {\displaystyle u=x+{\sqrt {d}}y} uma solução da equação de Pell, qualquer potência u n ; 2 n N {\displaystyle u^{n};2\leq n\in \mathbb {N} } é, também, uma solução da equação de Pell.

Demonstração

Usaremos o processo de indução. Assim, note que para n = 2 {\displaystyle n=2}

N ( u 2 ) = N ( ( x + d y ) ( x + d y ) ) = N ( x 2 + d y 2 + 2 d x y ) = ( x 2 + d y 2 + 2 d x y ) ( x 2 + y 2 2 d x y ) = ( x 2 + d y 2 ) 2 ( 2 d x y ) 2 = x 4 + 2 d x 2 y 2 + d 2 y 4 4 d x 2 y 2 = x 4 2 d x 2 y 2 + d 2 y 4 = ( x 2 d y 2 ) ( x 2 d y 2 ) = N ( u ) N ( u ) = 1 {\displaystyle {\begin{aligned}N(u^{2})&=N{\Big (}(x+{\sqrt {d}}y)(x+{\sqrt {d}}y){\Big )}\\&=N{\Big (}x^{2}+dy^{2}+2{\sqrt {d}}xy{\Big )}\\&=\left(x^{2}+dy^{2}+2{\sqrt {d}}xy\right)\left(x^{2}+y^{2}-2{\sqrt {d}}xy\right)\\&=\left(x^{2}+dy^{2}\right)^{2}-\left(2{\sqrt {d}}xy\right)^{2}\\&=x^{4}+2dx^{2}y^{2}+d^{2}y^{4}-4dx^{2}y^{2}\\&=x^{4}-2dx^{2}y^{2}+d^{2}y^{4}\\&=(x^{2}-dy^{2})(x^{2}-dy^{2})\\&=N(u)N(u)\\&=1\end{aligned}}}

Suponha que para todo n k , N ( u n ) = N ( u n 1 u ) = 1 {\displaystyle n\leq k,N(u^{n})=N(u^{n-1}\cdot u)=1} , então para n = k + 1 {\displaystyle n=k+1} , temos que

N ( u k + 1 ) = N ( u k u ) = N ( u k ) N ( u ) = 1 {\displaystyle {\begin{aligned}N(u^{k+1})&=N(u^{k}\cdot u)\\&=N(u^{k})N(u)\\&=1\end{aligned}}}

Assim, fica demonstrado que qualquer potência de uma solução é, ainda, uma solução da equação de Pell.

Proposição

Seja u = x + d y {\displaystyle u=x+{\sqrt {d}}y} a solução fundamental da equação de Pell e u k = x k + d y k {\displaystyle u^{k}=x_{k}+{\sqrt {d}}y_{k}} a k-ésima potência de u {\displaystyle u} , então

 
  
    
      
        
          
            
              
                
                  x
                  
                    k
                  
                
              
              
                
                =
                
                  x
                  
                    k
                    
                    1
                  
                
                x
                +
                d
                
                  y
                  
                    k
                    
                    1
                  
                
                y
              
            
            
              
                
                  y
                  
                    k
                  
                
              
              
                
                =
                
                  x
                  
                    k
                    
                    1
                  
                
                y
                +
                
                  y
                  
                    k
                    
                    1
                  
                
                x
              
            
          
        
      
    
    {\displaystyle {\begin{aligned}x_{k}&=x_{k-1}x+dy_{k-1}y\\y_{k}&=x_{k-1}y+y_{k-1}x\end{aligned}}}
  

Demonstração

Note que u k = u k 1 u {\displaystyle u^{k}=u^{k-1}u} , logo

u k = ( x k 1 + d y k 1 ) ( x + d y ) = x k 1 x + x k 1 d y + d y k 1 x + d y k 1 y x k + d y k = ( x k 1 x + d y k 1 y ) + d ( x k 1 y + y k 1 x ) {\displaystyle {\begin{aligned}u^{k}&=(x_{k-1}+{\sqrt {d}}y_{k-1})(x+{\sqrt {d}}y)\\&=x_{k-1}x+x_{k-1}{\sqrt {d}}y+{\sqrt {d}}y_{k-1}x+dy_{k-1}y\\\Rightarrow x_{k}+{\sqrt {d}}y_{k}&=(x_{k-1}x+dy_{k-1}y)+{\sqrt {d}}(x_{k-1}y+y_{k-1}x)\end{aligned}}}

Comparando os termos, temos que x k = x k 1 x + d y k 1 y  e  y k = x k 1 y + y k 1 x {\displaystyle x_{k}=x_{k-1}x+dy_{k-1}y{\text{ e }}y_{k}=x_{k-1}y+y_{k-1}x} .

Observação: note que u = x 1 + y 1 d {\displaystyle u=x_{1}+y_{1}{\sqrt {d}}} é a solução fundamental

O grupo dos inversíveis em Z [ d ] {\displaystyle \mathbb {Z} [{\sqrt {d}}]}

Seja a restrição,

N : Z [ d ] Z {\displaystyle N:\mathbb {Z} [{\sqrt {d}}]^{*}\rightarrow \mathbb {Z} ^{*}}

Onde Z [ d ] = { u Z [ d ] ; N ( u ) = 1 } {\displaystyle \mathbb {Z} [{\sqrt {d}}]^{*}=\{u\in \mathbb {Z} [{\sqrt {d}}];N(u)=1\}} e Z = { 1 , 1 } {\displaystyle \mathbb {Z} ^{*}=\{-1,1\}} (inversíveis em Z {\displaystyle \mathbb {Z} } ).

Note que a norma definida dessa forma é um homomorfismo.

Agora iremos verificar que G = ( Z [ d ] , ) {\displaystyle G=(\mathbb {Z} [{\sqrt {d}}]^{*},\cdot )} é um grupo:

  1. Associatividade: como Z [ d ] R {\displaystyle \mathbb {Z} [{\sqrt {d}}]\subset \mathbb {R} } e a operação é a multiplicação usual, garantimos a associatividade;
  2. Elemento neutro: Note que 1 u = u     u R {\displaystyle 1\cdot u=u\ \forall \ u\in \mathbb {R} } , logo, u Z [ d ] 1 u = u {\displaystyle u\in \mathbb {Z} [{\sqrt {d}}]^{*}\Rightarrow 1\cdot u=u} ;
  3. Inverso : Por construção, todo u Z [ d ] {\displaystyle u\in \mathbb {Z} [{\sqrt {d}}]^{*}} possui inverso.

Seja u Z [ d ] {\displaystyle u\in \mathbb {Z} [{\sqrt {d}}]^{*}} a solução fundamental da equação de Pell, então u k {\displaystyle u^{k}} ainda será uma solução da equação, logo, u {\displaystyle u} gera G {\displaystyle G} , portanto, G {\displaystyle G} é um grupo cíclico, logo, abeliano. Assim,

 
  
    
      
        G
        =
        {
        
          u
          
            0
          
        
        ,
        
          u
          
            1
          
        
        ,
        
        }
        =
        {
        1
        ,
        u
        ,
        
          u
          
            2
          
        
        ,
        .
        .
        .
        }
      
    
    {\displaystyle G=\{u^{0},u^{1},\dots \}=\{1,u,u^{2},...\}}
  



Aplicações

Da equação de Pell, temos que x 2 = 1 + d y 2 {\displaystyle x^{2}=1+dy^{2}} , logo

 
  
    
      
        
          
            
              x
              
                2
              
            
            
              y
              
                2
              
            
          
        
        =
        
          
            1
            
              y
              
                2
              
            
          
        
        +
        d
        
        
          
            x
            y
          
        
        =
        
          
            
              
                1
                
                  y
                  
                    2
                  
                
              
            
            +
            d
          
        
        .
      
    
    {\displaystyle {\frac {x^{2}}{y^{2}}}={\frac {1}{y^{2}}}+d\Longrightarrow {\frac {x}{y}}={\sqrt {{\frac {1}{y^{2}}}+d}}.}
  

Intuitivamente, podemos perceber que a razão x / y {\displaystyle x/y} nos dará boas aproximações para d {\displaystyle {\sqrt {d}}} quando 1 / y 2 {\displaystyle 1/y^{2}} for pequeno. Matematicamente, podemos formular essa ideia como

 
  
    
      
        
          lim
          
            y
            
            +
            
          
        
        
          
            
              
                1
                
                  y
                  
                    2
                  
                
              
            
            +
            d
          
        
        =
        
          
            d
          
        
        .
      
    
    {\displaystyle \lim _{y\rightarrow +\infty }{\sqrt {{\frac {1}{y^{2}}}+d}}={\sqrt {d}}.}
  

Com isso, podemos usar frações continuadas para escrever o número irracional d {\displaystyle {\sqrt {d}}} na forma d = [ a 0 , a 1 , , 2 a 0 ¯ ] {\displaystyle {\sqrt {d}}=[a_{0},{\overline {a_{1},\dots ,2a_{0}}}]}

Exemplo

Seja a equação de Pell

 
  
    
      
        
          x
          
            2
          
        
        
        2
        
          y
          
            2
          
        
        =
        1
      
    
    {\displaystyle x^{2}-2y^{2}=1}
  
.

A solução fundamental dessa equação é o par ( 3 , 2 ) {\displaystyle (3,2)} , pois, 3 2 2 ( 2 ) 2 = 9 8 = 1 {\displaystyle 3^{2}-2(2)^{2}=9-8=1} . Assim, u = 3 + 2 2 {\displaystyle u=3+2{\sqrt {2}}} , então u k = 3 x k 1 + 4 y k 1 + 2 ( x k 1 2 + y k 1 3 ) {\displaystyle u^{k}=3x_{k-1}+4y_{k-1}+{\sqrt {2}}(x_{k-1}2+y_{k-1}3)} , daí

 
  
    
      
        
          
            
              
                
                  u
                  
                    2
                  
                
              
              
                
                =
                3
                
                  x
                  
                    1
                  
                
                +
                4
                
                  y
                  
                    1
                  
                
                +
                (
                2
                
                  x
                  
                    1
                  
                
                +
                3
                
                  y
                  
                    1
                  
                
                )
                
                  
                    2
                  
                
                =
                17
                +
                12
                
                  
                    2
                  
                
              
            
            
              
                
                  u
                  
                    3
                  
                
              
              
                
                =
                3
                
                  x
                  
                    2
                  
                
                +
                4
                
                  y
                  
                    2
                  
                
                +
                (
                2
                
                  x
                  
                    2
                  
                
                +
                3
                
                  y
                  
                    2
                  
                
                )
                
                  
                    2
                  
                
                =
                99
                +
                70
                
                  
                    2
                  
                
              
            
            
              
                
                  u
                  
                    4
                  
                
              
              
                
                =
                3
                
                  x
                  
                    3
                  
                
                +
                4
                
                  y
                  
                    3
                  
                
                +
                (
                2
                
                  x
                  
                    3
                  
                
                +
                3
                
                  y
                  
                    3
                  
                
                )
                
                  
                    2
                  
                
                =
                577
                +
                408
                
                  
                    2
                  
                
              
            
          
        
      
    
    {\displaystyle {\begin{aligned}u^{2}&=3x_{1}+4y_{1}+(2x_{1}+3y_{1}){\sqrt {2}}=17+12{\sqrt {2}}\\u^{3}&=3x_{2}+4y_{2}+(2x_{2}+3y_{2}){\sqrt {2}}=99+70{\sqrt {2}}\\u^{4}&=3x_{3}+4y_{3}+(2x_{3}+3y_{3}){\sqrt {2}}=577+408{\sqrt {2}}\\\end{aligned}}}
  

 
  
    
      
        
      
    
    {\displaystyle \vdots }
  

Fazendo a razão entre os coeficientes das soluções, temos que

 
  
    
      
        
          
            3
            2
          
        
        =
        1
        ,
        5
        ;
        
        
          
            17
            12
          
        
        =
        1
        ,
        41
        
          
            6
            ¯
          
        
        ;
        
        
          
            99
            70
          
        
        =
        1
        ,
        41
        
          
            428571
            ¯
          
        
        ;
        
        
          
            577
            408
          
        
        =
        1.4142
        
          
            1568627450980392
            ¯
          
        
      
    
    {\displaystyle {\frac {3}{2}}=1,5;\quad {\frac {17}{12}}=1,41{\overline {6}};\quad {\frac {99}{70}}=1,41{\overline {428571}};\quad {\frac {577}{408}}=1.4142{\overline {1568627450980392}}}
  

Essas razões estão se aproximando de 2 = 1 , 414213562373 {\displaystyle {\sqrt {2}}=1,414213562373\dots }



[3] [4]

Referências

  1. «Pell Equation». Wolfram - MathWorld 
  2. Gondim, Rodrigo (2011). Aritmética em Retas e Cônicas (PDF). Sergipe: [s.n.] pp. 57 – 64. Consultado em 26 de novembro de 2019 
  3. Michel Waldschmidt (18 de fevereiro de 2016). «Pell's Equation» (PDF) 
  4. Bruce Ikenaga (2019). «The Pell-Fermat Equation» 
Ícone de esboço Este artigo sobre matemática é um esboço. Você pode ajudar a Wikipédia expandindo-o.
  • v
  • d
  • e
  • Portal da matemática