Primitive recursive function

In computability theory, a primitive recursive function is, roughly speaking, a function that can be computed by a computer program whose loops are all "for" loops (that is, an upper bound of the number of iterations of every loop is fixed before entering the loop). Primitive recursive functions form a strict subset of those general recursive functions that are also total functions.

The importance of primitive recursive functions lies in the fact that most computable functions that are studied in number theory (and more generally in mathematics) are primitive recursive. For example, addition and division, the factorial and exponential function, and the function which returns the nth prime are all primitive recursive. In fact, for showing that a computable function is primitive recursive, it suffices to show that its time complexity is bounded above by a primitive recursive function of the input size. It is hence not particularly easy to devise a computable function that is not primitive recursive; some examples are shown in section below.

The set of primitive recursive functions is known as PR in computational complexity theory.

Definition
A primitive recursive function takes a fixed number of arguments, each a natural number (nonnegative integer: {0, 1, 2, ...}), and returns a natural number. If it takes n arguments it is called n-ary.

The basic primitive recursive functions are given by these axioms:

More complex primitive recursive functions can be obtained by applying the operations given by these axioms:

The primitive recursive functions are the basic functions and those obtained from the basic functions by applying these operations a finite number of times.

Addition
A definition of the 2-ary function $$Add$$, to compute the sum of its arguments, can be obtained using the primitive recursion operator $$\rho$$. To this end, the well-known equations
 * $$\begin{array}{rcll}

0+y & = & y & \text{ and} \\ S(x)+y & = & S(x+y) &. \\ \end{array}$$ are "rephrased in primitive recursive function terminology": In the definition of $$\rho(g,h)$$, the first equation suggests to choose $$g = P_1^1$$ to obtain $$Add(0,y) = g(y) = y$$; the second equation suggests to choose $$h = S \circ P_2^3$$ to obtain $$Add(S(x),y) = h(x,Add(x,y),y) = (S \circ P_2^3)(x,Add(x,y),y) = S(Add(x,y))$$. Therefore, the addition function can be defined as $$Add = \rho(P_1^1,S \circ P_2^3)$$. As a computation example,
 * $$\begin{array}{lll}

& Add(1,7) \\ = & \rho(P_1^1,S \circ P_2^3) \; (S(0),7) & \text{ by Def. } Add, S \\ = & (S \circ P_2^3)(0,Add(0,7),7) & \text{ by case } \rho(g,h) \; (S(...),...) \\ = & S(Add(0,7)) & \text{ by Def. } \circ, P_2^3 \\ = & S( \; \rho(P_1^1,S \circ P_2^3) \; (0,7) \; ) & \text{ by Def. } Add \\ = & S(P_1^1(7)) & \text{ by case } \rho(g,h) \; (0,...) \\ = & S(7) & \text{ by Def. } P_1^1 \\ = & 8 & \text{ by Def. } S. \\ \end{array}$$

Doubling
Given $$Add$$, the 1-ary function $$Add \circ (P_1^1,P_1^1)$$ doubles its argument, $$(Add \circ (P_1^1,P_1^1))(x) = Add(x,x) = x+x$$.

Multiplication
In a similar way as addition, multiplication can be defined by $$Mul = \rho(C_0^1,Add \circ(P_2^3,P_3^3))$$. This reproduces the well-known multiplication equations:
 * $$\begin{array}{lll}

& Mul(0,y) \\ = & \rho(C_0^1,Add \circ(P_2^3,P_3^3)) \; (0,y) & \text{ by Def. } Mul \\ = & C_0^1(y) & \text{ by case } \rho(g,h) \; (0,...)\\ = & 0 & \text{ by Def. } C_0^1. \\ \end{array}$$ and
 * $$\begin{array}{lll}

