The Applications of
Probability to Cryptography

The underlying manuscript is held by the National Archives in the UK and can be accessed at www.nationalarchives.gov.uk using reference number HW 25/37. Readers are encouraged to obtain a copy.

The original work was under Crown copyright, which has now expired; the work is now in the public domain.

Provenance

Two Second World War research papers by Alan Turing were declassified recently. The papers, The Applications of Probability to Cryptography and its shorter companion Paper on Statistics of Repetitions, are available from from the National Archives in the UK at www.nationalarchives.gov.uk.

The released papers give the full text, along with figures and tables, and provide a fascinating insight into the preparation of the manuscripts, as well as the style of writing at a time when typographical errors were corrected by hand, and mathematical expression handwritten into spaces left in the text.

Working with the papers in their original format provides some challenges, so they have been typeset for easier reading and access. We recommend that the typeset versions are read with a copy of the original manuscript at hand.

This document contains the text and figures for The Applications of Probability to Cryptography, the companion paper is also available in typeset form from arXiv at www.arxiv.org/abs/1505.04715. These notes apply to both documents.

Separately, a journal article by Zabell provides an analysis of the papers and further background information.

The text

It is not our intent to cast Alan Turing’s manuscripts into a journal style article, but more to provide clearer access to his writing and, perhaps, to answer the questions “If Turing had have had access to typesetting software, what would his papers have looked like?”. Consequently no “house-style” copy-editing has been imposed. Occasional punctuation has been added to improve readability, some obvious errors corrected, and paragraph breaks added to ease the reading of long text blocks - and occasionally to give a better text flow. Turing uses typewriter underlining, single, and double quotes to indicate emphasis or style; these have been implemented using font format changes, double quotes are used as needed.

The manuscript has many typographical errors, deletions, and substitutions, all of which are indicated by over-typing, crossed out items, and handwritten pencil or ink annotations. These corrections have been implemented in this document to give the text that we presume Turing intended. Additionally, there are some hand written notes in the manuscript, which may or may not be by Turing; these are indicated by the use of footnotes.

British English spelling is used in the manuscript and this is retained, so words such as favour, neighbourhood, cancelling, etc. will be encountered. Turing appears to favour the spellings bigramme, trigramme, tetragramme, etc., although he is not always consistent; throughout this document the favoured rendering is used.

Turing’s wording is unchanged to give the full flavour of his original style. This means that “That is to say that we suppose for instance that ….. ” will be encountered - amongst others!

Both papers end abruptly, no summary or conclusion is offered, perhaps the papers are incomplete or pages are missing. To indicate the end of the manuscript we have marked the end of each paper with a printing sign - an infinity symbol between two horizontal bars.

In the section on a letter subtractor problem, reference is made to other methods to be discussed later in the paper. This does not happen - perhaps another indicator of an incomplete paper or missing pages.

Finally, Turing uses some forward page references that appear in the manuscript as see(p ), obviously intending to return and complete the reference. This also does not happen, so these references remain unresolved.

In short, we strive to represent Turing’s text as he wrote it.

Ciphertext, cleartext, etc.

In an attempt to capture the flavour of the time, ciphertext, cleartext, keys, etc. are displayed in a fixed pitch, bold, non-serif font to represent the typewriter, teletype, and telegraph machines that would have printed the original code, viz. CONDITIONS.

Mathematics

In the manuscript all mathematics is hand written in ink and pencil in spaces left between the typed text. Sometimes adequate space was left, other time not, and the handwriting spills into margins and adjacent lines, adding to the reading challenge. We have cast all mathematics into standard in-line or display formats as appropriate. We have used the mathcal font in places to capture the flavour of Turing’s handwriting, e.g. “the probability p” appears as “the probability $𝒫$”.

Turing uses no punctuation in his mathematics, this has been added to be consistent with modern practice; he also uses letters to reference equations - numbers are used in this document. In many places we have added parentheses to give clarity to an expression, and in some places where Turing is inconsistent in his uses of parentheses for a mathematical phrase (the expression for letter probability in the Vigenère in particular) we have chosen one format and been consistent in its use.

As Turing demonstrates a love of dense mathematics the algebraic multiplication symbol $\times$ has occasionally been used for readability, so all standard forms of multiplication will be encountered, viz., $ab,a\times b,a\cdot b$. Finally, convention suggests that the subject of a formula or expression sits on its own on the left hand side of the equals sign, with the subsidiary variables collected on the right hand side. Turing adheres to this convention as it suits him, his preference is retained.

In short, we strive to retain the elegance of Turing’s mathematics, whilst casting it into a modern format.

Figures and tables

All figures have been included with rearrangement of some items to improve clarity or document flow. Turing uses a variety of papers, styles, inks, pen, and pencil; these have all been represented in standard figure and table format.

Contents page

Turing provides a rudimentary Contents for The Applications of Probability to Cryptography, this has been reworked with some additions to make it more meaningful. Paper on Statistics of Repetitions, being much shorter, requires no Contents.

Editors

The editor can be contacted at: ian.taylor@maths.oxon.org.

Chapter 1 Introduction

1.1. Preamble

The theory of probability may be used in cryptography with most effect when the type of cipher used is already fully understood, and it only remains to find the actual keys. It is of rather less value when one is trying to diagnose the type of cipher, but if definite rival theories about the type of cipher are suggested it may be used to decide between them.

1.2. Meaning of probability and odds

I shall not attempt to give a systematic account of the theory of probability, but it may be worth while to define shortly probability and odds. The probability of an event on certain evidence is the proportion of cases in which that event may be expected to happen given that evidence. For instance if it is known the 20% of men live to the age of 70, then knowing of Hitler only Hitler is a man we can say that the probability of Hitler living to the age of 70 is 0.2. Suppose that we know that Hitler is now of age 52 the probability will be quite different, say 0.5, because 50% of men of 52 live to 70.

The odds of an event happening is the ratio $𝒫/(1-𝒫)$ where $𝒫$ is the probability of it happening. This terminology is connected with the common phraseology odds of 5:2 on meaning in our terminology that the odds are 5/2.

1.3. Probabilities based on part of the evidence

When the whole evidence about some event is taken into account it may be extremely difficult to estimate the probability of the event, even very approximately, and it may be better to form an estimate based on a part of the evidence, so that the probability may be more easily calculated. This happens in cryptography in a very obvious way. The whole evidence when we are trying to solve a cipher is the complete traffic, and the events in question are the different possible keys, and functions of the keys. Unless the traffic is very small indeed the theoretical answer to the problem “ What are the probabilities of the various keys? ” will be of the form “ The key $\mathrm{\dots }$ has a probability differing almost imperceptibly from 1 (certainty) and the other keys are virtually impossible”. But a direct attempt to determine these probabilities would obviously not be a practical method.

1.4. A priori probabilities

The evidence concerning the possibility of an event occurring usually divides into a part about which statistics are available, or some mathematical method can be applied, and a less definite part about which one can only use one’s judgement. Suppose for example that a new kind of traffic has turned up and that only three messages are available. Each message has the letter V in the 17th place and G in the 18th place. We want to know the probability that it is a general rule that we should find V and G in these places. We first have to decide how probable it is that a cipher would have such a rule, and as regards this one can probably only guess, and my guess would be about $1/5,000,000$. This judgement is not entirely a guess; some rather insecure mathematical reasoning has gone into it, something like this:-

The chance of there being a rule that two consecutive letters somewhere after the 10th should have certain fixed values seems to be about $1/500$ (this is a complete guess). The chance of the letters being the 17th and 18th is about $1/15$ (another guess, but not quite as much in the air). The probability of a letter being V or G is $1/676$ (hardly a guess at all, but expressing a judgement that there is no special virtue in the bigramme VG). Hence the chance is $1/(500\times 15\times 676)$ or about $1/5,000,000$. This is however all so vague, that it is more usual to make the judgment “$1/5,000,000$” without explanation.

The question as to what is the chance of having a rule of this kind might of course be resolved by statistics of some kind, but there is no point in having this very accurate, and of course the experience of the cryptographer itself forms a kind of statistics.

