Skip to content
Go back

How Is Compression Prediction?

Over the past few weeks, I have repeatedly encountered the same claim on Hacker News: compression is prediction. The recent discussion has approached it from both directions. Two 3Blue1Brown videos, Reinventing Entropy and But what is cross-entropy?, derive entropy and cross-entropy from the limits of source coding. An ngrok article follows the same mathematics through arithmetic coding and language models. Salvatore Sanfilippo asks how far the resulting identification between prediction and compression should be taken.

These explanations meet at one fact. A probabilistic model assigns a conditional probability to every possible continuation, and an entropy coder converts the probability assigned to the observed continuation into bits. For a sequence x1:nx_{1:n} and a model QQ, the resulting ideal payload length is

log2Q(x1:n)=i=1nlog2Q(xix<i)-\log_2 Q(x_{1:n}) = \sum_{i=1}^n -\log_2 Q(x_i\mid x_{<i})

up to the overhead introduced by the coding procedure. The quantity on the right is also the model’s cumulative logarithmic loss. In this setting, improving prediction under log-loss and reducing the encoded payload are the same optimization problem.

None of the underlying correspondence is new. Its foundations belong to classical information theory: Shannon connected probability to optimal code length, adaptive statistical compressors turned conditional estimates into codes long before modern language models, and the relation between learning and compression has been developed through minimum description length, MacKay’s treatment of information theory and inference, and work such as the Hutter Prize. Recent language-model results just instantiate this older correspondence at a new scale.

I have spent the last few years working on compression, information theory, and compressed data structures and wanted to give my two cents. I agree with the equivalence. What interests me is where it begins and where it ends. It describes the cost of encoding data under an agreed model, but a compression problem starts before that model can be applied and does not always end when the shortest bitstream has been produced. The encoder and decoder must agree on what kind of object is being represented, which alternatives remain possible, how the probability model is made available, and what the decoder must be able to do with the representation.

Throughout this article, compression means lossless compression unless stated otherwise. Even within that scope, compression can be defined before introducing a sequential model. A finite family of admissible objects gives a counting lower bound without identifying a next symbol. A fixed or data-dependent code can later be interpreted probabilistically, and a distribution over serialized objects can be factored into next-symbol conditionals. That reinterpretation does not choose the family of objects, pay for information unavailable to the decoder, or enforce operations such as random access.

The question is therefore not whether prediction and compression can be made mathematically equivalent. They can. The question is what must be fixed before the equivalence applies, which part of a complete representation its bit count measures, and what remains outside that measurement.

A note on level. This article is a bit technical, it assumes familiarity with undergraduate mathematics and elementary proof-style arguments, but no prior background in information theory is really required, although it may help.

Table of Contents

Open Table of Contents

Compression Before Probability

The ngrok article begins by distinguishing minification from what it calls “true” compression. A minifier removes comments, whitespace, and other parts of a source file that do not affect its execution. The resulting program is shorter, but the original source file cannot be reconstructed from it.

Whether this operation is lossless depends on what the representation is required to preserve. If the object is the original sequence of source bytes, minification is lossy. If the object is the program’s behaviour and the decoder may return any behaviourally equivalent program, a semantics-preserving minifier is lossless relative to that different contract. The transformation has not changed. The object being represented has.

This distinction precedes any probability model. Before asking how likely an object is, the encoder and decoder must agree on what counts as that object and when two decoded outputs count as equivalent. Only then does the length of a description become meaningful.

Once an individual object xx has been fixed, the most permissive effective descriptions are programs that produce it. After choosing a universal machine UU, the Kolmogorov complexity of a binary string xx is

KU(x)=min{p:U(p)=x}K_U(x) = \min \left\{ |p| : U(p)=x \right\}

Thus, KU(x)K_U(x) is the length of the shortest program that outputs xx. Any regularity that can be expressed algorithmically may shorten this description. A string containing a billion zeros has a long literal representation but a short program that prints one billion zeros. The definition does not require the string to have been sampled from a source, and it does not require one symbol to be predicted from the symbols preceding it.

The machine UU is part of the description language. Choosing a different universal machine changes which programs are available and therefore changes the exact value of the complexity. The invariance theorem bounds this dependence. For two fixed universal machines UU and VV, there is a constant cU,Vc_{U,V} such that

KU(x)KV(x)cU,V\left| K_U(x)-K_V(x) \right| \leq c_{U,V}

for every string xx. The constant may depend on the two machines, but not on xx. It accounts for the fixed program needed to simulate one description language in the other. The choice of machine is therefore part of the description language shared by whoever produces and interprets the program. This is the first instance of a recurring theme: the length of an object is meaningful only relative to information already fixed outside its description.

Kolmogorov complexity gives a limit on the effective description of an individual object, but it does not provide a general compression algorithm. The function KUK_U is not computable. No procedure can determine the length of the shortest program for every string, much less construct that program. A practical compressor must restrict the descriptions it is willing and able to consider.

One such restriction is that the object belongs to a finite family F\mathcal{F}. Once F\mathcal{F} has been fixed, a lossless representation must distinguish every member of that family from every other member. Consider a fixed-length encoding

C:F{0,1}C : \mathcal{F} \longrightarrow \{0,1\}^{\ell}

Lossless decoding requires CC to be injective. Since only 22^\ell binary strings of length \ell exist, injectivity implies

2F2^\ell \geq |\mathcal{F}|

and therefore

log2F\ell \geq \left\lceil \log_2|\mathcal{F}| \right\rceil

An agreed enumeration of F\mathcal{F} attains this bound by assigning each object an index and representing that index in binary. The quantity

log2F\log_2|\mathcal{F}|

