View source for Indução e Recursão Matemática

Este projeto descreve o uso da linguagem Maple para desenvolvimento da Indução Matemática e Recursão na disciplina de Fundamentos da Matemática para Computação II. Com uso da tradução livre do documento que foi disponibilizado do livro Matemática Discreta do Keneth  H. Rosen, descreveremos como a Indução pode ser facilmente entendida e utilizada com Maple .


==Introdução==
Neste capítulo descreveremos como o maple pode ser usado para ajudar na compreensão e construção de provas matemáticas. Capacidades computacionais podem não parecer particularmente relevantes para o estudo das provas, embora na realidade essas capacidades possam ser úteis em provas de várias maneiras. Neste capítulo, descrevemos como o maple pode ser útil para trabalhar com regras formais de inferência, descrevemos como pode ajudar a compreender provas construtivas e não construtivas. Além disso, mostramos como usar o maple para ajudar a desenvolver provas usando indução matemática, até mesmo mostrando sua utilidade para ambos o passo base e passo indutivo, na prova da fórmula de somatório. Ademais, mostraremos como o maple pode ser usado para computar termos de sequencias definidas recursivamente. Vamos também comparar a eficiência da geração de termos dessa sequencia via técnicas indutivas versus técnicas recursivas.

==Método de prova==

Embora o maple não possa receber  teoremas e resultados de provas para esses teoremas, pode receber expressões lógicas e simplificadas ou determinar características tais como: se uma expressão booleana pode ser satisfeita ou se é uma tautologia. Para trabalhar com expressões lógicas no maple, precisamos usar alguns dos recursos oferecidos pelo pacote de '''logic''' (um assunto abordado de maneira mais aprofundada no capítulo 9). 
Primeiramente examinaríamos os operadores lógicos: conjunção, disjunção, negação e implicação. Não existe (no Maple). Para estudar as condicionais, devemos trabalhar com os operadores booleanos inativos oferecidos pelo pacote '''logic'''.  Todos esses são iniciados com o caracter '''mexpr''' <math>\&</math> , por exemplo: usamos <math>\&and</math> invés de <math>and</math> e <math>\&not</math> invés de <math>not</math>. Em seguida, vão alguns exemplos do uso de operadores booleanos inativos:

[[File:imagem1.png]]

Nos preocupamos agora em determinar como o maple simplifica expressões booleanas caso estejam combinadas. Começamos com um simples exemplo de dupla negação: 

[[File:imagem2.png]]

Isto pode ser simplificado através do uso da função '''bsimp''' do maple.

[[File:imagem3.png]]

Agora, seguimos para um exemplo mais complexo, no qual o leitor pode confirmar a corretude da simplificação, construindo uma tabela verdade.

[[File:imagem4.png]]

O próximo exemplo ilustra a simplificação do Modus Ponens. Primeiro afirmamos que p implica q e  p é verdadeiro.

[[File:imagem5.png]]

Então simplificamos a expressão booleana,

[[File:imagem6.png]]

Determinando que q e p são verdadeiros, como nós já sabíamos que p era  verdadeiro, nós concluímos que q é verdadeiro.

A função '''bsimp''' é um simplificador geral para expressões booleanas construídas usando os operadores booleanos inativos.   A função computa uma expressão booleana simplificada equivalente ao seu argumento. Podemos também usar o maple para determinar se uma expressão é uma tautologia através do uso da função tautology oferecida pelo pacote logic.
[[File:imagem7.png]]

Agora mostramos como usar o maple para compreender algumas provas construtivas. Especificamente, vamos examinar como explorar a construção de uma lista sequencial  de números compostos.

[[File:imagem8.png]]

Enquanto o maple pode ser usado para gerar uma lista de n inteiros compostos consecutivos, gerados por prova, não é possível derivar a prova em si usando o maple. Devemos observar que este argumento não fornece o menor conjunto de n inteiros compostos consecutivos. Embora dado um inteiro positivo n, é possível usar o maple para encontrar a menor sequencia de n inteiros compostos consecutivos.
Através do maple abordaremos a prova não construtiva da existência de um número infinito de números primos. Por ser uma prova não construtiva, não podemos simplesmente criar um algoritmo para gerar um número primo maior, assumindo a existência de um número primo máximo. Embora a ideia chave da prova seja considerar a primalidade do inteiro $N! +1$, é possível que $N!+1$ seja por si só um número primo, porém mesmo que não seja seu maior fator primo deve ser maior que n. É possível encontrar o menor fator primo fatorando $N!+1$ diretamente, usando a rotina '''ifactor''' da biblioteca maple.