The remainder of the problem is then solved quite mathematically. Let us consider a large number of ciphers chosen at random. $N$ of them say. Of these $N/5,000,000$ of them will have the rule in question, and the remainder not. Now if we had three messages of each of the ciphers before us, we should find that for each of the ciphers with the rule, three messages have VG in the required place, but of the remaining $(4,999,999\times N)/5,000,000$ only a proportion $1/{676}^3$ will have them. Rejecting the ciphers which have not the required characteristics we are left with $N/5,000,000$ cases where the rule holds, and $(4,999,999\times N)/(5,000,000\times {676}^3)$ cases where it does not. This selection of ciphers is a random selection of ones which have all the known characteristics of the one in question, and therefore the odds in favour of the rule holding are:

	$$\frac{N}{5,000,000}:\frac{4,999,999\times N}{5,000,000\times {676}^3},$$
	$$i.e\mathit{\hspace{1em}}{676}^3:4,999,9999,$$
	$$\text{or about}\mathrm{\hspace{0.25em}60}:1\text{on}.$$

It should be noticed that the whole argument is to some extent fallacious, as it is assumed that there are only two possibilities, viz. that either VG must always occur in that position, or else that the letters in the 17th and 18th positions are wholly random. There are however many other possibilities worth consideration, e.g.

1. On the day in question we have VG in the position in question.

2. Or on another day we have some other fixed pair of letters.

3. Or in the positions 17, 18 we have to have one of the four combinations VG, RH, OM, IL and by chance VG has been chosen for all the three messages we have had.

4. Or the cipher is a simple substitution and VG is the substitute of some common bigramme, say TH.

The possibilities are of course endless, and it is therefore always necessary to bear in mind the possibility of there being other theories not yet suggested.

The a priori probability sometimes has to be estimated as above by some sort of guesswork, but often the situation is more satisfactory. Suppose for example that we know that a certain cipher is a simple substitution, the keys having no specially noticeable properties. Suppose also that we have 50 letters of such a message including five occurrences of P. We want to know how probable it it that P is the substitute of E. As before we have to answer two questions.

1. How likely is it that P would be the substitute of E neglecting the evidence of the five Es occurring in the message?

2. How likely are we to get 5 Ps?

   1. If P is not the substitute of E

   2. If P is the substitute of E.

I will not attempt to answer the second question for the present. The answer to the first is simply that the probability of a letter being the substitute of E is independent of what the letter is, and is therefore always $1/26$, in particular it is $1/26$ for the letter P. The only guesswork here is the judgement that the keys are chosen at random.

1.5. The Factor Principle