is the counting bound of the family. In the literature on succinct data structures, which studies representations whose space approaches information-theoretic lower bounds, the quantity log2F\log_2|\mathcal{F}| is sometimes called the worst-case entropy of the family. I will use counting bound because its derivation assumes neither a uniform distribution nor any sampling process. It only counts the alternatives that the representation must distinguish.

The family F\mathcal{F} is part of the information shared by the encoder and decoder. If the decoder knows only that the object belongs to a larger family G\mathcal{G}, then the representation must distinguish among the members of G\mathcal{G} instead. The lower bound becomes

log2G\log_2|\mathcal{G}|

A restriction from G\mathcal{G} to F\mathcal{F} saves bits only if the decoder already knows that restriction or if the representation communicates it. What counts as redundancy therefore depends on which alternatives have already been excluded.

Kolmogorov complexity and the counting bound answer different versions of the same preliminary question. The first considers the shortest effective description of one object. The second considers the number of bits needed to distinguish every object in a fixed finite family. Neither requires a probability distribution or a next-symbol predictor.

Possibilities Have Different Probabilities

The counting bound treats every admissible object symmetrically. To assign shorter descriptions to some objects, we need a rule that determines which objects receive them and which objects pay with longer descriptions. A probability distribution supplies that rule.

Let X\mathcal{X} be a finite set of possible objects. For each xXx\in\mathcal{X}, a source specifies a probability

P(x)=Pr(X=x)P(x)=\Pr(X=x)

where P(x)>0P(x)>0 and

xXP(x)=1\sum_{x\in\mathcal{X}}P(x)=1

A probability must be translated into a quantity measured in bits. If two independent outcomes occur with probabilities P(x)P(x) and P(y)P(y), their joint probability is the product P(x)P(y)P(x)P(y), while their bit costs should add. The logarithm performs this conversion. The information content of an outcome xx is

IP(x)=log2P(x)I_P(x) = -\log_2P(x)

An event with probability 2b2^{-b} has information content bb bits. More probable outcomes receive smaller values because fewer bits should be allocated to events that occur more often.

Before the source produces an outcome, its information content is not known. Its expected value is

H(P)=xXP(x)IP(x)=xXP(x)log2P(x)\begin{aligned} H(P) &= \sum_{x\in\mathcal{X}}P(x)I_P(x) \\ &= -\sum_{x\in\mathcal{X}}P(x)\log_2P(x) \end{aligned}

This is the Shannon entropy of the source. When PP is uniform, every object has information content

log2X\log_2|\mathcal{X}|

so Shannon entropy equals the counting bound of the family.

The values IP(x)I_P(x) are ideal bit costs derived from the source distribution. To obtain a representation, we need a binary code

C:X{0,1}C : \mathcal{X} \longrightarrow \{0,1\}^*

with codeword lengths

C(x)=C(x)\ell_C(x) = |C(x)|

The expected number of bits used by the code is

LP(C)=xXP(x)C(x)L_P(C) = \sum_{x\in\mathcal{X}} P(x)\ell_C(x)

The source suggests the ideal real-valued lengths log2P(x)-\log_2P(x), while an actual binary code must assign integer lengths. Before comparing their averages, we need to answer two questions. Which collections of integer lengths can form a decodable code? And, once such a code has been chosen, can its lengths themselves be interpreted probabilistically? Kraft–McMillan answers the first question and makes the second possible.

Assigning a different binary string to each object is sufficient when each codeword is presented in isolation. It is not sufficient when codewords are concatenated. A concatenated bitstream could admit two decompositions into codewords and therefore two possible source sequences. A code is uniquely decodable when every concatenation of its codewords has only one decomposition.

A prefix code guarantees this property by requiring that no codeword be a prefix of another. Its codewords can be placed at leaves of a binary tree. Each edge contributes one bit, and the depth of a leaf equals the length of its codeword.

Let mm be at least as large as the longest codeword. A codeword of length C(x)\ell_C(x) has

2mC(x)2^{m-\ell_C(x)}

descendants at depth mm. Since no codeword is a prefix of another, the sets of descendants belonging to distinct codewords are disjoint. The complete binary tree contains only 2m2^m nodes at depth mm, so

xX2mC(x)2m\sum_{x\in\mathcal{X}} 2^{m-\ell_C(x)} \leq 2^m

Dividing by 2m2^m gives

xX2C(x)1\sum_{x\in\mathcal{X}} 2^{-\ell_C(x)} \leq 1

This is the Kraft–McMillan inequality. The tree argument proves it for prefix codes. The same inequality is necessary for every uniquely decodable code, even when its codewords do not form leaves of a prefix tree. Conversely, any collection of non-negative integer lengths satisfying the inequality can be realized by a prefix code.

The diagram in Reinventing Entropy visualizes the prefix-code case of this inequality. A codeword of length \ell excludes all of its descendants and therefore occupies a fraction 22^{-\ell} of the available binary coding space. The video uses this geometry to explain why assigning one symbol a shorter codeword leaves less space for the others; Kraft–McMillan makes the same constraint precise and extends it to every uniquely decodable code.

The source probabilities satisfy a related identity. From the definition of information content,

2IP(x)=P(x)2^{-I_P(x)} = P(x)

and therefore

xX2IP(x)=xXP(x)=1\sum_{x\in\mathcal{X}} 2^{-I_P(x)} = \sum_{x\in\mathcal{X}} P(x) = 1

The ideal lengths IP(x)I_P(x) satisfy the same constraint as codeword lengths, with equality. They may nevertheless be fractional, so they do not necessarily specify a binary code.

So far, probability has produced ideal code lengths. The more important direction for our argument is the converse: a code chosen without any probabilistic model still determines probability weights through its lengths. Given a uniquely decodable code CC, define its Kraft sum

