# FRACTRAN

### From Seo Wiki - Search Engine Optimization and Programming Languages

**FRACTRAN** is a Turing-complete esoteric programming language invented by the mathematician John Conway. A FRACTRAN program is an ordered list of positive fractions together with an initial positive integer input *n*. The program is run by updating the integer *n* as follows:

- for the first fraction
*f*in the list for which*nf*is an integer, replace*n*by*nf* - repeat this rule until no fraction in the list produces an integer when multiplied by
*n*, then halt.

Conway gave an interesting formula for primes in FRACTRAN:

- <math>\left( \frac{17}{91}, \frac{78}{85}, \frac{19}{51}, \frac{23}{38}, \frac{29}{33}, \frac{77}{29}, \frac{95}{23}, \frac{77}{19}, \frac{1}{17}, \frac{11}{13}, \frac{13}{11}, \frac{15}{14}, \frac{15}{2}, \frac{55}{1} \right)</math>

Starting with *n*=2, this FRACTRAN program generates the following sequence of integers:

- 2, 15, 825, 725, 1925, 2275, 425, 390, 330, 290, 770, ... Template:OEIS

After 2, this sequence contains the following powers of 2:

- <math>2^2=4,\, 2^3=8,\, 2^5=32,\, 2^7=128,\, 2^{11}=2048,\, 2^{13}=8192,\, 2^{17}=131072,\, 2^{19}=524288,\, \dots</math> Template:OEIS

which are the prime powers of 2.

## Contents |

## Understanding a FRACTRAN program

A FRACTRAN program can be seen as a type of register machine where the registers are stored in prime exponents in the argument *n*.

Using Gödel numbering, a positive integer *n* can encode an arbitrary number of arbitrarily large positive integer variables.^{[note 1]} The value of each variable is encoded as the exponent of a prime number in the prime factorization of the integer. For example, the integer

- <math>60 = 2^2 \times 3^1 \times 5^1</math>

represents a register state in which one variable (which we will call v2) holds the value 2 and two other variables (v3 and v5) hold the value 1. All other variables hold the value 0.

A FRACTRAN program is an ordered list of positive fractions. Each fraction represents an instruction that tests one or more variables, represented by the prime factors of its denominator. For example:

- <math>f_1 = \frac{21}{20} = \frac{3 \times 7}{2^2 \times 5^1}</math>

tests v2 and v5. If <math>v_2 \ge 2</math> and <math>v_5 \ge 1</math>, then it subtracts 2 from v2 and 1 from v5 and adds 1 to v3 and 1 to v7. For example:

- <math>60 \cdot f_1 = 2^2 \times 3^1 \times 5^1 \cdot \frac{3 \times 7}{2^2 \times 5^1} = 3^2 \times 7^1</math>

Since the FRACTRAN program is just a list of fractions, these test-decrement-increment instructions are the only allowed instructions in the FRACTRAN language. In addition the following restrictions apply:

- Each time a variable is tested, it is also decremented.
- The same variable cannot be both decremented and incremented in a single instruction (otherwise the fraction representing that instruction would not be in its lowest terms). Therefore each FRACTRAN instruction consumes variables as it tests them.
- It is not possible for a FRACTRAN instruction to directly test if a variable is 0 (However, an indirect test can be implemented by creating a default instruction that is placed after other instructions that test a particular variable.).

## Creating simple programs

### Addition

The simplest FRACTRAN program is a single instruction such as

- <math>\left( \frac{3}{2} \right)</math>

This program can be represented as a (very simple) algorithm as follows:

FRACTRAN

InstructionCondition Action <math>\frac{3}{2}</math> v2 > 0 Subtract 1 from v2

Add 1 to v3v2 = 0 Stop

Given an initial input of the form <math>2^a 3^b</math>, this program will stop with a final output of the form <math> 3^{a + b} </math>. It therefore adds two integers together.

### Multiplication

We can create a "multiplier" by "looping" through the "adder". In order to do this we need to introduce states into our algorithm. This algorithm will take a number <math>2^a 3^b</math> and produce <math>5^{ab}</math>:

Current State Condition Action Next State A v7 > 0 Subtract 1 from v7

Add 1 to v3A v7 = 0 and

v2 > 0Subtract 1 from v2 B v7 = 0 and

v2 = 0 and

v3 > 0Subtract 1 from v3 A v7 = 0 and

v2 = 0 and

v3 = 0Stop B v3 > 0 Subtract 1 from v3

Add 1 to v5

Add 1 to v7B v3 = 0 None A

State B is a loop that adds v3 to v5 and also moves v3 to v7, and state A is an outer control loop that repeats the loop in state B v2 times. State A also restores the value of v3 from v7 after the loop in state B has completed.

We can implement states using new variables as state indicators. The state indicators for state B will be v11 and v13. Note that we require two state control indicators for one loop; a primary flag (v11) and a secondary flag (v13). Because each indicator is consumed whenever it is tested, we need an secondary indicator to say "continue in the current state"; this secondary indicator is swapped back to the primary indicator in the next instruction, and the loop continues.

Adding FRACTRAN state indicators and instructions to the multiplication algorithm table, we have:

FRACTRAN

InstructionCurrent State State

IndicatorsCondition Action Next State <math>\frac{3}{7}</math> A None v7 > 0 Subtract 1 from v7