& Mul(S(x),y) \\ = & \rho(C_0^1,Add \circ(P_2^3,P_3^3)) \; (S(x),y) & \text{ by Def. } Mul \\ = & (Add \circ(P_2^3,P_3^3)) \; (x,Mul(x,y),y) & \text{ by case } \rho(g,h) \; (S(...),...) \\ = & Add(Mul(x,y),y) & \text{ by Def. } \circ, P_2^3, P_3^3 \\ = & Mul(x,y)+y & \text{ by property of } Add. \\ \end{array}$$

Predecessor
The predecessor function acts as the "opposite" of the successor function and is recursively defined by the rules $$Pred(0) = 0$$ and $$Pred(S(n)) = n$$. A primitive recursive definition is $$Pred = \rho(C_0^0, P_1^2)$$. As a computation example,
 * $$\begin{array}{lll}

& Pred(8) \\ = & \rho(C_0^0, P_1^2) \; (S(7)) & \text{ by Def. } Pred, S \\ = & P_1^2(7,Pred(7)) & \text{ by case } \rho(g,h) \; (S(...),...) \\ = & 7 & \text{ by Def. } P_1^2. \\ \end{array}$$

Truncated subtraction
The limited subtraction function (also called "monus", and denoted "$$\stackrel.-$$") is definable from the predecessor function. It satisfies the equations
 * $$\begin{array}{rcll}

y \stackrel.- 0 & = & y & \text{and} \\ y \stackrel.- S(x) & = & Pred(y \stackrel.- x) &. \\ \end{array}$$ Since the recursion runs over the second argument, we begin with a primitive recursive definition of the reversed subtraction, $$RSub(y,x) = x \stackrel.- y$$. Its recursion then runs over the first argument, so its primitive recursive definition can be obtained, similar to addition, as $$RSub = \rho(P_1^1, Pred \circ P_2^3)$$. To get rid of the reversed argument order, then define $$Sub = RSub \circ (P_2^2,P_1^2)$$. As a computation example,
 * $$\begin{array}{lll}

& Sub(8,1) \\ = & (RSub \circ (P_2^2,P_1^2)) \; (8,1) & \text{ by Def. } Sub \\ = & RSub(1,8) & \text{ by Def. } \circ, P_2^2, P_1^2 \\ = & \rho(P_1^1, Pred \circ P_2^3) \; (S(0),8) & \text{ by Def. } RSub, S \\ = & (Pred \circ P_2^3) \; (0,RSub(0,8),8) & \text{ by case } \rho(g,h) \; (S(...),...) \\ = & Pred(RSub(0,8)) & \text{ by Def. } \circ, P_2^3 \\ = & Pred( \; \rho(P_1^1, Pred \circ P_2^3) \; (0,8) \; ) & \text{ by Def. } RSub \\ = & Pred(P_1^1(8)) & \text{ by case } \rho(g,h) \; (0,...) \\ = & Pred(8) & \text{ by Def. } P_1^1 \\ = & 7 & \text{ by property of } Pred. \\ \end{array}$$

Converting predicates to numeric functions
In some settings it is natural to consider primitive recursive functions that take as inputs tuples that mix numbers with truth values (that is $$t$$ for true and $$f$$ for false), or that produce truth values as outputs. This can be accomplished by identifying the truth values with numbers in any fixed manner. For example, it is common to identify the truth value $$t$$ with the number $$1$$ and the truth value $$f$$ with the number $$0$$. Once this identification has been made, the characteristic function of a set $$A$$, which always returns $$1$$ or $$0$$, can be viewed as a predicate that tells whether a number is in the set $$A$$. Such an identification of predicates with numeric functions will be assumed for the remainder of this article.

Predicate "Is zero"
As an example for a primitive recursive predicate, the 1-ary function $$IsZero$$ shall be defined such that $$IsZero(x) = 1$$ if $$x = 0$$, and $$IsZero(x) = 0$$, otherwise. This can be achieved by defining $$IsZero = \rho(C_1^0,C_0^2)$$. Then, $$IsZero(0) = \rho(C_1^0,C_0^2)(0) = C_1^0(0) = 1$$ and e.g. $$IsZero(8) = \rho(C_1^0,C_0^2)(S(7)) = C_0^2(7,IsZero(7)) = 0$$.