SC=xX2C(x)S_C = \sum_{x\in\mathcal{X}} 2^{-\ell_C(x)}

The Kraft–McMillan inequality gives SC1S_C\leq1. The quantities 2C(x)2^{-\ell_C(x)} may therefore sum to less than one, but normalization turns them into a probability distribution:

QC(x)=2C(x)SCQ_C(x) = \frac{2^{-\ell_C(x)}}{S_C}

Solving for the codeword length gives

C(x)=log2QC(x)log2SC\ell_C(x) = -\log_2Q_C(x)-\log_2S_C

The code lengths are therefore information contents under the induced distribution QCQ_C, shifted by the same non-negative amount log2SC-\log_2S_C. When the Kraft inequality is tight, so that SC=1S_C=1, the correspondence is exact:

C(x)=log2QC(x)\ell_C(x) = -\log_2Q_C(x)

A fixed-length code is the simplest instance. If every object receives the same length, then every value 2C(x)2^{-\ell_C(x)} is equal, and normalization produces the uniform distribution over the objects. A non-uniform code induces a non-uniform distribution in which shorter codewords correspond to more probable objects.

This does not mean that a probability distribution was required to define the code. The counting argument in the previous section produced a fixed-length code without assuming a source. It means that once a uniquely decodable code has been chosen, its lengths can always be given a probabilistic interpretation. Codes constrain probability assignments, and probability assignments suggest code lengths.

The induced distribution also proves that entropy lower-bounds expected code length. Averaging

C(x)=log2QC(x)log2SC\ell_C(x) = -\log_2Q_C(x)-\log_2S_C

under the source distribution gives

LP(C)=xXP(x)log2QC(x)log2SCL_P(C) = -\sum_{x\in\mathcal{X}} P(x)\log_2Q_C(x) -\log_2S_C

Subtracting the source entropy yields

LP(C)H(P)=xXP(x)log2P(x)QC(x)log2SC=DKL(PQC)log2SC\begin{aligned} L_P(C)-H(P) &= \sum_{x\in\mathcal{X}} P(x) \log_2 \frac{P(x)}{Q_C(x)} -\log_2S_C \\ &= D_{\mathrm{KL}}\left(P\mathbin\Vert Q_C\right) -\log_2S_C \end{aligned}

Both terms are non-negative. The first is the Kullback–Leibler divergence from the source distribution to the distribution induced by the code. The second is non-negative because SC1S_C\leq1. Therefore every uniquely decodable code satisfies

LP(C)H(P)L_P(C) \geq H(P)

The excess length has two sources. The divergence measures how poorly the code lengths match the source probabilities. The term log2SC-\log_2S_C measures unused coding capacity when the Kraft inequality is not tight.

The lower bound can be approached by rounding each ideal length upward:

(x)=log2P(x).\ell(x)=\left\lceil-\log_2P(x)\right\rceil.

Since 2(x)P(x)2^{-\ell(x)}\leq P(x), these lengths satisfy the Kraft–McMillan inequality and can therefore be realized by a prefix code. Moreover,

log2P(x)(x)<log2P(x)+1,-\log_2P(x)\leq\ell(x)<-\log_2P(x)+1,

so averaging under PP gives

H(P)LP(C)<H(P)+1.H(P)\leq L_P(C)<H(P)+1.

This is the one-symbol form of Shannon’s Source Coding Theorem: entropy lower-bounds every uniquely decodable code, while a prefix code can always remain within one bit of it.

For a block of rr independent outcomes drawn from PP, the same construction gives an expected block length smaller than rH(P)+1rH(P)+1. The excess per source symbol is therefore less than 1/r1/r and tends to zero as the block length grows. For dependent outcomes, the same statement applies to their joint block distribution, with rH(P)rH(P) replaced by the corresponding block entropy.

Up to this point, each value xXx\in\mathcal{X} has been treated as a complete object. The equivalence between code lengths and probability assignments is already present, but no next symbol has appeared. Prediction enters only after an object is represented as an ordered sequence and its probability is factored into conditional probabilities for the successive symbols.

When Compression Becomes Prediction

Now suppose that the object is an ordered sequence

x1:n=x1,,xnΣnx_{1:n}=x_1,\ldots,x_n\in\Sigma^n

over an alphabet Σ\Sigma. Write

x<i=x1,,xi1x_{<i}=x_1,\ldots,x_{i-1}

for the prefix before position ii. The probability of the complete sequence satisfies the chain rule

P(x1:n)=i=1nP(xix<i)P(x_{1:n}) = \prod_{i=1}^n P(x_i\mid x_{<i})

No independence assumption is involved. Each conditional distribution may depend on the complete prefix. Here, “next” means next in the agreed serialization. The complete sequence may already be available to the encoder, prediction refers to the conditional probabilities assigned along that serialization.

Applying log2-\log_2 turns the product into a sum:

log2P(x1:n)=i=1nlog2P(xix<i)-\log_2P(x_{1:n}) = \sum_{i=1}^n -\log_2P(x_i\mid x_{<i})

The information content of the complete sequence is therefore the sum of the conditional information contents of its symbols. Averaging over all possible sequences gives

H(X1:n)=i=1nH(XiX<i)H(X_{1:n}) = \sum_{i=1}^n H(X_i\mid X_{<i})

where

H(XiX<i)=E[log2P(XiX<i)]H(X_i\mid X_{<i}) = \mathbb{E}\left[-\log_2P(X_i\mid X_{<i})\right]

measures the information that remains at position ii after the prefix is known.

A compressor rarely knows the true conditional distributions. It instead uses a model QQ that returns, for every prefix, a distribution

Q(x<i)Q(\cdot\mid x_{<i})

over Σ\Sigma, with