[[File:imagem9.png]]

Podemos observar pelo resultado que, enquanto alguns desses números são primos, outros não são, a partir disso, podemos fazer a leitura do menor fator primo.
Para determinar o menor fator primo de cada um desses inteiros, podemos escrever a seguinte rotina:

[[File:imagem10.png]]


Ela usa o procedimento {\bf factorset} do pacote {\bf numtheory} para computar o conjunto de fatores do inteiro de entrada, e então simplesmente seleciona seu menor membro.

[[File:imagem11.png]]

Agora confrontamos nosso exemplo final do uso do maple explorando teoremas matemáticos. Neste caso vamos explorar a conjectura de Goldbach: que é, todo inteiro par maior que 4 pode ser expressado como a soma de dois primos.

[[File:imagem12.png]]

Agora, criamos um procedimento para examinar o Goldbach a conjectura mais automaticamente.

[[File:imagem13.png]]

==Indução Matemática==

O maple pode ser usado para auxiliar na elaboração de provas de várias afirmações matemáticas usando a indução matemática. De fato, com o maple como seu assistente, você pode conduzir inteiramente o processo de descoberta e averiguação de forma intuitiva. É provável que entre os primeiros exemplos de indução matemática, encontremos  a verificação da fórmula  <math>1+2+3+...+n = n(n+1)/2</math>,a soma dos primeiros n inteiros positivos. O maple se adequa para provar fórmulas como essa porque os passos envolvidos numa prova indutiva incluem manipulação simbólica. É possível gerar uma grande quantidade de dados numéricos para serem examinado.
[[File:imagem14.png]]

Através da geração de um grande conjunto de dados numéricos de pouca compreensão, eventualmente será possível determinar a fórmula acima. A saída mostra que $n^2$ é um pouco menos que o dobro da soma dos n primeiros inteiros para os valores de n testados. Partindo de um padrão, é possível perceber que a fórmula correta é uma função quadrática de n, resolva para encontrar os coeficientes e então teste se esse procedimento produz a fórmula correta.
Uma técnica útil para experimentos semelhantes é gerar listas de pares que consista da sequencia que esteja interessado  e de várias possibilidades que você elabore. Para investigar a hipótese de que a fórmula seja quadrática, devemos começar gerando uma lista de pares semelhantes ao seguinte: 

[[File:imagem15.png]]

Para explorar se a soma é uma função quadrática de n, podemos inserir um quadrático genérico em n e solucionar para encontrar os coeficientes a,b e c:

[[File:imagem16.png]]

Precisamos de três equações para serem solucionadas em busca dos três coeficientes:

[[File:imagem17.png]]

Agora, instruímos o maple para resolver essas equações e encontrar os três coeficientes.

[[File:imagem18.png]]

A nossa fórmula original torna-se:

[[File:imagem19.png]]

Neste ponto, as habilidades do maple permitem manipular expressões simbolicamente para ajudar a construir uma prova indutiva.  Segue um exemplo de como a prova interativa da fórmula acima pode ser conduzido no maple por indução matemática.
O termo geral da soma é:

[[File:imagem20.png]]

Enquanto o lado direito da fórmula é:

[[File:imagem21.png]]

Podemos usar o procedimento '''subs''' para verificar o passo base da indução; neste caso o  passo base é <math>n=1</math>

[[File:imagem22.png]]

Os resultados coincidem (concordam), então o passo base é estabelecido. Para o passo indutivo, usamos que a fórmula seja válida para n=k.

[[File:imagem23.png]]

Para somar k+1 termos, computamos:

[[File:imagem.png]]

Por fim, a fórmula para n=k+1  é:

[[File:imagem24.png]]

Os resultados coincidem, então o passo indutivo é verificado. A fórmula agora segue por indução matemática. Assim, podemos concluir que enquanto o maple ainda não é capaz de construir provas inteiramente por conta própria, é uma ferramenta muito efetiva para ser usada na construção de provas interativas.
Agora, vamos considerar um exemplo mais complicado. A fórmula do somatório.
<math>S = 1.1! + 2.2! +$...$+n.n!</math>
É bem menos óbvio que o exemplo anterior. Para descobri-lo, iniciaremos gerando alguns dados numéricos.

[[File:imagem25.png]]

Se um padrão não é imediatamente óbvio,podemos auxiliar nossa intuição gerando uma sequencia paralela.

[[File:imagem26.png]]

Observando isto um pouco, fica óbvio que estamos no caminho certo, vamos apenas fazer alguns pequeno ajustes.