Predicate "Less or equal"
Using the property $$x \leq y \iff x \stackrel.- y = 0$$, the 2-ary function $$Leq$$ can be defined by $$Leq = IsZero \circ Sub$$. Then $$Leq(x,y) = 1$$ if $$x \leq y$$, and $$Leq(x,y) = 0$$, otherwise. As a computation example,
 * $$\begin{array}{lll}

& Leq(8,3) \\ = & IsZero(Sub(8,3)) & \text{ by Def. } Leq \\ = & IsZero(5) & \text{ by property of } Sub \\ = & 0 & \text{ by property of } IsZero \\ \end{array}$$

Predicate "Greater or equal"
Once a definition of $$Leq$$ is obtained, the converse predicate can be defined as $$Geq = Leq \circ (P_2^2,P_1^2)$$. Then, $$Geq(x,y) = Leq(y,x)$$ is true (more precisely: has value 1) if, and only if, $$x \geq y$$.

If-then-else
The 3-ary if-then-else operator known from programming languages can be defined by $$\textit{If} = \rho(P_2^2,P_3^4)$$. Then, for arbitrary $$x$$,
 * $$\begin{array}{lll}

& \textit{If}(S(x),y,z) \\ = & \rho(P_2^2,P_3^4) \; (S(x),y,z) & \text{ by Def. } \textit{If} \\ = & P_3^4(x,\textit{If}(x,y,z),y,z) & \text{ by case } \rho(S(...),...) \\ = & y & \text{ by Def. } P_3^4 \\ \end{array}$$ and
 * $$\begin{array}{lll}

& \textit{If}(0,y,z) \\ = & \rho(P_2^2,P_3^4) \; (0,y,z) & \text{ by Def. } \textit{If} \\ = & P_2^2(y,z) & \text{ by case } \rho(0,...) \\ = & z & \text{ by Def. } P_2^2. \\ \end{array}$$. That is, $$\textit{If}(x,y,z)$$ returns the then-part, $$y$$, if the if-part, $$x$$, is true, and the else-part, $$z$$, otherwise.

Junctors
Based on the $$\textit{If}$$ function, it is easy to define logical junctors. For example, defining $$And = \textit{If} \circ (P_1^2,P_2^2,C_0^2)$$, one obtains $$And(x,y) = \textit{If}(x,y,0)$$, that is, $$And(x,y)$$ is true if, and only if, both $$x$$ and $$y$$ are true (logical conjunction of $$x$$ and $$y$$).

Similarly, $$Or = \textit{If} \circ (P_1^2,C_1^2,P_2^2)$$ and $$Not = \textit{If} \circ (P_1^1,C_0^1,C_1^1)$$ lead to appropriate definitions of disjunction and negation: $$Or(x,y) = \textit{If}(x,1,y)$$ and $$Not(x) = \textit{If}(x,0,1)$$.

Equality predicate
Using the above functions $$Leq$$, $$Geq$$ and $$And$$, the definition $$Eq = And \circ (Leq, Geq)$$ implements the equality predicate. In fact, $$Eq(x,y) = And( Leq(x,y), Geq(x,y) )$$ is true if, and only if, $$x$$ equals $$y$$.

Similarly, the definition $$Lt = Not \circ Geq$$ implements the predicate "less-than", and $$Gt = Not \circ Leq$$ implements "greater-than".

Other operations on natural numbers
Exponentiation and primality testing are primitive recursive. Given primitive recursive functions $$e$$, $$f$$, $$g$$, and $$h$$, a function that returns the value of $$g$$ when $$e \leq f$$ and the value of $$h$$ otherwise is primitive recursive.