aΣQ(ax<i)=1\sum_{a\in\Sigma}Q(a\mid x_{<i})=1

Returning only the most likely symbol is not sufficient. The encoder already knows which symbol occurs and needs a code length for that symbol, whether or not it was the model’s first choice. Two models may select the same most likely continuation while assigning different probabilities to the observed symbol.

The model defines a probability for the complete sequence:

Q(x1:n)=i=1nQ(xix<i)Q(x_{1:n}) = \prod_{i=1}^n Q(x_i\mid x_{<i})

Its cumulative logarithmic loss is

LQ(x1:n)=i=1nlog2Q(xix<i)=log2Q(x1:n)\begin{aligned} \mathcal{L}_Q(x_{1:n}) &= \sum_{i=1}^n -\log_2Q(x_i\mid x_{<i}) \\ &= -\log_2Q(x_{1:n}) \end{aligned}

A high probability assigned to the observed symbol produces a small loss. A low probability produces a large loss. A model used for lossless coding must assign positive probability to every symbol that may occur, since probability zero would give an infinite code length.

During generation, a model chooses or samples a symbol from this distribution. During compression, the actual symbol is already known. The distribution is used to determine how much of the code space that symbol receives. Prediction in this equivalence means assigning probabilities, not guessing one continuation and replacing the data with that guess.

The model still does not produce a bitstream. An entropy coder must convert its probability assignments into a decodable representation. Arithmetic coding begins with the interval [0,1)[0,1). At position ii, it partitions the current interval into adjacent subintervals whose widths are proportional to

Q(x<i)Q(\cdot\mid x_{<i})

and retains the subinterval assigned to the observed symbol xix_i.

If the current interval has width wi1w_{i-1}, the selected interval has width

wi=wi1Q(xix<i)w_i = w_{i-1}Q(x_i\mid x_{<i})

Starting from w0=1w_0=1, the final width is

wn=i=1nQ(xix<i)=Q(x1:n)w_n = \prod_{i=1}^n Q(x_i\mid x_{<i}) = Q(x_{1:n})

The emitted binary prefix identifies a dyadic interval. To decode the sequence unambiguously, that dyadic interval must be contained in the final arithmetic-coding interval, together with an agreed termination convention. The number of required bits is therefore

log2Q(x1:n)-\log_2Q(x_{1:n})

plus a bounded coding overhead in the ideal arithmetic-coding model.

Decoding repeats the same subdivisions. After recovering x<ix_{<i}, the decoder evaluates the same distribution Q(x<i)Q(\cdot\mid x_{<i}), partitions its interval in the same order, and determines which subinterval contains the encoded value. Encoder and decoder must begin from the same state, perform the same updates, use the same symbol ordering, and agree on where the sequence ends.

This proves one direction of the equivalence: a sequential probabilistic model can be converted into a lossless compressor whose ideal payload length is the model’s cumulative log-loss.

The converse also holds. A uniquely decodable compressor assigns lengths to complete strings, and those lengths induce probability weights. After normalization, the resulting distribution can be factored into next-symbol conditionals. Language Modeling Is Compression also constructs predictors directly from changes in compressed length when candidate symbols are appended to a prefix.

The equivalence is therefore not limited to compressors explicitly implemented as a probability model followed by arithmetic coding. At a mathematical level, codes and probability assignments can be translated into one another. This generality also limits what the equivalence tells us. It applies after the possible objects, their serialization, and the information available to the decoder have been fixed. It does not determine any of them.

The log-loss measures the data encoded under QQ. It does not yet account for how QQ was chosen or how the decoder obtains it.

The Source Is Unknown

The identity above assumes that the model QQ is already available. In practice, however, the source distribution is usually unknown and must be estimated from the same sequence we want to encode. This raises a question: does empirical entropy describe a genuine property of the sequence, or only the probabilistic model we have fitted to it?

We will reach the answer in two independent ways. First, we will find the zero-order distribution that assigns the observed sequence the smallest log-loss. Then we will discard the generative interpretation, fix only the observed symbol counts, and count how many sequences remain possible. The two description lengths will differ by only a logarithmic number of bits.

Fitting a Zero-Order Model

Let

S=s1,,snΣn,n1S=s_1,\ldots,s_n\in\Sigma^n,\qquad n\geq 1

be the observed sequence. Consider first the family of zero-order models, which use the same distribution at every position and ignore the preceding symbols. Let q(a)q(a) be the probability assigned to a symbol aΣa\in\Sigma. The model assigns the complete sequence the probability

q(S)=i=1nq(si)q(S) = \prod_{i=1}^n q(s_i)

If nan_a denotes the number of occurrences of aa in SS, equal factors can be collected:

q(S)=aΣq(a)naq(S) = \prod_{a\in\Sigma} q(a)^{n_a}

The corresponding logarithmic loss is

log2q(S)=aΣnalog21q(a)-\log_2q(S) = \sum_{a\in\Sigma} n_a\log_2\frac{1}{q(a)}

Once the sequence has been observed, the counts determine which distribution in this model family assigns it the smallest loss. Define the empirical distribution

P^S(a)=nan\widehat{P}_S(a) = \frac{n_a}{n}

For any distribution qq that assigns positive probability to every symbol occurring in SS,

log2q(S)=aΣnalog21q(a)=naΣP^S(a)log21q(a)\begin{aligned} -\log_2q(S) &= \sum_{a\in\Sigma} n_a\log_2\frac{1}{q(a)} \\ &= n\sum_{a\in\Sigma} \widehat{P}_S(a) \log_2\frac{1}{q(a)} \end{aligned}

Under the empirical distribution, the probability assigned to the complete sequence is

P^S(S):=i=1nP^S(si)\widehat{P}_S(S) := \prod_{i=1}^n \widehat{P}_S(s_i)

