Рекурсивные функции — различия между версиями

Материал из Викиконспекты
Перейти к: навигация, поиск
(Работа со списками фиксированной длины)
(Теорема о рекурсии)
Строка 126: Строка 126:
  
 
=== Теорема о примитивной рекурсивности вычислимых функций ===
 
=== Теорема о примитивной рекурсивности вычислимых функций ===
==Теорема о рекурсии==
 
 
 
{{Теорема
 
{{Теорема
|id=th1
+
|statement= Если для  вычислимой функции <tex> F </tex> существует примитивно рекурсивная функция <tex> T </tex>, такая что для любого входа <tex> x </tex> максимальное количество  за которое будет посчитана <tex> F(x) </tex>  на MT равно  <tex> T(x) </tex>, то <tex> F </tex> примитивно рекурсивная функция.  
|about=О рекурсии
 
|statement= Пусть <tex>V(n, x)</tex> {{---}} вычислимая функция. Тогда найдётся такая вычислимая <tex>p</tex>, что <tex>\forall y</tex> <tex>p(y) = V(p, y)</tex>.
 
 
|proof=
 
|proof=
Приведем конструктивное доказательство теоремы.
+
Каждому состоянию MT поставим в соответствие список из четырех чисел [L,R,S,C],где
Пусть есть вычислимая <tex>V(x,y)</tex>. Будем поэтапно строить функцию <tex>p(y)</tex>. <br> Предположим, что у нас в распоряжении есть функция <tex>getSrc()</tex>, которая вернет код <tex>p(y)</tex>. Тогда саму <tex>p(y)</tex> можно переписать так:
+
<tex> L </tex> - состояние MT слева от головки ленты, представлено в виде числа в системы счисления с основанием равным алфавиту MT, число записано слева направо. Пробелу соответствует ноль, чтобы число было конечным.
 
+
<tex> R </tex> - состояние MT справа от головки, представлено аналогично <tex> L </tex> только число записано справа налево.   
<code><font size = "3em">
+
<tex> S </tex> - номер текущего состояния
p(y){
+
<tex> C </tex> - символ на который указывает головка ленты.
      V(x,y) {...}
+
Тогда всем переходам соответствует функция <tex> f </tex> принимающая состояние МТ и возвращающая новое состояние.
+
Покажем что она примитивно рекурсивная . При применении перехода в <tex> C </tex> записывается новый символ,затем из-за сдвига головки  а  <tex> L </tex> и <tex> R </tex>, в конец добавляется новая цифра или удаляется старая, затем в <tex> C </tex> записываетcя символ после сдвига, и в конце перехода в <tex> S </tex> записывается новое состояние автомата. Операции добавления в конец цифры или удаления последней цифры легко выражаются через деление и умножения, следовательно они примитивно рекурсивные. Все остальные операции являются простыми операциями над списками, а значит они тоже примитивно рекурсивные. Из этого следует что применения перехода {{---}} примитивно рекурсивная функция. В силу того что нужный переход можно выбрать используя конечное число функций if следует что и <tex> f </tex> тоже примитивно рекурсивная функция.
      main() {
 
          return V(getSrc(), y)
 
      }
 
 
 
      string getSrc() {...}
 
  }
 
</font></code>
 
Теперь нужно определить функцию <tex>getSrc()</tex>. Предположим, что внутри <tex>p(y)</tex> мы можем определить функцию <tex>getOtherSrc()</tex>, состоящую из одного оператора <tex>return</tex>, которая вернет весь предшествующий ей код. Тогда <tex>p(y)</tex> перепишется так.
 
<code><font size = "3em">
 
  p(y){
 
      V(x,y) {...}
 
 
      main() {
 
          return V(getSrc(), y)
 
      }
 
 
 
      string getSrc() {
 
          string src = getOtherSrc();
 
          return src + "string getOtherSrc() {" + "\n" + "return" + src + "\n" + "}";
 
      }
 
 
 
      string getOtherSrc() {...}
 
  }
 
</font></code>
 
 
 
Теперь <tex>getOtherSrc()</tex> определяется очевидным образом, и мы получаем '''итоговую версию''' функции <tex>p(y)</tex>
 
<code><font size = "3em">
 
  p(y){
 
      V(x,y) {...}
 
 
      main() {
 
          return V(getSrc(), y)
 
      }
 
 
 
      string getSrc() {
 
          string src = getOtherSrc();
 
          return src + "string getOtherSrc() {" + "\n" + "return" + src + "\n" + "}";
 
      }
 
 
 
      string getOtherSrc() {
 
          return "  p(y){            // Возвращаем весь предыдущий код
 
                    V(x,y) {...}
 
   
 
                      main() {
 
                          return V(getSrc(), y)
 
                      }
 
 
 
                      string getSrc() {
 
                          string src = getOtherSrc();
 
                          return src + "string getOtherSrc() {" + "\n" + "return" + src + "\n" + "}";
 
                  }";
 
      }
 
  }
 
</font></code>
 
 
 
}}
 
Если говорить неформально, теорема о рекурсии утверждает, что внутри программы можно использовать ее код. Это упрощает доказательство некоторых теорем.
 
 
 
Приведем так же альтернативую формулировку теоремы и альтернативное (неконструктивное) доказательство.
 
 
 
{{Теорема
 
 
 
|about=о неподвижной точке, Клини
 
|statement= Пусть <tex>U</tex> {{---}} [[Диагональный_метод|универсальная функция]], <tex>h</tex> {{---}} всюду определённая [[Вычислимые_функции|вычислимая функция]]. Тогда найдется такое <tex>n</tex>, что <tex>U_n=U_{h(n)}</tex>.
 
Другими словами: нельзя найти алгоритма, преобразующего про-
 
граммы, который бы по каждой программе давал другую (не эквива-
 
лентную ей).
 
|proof=
 
Начнём с доказательства леммы.
 
{{Лемма
 
|statement= Пусть на натуральных числах задано отношение эквивалентности <tex>\equiv</tex>. Тогда следующие два утверждения не могут быть выполнены одновременно: <br>
 
# Пусть <tex>f</tex> {{---}} вычислимая функция. Тогда существует всюду определённое вычислимое <tex>\equiv</tex> {{---}} продолжение <tex>g</tex> функции <tex>f</tex>, то есть такая <tex>g</tex>, что <tex>D(g)=N</tex> и <tex>\forall x</tex> такого, что <tex>f(x) \ne \perp</tex>, выполнено <tex>f(x) \equiv g(x)</tex>.
 
# Найдётся такая всюду определенная вычислимая <tex>h</tex>, что <tex>\forall n </tex> выполнено <tex>h(n) \not\equiv n</tex>.
 
|proof=
 
Приведем доказательство от противного. Пусть оба утверждения выполнены. <br>
 
Определим функцию <tex>f</tex> так: <tex>f(x)=U(x,x)</tex>. Заметим, что никакая всюду вычислимая функция не отличается от <tex>f</tex> всюду. <br> Согласно первому утверждению найдётся всюду определённое вычислимое <tex>\equiv</tex> {{---}} продолжение <tex>g</tex> функции <tex>f</tex>. <br> Определим функцию <tex>t</tex> так: <tex>t(x)=h(g(x))</tex>, где <tex>h</tex> {{---}} функция из второго утверждения. <br >Если <tex>f(x) \ne \perp</tex>, то <tex>f(x)=g(x) \ne h(g(x))=t(x)</tex>, то есть <tex>f(x) \ne t(x)</tex>. Если <tex>f(x)= \perp</tex>, то <tex>f(x) \ne t(x)</tex>, так как <tex>t</tex> всюду определена. Значит, <tex>f</tex> всюду отлична от <tex>t</tex>, получили противоречие.
 
}}
 
Теперь определим отношение <tex>\equiv</tex> так: <tex>x \equiv y \Leftrightarrow U_x = U_y</tex>. Покажем, что для него выполнено первое утверждение леммы. <br> Для заданной <tex>f</tex> определим <tex>V(n,x) = U(f(n), x)</tex>. <br> Так как <tex>U</tex> {{---}} универсальная функция, то найдётся такая всюду определенная вычислимая функция <tex>s</tex>, что <tex>V(n,x) = U(s(n), x)</tex>. <br> Тогда  <tex>\forall x </tex> и <tex> n </tex> будет выполнено <tex>U(f(n), x) = U(s(n), x)</tex>. Значит, <tex>\forall n </tex> <tex> s(n) \equiv f(n)</tex>, то есть <tex>s</tex> {{---}} всюду определенное <tex>\equiv</tex> {{---}} продолжение <tex>f</tex>.
 
Значит, для нашего отношения эквивалентности второе утверждение леммы не верно, то есть для любого вычислимого всюду определенного <tex>h</tex> <tex> \exists n</tex> такое, что <tex>U_{h(n)} = U_n</tex>.
 
 
}}
 
}}