Nearly all applications of probability to cryptography depend on the factor principle (or Bayes’ Theorem). This principle may first be illustrated by a simple example. Suppose that one man in five dies of heart failure, and that of the men who die of heart failure two in three die in their beds, but of the men who die from other causes only one in four dies in their beds. (My facts are no doubt hopelessly inaccurate). Now suppose we know that a certain man died in his bed. What is the probability that he died of heart failure? Of all numbering N say we find that

	$N\times$	$(1/5)\times (2/3)$		die in their beds of heart failure
	$N\times$	$(1/5\times (1/3)$		$\mathrm{\dots }\mathit{\text{ elsewhere }}\mathit{\hspace{1em}}\mathrm{\dots }\mathrm{\dots }\mathrm{\dots }\mathrm{\dots }\mathrm{\dots }\mathrm{\dots }$
	$N\times$	$(4/5)\times (1/4)$		die in their beds from other causes
	$N\times$	$(4/5\times (3/4)$		$\mathrm{\dots }\mathit{\text{ elsewhere }}\mathit{\hspace{1em}}\mathrm{\dots }\mathrm{\dots }\mathrm{\dots }\mathrm{\dots }\mathrm{\dots }\mathrm{\dots }$

Now as our man died in his bed we do not need to consider the cases of men who did not die in their beds, and these consist of

	$N\times$	$(1/5)\times (2/3)$		cases of heart failure and
	$N\times$	$(4/5)\times (1/4)$		from other causes

and therefore the odds are $1\times (2/3):4\times (1/4)$in favour of heart failure. If this had been done algebraically the result would have been

$$\begin{array}{c}\mathit{\text{A posteriori odds of the theory}}\hfill \\ \hfill =\mathit{\text{A priori odds of the theory}}\\ \hfill \times \frac{\mathit{\text{Probability of the data being fulfilled if the theory is true}}}{\mathit{\text{Probability of the data being fulfilled if the theory is false}}}.\end{array}$$

In this the theory is that the man died of heart failure, and the data is that he died in his bed.

The general formula above will be described as the factor principle, the ratio

$$\frac{\mathit{\text{Probability of the data if the theory is true}}}{\mathit{\text{Probability of the data if the theory is false}}},$$

is called the factor for the theory on account of the data.

1.6. Decibanage

Usually when we are estimating the probability of a theory there will be several independent pieces of evidence e.g. following our last example, where we want to know whether a certain man died of heart failure or not, we may know

1. He died in his bed

2. His father died of heart failure

3. His bedroom was on the ground floor

and also have statistics telling us

1. 2/3 of men who die of heart failure die in their beds

2. 2/5 ……………………………have fathers who died of heart failure

3. 1/2 ……………………………have bedroom on the ground floor

4. 1/4 of men who died from other causes die in their beds

5. 1/6 ……………………………have fathers who died of heart failure

6. 1/20 of men who die of other cause have their bedrooms on the ground floor

Let us suppose that the three pieces of evidence are independent of one another if we know that he died of heart failure, and also if we know that he did not die of heart failure. That is to say that we suppose for instance that knowing that he slept on the ground floor does not make it any more likely that he died in his bed if we knew all along that he died of heart failure. When we make these assumptions the probability of a man who died of heart failure satisfying all three conditions is obtained simply by multiplication, and is $(2/3)\times (2/5)\times (1/2)$ and likewise for those who died from other causes the probability is $(1/4)\times (1/6)\times (1/20)$, and the factor in favour of the heart theory failure is

$$\frac{(2/3)\times (2/5)\times (1/2)}{(1/4)\times (1/6)\times (1/20)}.$$

We may regard this as the product of three factors $(2/3)/(1/4)$ and $(2/5)/(1/6)$ and $(1/2)/(1/20)$ arising from from the three independent pieces of evidence. Products like this arise very frequently, and sometimes one will get products involving thousands of factors, and large groups of these factors may be equal. We naturally therefore work in terms of the logarithms of the factors. The logarithm of the factor, taken to the base ${10}^{1/10}$ is called decibanage in favour of the theory. A deciban is a unit of evidence; a piece of evidence is worth a deciban if it increase the odds of the theory in the ratio ${10}^{1/10}:1$. The deciban is used as a more convenient unit that the ban. The terminology was introduced in honor of the famous town of Banbury.

Using this terminology we might say that the fact that our man died in bed scores 4.3 decibans in favour of the heart failure theory $(10\log (8/3)=4.3)$. We score a further 3.8 decibans for his father dying of heart failure, and 10 for his having his bedroom on the ground floor, totalling 18.1 decibans. We then bring in the a priori odds 1/4 or ${10}^{-6/10}$ and the result is the the odds are ${10}^{12.1/10}$, or as we may say “12.1 deciban up on evens”. This means about 16:1 on.

Chapter 2 Straightforward Cryptographic Problems

2.1. Vigenère

The factor principle can be applied to the solutions of a Vigenère problem with great effect. I will assume here that the period of the cipher has already been determined. Probability theory may be applied to this part of the problem also, but that is not so elementary. Suppose our cipher, written out in its correct period is

D K Q H S H Z N M P

R C V X U H T E A Q

X H P U E P P S B K

T W U J A G D Y O J

T H W C Y D Z H G A

P Z K O X O E Y A E

B O K B U B P I K R

W W A C E J P H L P

T U Z Y F H L R Y C

Figure 1. Vigenère problem.

(It is only by chance that it makes a rectangular array.)

Let us try to find the key for the first column, and for the moment let us only take into account the evidence afforded by the first letter D. Let us first consider the key B. The factor principle tells us

$$\begin{array}{c}\mathit{\text{Odds in favour of key }}\text{𝙱}=\mathit{\text{ A priori odds in favour of key }}\text{𝙱}\hfill \\ \hfill \times \frac{\mathit{\text{Probability of getting }}\text{𝙳}\mathit{\text{ in cipher if key is }}\text{𝙱}}{\mathit{\text{Probability of getting key }}\text{𝙳}\mathit{\text{ in cipher if key is not }}\text{𝙱}}\end{array}$$

Now the a priori odds in favour of key B may be taken as 1/25. The probability of getting D in the cipher with the key B is just the probability of getting C in the clear which (using the count on 1000 letters in Fig 2) is 0.021. If however the key is not B we can have any letter other the C in the clear, and the probability is (1 - 0.021)/25. Using the evidence of the D then the odds in favour of the key B are

$$\frac{1}{25}\times \left(\frac{25\times 0.021}{1-0.021}\right).$$

We may then consider the effect of the next letter in the column R which gives a further factor of (25 x 0.064)/(1 - 0.064). We are here assuming that the evidence of the R is independent of the evidence of the D. This is not quite correct, but is a useful approximation; a more accurate method of calculation will be given later. Let us write ${𝒫}_{\alpha }$ for the frequency of the letter $\alpha$ in plain language. Then our final estimate for the odds in favour of key B is

$$\frac{1}{25}\prod _i\frac{25{𝒫}_{{\alpha }_i-1}}{1-{𝒫}_{{\alpha }_i-1}}.$$

where ${\alpha }_1,{\alpha }_2,\mathrm{\dots }$ is the series of letters in the 1st column, and we use the letters and numbers interchangeably, A meaning 1, B meaning 2, $\mathrm{\dots }$, Z meaning 26 or 0. More generally for key $\beta$ the odds are

$$\frac{1}{25}\prod _i\frac{25{𝒫}_{{\alpha }_i-\beta +1}}{1-{𝒫}_{{\alpha }_i-\beta +1}}.$$

The value of this can be calculated by having a table of the decibanage corresponding to the factors $25{𝒫}_{\alpha }/(1-{𝒫}_{\alpha })$. One then decodes the column with the various possible keys, looks up the decibanage, and adds them up.

The most convenient form for doing this is a table of values of $20{\log }_{10}[25{𝒫}_{\alpha }/(1-{𝒫}_{\alpha })]$, taken to the nearest integer, or as we may say, the values of the score in half decibans. One may also have columns showing multiples of these, and the table made of double height (Fig 2.3). For the first column with key B the decoded column is CQWS••OAV, and we score -5 for C, -26 for Q, -5 for W, 17 for the three letters S, 5 for O, 7 for A and -10 V, totalling -17. These calculations can be done very quickly by the use of the transparent gadget Fig 2.4 , in which squares are ringed in pencil to show the number of letters occurring in the column.

A	84	J	2	S	73
B	23	K	5	T	81
C	21	L	38	U	19
D	46	M	34	V	11
E	116	N	66	W	21
F	20	O	66	X	16
G	25	P	15	Y	24
H	49	Q	2	Z	3
I	76	R	64

Figure 2. Count on 1000 letters.

(English text)

The value for X has been taken more of less at random as a compromise

between real language & telegraphese. Also I added to each entry (see p ).

Figure 2.3. Table for scoring a Vigenère.

In units of half a deciban.

Editor -This table is not referenced in the manuscript, does not have a page number, and nothing appear to indicate its purpose. It is presented here for completeness.

Figure 2.4. Apparatus for scoring a Vigenère.

Pencil marks arranged for 1st wheel of Fig. 1.

The gadget may be placed over Fig 2.3 in various positions corresponding to the various keys. The score is obtained by adding up the numbers showing through the various squares. In Fig 2.5 the alphabet has been written in a vertical below the cipher text of Fig 1, each letter representing a possible key. The score for each key has been written opposite the key, and under the relevant column. An X denotes a bad score, not worth adding up. Usually these will be -15 or worse. It will be seen that for the first column P, having a score of 43 is extremely likely to be right, especially as there is no other score better than 8. If we neglect this latter fact the odds for the key are $(1/25){10}^{2.15}$ i.e. about 5:1 on. The effect of decoding this column with key P has been shown underneath.

For the second column the best key is O, but is by no means so certain as the first column. The decode for this column is also shown, and provides very satisfactory combinations with the first column, confirming both keys. (This confirmation could also be based on probability theory, given a table of bigramme frequencies). In the third column I and C are best although D would be very possible, and in the fourth column Q and U are best.

Figure 2.5. Scoring and solving a Vigenère.

Writing down the possible decodes we see that the first line must read OWING and this makes the other lines read CONDI, ITHAS, EIMPO, ETOIM, ALCUL, MACHI, HISIS, EGRET. By filling in the word CONDITIONS the whole can now be decoded.

A more accurate argument would run as follows. For the first column, instead of setting up as rival theories the two possibilities that B is the key and that B is not we can set up 26 rival theories that the key is A or B or … Z, and we may apply the factor principle in the form:-

$$\begin{array}{c}\frac{\mathit{\text{A posteriori probability of key A}}}{\mathit{\text{A priori probability of key }}\text{𝙰}\times \mathit{\text{Probability of getting the given column with key }}\text{𝙰}},\hfill \\ \hfill =\frac{\mathit{\text{A posteriori probability of key }}\text{𝙱}}{\mathit{\text{A priori probability of key }}\text{𝙱}\times \mathit{\text{Probability of getting the given column with key }}\text{𝙱}},\end{array}$$

= etc.

The argument to justify this form of factor principle is really the same as for the original form. Let $q_{\beta }$ be the a priori probability of key $\beta$. Then out of N cases we have $N{q}_{\beta }$ cases of key $\beta$. Let $𝒫(\beta ,C)$ be the probability of getting the column C with key $\beta$, then we have rejected the cases where we get columns other than C we find that there are $N{q}_{\beta }𝒫(\beta ,C)$ cases of key $\beta$ i.e. the a posteriori probability of key $\beta$ is $KN{q}_{\beta }𝒫(\beta ,C)$, where $K$ is independent of $\beta$.

We have therefore to calculate the probability of getting the column C with key $\beta$ and this is simply ${\prod }_i{𝒫}_{\left({\alpha }_i-\beta +1\right)}$, i.e. the product of the frequencies of the decode letters which we get if the key is $\beta$.

Since the a priori probabilities of the keys are all equal we may say that the a posteriori probabilities are in the ratio ${\prod }_i{𝒫}_{{\alpha }_i-\beta +1}$ i.e. in the ratio ${\prod }_{i}26{𝒫}_{{\alpha }_i-\beta +1}$ which is more convenient for calculation. The final value for the probability is then

$$\frac{\prod _{i}26{𝒫}_{{\alpha }_i-\beta +1}}{\sum _{\beta }\prod _{i}26{𝒫}_{{\alpha }_i-\beta +1}}.$$

The calculation of the product ${\prod }_{i}26{𝒫}_{{\alpha }_i-\beta +1}$ may be done by the method recommended before for

$$\prod _i\frac{25{𝒫}_{{\alpha }_i-\beta +1}}{1-{𝒫}_{{\alpha }_i-\beta +1}}.$$

$($The table in Fig 3 was in fact made up for ${\prod }_{i}26{𝒫}_{{\alpha }_i-\beta +1}$. The differences between the two tables would of course be rather slight$)$. The new result is more accurate than the old because of the independence assumption in the original result.

If we only want to know the ratios of the probabilities of the various keys there is no need to calculate the denominator ${\sum }_{\beta }{\prod }_{i}26{𝒫}_{{\alpha }_i-\beta +1}$. This denominator has however another importance: it gives us some evidence about other assumptions, such as that the cipher is Vigenère, and that the period is 10. This aspect will be dealt with later (p. ).

2.2. A letter subtractor problem

A substitution with the period $91\times 95\times 99$ is obtained by superimposing three substitutions of periods 91, 95, and 99, each substitution being a Vigenère composed of slides of 0, 1, 2, 3, 4, 5, 6, 7, 8, or 9. The three substitutions are known in detail, but we do not know for any given message at what point in the complete substitution to begin. For many messages however we can provide a more or less probable crib. How can we test the probability of a crib before attempting to solve it? It may be assumed that approximately equal numbers of slides 0, 1, …,  9 occur in each substitution.

The principle of the calculation is that owing to the way in which the substitution is built up, not all slides are equally frequent, e.g. a slide of 25 can only be the sum of slides of 9, 8, and 8, or 9, 9, and 7 whilst a slide of 15 can be any of the following

9,6,0	8,7,0	7,7,1	6,6,3
9,5,1	8,6,1	7,6,2	6,5,4
9,4,2	8,5,2	7,5,3
9,3,3	8,4,3	7,4,4

A crib will therefore, other things being equal, be more likely if it requires a slide of 15 than if it requires a slide of 25. The problem is to make the best use of this principle, by determining the probability of the crib with reasonable accuracy, but without spending long over it.

We have to find the probability of getting a given slide. To do this we can apply several methods.

1. We can produce a long stretch of key by addition and take a count of the resulting slides. This is obviously a very general method, and requires no special mathematical technique. It may be rather laborious, but by interpreting a small count with common sense one can probably get quite good results.

2. There are 1000 possible combinations of slides all equally likely viz. 000, 001, …, 999. We can add up the digits in these and take the remainder on division by 26, and then count the number of combinations giving each of the possible remainders.

3. We can make use of a trick which might appear to be rather special, but is really applicable to a multitude of problems. Consider the expression

   $$f(x)={\left(1+x+x^2+\mathrm{\cdots }+x^9\right)}^3.$$

   For each possible way of expressing a number n as the sum of three numbers 0, …, 9, say $n=m_1+m_2+m_3$, there is a term $x^{m_1}{x}^{m_2}{x}^{m_3}$ in $f(x)$, $x^{m_1}$ coming out of the first factor, $x^{m_2}$ out of the second, and $x^{m_3}$ out of the third. Hence the number of ways of expressing n in the form $n=m_1+m_2+m_3$, is the coefficient of $x^n$ in $f(x)$ i.e. in

   $$\frac{{\left(1-x^{10}\right)}^3}{{\left(1-x\right)}^3},$$

   or in

   $$\left(1-3x^{10}+3x^{20}-x^{30}\right){\left(1-x\right)}^{-3}.$$

   Expanding ${\left(1-x\right)}^{-3}$ by the binomial theorem

   $$\begin{array}{c}{\left(1-x\right)}^{-3}=1+3x+6x^2+10x^3+15x^4+21x^5+28x^6+36x^7\hfill \\ \hfill +45x^8+55x^9+66x^{10}+78x^{11}+91x^{12}+105x^{13}\\ \hfill +120x^{14}+136x^{15}+153x^{16}+171x^{17}+190x^{18}\\ \hfill +210x^{19}+231x^{20}+253x^{21}+276x^{22}+300x^{23}\\ \hfill +325x^{24}+351x^{25}+378x^{26}+406x^{27}+435x^{28}+\mathrm{\dots }.\end{array}$$

   Now multiply by $1-3x^{10}+3x^{20}-x^{30}$ and we get

   $$\begin{array}{c}f(x)=1+3x+6x^2+10x^3+15x^4+21x^5+28x^6+36x^7\hfill \\ \hfill +45x^8+55x^9+63x^{10}+69x^{11}+73x^{12}+75x^{13}\\ \hfill +75x^{14}+73x^{15}+69x^{16}+63x^{17}+55x^{18}\\ \hfill +45x^{19}+36x^{20}+28x^{21}+21x^{22}+15x^{23}\\ \hfill +10x^{24}+6x^{25}+3x^{26}+x^{27}\end{array}$$

   This means to say that the chances of getting totals 0, 1, 2, …  are in the ratio 1, 3, 6, 10, …  The chances of getting remainders of 0, 1, 2, …  on division by 26 are in the ration 4, 4, 6, 10, 15, …  To get true probabilities these must be divided by their total which is conveniently 1000.

4. There are two other methods, both connected with the last method but not relying so much on the special features of the problem. They will be discussed later.

Suppose then that the probabilities have been calculated by one method or the other (as in fact we have done under (c)). We can then estimate the values of cribs. Let us suppose that a possible crib for a message beginning MVHWUSXOWBVMMK was AMBASSADOR so that the slides were 12, 9, 6, 22, 2, 0, 23, 11, 14. The slide of 12 gives us some slight evidence in favour of the crib being right for slides of 12 occur with frequency 0.073 with right cribs, whilst with wrong cribs they occur with frequency only 1/26. The factor in favour of the crib is therefore $26\times 0.073$ or about 1.9. A similar calculation may be made for each of the slides, but of course the work may be greatly speeded up by having the values of the factors 26 $C_s/1000$ in half decibans tabulated: here $C_s$ is the coefficient of $x^s$ in the above polynomial $f(x)$. The table is given below (Fig 6)

1	0	-20
2	25	-16
3	24	-12
4	23	-8
5	22	-6
6	21	-3
7	20	-1
8	19	1
9	18	3
10	17	4
11	16	5
12	15	6
13	14	6

Figure 6. Scores in half decibans of the various slides.

Evaluating this crib by means of this table we score

$$6+3-3-6-16-20-8+5+6(=-33),$$

i.e. the crib is worse by a factor of ${10}^{-33/20}$ than it was before e.g. if the a priori odds of the crib were 2:1 against it becomes 98:1 against. This crib was in fact made up at random i.e. the letters of the cipher text were chosen at random.

Now let us take one made up correctly, i.e. really enciphered by the method in question, but with a random chosen key.

N	Y	X	L	N	X	I	Q	H	H
A	M	B	A	S	S	A	D	O	R
13	12	22	11	21	5	8	13	19	16

(slides)

This scores 15 so that if it were originally 2:1 against, it now becomes nearly 3:1 on.

Having decided on a crib the natural way to test it is to have a catalogue of the positions in which a given series of slides is obtained if the 91 period component is omitted. We make 91 different hypotheses as to this third component, draw an inference as to what is the part of the slide arising from the components of periods 95 and 99 combined. This we look up in the catalogue. This process is fairly lengthy, and as the scoring of the crib takes only a minute it is certainly worth doing.

2.3. Theory of repeats

Suppose we have a cipher in which there are several very long series of substitutions which can be used for enciphering a message, but that one may sometimes get two messages enciphered with the same series of substitutions (or possibly, the series of substitutions for one message being those for another with some at the beginning omitted). In such a case let us say that the messages fit, or that they fit at such and such a distance, the distance being the number of substitutions which have to be omitted from the one series to obtain the other series. One will frequently want to know whether two messages fit or not, and we may find some evidence about this by examining the repeats between them.

By the repeats between them I mean this. One writes out the cipher texts of the two messages with the letters which are thought to have been enciphered with the same substitution under one another. One then writes under these messages a series of letters O and X, an O being written where the cipher texts differ and an X where they agree. The series of letters O and X will begin where the second message begins and end where the first to end ends. This series of letters O and X may be called the repetition figure. It may be completed by adding at the ends an indication of how many letters there are which do not overlap, and which message they belong to.

As an example:

GFRLIKQGVBMILAFIXMMOROGBYSKYXDAZCHMUMRKBZLDLDDOHCMVTIPRSD

VLOVDYQCEJSOPYGBMBKYXDAZNBFIOPTFCXDOD

${}^8$ XOOOOOOOOOOXOOXXOOXXXXXXOOOOOOOOOOXOX${}^{\mathrm{11}}$

On the whole one expects that a fit is more likely to be right the more letters X there are in the repetition figure, and that long series of letters X are especially desirable. This is because it would not be very unusual for two fairly common words to lie directly under one another when the clear texts are written out, thus

THEMAINCONVOYWILLARRIVE …

ALLCONVOYSMUSTREPORT …

XOOXXXXXXOOOOOXOOOO …

If the corresponding cipher texts really fit, i.e. if the letters in the same column are enciphered with the same substitution, then the condition for an X in the repetition figure of the cipher texts is that there be an X in the repetition figure of the corresponding clear text. Now series of several consecutive letters X can occur quite easily as above by two identical words coming under one another, or by such combinations as

ITISEASIERTOTEACHTHANALGEBRA …

THERAINWASSUCHTHATHECOULD …

OOOOOOOOOOOOXXXXXOOOOOOOO …

if the messages really fit, but if not they can only occur by complete coincidence. One therefore tends to believe that there is a fit when one gets such series of letters X. As regards single cases of X the value of them is not so clear, but one can see that if ${𝒫}_{\alpha }$ is the frequency of letters $\alpha$ in plain language then the frequency of letters X as a whole in comparison of plain language with plain language is ${\sum }_{\alpha }{𝒫}_{\alpha }^{\mathrm{\hspace{0.33em}2}}$, whilst for wrong fits of cipher text it is 1/26 which is necessarily less. Given a sufficiently long repetition figure one should therefore be able to tell whether it is a fit or not simply by counting the letters X and O.

So much is well known. The real point of this section is to show these ideas can be developed into an accurate method of estimating the probabilities of fits.

2.3.1. Simple form of theory

The complete theory takes account of the various possible lengths of repeat. As this theory is somewhat complicated it will be as well to give first two simplified forms of the theory. In both cases the simplification arises by neglecting a part of the evidence. In the first simplified form of theory we neglect all evidence except the number of letters X and the number of letters O. In the other simplified form the evidence is the number of series of (say) four consecutive letters X in a repetition figure.

When our evidence is just the number of times X occurs in the repetition figure, (n let us say) and the length of the repetition figure (N say), then the factor in favour of the fit is

$$\frac{\mathit{\text{Probability of a right repetition figure of length N and n occurrences of }}\text{𝚇}}{\mathit{\text{Probability of a wrong repetition figure of length N having n occurrences of }}\text{𝚇}}.$$

As an approximation we may assume that the numerator of this expression has the same value as if the right repetition figures were produced letter by letter by independent random choices, with a certain fixed probability of getting an X at each stage. This probability will have to be $\beta ={\sum }_{\alpha }{𝒫}_{\alpha }^{\mathrm{\hspace{0.33em}2}}$. The numerator is then

$$\begin{array}{c}\left(\mathit{\text{Number of repetition patterns with length N and n occurrences of }}\text{𝚇}\right)\hfill \\ \hfill \times \left(\mathit{\text{Probability of getting a given such repetition pattern by the process just mentioned}}\right),\end{array}$$

which we may write as $R(N,n)Q(N,n)$. Now let us denote by $y_i$ the $i$th symbol of the given repetition pattern and put ${\tau }_x=\beta$ and ${\tau }_0=1-\beta$. Then $Q(N,n)$, the probability of getting the repetition pattern is ${\prod }_{i=1}^N{\tau }_{y_i}$ which simplifies to ${\beta }^n{(1-\beta )}^{N-n}$. We may do a similar calculation for the denominator, but here we must take $\beta =1/26$ since all letters occur equally frequently in the cipher. The denominator is then

$$R(N,n){\left(\frac{1}{26}\right)}^n{\left(\frac{25}{26}\right)}^{N-n}.$$

In dividing to find the factor for the fit $R(N,n)$ cancels out, leaving

$${\left(26\beta \right)}^n{\left(\frac{26}{25}\left(1-\beta \right)\right)}^{N-n}.$$

In other words we score a factor of $26\beta$ for an X and a factor of $(26/25)(1-\beta )$ for an O. More convenient is to regard it as $10{\log }_{10}\left[(25\beta )/(1-\beta )\right]$ decibans for an X and $10{\log }_{10}[(26/25)/(1-\beta )]$ per unit length of repetition figure (per unit overlap).

An alternative argument, leading to the same result, runs as follows. Having decided to neglect all evidence except the overlap and the number of repeats we pretend that nothing else matters, i.e. that the form of the figure is irrelevant. In this case we can regard each letter of the repetition figure as independent evidence about the fit. If we get an X the factor for the fit is

$$\frac{\mathit{\text{Probability of getting an }}\text{𝚇}\mathit{\text{ if the fit is right}}}{\mathit{\text{Probability of getting an }}\text{𝚇}\mathit{\text{ if the fit is wrong}}},$$

i.e. $\beta /(1/26)$. Similarly the factor for an O is $(1-\beta )/(25/26)$.

In either form of argument it is unnecessary to calculate the number $R(N,n)$. In this particular case there is no particular difficulty about about it: it is the binomial coefficient. In some similar problems this cancelling out is a great boon, as we might not be able to find any simple form for the factor which cancels. The cancelling out is a normal feature of this kind of problem, and it seems quite natural that it should happen when we think of the second form of argument in which we think of the evidence as consisting of a number of independent parts.

The device of assuming, as we have done here, that the evidence which is not available is irrelevant can often be used and usually leads to good results. It is of course not supposed that the evidence really is irrelevant, but only that the error resulting from the assumption when used in this kind of way is likely to be small.

2.3.2. Second simplified form of theory

In the second simplified form of theory we take as our evidence that a particular part of the repetition figure is OXXXXO (say, or alternatively OXXXXXO say). The factor is then

$$\frac{\mathit{\text{Frequency of }}\text{𝙾𝚇𝚇𝚇𝚇𝙾}\mathit{\text{ in right repetition figures}}}{\mathit{\text{Frequency of }}\text{𝙾𝚇𝚇𝚇𝚇𝙾}\mathit{\text{ in wrong repetition figures}}}.$$

The denominator is

$${\left(\frac{1}{26}\right)}^4{\left(\frac{25}{26}\right)}^2,$$

and the numerator may be estimated by taking a sample of language hexagrams and counting the number of pairs that have the repetition figure OXXXXO. The expectation of the number of such pairs is the sum for all pairs of the probabilities of those pairs having the desired repetition figure i.e. is the number of such pairs (viz $N(N-1)/2$ where $N$ is the size of the sample) multiplied by the frequency of OXXXXO repetition figures. This frequency may therefore be obtained by division if we equate the expected number of these repetition figures to the actual number.

2.3.3. General form of theory

It is not of course possible to have statistics of every conceivable repetition figure. We must make some assumptions to reduce the variety that need to be considered. The following assumption is theoretically very convenient, and also appears to be a very good approximation.

The probability of repeats at two points known to be separated by a point where there is known to be no repeat are independent.

We may also assume that the probability of a repeat is independent of anything but the repetition figure in this neighbourhood. (We may however as a refinement produce different positions in a message). We can therefore think of repetition figures as being produced by selecting the symbols of the figure consecutively, the probability of getting an X at each stage being determined by the repetition figure from the point in question back as far as the last O. Sometimes this will take us back as far as the beginning of the message, and will include the number telling us how many more letters there are which do not repeat at all. We need in practice only distinguish two cases, where this number is 0 and when it is more. We may also neglect the question as to which message occurs first. We therefore have to distinguish the following cases

O	$a_0$	some	$b_0$	none	$c_0$
OX	$a_1$	some X	$b_1$	none X	$c_1$
OXX	$a_2$	some XX	$b_2$	none XX	$c_2$
OXXX	$a_3$	some XXX	$b_3$	none XXX	$c_3$
…		…		…

The entries $a_0,a_1,b_0,etc.$ opposite the repetition figures are the notations we are adopting for the probability of getting another X following such a figure. Strictly speaking we should also bring in a notation for the probability of the message coming to an end after any given repetition figure. As the repeats at the end of a comparison do not appear to behave very differently from those in the main part of the message I shall neglect this complication by assuming that the probability of getting an O added to the probability of getting an X is 1, and that afterwards one cuts off the end of the series arbitrarily.

Let us calculate the factor for the repeat figure

Underneath each symbol has been written the probability that one would get that symbol, knowing the ones which precede, both for the case of a right and of a wrong repetition figure. The factor for the fit is the product of the first row divided by the product of the second. It is convenient to split this up as indicated by the vertical lines into the product of

		${\displaystyle \frac{c_0{c}_1{c}_2{c}_3\left(1-c_4\right)}{{(1/26)}^4\times \left(25/26\right)}},$
		${\displaystyle \frac{1-a_0}{\left(25/26\right)}},\text{- occurring three times},$
		${\displaystyle \frac{a_0(1-a_1)}{\left(1/26\right)\times \left(25/26\right)}},$
		${\displaystyle \frac{a_0{a}_1{a}_2\left(1-a_3\right)}{{(1/26)}^3\times \left(25/26\right)}},$
		${\displaystyle \frac{a_0{a}_1}{{(1/26)}^2}},$

and this product may be put into the form of the product of

$$\frac{c_0{c}_1{c}_2{c}_3\left(1-c_4\right)}{{\left(1/26\right)}^4\times \left(25/26\right)}\times {\left(\frac{1-a_0}{\left(25/26\right)}\right)}^{-5},$$

- which we call the factor for an

initial tetragramme repeat level,

$$\frac{a_0(1-a_1)}{\left(1/26\right)\times \left(25/26\right)}\times {\left(\frac{1-a_0}{\left(25/26\right)}\right)}^{-2},$$

- the factor for a single repeat,

$$\frac{a_0{a}_1{a}_2\left(1-a_3\right)}{{\left(1/26\right)}^3\times \left(25/26\right)}\times {\left(\frac{1-a_0}{\left(25/26\right)}\right)}^{-4},$$

- the factor for a trigramme,

$$\frac{1-a_0}{1-a_2},$$

- the correction for a final bigramme,

$${\left(\frac{1-a_0}{\left(25/26\right)}\right)}^{16},$$

- the factor for an overlap of 16,

$$\frac{a_0{a}_1\left(1-a_2\right)}{{\left(1/26\right)}^2\times \left(25/26\right)}\times {\left(\frac{1-a_0}{\left(25/26\right)}\right)}^{-3},$$

- the factor for a trigramme.

We shall neglect the correction for a final bigramme (or whatever it may be). It is in any case rather small, and vanishes if the repetition figure ends with O; also with our conventions the whole question of the ends of repetition figures has been left rather in doubt.

Now let us put

	$$a_0{a}_1\mathrm{\dots }{a}_r(1-a_{r+1})=k_r,$$
	$$b_0{b}_1\mathrm{\dots }{b}_r(1-b_{r+1})=j_r,$$
	$$c_0{c}_1\mathrm{\dots }{c}_r(1-c_{r+1})=i_r.$$

The values of the $i_r$ can be obtained as follows. We take a number of plain language messages and leave out two or three words at the beginning. Then combine the messages to form one long message; this message may be made to eat its own tail i.e. it may be written round a circle. If the message were compared with itself in every possible position, except level, we should expect to get repetition figures which when divided up as shown by vertical lines after each O, containing $(N(N-1)/2)k_r(=N_r)$ parts which consist of r symbols O, or as we may say $N_r$ actual r-gramme repeats, where h is the probability of an O .

The values of $N_r$ can be calculated given the apparent number of r-gramme repeats $M_r$ for each $r$. This apparent number of r-gramme repeats is the number of series of r consecutive symbols X in the repetition figures regardless of what precedes or follows the series.

By considering the ways in which an actual repeat can give rise to the apparent repeat of various lengths we see that

$$M_r=N_r+2N_{r+1}+3N_{r+2}+\mathrm{\dots },$$

and therefore

$$M_r-M_{r+1}=N_r+N_{r+1}+N_{r+2}+\mathrm{\dots },$$

and

$$\left(M_r-M_{r+1}\right)-\left(M_{r+1}-M_{r+2}\right)=N_r.$$

The calculation of $j_r$ may perhaps best be done by comparing the beginners of a number of messages with the long circular message, and the values of $i_r$ by comparing the beginners among themselves. A similar technique of actual and apparent numbers of repeats can be used. I shall not go into this in detail. The formulae required may now be assembled.

	${\mu }_r$	$=\text{decibanage for an }\text{𝑟}\text{-gramme repeat},$
	$\gamma$	$=\text{negative decibanage for unit overlap},$
	$S_{\beta ,r}$	$=\text{number of occurrences in the statistics of the }\text{𝑟}\text{-gramme }\beta ,$
	$N$	$=\text{total number of letters in the statistics}.$

Then if

	$M_r$	$={\displaystyle \sum _{\beta }}{\displaystyle \frac{S_{\beta ,r}\left(S_{\beta ,r}-1\right)}{2}},$
	$N_r$	$=M_r-2M_{r+1}+M_{r+2},$
	$L$	$={\displaystyle \frac{N(N-1)}{2}},$
	$k_r$	$={\displaystyle \frac{N_r}{Lh}}.$

$h$ may be calculated as follows. From the identity

$$\left(1-a_0\right)+a_0\left(1-a_1\right)+a_0{a}_1\left(1-a_2\right)+\mathrm{\dots }=1,$$

we get

$$k_0+k_1+k_2+\mathrm{\dots }=1,$$

$$\therefore \mathit{\hspace{1em}\hspace{1em}}\frac{L-M_1}{Lh}=1,$$

$$\left(1-a_0\right)=k_0=\frac{N_0}{L-M_1}=\frac{L-2M_1+M_2}{L-M_1},$$

$${\mu }_r=10{\log }_{10}\left(\frac{{26}^{r+1}{k}_r}{25}\right)+\left(r+1\right)\nu ,$$

$$\nu =-10{\log }_{10}\left(\frac{26(1-a_0)}{25}\right).$$

2.4. Transposition ciphers

2.4.1. A probability problem

In making calculations about substitution ciphers we have often found it useful to treat the plain language as if it were produced by independent choices for the letters, using certain fixed frequencies with which the letters are chosen. Our method for Vigenère and one of the simplified forms of repeat theory could be based on this sort of assumption. With a transposition cipher however such an assumption would be useless or worse than useless, for it would result in the conclusion that all transpositions were equally likely. We have therefore to take a slightly less crude assumption, and the one which suggests itself is that the letters forming the plain language are chosen consecutively, the probability of getting a particular letter depending only on what the letter is and what the preceding letter was. It is easily verified the if ${𝒫}_{\alpha \beta }$ is the proportion of bigrammes $\alpha \beta$ in plain language and ${𝒫}_{\alpha }$ the frequency of the letter $\alpha$ then the probability $q_{\alpha \beta }$ of a letter $\beta$ following an $\alpha$ is ${𝒫}_{\alpha \beta }/{𝒫}_{\alpha }$. The probability of a piece of plain language of length $L$ letters saying ${\alpha }_1{\alpha }_2\mathrm{\dots }{\alpha }_L$ is then

$${𝒫}_{{\alpha }_1}\times q_{{\alpha }_1{\alpha }_2}\times q_{{\alpha }_2{\alpha }_3}\times q_{{\alpha }_3{\alpha }_4}\times \mathrm{\dots }\times q_{{\alpha }_{(L-1)}{\alpha }_L},$$

which may also be written as

$$𝒥({\alpha }_1,\mathrm{\dots },{\alpha }_L).$$

We may also calculate the probability of a given piece of plain language having certain given letters in given places, the remainder of the message being unspecified. The probability is given by

$$\sum ({\xi }_1,\mathrm{\dots },{\xi }_L\text{ consistent with data })𝒥({\xi }_1,\mathrm{\dots },{\xi }_L),$$

and if the data is that the known letters are

$$\underset{n_1\text{ dots }}{\mathrm{\cdots }}{\beta }_1\mathit{\hspace{1em}}\underset{n_2\text{ dots }}{\mathrm{\cdots }}{\beta }_2\mathit{\hspace{1em}}\mathrm{\cdots }\mathit{\hspace{1em}}\mathrm{\cdots }{\beta }_{r-1}\mathit{\hspace{1em}}\underset{n_r\text{ dots }}{\mathrm{\cdots }}{\beta }_r\mathit{\hspace{1em}}\mathrm{\cdots },$$

(2.1)

it is approximately

$$\prod _r{𝒫}_{{\beta }_r}\cdot \prod _{n_{r+1}=0}\frac{{𝒫}_{{\beta }_r{\beta }_{r+1}}}{{𝒫}_{\beta r}{𝒫}_{{\beta }_{r+1}}}.$$

(2.2)

A more or less rigorous deduction of this approximation from the assumptions above is given at the end of the section. For the present let us see how it can be applied. If we have two theories about the transposition of which the one requires the above pattern of letters, and the other brings the same letters in to positions in which no two of them are consecutive, then the factor in favour of the first as compared with the second is

$$\prod _{n_{r+1}=0}\frac{{𝒫}_{{\beta }_r{\beta }_{r+1}}}{{𝒫}_{\beta r}{𝒫}_{{\beta }_{r+1}}}.$$

We can apply this straightforwardly to the case of a simple transposition by columns. The following text is known to be a simple transposition of a certain type of German text with a key length of not more than 15.

S A T P T W S F A S T A U T E E A I E U F H W T J T D D G C

N L T S E F C U I E B O E Y Q H G T J T E E F I E O R T A R

U R N L N N N N A I E O T U S H L E S B F B R N D X G N J H

U A N W R

To solve this transposition, we may try comparing the first six letters of S A T P T W which we know form part of one column with each other series of six letters in the message, for we know that one such comparison will give entirely bigrammes occurring in the decode. We may try first

S F

A A

T S

P T

T A

W U

The factor for a transposition which brings these letters together, as compared with one which leaves them apart is

$$\frac{{𝒫}_{SF}}{{𝒫}_S{𝒫}_F}\times \frac{{𝒫}_{AA}}{{𝒫}_A{𝒫}_A}\times \mathrm{\dots }\times \frac{{𝒫}_{WU}}{{𝒫}_W{𝒫}_U}.$$

By using a table of values of

$$20{\log }_{10}\left(\frac{{𝒫}_{\alpha \beta }}{{𝒫}_{\alpha }{𝒫}_{\beta }}\right),$$

made up for the type of traffic in question, and given to the nearest integer (table of values of ${𝒫}_{\alpha \beta }/({𝒫}_{\alpha }{𝒫}_{\beta })$ expressed in half-decibans) we get the product by addition. Such a table is shown in Fig 6. The scores for this particular columns are SF -7, AA -7, TS -2, PT -10, TA -3, WU -13, totalling -36. If we consider this combination as a priori about 100:1 against (there are 95 letters in the message) it is a posteriori about 3000:1 against.

Figure 2.6. Exclusive bigramme scores in half decibans, i.e. $20{\log }_{10}\left(\frac{{𝒫}_{\alpha \beta }}{{𝒫}_{\alpha }{𝒫}_{\beta }}\right)$, for a certain kind of German traffic.

Figure 2.7. Scoring the matching of columns in a simple transposition. Correct matchings noted.

Similar scoring may be done for every possible comparison of S A T P T W with six consecutive letters of the message. The comparison may be made both with S A T P T W as earlier and as later column; one may also use the last six letters of the message H U A N W R.

The results of doing this are shown in Fig 7. The message has been written out vertically. The first columns of figures after the message gives the score for S A T P T W as earlier column, entered against the first letter of the later column, e.g. the -36 as calculated above gets entered against the F of F A S T A U. The second column after the message consists of the scores for H U A N W R as first column [and the column before the message gives scores for H U A N W R as second column]. One of these columns has been worked out in detail but in the other two crosses have been put in where the scores are very bad.

The scores which eventually turned to to be right are ringed. The fourth comparison, which did not have to be done scored very badly viz. -27. Amongst the good scores which were wrong there was one score of 37. It was not difficult to see that this one was wrong as most of the score came from W O with requires Z to precede it, and there was no Z in the message. Apart from this fact the comparison was about evens, although if we take into account the fact that there was no better score it would be better. [We have already had a case of this kind of thing in connection with Vigenère; if the various positions are a priori equally likely and the factors are $f_1,f_2,\mathrm{\dots },f_N$ then the value $f_r/\sum f_i$ for the probability of the $r$th alternative is better than $\left(f_r/N\right)/\left(1+f_r/N\right)$].

2.4.2. The Probability Formula

(Semi) Rigorous deduction of the formula (2.2) on page 2.2. (This is something of a digression).

The probability of a piece of plain language coinciding where necessary with the data (2.1) on page 2.1 is

$${𝒫}_{{\beta }_1}{𝒯}_{n_2,{\beta }_1{\beta }_2}{𝒯}_{n_3,{\beta }_2{\beta }_3}\mathrm{\dots }{𝒯}_{n_m,{\beta }_{m-1}{\beta }_m},$$

where

$${𝒯}_{n,\alpha \beta }\mathit{\hspace{1em}}\text{is}\mathit{\hspace{1em}}\sum _{{\eta }_1{\eta }_1\mathrm{\dots }{\eta }_n}{q}_{\alpha {\eta }_1}{q}_{{\eta }_1{\eta }_2}\mathrm{\dots }{q}_{{\eta }_n\beta },$$

since

$$\sum _{{\eta }_1\mathrm{\dots }{\eta }_{n_1}}{𝒫}_{{\eta }_1}{q}_{{\eta }_1{\eta }_2}\mathrm{\dots }{q}_{{\eta }_{n_1}{\beta }_1}={𝒫}_{{\beta }_1}.$$

We can put

$${𝒯}_{n,\alpha \beta }={\left({𝒬}^{n+1}\right)}_{\alpha \beta },$$

where $𝒬$ is the matrix whose $\alpha \beta$ coefficient is $q_{\alpha \beta }$. The formula (2.2) on page 2.2 would then be accurate if we could say that for $n>0,$

$${\left({𝒬}^{n+1}\right)}_{\alpha \beta }={𝒫}_{\beta }.$$

This is not true, but it is true that except for very special values for $q_{\alpha \beta }$,

$${\left({𝒬}^n\right)}_{\alpha \beta }\to {𝒫}_{\beta },\text{ as }n\to \mathrm{\infty },$$

and the convergence is rather rapid.

To prove this I shall assume that the eigenvalues of $𝒬$ are all different in modulus. In this case we can find a matrix $𝒰$ with unit determinant, such that ${𝒰}^{-1}𝒬𝒰$ is in the diagonal form

since $𝒬𝒰=𝒰ℳ$ we have

$$\sum _{\gamma }{q}_{\alpha \gamma }{u}_{\gamma \beta }=\sum _{\xi }{u}_{\alpha \xi }{m}_{\xi \beta },$$

i.e.

$$\sum _{\gamma }{q}_{\alpha \gamma }{u}_{\gamma \beta }={\mu }_{\beta }{\mu }_{\alpha \beta }.$$

That is, for each $\beta ,u_{\alpha \beta }$ provides a solution of

$$\sum _{\gamma }{q}_{\alpha \gamma }{l}_{\gamma }=\mu l_{\alpha },$$

(2.3)

with $\mu ={\mu }_{\beta }$. Conversely if we have any solution of (2.3) then $\mu ={\mu }_{\vartheta },l_{\alpha }=k{u}_{\alpha \vartheta }$ for some $k,\vartheta$ and all $\alpha$, for as $𝒰$ is non singular we can find numbers $c_{\gamma }$ such that

$$l_{\alpha }=\sum _{\gamma }{u}_{\alpha \gamma }{c}_{\gamma }\text{ for all }\alpha ,$$

and then substituting in (2.3) we get

$$\sum _{\gamma ,\delta }{q}_{\alpha \gamma }{u}_{\gamma \delta }{c}_{\delta }=\mu \sum _{\delta }{u}_{\alpha \delta }{c}_{\delta },$$

i.e.

$$\sum _{\delta }\left({\mu }_{\delta }{c}_{\delta }-\mu c_{\delta }\right)u_{\alpha \delta }=0.$$

Which, since $𝒰$ is nonsingular implies $\mu ={\mu }_{\delta }\text{ or }{c}_{\delta }=0\text{ for all }\delta$.

As the series ${\mu }_1,\mathrm{\dots },{\mu }_{26}$ are all different there is only one value $\vartheta$ of $\delta$ for which $\mu ={\mu }_{\delta }$ and so $l_{\alpha }=c_{\vartheta }{u}_{\alpha \vartheta }$ for all $\alpha$. Now putting $l_{\alpha }=1$ for all $\alpha$ we see that one member of the series ${\mu }_1,\mathrm{\dots },{\mu }_{26}$ is 1, for (2.3) is certainly satisfied.

I shall prove that the remaining eigenvalues satisfy $|\mu |\le 1$. We first prove that if $\mu \ne 1$ then $\sum p_{\alpha }{l}_{\alpha }=0$. This follows by multiplying (2.3) on each side by ${𝒫}_{\alpha }$ and summing. Since

$$q_{\alpha \beta }=\frac{{𝒫}_{\alpha \beta }}{{𝒫}_{\alpha }}\text{ and }\sum _{\alpha }{𝒫}_{\alpha \beta }={𝒫}_{\beta },$$

we get

$$\sum _{\alpha \gamma }{q}_{\alpha \gamma }{l}_{\gamma }=\sum p_{\gamma }{l}_{\gamma }=\mu \sum p_{\alpha }{l}_{\alpha },$$

which implies

$$\mu =1\text{ or }\sum p_{\alpha }{l}_{\alpha }=0.$$

Next we show that each $\mu$ for which $|\mu |>1$ is real and positive. Let $l_{\alpha }$ satisfy (2.3) with $|\mu |>1$; then the eigenvalue for ${\overline{l}}_{\alpha }$ is $\overline{\mu }$ and so

$$\sum _{\beta }{\left({𝒬}^r\right)}_{\alpha \beta }\left(1+\epsilon \left(l_{\beta }+{\overline{l}}_{\beta }\right)\right)=1+2\epsilon \mathrm{\Re }{\mu }^r{l}_{\alpha }.$$

If $\epsilon >0$ has been chosen so small that $\mathrm{\Re }\epsilon l_{\beta }>-1/2$ for all $\beta$ then the L.H.S. is positive for the coefficients in the matrix are positive, whereas the R.H.S. is negative for suitably chosen $\mu$, unless $l_{\alpha }=0$. If now $\mu >1$ we may take it that $l_{\alpha }$ is real for each $\alpha$. As it must satisfy $\sum p_{\alpha }{l}_{\alpha }=0$ it is negative for some $\alpha$, but then

$$\sum _{\beta }{\left({𝒬}^r\right)}_{\alpha \beta }\left(1+\epsilon l_{\beta }\right)=1+\epsilon {\mu }^r{l}_{\alpha },$$

and if $\epsilon$ is chosen so that $1+\epsilon l_{\beta }>0$ for all $\beta$ the L.H.S is positive whereas the R.H.S is negative for sufficiently large $r$.

All the eigenvalues therefore satisfy $|\mu |\le 1$ as the eigenvalues are all different in modulus this means that except for one value of $\mu$. Then as $r\to \mathrm{\infty },{ℳ}^r$ tends to a matrix which has only one element different from 0, and that a 1 on the diagonal, say in position $\sigma \sigma$.

Calling this matrix ${𝒳}_{\sigma }$ the series of matrices ${𝒬}^r$ tends to the matrix ${𝒰}^{-1}{𝒳}_{\sigma }𝒰$. This matrix is the one and only one $𝒴$ which satisfies $𝒴𝒬=𝒴,{𝒴}^{\text{2}}=𝒴,𝒴\ne \text{0}$ and is therefore the one whose $\alpha \beta$ coefficient is ${𝒫}_{\beta }$.

2.4.3. Another probability problem

There is another probability problem that arises in connection with simple transpositions. With a message of length $L$, and a key length of $K$ what is the probability that the $m$th letter will be at the bottom of a column? Let $D$ be the length of the short columns i.e. $D=\left[L/K\right]$, and let $E=L-DK$. Then if the $m$th letter is at the bottom of the $w$th column we must have

$$\frac{m}{D+1}\le w\le \frac{m}{D},$$

and there will be $(D+1)w-m$ short and $m-Dw$ long columns among these first $w$ columns. There are

$$\left(\begin{array}{c}n\\ k\end{array}\right)=C(n,k)=\frac{P(n,k)}{P(k,k)}=\frac{n!}{(n-k)!k!}$$

$$\left(\begin{array}{c}w\\ m-Dw\end{array}\right)\left(\begin{array}{c}K-w\\ E-m+Dw\end{array}\right)$$

ways in which the short and long columns can be arranged consistently with this, and altogether $\left(\begin{array}{c}K\\ E\end{array}\right)$ ways in which the columns can be arranged, so that the probability of the $m$ the letter being at the bottom of a column is

$$\sum _{\left(m/D+1\right)\le w\le \left(m/D\right)}\left(\begin{array}{c}w\\ m-Dw\end{array}\right)\left(\begin{array}{c}K-w\\ E-m+Dw\end{array}\right)/\left(\begin{array}{c}K\\ E\end{array}\right).$$

There will normally be very few terms in the sum.

Let us take the case of the message of length 133 and consider the 45th letter, assuming the key length is between 10 and 20 (inclusive). $L_O{}^B.\mathit{\hspace{1em}}L=133,m=45.$

		$K=10,D=13,E=3,{\displaystyle \frac{m}{D+1}}=3+,{\displaystyle \frac{m}{D}}=3+$		no terms
		$K=11,D=12,E=1,{\displaystyle \frac{m}{D+1}}=3+,{\displaystyle \frac{m}{D}}=3+$		no terms
		$K=12,D=11,E=1,{\displaystyle \frac{m}{D+1}}=3+,{\displaystyle \frac{m}{D}}=4+$
		$\mathrm{\hspace{1em}}\mathit{\text{only terms }}w=4\mathit{\text{, }}m-Dw=1\mathit{\text{ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}probability is:}}$		$\left({\displaystyle \begin{array}{c}4\\ 1\end{array}}\right)\left({\displaystyle \begin{array}{c}8\\ 0\end{array}}\right)/\left({\displaystyle \begin{array}{c}12\\ 1\end{array}}\right)={\displaystyle \frac{4}{12}}$
		$K=13,D=10,E=3,{\displaystyle \frac{m}{D+1}}=4+,{\displaystyle \frac{m}{D}}=4+$		no terms
		$K=14,D=9,E=7,{\displaystyle \frac{m}{D+1}}=4+,{\displaystyle \frac{m}{D}}=5$
		$\mathrm{\hspace{1em}}\mathit{\text{only terms }}w=5\mathit{\text{, }}m-Dw=0\mathit{\text{ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}probability is:}}$		$\left({\displaystyle \begin{array}{c}5\\ 0\end{array}}\right)\left({\displaystyle \begin{array}{c}9\\ 7\end{array}}\right)/\left({\displaystyle \begin{array}{c}14\\ 7\end{array}}\right)={\displaystyle \frac{3}{286}}=0.0105,$
		$K=15,D=8,E=13,{\displaystyle \frac{m}{D+1}}=5,{\displaystyle \frac{m}{D}}=5+$
		$\mathrm{\hspace{1em}}\mathit{\text{only terms }}w=5\mathit{\text{, }}m-Dw=5\mathit{\text{ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}probability is:}}$		$\left({\displaystyle \begin{array}{c}5\\ 5\end{array}}\right)\left({\displaystyle \begin{array}{c}10\\ 8\end{array}}\right)/\left({\displaystyle \begin{array}{c}15\\ 13\end{array}}\right)={\displaystyle \frac{3}{7}}=0.428,$
		$K=16,D=8,E=5,{\displaystyle \frac{m}{D+1}}=5+,{\displaystyle \frac{m}{D}}=5+$
		$\mathrm{\hspace{1em}}\mathit{\text{only terms }}w=5\mathit{\text{, }}m-Dw=5\mathit{\text{ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}probability is:}}$		$\left({\displaystyle \begin{array}{c}5\\ 5\end{array}}\right)\left({\displaystyle \begin{array}{c}11\\ 0\end{array}}\right)/\left({\displaystyle \begin{array}{c}16\\ 5\end{array}}\right)={\displaystyle \frac{1}{4368}}=0.000229,$
		$K=17,D=7,E=14,{\displaystyle \frac{m}{D+1}}=5+,{\displaystyle \frac{m}{D}}=6+$
		$\mathrm{\hspace{1em}}\mathit{\text{only terms }}w=6\mathit{\text{, }}m-Dw=3\mathit{\text{ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}probability is:}}$		$\left({\displaystyle \begin{array}{c}6\\ 3\end{array}}\right)\left({\displaystyle \begin{array}{c}11\\ 11\end{array}}\right)/\left({\displaystyle \begin{array}{c}17\\ 14\end{array}}\right)={\displaystyle \frac{1}{34}}=0.0307$
				$\mathrm{\hspace{1em}\hspace{1em}}\mathit{\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}}(Editor-1/34\text{ is }\mathrm{0.0294.})$
		$K=18,D=7,E=7,{\displaystyle \frac{m}{D+1}}=5+,{\displaystyle \frac{m}{D}}=6+$
		$\mathrm{\hspace{1em}}\mathit{\text{only terms }}w=6\mathit{\text{, }}m-Dw=3\mathit{\text{ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}probability is:}}$		$\left({\displaystyle \begin{array}{c}6\\ 3\end{array}}\right)\left({\displaystyle \begin{array}{c}12\\ 4\end{array}}\right)/\left({\displaystyle \begin{array}{c}18\\ 7\end{array}}\right)={\displaystyle \frac{4950}{15912}}=0.311,$
		$K=19,D=7,E=0,\mathit{\text{ probability is:}}$		$=0$
		$K=20,D=6,E=13,{\displaystyle \frac{m}{D+1}}=6+,{\displaystyle \frac{m}{D}}=7+$
		$\mathrm{\hspace{1em}}\mathit{\text{only terms }}w=7\mathit{\text{, }}m-Dw=3\mathit{\text{ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}probability is:}}$		$\left({\displaystyle \begin{array}{c}7\\ 3\end{array}}\right)\left({\displaystyle \begin{array}{c}13\\ 4\end{array}}\right)/\left({\displaystyle \begin{array}{c}20\\ 7\end{array}}\right)={\displaystyle \frac{35\times 143}{15504}}=\mathrm{0.323.}$

 

$\mathrm{\infty }$