Its logarithmic loss is

log2P^S(S)=naΣP^S(a)log21P^S(a)-\log_2\widehat{P}_S(S) = n\sum_{a\in\Sigma} \widehat{P}_S(a) \log_2 \frac{1}{\widehat{P}_S(a)}

Subtracting the two quantities gives

log2q(S)+log2P^S(S)=naΣP^S(a)log2P^S(a)q(a)=nDKL(P^Sq)0\begin{aligned} -\log_2q(S) + \log_2\widehat{P}_S(S) &= n\sum_{a\in\Sigma} \widehat{P}_S(a) \log_2 \frac{\widehat{P}_S(a)}{q(a)} \\ &= nD_{\mathrm{KL}}\left(\widehat{P}_S \mathbin\Vert q\right) \\ &\geq 0 \end{aligned}

Terms with P^S(a)=0\widehat{P}_S(a)=0 contribute zero. Since the divergence is non-negative, no zero-order distribution assigns the observed sequence a smaller log-loss than P^S\widehat{P}_S. Equality holds exactly when

q=P^Sq=\widehat{P}_S

The empirical frequencies are therefore the maximum-likelihood estimate within the family of zero-order models. Equivalently, they minimize the in-sample logarithmic loss over that family.

The resulting cost per symbol is the zero-order empirical entropy:

H0(S)=aΣnanlog2nna\mathcal{H}_0(S) = \sum_{a\in\Sigma} \frac{n_a}{n} \log_2\frac{n}{n_a}

with the convention that terms for which na=0n_a=0 contribute zero. Multiplying by nn gives

nH0(S)=aΣnalog2nna=log2P^S(S)n\mathcal{H}_0(S) = \sum_{a\in\Sigma} n_a\log_2\frac{n}{n_a} = -\log_2\widehat{P}_S(S)

Unlike Shannon entropy, H0(S)\mathcal{H}_0(S) is not defined from a source distribution that exists independently of the data. It is a property of the individual sequence SS, obtained by fitting a zero-order model to its observed symbol frequencies. It does not assert that the sequence was generated by independent draws from that distribution.

Counting Sequences Instead

We now ask whether essentially the same cost can be recovered without treating the empirical frequencies as a source distribution. The binary case makes the counting argument explicit. Consider a binary string BB of length nn containing exactly mm ones. If the decoder knows nn and mm, then BB belongs to the family

Bn,m={B{0,1}n:B contains exactly m ones}\mathcal{B}_{n,m} = \left\{ B\in\{0,1\}^n : B\text{ contains exactly }m\text{ ones} \right\}

A member of this family is determined by choosing which mm positions contain a one, so

Bn,m=(nm)|\mathcal{B}_{n,m}| = \binom{n}{m}

The counting bound from the first section says that identifying an arbitrary member of this family requires

log2(nm)\log_2\binom{n}{m}

bits, up to integer rounding.

For 0<m<n0<m<n, set

p=mnp=\frac{m}{n}

The endpoint cases m=0m=0 and m=nm=n contain only one binary string and have both counting bound and empirical entropy equal to zero. Under the zero-order model that assigns probability pp to a one, every member of Bn,m\mathcal{B}_{n,m} receives the same probability:

pm(1p)nm=(mn)m(nmn)nm=2nH0(B)\begin{aligned} p^m(1-p)^{n-m} &= \left(\frac{m}{n}\right)^m \left(\frac{n-m}{n}\right)^{n-m} \\ &= 2^{-n\mathcal{H}_0(B)} \end{aligned}

The total probability assigned to the family is therefore

(nm)2nH0(B)\binom{n}{m} 2^{-n\mathcal{H}_0(B)}

Since this probability cannot exceed one,

(nm)2nH0(B)1\binom{n}{m} 2^{-n\mathcal{H}_0(B)} \leq 1

and hence

log2(nm)nH0(B)\log_2\binom{n}{m} \leq n\mathcal{H}_0(B)

For the reverse bound, consider the number of ones produced by the fitted Bernoulli model. This count can take only the n+1n+1 values from 00 to nn. When p=m/np=m/n, the count mm is a mode of the resulting binomial distribution. Its probability is therefore at least the average probability of the possible counts:

(nm)pm(1p)nm1n+1\binom{n}{m} p^m(1-p)^{n-m} \geq \frac{1}{n+1}

Substituting the expression in terms of empirical entropy gives

(nm)2nH0(B)1n+1\binom{n}{m} 2^{-n\mathcal{H}_0(B)} \geq \frac{1}{n+1}

Taking logarithms yields

nH0(B)log2(n+1)log2(nm)n\mathcal{H}_0(B) - \log_2(n+1) \leq \log_2\binom{n}{m}

Together,

nH0(B)log2(n+1)log2(nm)nH0(B)n\mathcal{H}_0(B) - \log_2(n+1) \leq \log_2\binom{n}{m} \leq n\mathcal{H}_0(B)

The two routes therefore differ by at most log2(n+1)\log_2(n+1) bits. Dividing by nn, their difference per symbol is at most

log2(n+1)n,\frac{\log_2(n+1)}{n},

which tends to zero. The counting argument fixes the number of ones and distinguishes the binary strings that remain possible. The probabilistic argument fits that same frequency and evaluates the observed string under it. They begin from different assumptions but asymptotically assign the same cost per symbol.

If nn is known but mm is not, the value of mm must also be represented. There are n+1n+1 possible values, so a fixed-width representation uses

log2(n+1)\left\lceil \log_2(n+1) \right\rceil

bits. Once this cost is included, the counting and probabilistic descriptions agree within the same lower-order term.

The same relation holds over a general alphabet. Let the composition of SS be the vector

