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 $$\&$$, por exemplo: usamos $$\&and$$ invés de $$and$$ e $$\&not$$ invés de $$not$$. Em seguida, vão alguns exemplos do uso de operadores booleanos inativos:



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:



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



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



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



Então simplificamos a expressão booleana,



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.



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.



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.



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:



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



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.



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



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  $$1+2+3+...+n = n(n+1)/2$$ ,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.



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:



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:



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



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



A nossa fórmula original torna-se:



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 é:



Enquanto o lado direito da fórmula é:



Podemos usar o procedimento subs para verificar o passo base da indução; neste caso o passo base é $$n=1$$



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$$.



Para somar k+1 termos, computamos:



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



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. $$S = 1.1! + 2.2! +...+n.n!$$ É bem menos óbvio que o exemplo anterior. Para descobri-lo, iniciaremos gerando alguns dados numéricos.



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



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



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

A prova indutiva pode ser conduzida assim como foi feito no primeiro exemplo.



O passo base é:



O passo indutivo é:



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 ->), cada um desses métodos envolve meios interativos e recursivos de definição. Iniciamos nosso estudo, usando a função -> do maple. Se nós quiséssemos definir uma função polinomial $$A(n)= 3n^3 + 41n^2- 3n + 101$$ nós usaríamos o seguinte comando:



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



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.



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



Enquanto a notação $$->$$ 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.



Esse método processual engloba ambos os casos, o caso base (quando $$n <=2$$) 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.



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