Версия 00:03, 20 января 2013

Эта статья находится в разработке!

Все рассматриваемые здесь функции действуют из подмножества [math] \mathbb {N}^t [/math] в [math] \mathbb {N} [/math], где [math] t [/math] - любое натуральное число.Также будем считать что [math] 0[/math] натуральное число.

Примитивно рекурсивные функции

Основные определения

Рассмотрим следующие правила преобразования функций:

  • Рассмотрим [math] k [/math]-местную функцию [math] f(x_1,\ldots,x_k) [/math] и [math] k [/math] [math]n [/math]-местных функций [math] g_i(x_1,x_2,\ldots,x_n) [/math]. Тогда после преобразования у нас появится [math] n [/math] - местная функция [math] F = f(g_1(x_1,\ldots,x_n),\ldots, g_k(x_1,\ldots,x_n)) [/math].
Это правило называется правилом подстановки
  • Рассмотрим [math] k [/math]-местную функцию [math] f [/math] и [math] k + 2 [/math]-местную функцию [math] h [/math]. Тогда после преобразования у нас будет [math] k+1 [/math] -местная функция [math] g [/math], которая определена следующим образом:
[math]g(x_1,\ldots,x_n,0)=f(x_1,\ldots,x_n)[/math]
[math]g(x_1,\ldots,x_n,y+1)=h(x_1,\ldots,x_n,y,g(x_1,\ldots, x_n,y))[/math]
Это правило называется правилом рекурсии,при этом будем говорить что рекурсия запускается по аргументу [math] y [/math].