[[File:imagens.png]]

Desta evidencia, devemos inferir  a conjectura que a fórmula para nosso somatório é: <math>S=(n+1)!-1</math>

A prova indutiva pode ser conduzida assim como foi feito no primeiro exemplo.
[[File:imagem28.png]]

O passo base é:

[[File:entrenseind.png]]

O passo indutivo é:

[[File:imagem29.png]]

Usando um pouco de manipulação algébrica, vemos que as duas últimas fórmulas são iguais. Isso completa a prova via indução matemática. Concluímos que nosso palpite sobre a fórmula está correto.

==Definições Recursiva e Interativa==
As funções do maple podem ser definidas tanto processualmente (usando a função proc) como explicitamente (usando a notação <math>-></math>), cada um desses métodos envolve meios interativos e recursivos de definição. Iniciamos nosso estudo, usando a função <math>-></math> do maple. Se nós quiséssemos definir uma função polinomial <math>A(n)= 3n^3 + 41n^2- 3n + 101</math> nós usaríamos o seguinte comando:

[[File:imagem30.png]]

Agora, se quiséssemos definir uma função recursivamente, digamos: 
<math>b(n) =b(n-1)^2 +  2b(n-1) +6</math>, com a condição inicial <math>b(0)=2</math>, então informaríamos (inserir):

[[File:imagem31.png]]

Se quiséssemos ver uma sequencia de valores para a função b, podemos usar a função seq, para exibir as saídas de um dado intervalo de entradas.

[[File:imagem32.png]]

Agora, criaremos uma função similar a '''b''', chamada '''f1''', que encontrará os números Fibonacci.

[[File:image33.png]]

Enquanto a notação <math>-></math>  para funções é conveniente e intuitiva, ela não oferece todas as facilidades para melhoria da eficiência que estão disponíveis no uso do comando '''proc'''. Para forçar o maple a calcular esses valores de forma eficiente, usamos a opção '''remember''' para definições de procedimentos afetados pelo uso do '''proc'''. Esta opção exige que o maple lembre de qualquer valor para procedimento que já tenha sido computado através do armazenamento destes em uma tabela Fibonacci.

[[File:imagem34.png]]

Esse método processual engloba ambos os casos, o caso base (quando <math>n <=2</math>) e os casos indutivos (como na condição '''else'''). Adicionalmente, o procedimento tem a indicação da função remember, forçando o maple a manter controle sobre quais valores da função já foram encontrados, para que estes sejam diretamente verificados ao invés de serem novamente computados.

[[File:imagem35.png]]

Agora, para ilustrar a diferença em complexidade computacional, compararemos os métodos processual e “->” usando a função {\bf time} do maple.

[[File:imagem36.png]]

Então , fica claro que a função {\bf remember} pode fazer uma enorme diferença em complexidade de tempo. Outra maneira de melhorar a eficiência de uma função definida recursivamente é reescrevê-la para evitar o uso de recursão. Ao invés disso, reestruturamos para que use um algoritmo iterativo. Na construção de um algoritmo iterativo, os componentes chave consistem em criar uma forma de loop (um {\bf for} ou {\bf while}) que computará valores começando do menor para o maior. Este método de programação é chamado de {\bf bottom up}: onde os menores valores de uma sequencia são computadose então usados para valores maiores.

[[File:imagem37.png]]

Faça um contraste entre esse procedimento e o procedimento recursivo {\bf f2} que foi definido anteriormente.

[[File:imagem38.png]]

Ambos os casos base e passo recursivo são explicitamente definidos no corpo do procedimento. O algoritmo primeiro tenta computar diretamente o valor verdadeiro  e pede valores de subcasos conforme é solicitado. Este método de programação é conhecido como top-down, pois os valores maiores são computados através da quebra da entrada em partes menores e da combinação dos resultados. Similar a atravessar uma árvore binária. Perceba que o procedimento recursivo com a opção remember e o procedimento iterativo tem praticamente a mesma performance. Para os primeiros vinte números Fibonacci, obtemos:

[[File:imagem39.png]]

Que é comparável às vezes obtidas para '''f2'''. Note que a implementação puramente recursiva '''f1''' não possibilita o uso para f2(100). De fato, um bom exercício é mostrar que para fazê-lo, '''f2''' precisaria de aproximadamente

[[File:imagem40.png]]

Chamadas para operar todos os subcasos que surgirem. Mesmo a um bilhão de subcasos por segundo, isso exigiria mais de seis mil anos para sua conclusão.

==Computações e Explorações==