(na)aΣ(n_a)_{a\in\Sigma}

and consider its type class

T(na)={TΣn:T contains exactly na occurrences of each aΣ}\mathcal{T}(n_a) = \left\{ T\in\Sigma^n : T\text{ contains exactly }n_a\text{ occurrences of each }a\in\Sigma \right\}

where aΣna=n\sum_{a\in\Sigma}n_a=n. A sequence in this class is obtained by choosing which positions contain each symbol, so

T(na)=n!aΣna!|\mathcal{T}(n_a)| = \frac{n!}{\prod_{a\in\Sigma}n_a!}

Every sequence in this type class receives the same probability under the empirical distribution:

P^S(T)=aΣna>0(nan)na=2nH0(S)\begin{aligned} \widehat{P}_S(T) &= \prod_{\substack{a\in\Sigma \\ n_a>0}} \left(\frac{n_a}{n}\right)^{n_a} \\ &= 2^{-n\mathcal{H}_0(S)} \end{aligned}

The total probability assigned to the type class is

T(na)2nH0(S)|\mathcal{T}(n_a)| 2^{-n\mathcal{H}_0(S)}

Since this probability cannot exceed one,

T(na)2nH0(S)|\mathcal{T}(n_a)| \leq 2^{n\mathcal{H}_0(S)}

and therefore

log2n!aΣna!nH0(S)\log_2 \frac{n!}{\prod_{a\in\Sigma}n_a!} \leq n\mathcal{H}_0(S)

For the reverse direction, there are at most

(n+1)Σ(n+1)^{|\Sigma|}

possible type vectors. Under the multinomial distribution P^S\widehat{P}_S, the observed count vector (na)aΣ(n_a)_{a\in\Sigma} is a mode. Its probability must therefore be at least the reciprocal of the number of possible types:

T(na)2nH0(S)1(n+1)Σ|\mathcal{T}(n_a)|2^{-n\mathcal{H}_0(S)} \geq \frac{1}{(n+1)^{|\Sigma|}}

Rearranging gives

T(na)2nH0(S)(n+1)Σ|\mathcal{T}(n_a)| \geq \frac{2^{n\mathcal{H}_0(S)}}{(n+1)^{|\Sigma|}}

Taking logarithms produces

nH0(S)Σlog2(n+1)log2n!aΣna!nH0(S)n\mathcal{H}_0(S) - |\Sigma|\log_2(n+1) \leq \log_2 \frac{n!}{\prod_{a\in\Sigma}n_a!} \leq n\mathcal{H}_0(S)

This is the general form of the method of types. Thus the counting description and the fitted zero-order log-loss differ by at most Σlog2(n+1)|\Sigma|\log_2(n+1) bits. For a fixed alphabet, the difference per symbol is at most

Σlog2(n+1)n,\frac{|\Sigma|\log_2(n+1)}{n},

and therefore vanishes as nn grows. This is the precise sense in which zero-order empirical entropy has both a probabilistic and a combinatorial interpretation.

Adding Context

Zero-order empirical entropy ignores order: rearranging the symbols leaves every count nan_a unchanged. To let the preceding symbols matter, fix a context length 0k<n0\leq k<n. For each context ωΣk\omega\in\Sigma^k, let SωS_\omega be the sequence obtained by collecting, in their original order, all symbols that follow occurrences of ω\omega.

Each SωS_\omega can be fitted independently with its own zero-order model. The resulting kk-th order empirical entropy is

Hk(S)=1nωΣkSω>0SωH0(Sω).\mathcal{H}_k(S) = \frac{1}{n} \sum_{\substack{\omega\in\Sigma^k \\ |S_\omega|>0}} |S_\omega| \mathcal{H}_0(S_\omega).

In other words, positions are grouped by their preceding context, and the zero-order analysis is repeated inside each group.

To express this explicitly in terms of counts, for aΣa\in\Sigma and ωΣk\omega\in\Sigma^k define

nωa={i{k+1,,n}:sik:i1=ω, si=a}n_{\omega a} = \left| \left\{ i\in\{k+1,\ldots,n\} : s_{i-k:i-1}=\omega,\ s_i=a \right\} \right|

and let

nω=aΣnωan_\omega = \sum_{a\in\Sigma} n_{\omega a}

be the number of symbols observed after ω\omega.

Within this group, the empirical conditional distribution is

P^S(aω)=nωanω\widehat{P}_S(a\mid\omega) = \frac{n_{\omega a}}{n_\omega}

for every observed context with nω>0n_\omega>0. The same maximum-likelihood argument used in the zero-order case applies independently to each context. The smallest log-loss obtained by assigning one distribution to the symbols following ω\omega is

aΣnωa>0nωalog2nωnωa\sum_{\substack{a\in\Sigma \\ n_{\omega a}>0}} n_{\omega a} \log_2 \frac{n_\omega}{n_{\omega a}}

Summing over the observed contexts gives

nHk(S)=ωΣknω>0aΣnωa>0nωalog2nωnωan\mathcal{H}_k(S) = \sum_{\substack{\omega\in\Sigma^k \\ n_\omega>0}} \sum_{\substack{a\in\Sigma \\ n_{\omega a}>0}} n_{\omega a} \log_2 \frac{n_\omega}{n_{\omega a}}

which matches the kk-th order empirical entropy under the boundary convention that only positions with a complete length-kk context contribute to the sum.

The first kk symbols have no complete length-kk context. They may be encoded separately as one block in

log2Σk=klog2Σ\left\lceil \log_2|\Sigma|^k \right\rceil = \left\lceil k\log_2|\Sigma| \right\rceil

bits, or handled through an agreed boundary convention. When k=0k=0, there is only the empty context and the definition reduces to H0(S)\mathcal{H}_0(S).