Определение:
Примитивно рекурсивными называют функции, которые можно получить с помощью правил подстановки и рекурсии из константной функции [math] \textbf 0 [/math], функции [math] I(x) = x + 1, [/math] и набора функций [math] P_{n,k}(x_1,\ldots,x_n) = x_k,[/math] где [math] k \le n [/math].

Заметим, что если [math] f [/math][math] n [/math]-местная примитивно рекурсивная функция, то она определена на всем множестве [math] \mathbb {N}^{n} [/math], так как [math] f [/math] получается путем правил преобразования из всюду определенных функций, и правила преобразование не портят всюду определенность. Говоря неформальным языком, рекурсивные функции напоминают программы, у которых при любых входных данных все циклы и рекурсий завершатся за конечное время.

Благодаря проекторам мы можем делать следующие преобразования:

  • В правиле подстановки можно использовать функции с разным числом аргументов. Например, подстановка [math] F(x,y) = f(g(y),h(x,x,y)) [/math] эквивалентна [math] F(x,y,z) = f(g(P_{2,2}(x,y)),h(P_{2,1}(x,y),P_{2,1}(x,y),P_{2,2}(x,y))) [/math], но если [math] F [/math] не константная функция то все подставляемые функции должны иметь хотя бы один аргумент.
  • В рекурсии не обязательно вести индукцию по последнему аргументу. Следует из того что мы можем с помощью проекторов поставить требуемый аргумент на последнее место.