Add 1 to v3A <math>\frac{11}{2}</math> v7 = 0 and

v2 > 0Subtract 1 from v2 B <math>\frac{1}{3}</math> v7 = 0 and

v2 = 0 and

v3 > 0Subtract 1 from v3 A v7 = 0 and

v2 = 0 and

v3 = 0Stop <math>\frac{5 \cdot 7 \cdot 13}{3 \cdot 11}, \frac{11}{13}</math> B v11, v13 v3 > 0 Subtract 1 from v3

Add 1 to v5

Add 1 to v7B <math>\frac{1}{11}</math> v3 = 0 None A

When we write out the FRACTRAN instructions, we must put the state A instructions last, because state A has no state indicators - it is the default state if no state indicators are set. So as a FRACTRAN program, the multiplier becomes:

- <math>\left( \frac{455}{33}, \frac{11}{13}, \frac{1}{11}, \frac{3}{7}, \frac{11}{2}, \frac{1}{3} \right)</math>

With input 2^{a}3^{b} this program produces output 5^{ab}. ^{[note 2]}

### Subtraction and division

In a similar way, we can create a FRACTRAN "subtracter", and repeated subtractions allow us to create a "quotient and remainder" algorithm as follows:

FRACTRAN

InstructionCurrent State State

IndicatorsCondition Action Next State <math>\frac{7 \cdot 13}{2 \cdot 3 \cdot 11}, \frac{11}{13}</math> A v11, v13 v2 > 0 and

v3 > 0Subtract 1 from v2

Subtract 1 from v3

Add 1 to v7A <math>\frac{1}{3 \cdot 11}</math> v2 = 0 and

v3 > 0Subtract 1 from v3 X <math>\frac{5 \cdot 17}{11}</math> v3 = 0 Add 1 to v5 B <math>\frac{3 \cdot 19}{7 \cdot 17}, \frac{17}{19}</math> B v17, v19 v7 > 0 Subtract 1 from v7

Add 1 to v3B <math>\frac{11}{17}</math> v7 = 0 None A <math>\frac{1}{3}</math> X v3 > 0 Subtract 1 from v3 X v3 = 0 Stop

Writing out the FRACTRAN program, we have:

- <math>\left( \frac{91}{66}, \frac{11}{13}, \frac{1}{33}, \frac{85}{11}, \frac{57}{119}, \frac{17}{19}, \frac{11}{17}, \frac{1}{3} \right)</math>

and input 2^{n}3^{d}11 produces output 5^{q}7^{r} where *n* = *qd* + *r* and 0 ≤ *r* < *d*.

## Conway's prime algorithm

Conway's prime generating algorithm above is essentially a quotient and remainder algorithm within two loops. Given input of the form <math>2^n 7^m</math> where 0 ≤ *m* < *n*, the algorithm tries to divide *n*+1 by each number from *n* down to 1, until it finds the largest number *k* that is a divisor of *n*+1. It then returns 2^{n+1} 7^{k-1} and repeats. The only times that the sequence of state numbers generated by the algorithm produces a power of 2 is when *k* is 1 (so that the exponent of 7 is 0), which only occurs if the exponent of 2 is a prime. A step-by-step explanation of Conway's algorithm can be found in Havil(2007).

## Other examples

The following FRACTRAN program:

- <math>\left( \frac{33}{20}, \frac{5}{11}, \frac{13}{10}, \frac{1}{5}, \frac{2}{3}, \frac{10}{7}, \frac{7}{2} \right)</math>

calculates the Hamming weight H(*a*) of the binary expansion of *a* i.e. the number of 1s in the binary expansion of *a*.^{[1]} Given input 2^{a}, its output is 13^{H(a)}. The program can be analysed as follows:

FRACTRAN

InstructionCurrent State State

IndicatorsCondition Action Next State <math>\frac{3 \cdot 11}{2^2 \cdot 5}, \frac{5}{11}</math> A v5, v11 v2 > 1 Subtract 2 from v2

Add 1 to v3A <math>\frac{13}{2 \cdot 5}</math> v2 = 1 Subtract 1 from v2

Add 1 to v13B <math>\frac{1}{5}</math> v2 = 0 None B <math>\frac{2}{3}</math> B None v3 > 0 Subtract 1 from v3

Add 1 to v2B <math>\frac{2 \cdot 5}{7}</math> v3 = 0 and

v7 > 0Subtract 1 from v7

Add 1 to v2A <math>\frac{7}{2}</math> v3 = 0 and

v7 = 0 and

v2 > 0Subtract 1 from v2

add 1 to v7B v2 = 0 and

v3 = 0 and

v7 = 0Stop

## See also

## Notes

- ↑ Gödel numbering cannot be directly used for negative integers, floating point numbers or text strings, although conventions could be adopted to represent these data types indirectly. Proposed extensions to FRACTRAN include FRACTRAN++ and Bag.
- ↑ A similar multiplier algorithm is described at the Esolang FRACTRAN page.

## References

- Conway, John H.; Guy, Richard K. (1996).
*The Book of Numbers*. Springer-Verlag New York, Inc.. ISBN 038797993X. - Havil, Julian (2007).
*Nonplussed!*. Princeton University Press. ISBN 0691120560.