Nesta seção do material, exploraremos o modo como o maple pode ser usado para resolver as questões 4,5 e 8 da seção “computações e explorações” do livro.
1 -  Quantos pares de números primos podem ser encontrados?
Para determinar quantos pares de números primos existem, usaremos o pacote “numtheory” do maple, que contém as funções '''nextprime''', '''prevprime''' e '''ithprime'''

[[File:imagem41.png]]

Agora, após formada uma lista de primos, queremos extrair quaisquer pares de primos que ocorram nessa lista.

[[File:imagem42.png]]

Ao invés de dar como resultados os pares de números primos, o número de sequências de primos  pode indicar um padrão, então construímos uma lista dos {\bf i}’s que formarem pares.

[[File:imagem43.png]]

Parece não haver um padrão óbvio ocorrendo.
2- Determine quais números Fibonacci são divisíveis por 5, quais são divisíveis por 7 e quais são divisíveis por 111. Prove que suas conjecturas estão corretas.
Primeiro vamos gerar dados para trabalhá-los (manipulá-los).

[[File:imagem44.png]]

Queremos determinar os índices n para os quais o enésimo número Fibonacci é divisível por 5. Uma maneira de fazer isso é construindo uma lista, através de testes com os dados acima e adicionando à lista somente os índices n para os quais o teste retorne {\bf verdadeiro}.

[[File:imagem45.png]]

Isso constrói uma lista indicando quais entre os primeiros 500 números Fibonacci são múltiplos de 5. Os dados indicam que o enésimo número Fibonacci <math>F_n</math> é divisível por 5, somente se n é divisível por 5. Para obter evidências para a conversão, devemos testar se <math>F_5n</math> é divisível por 5, para tantos n quanto forem possíveis. Para que nosso teste seja conciso e ainda permita testar um grande intervalo(série) de valores, vamos implementá-lo de maneira que nenhum resultado seja produzido, a menos que seja encontrado um contra exemplo.

[[File:imagem46.png]]

Assim, não há contra exemplo entre os primeiros 5000 números Fibonacci. Você pode testar com valores maiores que 1000, para obter novas evidências. Outra abordagem ligeiramente diferentepode ser usada para encontrar os números Fibonacci divisíveis por um dado inteiro, neste caso, o número 7. Nós simplesmente construímos o teste de divisibilidade no comando para gerar dados.

[[File:imagem47.png]]

Podemos agora selecionar os índices dos pares cujo segundo membro seja igual a 0.

[[File:imagem48.png]]

Podemos perceber um padrão nesses dados, como o seguinte:

[[File:imagem49.png]]

Você pode tentar averiguar se esse padrão persiste, substituindo 500 na definição de '''fib_list''' por números muito maiores. (O teste da divisibilidade por 111, nós deixamos pra você).
2 – A notória conjectura <math>3x + 1</math> (também conhecida como conjectura de Collatz e por vários outros nomes) afirma que: independente de qual inteiro x você escolha para iniciar, em iteração com a função f(x), onde f(x) = x/2, se x é par e f(x) = 3x+1 se x é ímpar, sempre produz o inteiro 1. Cheque essa conjectura para tantos inteiros positivos possíveis.
Para inicar, precisamos definir a função, a qual examinaremos.

[[File:imagem50.png]]

Agora escrevemos uma função que fará a iteração da função Collatz até que o valor obtido seja igual a 1, nós incluímos uma variável “count” por dois motivos: Primeiro, queremos ter uma idéia de quanto tempo leva para que as iterações estabilizem; Segundo, já que não sabemos ao certo se as iterações vão estabilizar para um dado valor de entrada “seed”, nos codificamos um limite superior para o número de iterações a serem computadas.

[[File:imagem51.png]]

Para averiguar a conjectura para os 1000 primeiros inteiros, podemos usar a função {\bf IC} como no exemplo a seguir:

[[File:imagem52.png]]


Perceba que, o fato de a função ter eventualmente parado é a averiguação que buscamos.

===Exercícios e Projetos===

==Conclusão==

A elaboração deste trabalho nos proporcionou um novo conhecimento sobre a liguagem Maple e como podemos aplica-la no entendimento da Indução e Recursão na Matemática em Fundamentos da Matemática para a Computação II. Com isso podemos perceber o quanto é importante o aprendizado contínuo na área da Tecnologia da Informação para nos mantermos atualizados e em crescente aprimoração seja acadêmica ou profissional .

==Referências==
O desenvolvimento desse trabalho é proviniente da página Online Learning Center do \cite{boulic:91}, sendo a tradução do mesmo.