Longer contexts divide the observed positions into smaller groups. Refining a group cannot increase the minimum fitted log-loss when both models are evaluated on the same set of positions, since the refined model can always reuse the distribution of the original group. More context can therefore reduce the empirical data term.

This reduction eventually exposes a limitation of the measure. If every observed context is followed by only one distinct symbol, then every SωS_\omega is constant and

H0(Sω)=0\mathcal{H}_0(S_\omega)=0

for all observed contexts. Consequently,

Hk(S)=0\mathcal{H}_k(S)=0

This always occurs at k=n1k=n-1, where only one position has a complete context, and it may occur much earlier when sufficiently long contexts determine their observed continuations.

The value Hk(S)=0\mathcal{H}_k(S)=0 does not give a zero-length lossless representation of SS. It gives a zero data term after the fitted context distributions are available. A decoder can assign probability one to an observed continuation only if it already knows which continuation followed that context. As the empirical loss decreases, information may have moved from the encoded sequence into the fitted model. The complete description must account for how the decoder obtains that model.

The Model Is Part of the Message

That missing information is the model itself.

Let MM contain everything the decoder needs to reproduce the probabilities used by the encoder. A two-part description has length

L(M,S)=L(M)+L(SM)L(M,S) = L(M)+L(S\mid M)

The term L(M)L(M) describes the model. The term L(SM)L(S\mid M) encodes the sequence using that model. In a sufficiently expressive model family, minimizing only the second term rewards memorization. A more complex model improves the complete description only when the reduction in L(SM)L(S\mid M) exceeds the additional cost of describing MM. This is the basic two-part form of the minimum description length principle.

For the order-kk model above, this cost is concrete. A naïve dense table contains σk\sigma^k contexts and σ\sigma continuation counts for each context. Storing every count in a fixed-width field gives an upper bound of

σk+1log2(n+1)\sigma^{k+1}\left\lceil\log_2(n+1)\right\rceil

bits. Sparse representations reduce this cost, but they must still identify which contexts and continuations occur.

How the model is charged depends on what the encoder and decoder already share. There are three relevant cases: the model may be fixed in advance, transmitted with the data, or reconstructed online while decoding.

If the model is fixed by a file format, built into the decoder, or otherwise shared in advance, it does not belong to the conditional description length of an individual message. It remains part of the system that makes that description meaningful, but charging its complete size to every message would also be misleading.

If it is fitted offline and not already available to the decoder, a description of its tables, parameters, or weights must accompany the encoded data. Language Modeling Is Compression calls the ratio obtained without parameter size the raw compression rate. Its adjusted compression rate adds the parameter size to the compressed output. A larger model may obtain a lower log-loss while producing a worse adjusted rate when it is amortized over too little data.

In prequential coding, encoder and decoder begin from the same initial state. Let Qi1Q_{i-1} be the model available after the prefix s<is_{<i} has been processed. The ideal prequential length is

Lpreq(S)=i=1nlog2Qi1(sis<i)L_{\mathrm{preq}}(S) = \sum_{i=1}^n -\log_2 Q_{i-1}(s_i\mid s_{<i})

After decoding sis_i, the decoder performs the same update as the encoder and reconstructs QiQ_i. The final parameters therefore need not be transmitted.

The model cost is paid through the online log-loss rather than through a separate description of the final parameters. Before sufficient data have been observed, the current model will typically predict less well than a model fitted to the complete sequence. More generally, the difference appears as prequential regret relative to that offline fit. The initialization, update rule, training procedure, numerical conventions, and any randomness affecting them must be shared or described. Any unshared information needed to reproduce the learning procedure must be added to the prequential length.

The data term also separates uncertainty in the source from mismatch in the model. Suppose objects are generated according to a distribution PP, while the compressor assigns probabilities according to QQ. If Q(x)>0Q(x)>0 whenever P(x)>0P(x)>0, the expected ideal data length is the cross-entropy

H(P,Q)=xP(x)log2Q(x)H(P,Q) = -\sum_x P(x)\log_2 Q(x)

Subtracting the source entropy gives

H(P,Q)H(P)=xP(x)log2P(x)Q(x)=DKL(PQ)\begin{aligned} H(P,Q)-H(P) &= \sum_x P(x)\log_2\frac{P(x)}{Q(x)} \\ &= D_{\mathrm{KL}}(P\mathbin\Vert Q) \end{aligned}

and therefore

H(P,Q)=H(P)+DKL(PQ)H(P,Q) = H(P)+D_{\mathrm{KL}}(P\mathbin\Vert Q)

For sequential distributions, relative entropy decomposes across positions:

DKL(P1:nQ1:n)=i=1nEX<iP[DKL(P(X<i)Q(X<i))]D_{\mathrm{KL}}\left(P_{1:n}\mathbin\Vert Q_{1:n}\right) = \sum_{i=1}^n \mathbb{E}_{X_{<i}\sim P} \left[ D_{\mathrm{KL}} \left( P(\cdot\mid X_{<i}) \mathbin\Vert Q(\cdot\mid X_{<i}) \right) \right]

Each term is the expected number of additional bits paid at one position because the model’s conditional distribution differs from the source distribution.

This distinction matters when discussing whether a better predictor “reduces entropy.” For a fixed source PP, improving QQ under expected logarithmic loss means reducing the cross-entropy H(P,Q)H(P,Q), equivalently reducing the mismatch term DKL(PQ)D_{\mathrm{KL}}(P\mathbin\Vert Q). It does not change H(P)H(P). The ngrok article correctly associates better probability estimates with shorter encodings, but its final use of entropy merges these two quantities.