В дальнейшем вместо [math] P_{n,k}(x_1,\ldots,x_k) [/math] будем писать просто [math] x_k [/math], подразумевая требуемое нам [math] n [/math].

Арифметические операции на примитивно рекурсивных функциях

n -местный ноль

[math] \textbf 0 [/math] - функция нуля аргументов.

Выразим сначала [math] \textbf 0^1 [/math]

[math] \textbf 0^{1}(0) = \textbf 0 [/math]

[math] \textbf 0^{1}(y+1) = h(y,\textbf 0^{1}(y)) [/math], где [math] h(x,y) = y [/math]

Теперь выразим [math] \textbf 0^n [/math]

[math] \textbf 0^{n}(x_1,\ldots,x_{n-1},0) = \textbf 0^{n-1} [/math]

[math] \textbf 0^{n}(x_1,\ldots,x_{n-1},y+1) = h(x_1,\ldots,x_{n-1},\textbf 0^{n}(y)) [/math], где [math] h(x_1,\ldots, x_n,y) = y [/math]

Константа [math] \textbf M [/math] равна [math] I(\textbf{M-1}) [/math]

[math] \textbf M^n [/math] - n местная константа, получается аналогичным к [math] \textbf 0^n [/math] образом.

Сложения

[math] sum(x,0) = x [/math]

[math] sum(x,y+1) = h(x,y,sum(x,y)) [/math] , где [math] h(x,y,z)=I(z) [/math]

Умножения

[math] prod(x,0) = \textbf 0^1(x) [/math]

[math] prod(x,y+1) = h(x,y,prod(x,y)) [/math], где [math] h(x,y,z)=sum(x,z) [/math]

Вычитания

Если [math] x \lt y [/math], то [math] sub(x,y) = 0 [/math] , иначе [math] sub(x,y) = x - y [/math].

Рассмотрим сначала вычитания единицы [math] sub_{1}(x) = x - 1 [/math]

[math] sub_1(0) = \textbf 0 [/math]

[math] sub_1(x+1) = h(x,sub_1(x)) [/math], где [math] h(x,y) = y [/math]

Теперь рассмотрим [math] sub(x,y) [/math]

[math] sub(x,0) = x [/math]

[math] sub(x,y+1) = h(x,y,sub(x,y)) [/math], где [math] h(x,y,z) =sub_1(z) [/math]

Операции сравнения

[math] eq(x,y) = 1 [/math] если [math] x = y [/math], иначе [math] eq(x,y) = 0 [/math]

[math] le(x,y) = 1 [/math] если [math] x \le y [/math], иначе [math] lq(x,y) = 0 [/math] [math] lower(x,y) = 1 [/math] если [math] x \lt y [/math], иначе [math] lower(x,y) = 0 [/math]

Сначала выразим [math] eq_{0}(x) = eq(x,0) [/math]

[math] eq_0(0) =I(\textbf 0) [/math]

[math] eq_0(y+1) = h(y,eq(y)) [/math] , где [math] h(y,eq(y)) = \textbf 0^2(x,y) [/math]

Теперь все остальные функции

[math] le(x,y) = eq_0(sub(x,y)) [/math]

[math] eq(x,y) = mul(le(x,y),le(y,x)) [/math]

[math] lower(x,y) = mul(le(x,y),le(I(x),y)) [/math]

IF

[math] if(0,x,y) = y [/math]

[math] if(c+1,x,y) = h(c,x,y,if(c,x,y)) [/math] , где [math] h(c,x,y,d) = x [/math]

Деление