Operations on integers and rational numbers
By using Gödel numberings, the primitive recursive functions can be extended to operate on other objects such as integers and rational numbers. If integers are encoded by Gödel numbers in a standard way, the arithmetic operations including addition, subtraction, and multiplication are all primitive recursive. Similarly, if the rationals are represented by Gödel numbers then the field operations are all primitive recursive.

Some common primitive recursive functions
The following examples and definitions are from Kleene (1952) pp. 223–231. Many appear with proofs. Most also appear with similar names, either as proofs or as examples, in Boolos-Burgess-Jeffrey 2002 pp. 63–70; they add the logarithm lo(x, y) or lg(x, y) depending on the exact derivation.

In the following the mark " ' ", e.g. a', is the primitive mark meaning "the successor of", usually thought of as " +1", e.g. a +1 =def a'. The functions 16–20 and #G are of particular interest with respect to converting primitive recursive predicates to, and extracting them from, their "arithmetical" form expressed as Gödel numbers.


 * Addition: a+b
 * Multiplication: a×b
 * Exponentiation: ab
 * Factorial a! : 0! = 1, a'! = a!×a'
 * pred(a): (Predecessor or decrement): If a > 0 then a−1 else 0
 * Proper subtraction a ∸ b: If a ≥ b then a−b else 0
 * Minimum(a1, ... an)
 * Maximum(a1, ... an)
 * Absolute difference: | a−b | =def (a ∸ b) + (b ∸ a)
 * ~sg(a): NOT[signum(a)]: If a=0 then 1 else 0
 * sg(a): signum(a): If a=0 then 0 else 1
 * a | b: (a divides b): If b=k×a for some k then 0 else 1
 * Remainder(a, b): the leftover if b does not divide a "evenly". Also called MOD(a, b)
 * a = b: sg | a − b | (Kleene's convention was to represent true by 0 and false by 1; presently, especially in computers, the most common convention is the reverse, namely to represent true by 1 and false by 0, which amounts to changing sg into ~sg here and in the next item)
 * a < b: sg( a' ∸ b )
 * Pr(a): a is a prime number Pr(a) =def a>1 & NOT(Exists c)1 1 then the greatest x such that bx | a else 0


 * In the following, the abbreviation x =def x1, ... xn; subscripts may be applied if the meaning requires.


 * #A: A function φ definable explicitly from functions Ψ and constants q1, ... qn is primitive recursive in Ψ.
 * #B: The finite sum Σy<z ψ(x, y) and product Πy<zψ(x, y) are primitive recursive in ψ.
 * #C: A predicate P obtained by substituting functions χ1,..., χm for the respective variables of a predicate Q is primitive recursive in χ1,..., χm, Q.
 * #D: The following predicates are primitive recursive in Q and R:
 * NOT_Q(x).
 * Q OR R: Q(x) V R(x),
 * Q AND R: Q(x) & R(x),
 * Q IMPLIES R: Q(x) → R(x)
 * Q is equivalent to R: Q(x) ≡ R(x)


 * #E: The following predicates are primitive recursive in the predicate R:
 * (Ey)y<z R(x, y) where (Ey)y<z denotes "there exists at least one y that is less than z such that"
 * (y)y<z R(x, y) where (y)y<z denotes "for all y less than z it is true that"
 * μyy<z R(x, y). The operator μyy<z R(x, y) is a bounded form of the so-called minimization- or mu-operator: Defined as "the least value of y less than z such that R(x, y) is true; or z if there is no such value."


 * #F: Definition by cases: The function defined thus, where Q1, ..., Qm are mutually exclusive predicates (or "ψ(x) shall have the value given by the first clause that applies), is primitive recursive in φ1, ..., Q1, ... Qm:
 * φ(x) =
 * φ1(x) if Q1(x) is true,
 * φm(x) if Qm(x) is true
 * φm+1(x) otherwise
 * φm+1(x) otherwise


 * #G: If φ satisfies the equation:
 * φ(y,x) = χ(y, COURSE-φ(y; x2, ... xn ), x2, ... xn then φ is primitive recursive in χ. The value COURSE-φ(y; x2 to n ) of the course-of-values function encodes the sequence of values φ(0,x2 to n), ..., φ(y-1,x2 to n) of the original function.

Relationship to recursive functions
The broader class of partial recursive functions is defined by introducing an unbounded search operator. The use of this operator may result in a partial function, that is, a relation with at most one value for each argument, but does not necessarily have any value for any argument (see domain). An equivalent definition states that a partial recursive function is one that can be computed by a Turing machine. A total recursive function is a partial recursive function that is defined for every input.

Every primitive recursive function is total recursive, but not all total recursive functions are primitive recursive. The Ackermann function A(m,n) is a well-known example of a total recursive function (in fact, provable total), that is not primitive recursive. There is a characterization of the primitive recursive functions as a subset of the total recursive functions using the Ackermann function. This characterization states that a function is primitive recursive if and only if there is a natural number m such that the function can be computed by a Turing machine that always halts within A(m,n) or fewer steps, where n is the sum of the arguments of the primitive recursive function.

An important property of the primitive recursive functions is that they are a recursively enumerable subset of the set of all total recursive functions (which is not itself recursively enumerable). This means that there is a single computable function f(m,n) that enumerates the primitive recursive functions, namely: f can be explicitly constructed by iteratively repeating all possible ways of creating primitive recursive functions. Thus, it is provably total. One can use a diagonalization argument to show that f is not recursive primitive in itself: had it been such, so would be h(n) = f(n,n)+1. But if this equals some primitive recursive function, there is an m such that h(n) = f(m,n) for all n, and then h(m) = f(m,m), leading to contradiction.
 * For every primitive recursive function g, there is an m such that g(n) = f(m,n) for all n, and
 * For every m, the function h(n) = f(m,n) is primitive recursive.

However, the set of primitive recursive functions is not the largest recursively enumerable subset of the set of all total recursive functions. For example, the set of provably total functions (in Peano arithmetic) is also recursively enumerable, as one can enumerate all the proofs of the theory. While all primitive recursive functions are provably total, the converse is not true.

Limitations
Primitive recursive functions tend to correspond very closely with our intuition of what a computable function must be. Certainly the initial functions are intuitively computable (in their very simplicity), and the two operations by which one can create new primitive recursive functions are also very straightforward. However, the set of primitive recursive functions does not include every possible total computable function—this can be seen with a variant of Cantor's diagonal argument. This argument provides a total computable function that is not primitive recursive. A sketch of the proof is as follows:

This argument can be applied to any class of computable (total) functions that can be enumerated in this way, as explained in the article Machine that always halts. Note however that the partial computable functions (those that need not be defined for all arguments) can be explicitly enumerated, for instance by enumerating Turing machine encodings.

Other examples of total recursive but not primitive recursive functions are known:
 * The function that takes m to Ackermann(m,m) is a unary total recursive function that is not primitive recursive.
 * The Paris–Harrington theorem involves a total recursive function that is not primitive recursive.
 * The Sudan function
 * The Goodstein function

Constant functions
Instead of $$C_n^k$$, alternative definitions use just one 0-ary zero function $$C_0^0$$ as a primitive function that always returns zero, and built the constant functions from the zero function, the successor function and the composition operator.

Weak primitive recursion
The 1-place predecessor function is primitive recursive, see section. Fischer, Fischer & Beigel removed the implicit predecessor from the recursion rule, replacing it by the weaker rule

\begin{array}{lcl} f (0, x_1, \ldots, x_k) & = & g (x_1, \ldots, x_k) & \text{and} \\ f (S(y), x_1, \ldots, x_k) & = & h (S(y), f (y, x_1, \ldots, x_k), x_1, \ldots, x_k) \end{array} $$ They proved that the predecessor function still could be defined, and hence that "weak" primitive recursion also defines the primitive recursive functions.

Iterative functions
Weakening this even further by using functions $$h$$ of arity k+1, removing $$y$$ and $$S(y)$$ from the arguments of $$h$$ completely, we get the iteration rule:

$$\begin{array}{lcll} f(0,x_1,\ldots,x_k) & = & g(x_1,\ldots,x_k) &\textrm{and} \\ f(S(y),x_1,\ldots,x_k) & = & h(f(y,x_1,\ldots,x_k),x_1,\ldots,x_k) \end{array}$$

The class of iterative functions is defined the same way as the class of primitive recursive functions except with this weaker rule. These are conjectured to be a proper subset of the primitive recursive functions.

Additional primitive recursive forms
Some additional forms of recursion also define functions that are in fact primitive recursive. Definitions in these forms may be easier to find or more natural for reading or writing. Course-of-values recursion defines primitive recursive functions. Some forms of mutual recursion also define primitive recursive functions.

The functions that can be programmed in the LOOP programming language are exactly the primitive recursive functions. This gives a different characterization of the power of these functions. The main limitation of the LOOP language, compared to a Turing-complete language, is that in the LOOP language the number of times that each loop will run is specified before the loop begins to run.

Computer language definition
An example of a primitive recursive programming language is one that contains basic arithmetic operators (e.g. + and −, or ADD and SUBTRACT), conditionals and comparison (IF-THEN, EQUALS, LESS-THAN), and bounded loops, such as the basic for loop, where there is a known or calculable upper bound to all loops (FOR i FROM 1 TO n, with neither i nor n modifiable by the loop body). No control structures of greater generality, such as while loops or IF-THEN plus GOTO, are admitted in a primitive recursive language.

The LOOP language, introduced in a 1967 paper by Albert R. Meyer and Dennis M. Ritchie, is such a language. Its computing power coincides with the primitive recursive functions. A variant of the LOOP language is Douglas Hofstadter's BlooP in Gödel, Escher, Bach. Adding unbounded loops (WHILE, GOTO) makes the language general recursive and Turing-complete, as are all real-world computer programming languages.

The definition of primitive recursive functions implies that their computation halts on every input (after a finite number of steps). On the other hand, the halting problem is undecidable for general recursive functions.

Finitism and consistency results
The primitive recursive functions are closely related to mathematical finitism, and are used in several contexts in mathematical logic where a particularly constructive system is desired. Primitive recursive arithmetic (PRA), a formal axiom system for the natural numbers and the primitive recursive functions on them, is often used for this purpose.

PRA is much weaker than Peano arithmetic, which is not a finitistic system. Nevertheless, many results in number theory and in proof theory can be proved in PRA. For example, Gödel's incompleteness theorem can be formalized into PRA, giving the following theorem:
 * If T is a theory of arithmetic satisfying certain hypotheses, with Gödel sentence GT, then PRA proves the implication Con(T)→GT.

Similarly, many of the syntactic results in proof theory can be proved in PRA, which implies that there are primitive recursive functions that carry out the corresponding syntactic transformations of proofs.

In proof theory and set theory, there is an interest in finitistic consistency proofs, that is, consistency proofs that themselves are finitistically acceptable. Such a proof establishes that the consistency of a theory T implies the consistency of a theory S by producing a primitive recursive function that can transform any proof of an inconsistency from S into a proof of an inconsistency from T. One sufficient condition for a consistency proof to be finitistic is the ability to formalize it in PRA. For example, many consistency results in set theory that are obtained by forcing can be recast as syntactic proofs that can be formalized in PRA.

History
Recursive definitions had been used more or less formally in mathematics before, but the construction of primitive recursion is traced back to Richard Dedekind's theorem 126 of his Was sind und was sollen die Zahlen? (1888). This work was the first to give a proof that a certain recursive construction defines a unique function.

Primitive recursive arithmetic was first proposed by Thoralf Skolem in 1923.

The current terminology was coined by Rózsa Péter (1934) after Ackermann had proved in 1928 that the function which today is named after him was not primitive recursive, an event which prompted the need to rename what until then were simply called recursive functions.