The distribution used for evaluation also matters. If QSQ_S is fitted to an observed sequence SS, minimizing its in-sample log-loss does not imply that it will minimize log-loss on future data. For a future distribution RR,

EXR[log2QS(X)]=H(R)+DKL(RQS).\mathbb{E}_{X\sim R}\left[-\log_2 Q_S(X)\right] = H(R)+D_{\mathrm{KL}}(R\mathbin\Vert Q_S).

A model may therefore compress its training sequence more tightly while predicting future observations worse, either because it has overfit the sample or because the source has changed. Compression length and prediction log-loss remain the same quantity when evaluated on the same data, generalization between different samples or distributions is a separate question.

The Shortest Bitstream May Be the Wrong Representation

Every code considered so far has been judged by one operation: reconstructing the complete object. The two-part description length says nothing about what can be done with the encoded data before that reconstruction is complete.

A bit-packed vector and an entropy-coded stream make this difference concrete. Nothing in this factorization requires the sequence order to represent physical time. Consider a vector

I have implemented both bit-packed and entropy-coded representations of a vector in Rust, you can have a look here

A=(a0,,an1)A=(a_0,\ldots,a_{n-1})

whose values belong to {0,,u1}\{0,\ldots,u-1\}. Let

b=log2ub=\lceil\log_2 u\rceil

A bit-packed representation assigns exactly bb consecutive bits to each value, using nbnb bits for the payload apart from alignment and metadata. The representation of aia_i begins at bit position ibib. If the storage-word width is wbw\geq b, recovering aia_i requires reading at most two adjacent words, shifting their contents, and applying a mask. The addresses and shifts are computed directly from ii, so access takes O(1)O(1) time.

This representation does not exploit differences in frequency. Every value receives the same number of bits. If the values follow a non-uniform distribution, or if their probabilities depend on earlier values, an entropy coder may produce a shorter stream:

L(AM)i=0n1log2Q(aia<i)L(A\mid M) \approx \sum_{i=0}^{n-1} -\log_2 Q(a_i\mid a_{<i})

The shorter stream provides a different access contract. In an ordinary arithmetic-coded stream, the decoding state at position ii depends on the symbols that precede it. If the model also uses their context, its next distribution depends on the same prefix. Recovering aia_i requires decoding from the beginning of the stream or from an earlier checkpoint whose coding and model states have been stored.

The bit-packed vector may occupy more space while answering vector[i] directly. Its compression comes from restricting the possible value at each position to an alphabet of size uu, rather than from predicting which value will occur. When all values remain equally plausible, fixed-width packing uses the information supplied by that restriction without requiring a non-uniform model.

A compressed archive needs an encoding EE and a decoder DD satisfying

D(E(A))=AD(E(A))=A

Once the whole vector can be reconstructed, the encoding has fulfilled its contract. A compressed representation may be required to satisfy a stronger condition. It must support an access algorithm such that

Access(R(A),i)=ai\operatorname{Access}(R(A),i)=a_i

without first reconstructing all of AA.

Checkpoints can give an entropy-coded stream faster access, but each checkpoint occupies space. Smaller blocks reduce the amount of decoding required for an access and increase the number of stored states. Larger blocks save metadata and increase access time. The objective is no longer to minimize the bitstream without qualification. It becomes

minRR(A)subject toTAccess(R)t\min_R |R(A)| \qquad \text{subject to} \qquad T_{\operatorname{Access}}(R)\leq t

for a chosen access-time bound tt.

The operation need not be random access, and the object need not be a vector. The same issue arises whenever compressed data must be searched, traversed, compared, or partially decoded. The required operations constrain which short descriptions are useful and how much auxiliary information they need.

Bitstream length alone answers only the bare reconstruction problem. A compressed representation must also encode enough structure for its required operations. The decoder contract now has three explicit parts: the objects it must distinguish, the information it already shares with the encoder, and the operations it must perform without full reconstruction. Prediction determines conditional code lengths inside this contract. It does not determine the contract itself.

So, Is Compression Prediction?

So, is compression prediction? Once the coding problem has been fixed, yes, but in a precise sense. The equivalence concerns the ideal length of the encoded payload under an agreed model.

The qualification matters because a description-length problem can be posed before a next-symbol predictor exists. Kolmogorov complexity does so for an individual string, while the counting bound does so for a finite family of objects.

What I felt the other explanations left implicit was the compression problem that must be fixed before the equivalence becomes meaningful. The encoder and decoder need an agreed family of objects, a serialization, a boundary between transmitted and shared information, and a decoding contract. None of these choices is determined by next-symbol prediction.

Once those choices have been made, the equivalence is exact at the level of code lengths. A sequential probability model assigns ideal payload lengths through logarithmic loss. Conversely, every uniquely decodable code over a fixed object family induces the distribution

QC(x)=2C(x)SC,Q_C(x)=\frac{2^{-\ell_C(x)}}{S_C},

and its codeword lengths satisfy

C(x)=log2QC(x)log2SC.\ell_C(x)=-\log_2Q_C(x)-\log_2S_C.

After the objects have been serialized, the induced distribution can in turn be factored into next-symbol conditionals.

That generality is also the limit of the slogan Compression is Prediction. Recasting a representation probabilistically does not explain why its objects were chosen, whether its model must be transmitted, or which operations it supports. A model with lower log-loss can produce a larger complete file after its parameters are included. An empirical entropy of zero can still leave the decoder without the model needed to reconstruct the sequence.

Compression is therefore prediction after the coding problem has been fixed, and only at the level measured by the induced code lengths. For a shared sequential model, cumulative log-loss gives the ideal payload length up to coding overhead. It does not define what must be represented, what the decoder already knows, or what the representation must allow the decoder to do.


Share this post on:

Next Post
Who Owns the Memory? Part 3: How Big Is your Type?