[math] divide(x,y) = \lfloor {\frac{x}{y}} \rfloor [/math], если [math] y \gt 0 [/math]. Если же [math] y = 0 [/math], то [math] divide(x,0) [/math] и все связанные с делением функции равны каким то ,не интересными для нас, числами.

Сначала определим [math] divmax(x,y) [/math] — функция равна максимальному числу меньшему [math] x [/math] и которое нацело делится на [math] y [/math].

[math] divmax(0,y) =\textbf 0^{1} [/math]

[math] divmax(x+1,y) = h(x,y,divmax(x,y)) [/math], где [math] h(x,y,z) = if(eq(sub(I(x),z),y),I(x),z) [/math],

или не формально если [math] x+1 - y = z [/math] то [math] h(x,y,z) = x+1 [/math], иначе [math] h(x,y,z) = y [/math]

Теперь само деления

[math] divide(0,y) =\textbf 0^{1} [/math]

[math] divide(x,y) = h(x,y,divide(x,y)) [/math], где [math] h(x,y,z) = sum(z,eq(I(x),divmax(I(x),y))) [/math]

или не формально если [math] x+1~\vdots z [/math] то [math] h(x,y,z) = z+1 [/math], иначе [math] h(x,y,z) = z [/math]

Остаток от деления выражается так:

[math] mod(x,y) = sub(x,mul(y,divide(x,y))) [/math]

Работа со списками фиксированной длины

С помощью описанных выше арифметических операций можно выразить проверку на простоту числа и поиск [math] n [/math] - того простого числа. Рассмотрим список из натуральны чисел [math] [x_1,\ldots,x_n] [/math], тогда ему в соответствия можно поставить число [math] p_1^{x_1+1} \cdot p_2^{x_2+1} \cdot \ldots \cdot p_n^{x_n+1} [/math], где [math] p_i - i[/math]-тое простое число. Как видно из представления,создания списка, взятие [math] i [/math] - того элемента и остальные операции являются простыми арифметическими операциями, а следовательно примитивно рекурсивными. Поэтому будем считать что у примитивно рекурсивной функций аргументы и результат могут быть списками из натуральных чисел.

Теорема о примитивной рекурсивности вычислимых функций

Теорема:
Если для вычислимой функции [math] F [/math] существует примитивно рекурсивная функция [math] T [/math], такая что для любого входа [math] x [/math] максимальное количество за которое будет посчитана [math] F(x) [/math] на MT равно [math] T(x) [/math], то [math] F [/math] примитивно рекурсивная функция.
Доказательство:
[math]\triangleright[/math]

Каждому состоянию MT поставим в соответствие список из четырех чисел [L,R,S,C],где [math] L [/math] - состояние MT слева от головки ленты, представлено в виде числа в системы счисления с основанием равным алфавиту MT, число записано слева направо. Пробелу соответствует ноль, чтобы число было конечным. [math] R [/math] - состояние MT справа от головки, представлено аналогично [math] L [/math] только число записано справа налево. [math] S [/math] - номер текущего состояния [math] C [/math] - символ на который указывает головка ленты. Тогда всем переходам соответствует функция [math] f [/math] принимающая состояние МТ и возвращающая новое состояние.

Покажем что она примитивно рекурсивная . При применении перехода в [math] C [/math] записывается новый символ,затем из-за сдвига головки а [math] L [/math] и [math] R [/math], в конец добавляется новая цифра или удаляется старая, затем в [math] C [/math] записываетcя символ после сдвига, и в конце перехода в [math] S [/math] записывается новое состояние автомата. Операции добавления в конец цифры или удаления последней цифры легко выражаются через деление и умножения, следовательно они примитивно рекурсивные. Все остальные операции являются простыми операциями над списками, а значит они тоже примитивно рекурсивные. Из этого следует что применения перехода — примитивно рекурсивная функция. В силу того что нужный переход можно выбрать используя конечное число функций if следует что и [math] f [/math] тоже примитивно рекурсивная функция.
[math]\triangleleft[/math]