Então, fica claro que a função 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 for ou while) que computará valores começando do menor para o maior. Este método de programação é chamado de bottom up: onde os menores valores de uma sequencia são computadose então usados para valores maiores.



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



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:



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



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



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



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 i’s que formarem pares.



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).



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 verdadeiro.



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 $$F_n$$ é divisível por 5, somente se n é divisível por 5. Para obter evidências para a conversão, devemos testar se $$F_{5n}$$ é 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.



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 diferente pode 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.



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



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



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 $$3x + 1$$ (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.



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.



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



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

Exercícios e Projetos
1.Use Maple para encontrar e provar a fórmula do soma do primeiros $$k^N$$ (K elevado a enésima potência) inteiros positivos para N = 4, 5, 6, 7, 8, 9 e 10.

2.Use Maple para estudar a função McCarthy 911.

3.Escreva um procedimento no Maple para encontrar a menr (isto é, a primeira) sequência consecutiva de N inteiros positivos compostos, para um inteiro N positivo arbitrário.

4.Use Maple para desenvolver um procedimento para gerar números Ulam. Faça e estude numericamente as conjecturas sobre a distribuição desses números.

5.Escreva um procedimento no Maple que receba um inteiro K como entrada e determine se é ou não o produto dos primeiros K primos mais 1, e se é primo ou não, através da fatoração deste número.

6.Outra maneira de mostrar que existem infinitos primos é assumir que existem apenas N primos $$p_1, p_2, ..., p_n$$. Mas isso é uma contradição já que $$p_1, p_2, ..., p_{n+1}$$ tem ao menos um fator primo que não é divisível por $$p_i$$, i = 1, 2, ..., N. Encontre o menor fator primo de $$2, 3, ..., p_{n+1}$$ para todos os inteiros positivos N que não excedam 200. Para os quais N é este número primo.

7.O número Lucas satisfaz a recorrência $$L_{n} = L{n-1} + L_{n-2}$$ e as condições iniciais $$L_0 = 2$$ e $$L_1 = 1$$. Use o Maple para obter evidências para as conjecturas sobre a divisibilidade dos números Lucas por outros divisores inteiros diferentes.

8.A sequência $$a_1, a_2, a_3, ...$$ é chamada de periódica se existirem inteiros positivos N e p, para os quais $$a_n = a_{n+p}$$, para todo $$n \ge N$$. O menor inteiro p, para o qual isso é verdadeiro é chamado de período da sequência $$a_1, a_2, a_3, ...$$ se diz que é o módulo periódico m, para um inteiro positivo M, se a sequência $$a_1~ mod~ m$$ ,$$ a_2~ mod~ m$$, $$a_3~ mod~ m$$, ... é periódica. Use maple para determinar se a sequência Fibonacci é modulo periódico m, para vários inteiros m e, se for, encontre o período.Você pode, através da inspeção de valores de m diferentes o suficiente, fazer alguma conjectura a respeito da relação entre m e o período? Faça o mesmo para sequências que julgar interessante.

Extra
Nesta sessão serão tratados algumas demonstrações sobre indução e recursão contidas na pagina 240, cujo Extra 1, 2, 3 e 4 são implentações na linguagem de C++ como forma de uma abordagem alternativa como sugestão do Professor Umberto Rivieccio da disciplina de FMC II.


 * Exemplo 1

Use o princípio da indução matemática para provar que $$1 + 5 + 5^2 +5^3 +...+5n = \frac {5^{n+1}-1}{4} $$ para todo $$  n \ge 0$$

Seja p(n) a afirmação $$1+5^2+5^3+...+5n = \frac {5^{n+1} -1}{4}$$

Passo base: p(0) : $$1= \frac {5^{0+1}-1}{4}$$, (perceba que a soma no lado esquerdo de p(0) inicia e termina com o primeiro termo 1, consequentemente é apenas o primeiro termo 1). $$p(0)$$ é verdadeiro pois ambos os lados são iguais a 1.

Passo indutivo: $$P(k) \to p(k+1)$$: suponha para qualquer k, $$p(k)$$ é verdadeiro, isto é, $$1+5+5^2+5^3+...+ 5k = \frac {5^{k+1} -1}{4}$$ Precisamos mostrar que a próxima afirmação $$p(k+1)$$, é verdadeira: $$1+5+5^2 +5^3 +...+5^{k+1} = \frac {5^{k+2} -1}{4}$$ Para fazer isso, iniciamos com $$p(k)$$ e adicionamos o próximo termo, $$5^{k+1}$$, em ambos os lados, depois mostramos que essa é a afirmação $$p(k+1)$$.

$$1+5^2+5^3+...+5k = \frac {5^{k+1} -1}{4}$$

$$+5^{k+1} -- 5^{k+1}$$

$$1+5^2+5^3+...+5k+5^{k+1} = \frac {5^{k+1} -1 +4.5^{k+1}}{4}$$

$$=\frac {(1+4)5^{k+1} -1}{4} $$

$$=\frac {5.5^{k+1} -1}{4}$$

$$= \frac {5^{k+2} -1}{4}$$

Isto é, $$p(k+1)$$. Portanto uma afirmação verdadeira p(k, é seguida por outra afirmação verdadeira $$p(k+1)$$, por isso $$p(k) \to p(k+1)$$ é verdadeira. Assim, pelo princípio da indução matemática ,$$ p(n)$$ é verdadeiro para todo $$n\ge 0$$. Nota: Alternativamente, a prova de $$p(k) \to p(k+1)$$ pode ser escrita dessa forma. Começamos escrevendo o lado esquerdo da equação que é $$p(k+1)$$. Então mostramos que pode ser reescrita para dar o lado direito de $$p(k+1)$$. Perceba que a suposição que $$p(k)$$ é verdadeira está sendo usada em substituição do primeiro passo.

$$+5+5^2+5^3 + ...+ 5k +5^{k+1} = \frac {5^{k+1} -1}{4} + 5^{k+1}$$

$$ = \frac {5^{k+1} -1} + 4.5^{k+1}{4}$$

$$= \frac {(1+4).5^{k+1} -1}{4}$$

$$=\frac {5.5^{k+1}-1}{4}$$

$$=\frac {5^{k+2}-1}{4}$$


 * Exemplo 2

Prove que :

$$\sum_{k=1}^n (2k+3)= n(n+4)$$ para todo $$n\ge 1$$

Seja $$p(n)= 5+7+9+...+(2n+3)= n.(n+4)$$

Passo base: $$p(1)$$ afirma que $$2.1+3=1(1+4)$$, é verdadeiro, já que ambos os lados são iguais a 5.

Passo indutivo: $$p(k) \to p(k+1)$$: suponha que $$p(k)$$ é verdadeiro, isto é: $$5+7+9+...+(2k+3 )=k(k+4)$$ Adicionamos $$(2k+3)$$ a ambos os lados e simplificamos :

$$5+7+9+...+(2k+3)+(2k+5) = k.(k+4)+(2k+5)$$

$$= k^2 + 6k + 5$$

$$= (k+1).(k+5)$$ Que é $$p(k+1)$$, assim, $$p(k) \to p(k+1)$$ é verdadeiro. Portanto, pelo princípio da indução matemática,$$ 5+7+9+...+(2n+3)= n(n+4)$$ é verdadeiro para todo $$n \ge 1$$.


 * Exemplo 3

Encontre uma fórmula para $$(1 – \frac {1}{2}^2)$$ (1- 1/3²) (1 – ¼²)...(1 – 1/n²) Para n>=2, use o princípio de indução matemática para provar que sua fórmula está correta. Primeiro precisamos de um palpite para a fórmula do produto. Usando n= 2,3,4,5 obtemos:

(1-1/2²) = ¾ (1-1/2² ) (1-1/3²)= 3/4 . 8/9 = 2/3 (1-1/2²) (1-1/3²) (1-1/4²) = 3/4 . 8/9 . 15/16 = 5/8 (1-1/2²) (1-1/3²) (1-1/4²) (1-1/5²) = 3/4 . 8/9 . 15/16 . 24/25 = 3/5 Assim, os produtos são 3/4, 2/3, 5/8, 3/5. Podemos reescrever essas frações como 3/4, 4/6, 5/8, 6/10. Isso sugere n+1/2n como a forma geral da soma.

Seja p(n) : (1-1/2²) (1-1/3²) (1-1/4²)...(1-1/n²) = n+1/2n Agora, tentamos mostrar que p(n) é verdadeiro para todo n>=2.

Passo base: p(2) é verdadeiro. P(2) afirma que (1-1/2²) = 2+1/2.2, que é verdade pois, ambos os lados são iguais a ¾. Passo indutivo: P(k)p(k+1): suponha que p(k) é verdadeiro para algum k. portanto; (1-1/2²) (1-1/3²) (1-1/4²)...(1-1/k²) = k+1/2k Multiplica-se ambos os lados da equação por (1-1/(k+1)²) para obter (1-1/2²) (1-1/3²) (1-1/4²) ...(1-1/k²) (1-1/(k+1)²) = k+1/2k (1-1/(k+1)²)

= k+1/2k ((k+1)² - 1/ (k+1)²) = k+1/2k. k(k+2)/(k+1)² = k+2/2(k+1) Que é p(k+1). Portanto, pelo princípio da indução matemática, (1-1/2²) (1-1/3²) (1-1/4²)...(1-1/n²) = n+1/2n é verdadeiro para todo n >=2.

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
Discrete Mathematics and Its Applications. Disponível em:  Acesso em: 12 de Maio de 2016.