How to use this course
Leonid Kantorovich is the only person to have received the Nobel Memorial Prize in Economic Sciences for work that began as pure mathematics. Economists know him for linear programming and for the idea that an optimal plan carries its own prices. Analysts know him for Kantorovich spaces, the Newton–Kantorovich theorem and the Kantorovich inequality. Probabilists and geometers know the Kantorovich duality and the Kantorovich–Rubinstein distance, which sit at the centre of modern optimal transport. This course presents all three, alongside a career shaped by the Soviet system: early fame, two decades of official suspicion of his economics, then late honours.
- Linear algebra and multivariable calculus.
- Basic real analysis: normed spaces, completeness, compactness.
- For Lessons 05, 08 and 09, some measure theory (the Lebesgue integral and Borel measures) is helpful.
- No economics is assumed.
Each technical lesson opens with a paragraph for the general reader, then gives precise statements, proof sketches and exercises with folded solutions. Lessons 07 to 10 can be read independently of 04 to 06 by a reader mainly interested in optimization and transport.
Timeline
| 1912 | Born 19 January in St Petersburg, the youngest child of Vitaly Kantorovich, a physician, and Paulina Zaks. |
| 1922 | Death of his father. |
| 1926 | Enters Leningrad State University, aged 14. |
| 1929–30 | Research on analytic and projective sets, solving problems of Nikolai Luzin; lecture at the First All-Union Mathematical Congress in Kharkov (1930). |
| 1930 | Graduates from Leningrad State University; begins teaching. |
| 1934 | Full professor, aged 22. |
| 1935 | Doctor of Science when the USSR restores academic degrees. First papers on semi-ordered spaces; meets John von Neumann and others at the Moscow Topological Conference. |
| 1936 | Approximate Methods of Higher Analysis, with Vladimir Krylov. |
| 1938 | Consults for the Plywood Trust in Leningrad. Marries Natalya Ilyina, a physician. |
| 1939 | Mathematical Methods of Organizing and Planning Production, the founding text of linear programming. |
| 1941–42 | During the siege of Leningrad computes safe vehicle spacing for the ice road across Lake Ladoga. |
| 1942 | On the translocation of masses, in the Doklady of the USSR Academy of Sciences. Writes The Best Use of Economic Resources, unpublished until 1959. |
| 1944 | Returns to Leningrad; heads the Department of Approximate Methods at the Academy’s Mathematical Institute. |
| 1948 | Functional analysis and applied mathematics in Uspekhi. Assigned to computations for the Soviet atomic project. |
| 1949 | Stalin (State) Prize. |
| 1951 | Rational Cutting of Industrial Materials, with Victor Zalgaller. |
| 1958 | Corresponding Member of the USSR Academy of Sciences. Kantorovich–Rubinstein paper on the transport metric. The 1939 booklet appears in English (1960, Management Science). |
| 1959 | The Best Use of Economic Resources; Functional Analysis in Normed Spaces, with Gleb Akilov. |
| 1960 | Moves to Novosibirsk and the Siberian Branch of the Academy. |
| 1964 | Full Member (Academician) of the USSR Academy of Sciences. |
| 1965 | Lenin Prize, with Vasily Nemchinov and Viktor Novozhilov. |
| 1971 | Moves to Moscow to head research at the Institute of National Economic Management. |
| 1975 | Nobel Memorial Prize in Economic Sciences, with Tjalling Koopmans. |
| 1986 | Dies in Moscow on 7 April, aged 74; buried at Novodevichy Cemetery. |
A Petrograd childhood, 1912–1926
Leonid Vitalyevich Kantorovich was born in St Petersburg on 19 January 1912 into a Jewish family. His father, Vitaly Moiseevich Kantorovich, was a doctor; his mother was Paulina Grigoryevna Zaks. He was the youngest of five children. Two of his elder brothers became physicians.
His childhood coincided with war, revolution and civil war. In his Nobel autobiography he recalled the February and October Revolutions of 1917 and a year spent in Belorussia during the Civil War. His father died in 1922, and his mother brought him up. His interest in science appeared around 1920, and it developed quickly enough that in 1926, at the age of fourteen, he entered the mathematics department of Leningrad State University.
Professor at twenty-two, 1926–1937
At the university his teachers included Vladimir Smirnov, Grigorii Fichtenholz and Boris Delaunay. Among his friends were Sergei Sobolev, Solomon Mikhlin, Isidor Natanson and the Faddeevs. He also attended the history lectures of Evgeny Tarle, which first drew him to political economy. He began research in his second year in Fichtenholz’s seminar, on the descriptive set theory then being created in Moscow by Nikolai Luzin. In 1929–30 he solved several of Luzin’s problems on analytic and projective sets and reported them at the First All-Union Mathematical Congress in Kharkov in 1930, where he met Sergei Bernstein, Pavel Alexandrov, Andrey Kolmogorov and Aleksandr Gelfond, and foreign guests including Jacques Hadamard and Paul Montel.
He graduated in 1930, aged eighteen. The Soviet Union had abolished academic degrees, so he began teaching at once: at the Naval Engineering School, at the Leningrad Institute of Industrial Construction Engineering and at the university. His applied work of these years produced new methods for conformal mapping and for variational problems, collected in Approximate Methods of Higher Analysis (1936) with Vladimir Krylov. He became a full professor in 1934. When degrees were restored in 1935 he received the doctorate of science.
From 1935 he developed a theory of partially ordered function spaces (Lesson 05). At the Moscow Topological Conference that year he met John von Neumann, Garrett Birkhoff, Albert Tucker and Maurice Fréchet. Torsten Carleman invited a memoir on functional equations for Acta Mathematica. In 1937 he won first prize in a Leningrad competition for young researchers.
The plywood trust and the war, 1938–1944
In 1938 the laboratory of the Leningrad Plywood Trust asked the young professor for help. Several machines could each process several kinds of raw material at different rates, and the trust wanted to distribute the work so as to maximize output in a fixed assortment. Mathematically this meant maximizing a linear function over a convex polytope with an enormous number of vertices. Kantorovich found a general method, which he called the method of resolving multipliers, and saw at once that the same structure appeared in many planning problems: allocating crops to fields, cutting sheets with least waste, routing freight. In July 1939 Leningrad University published his booklet Mathematical Methods of Organizing and Planning Production. It contained the essential ideas of linear programming, eight years before George Dantzig’s simplex method.
The booklet was hardly noticed outside Leningrad, and inside the Soviet Union its implications were unwelcome. The resolving multipliers behaved like prices, and the claim that an optimal plan determines prices sounded too close to marginalist economics, which Soviet ideology rejected. Official criticism followed, and his economic writing stalled for most of two decades.
The road of life
In June 1941 Germany invaded the Soviet Union, and by September Leningrad was besieged. Kantorovich, now a professor at the Higher Engineering-Technical School of the Navy, was given the rank of major. The only supply route into the city in winter ran across the ice of Lake Ladoga. He calculated the safe distance between vehicles as a function of ice thickness and air temperature, and in December 1941 and January 1942 he walked the ice between the lorries to check that the spacing held. He received the Order of the Patriotic War and the medal for the Defence of Leningrad. His school was later evacuated to Yaroslavl, where he wrote a probability textbook for military engineers.
In 1942 two of his most lasting works appeared or were written. The Academy’s Doklady published the four-page note On the translocation of masses (Lesson 08), and he completed the manuscript of The Best Use of Economic Resources, which would wait seventeen years for publication.
Approximation and variational methods
For the general reader. Many problems of engineering ask for a function that makes some energy as small as possible. Exact solutions are rare. In the early 1930s Kantorovich found ways to reduce such problems in two variables to ordinary differential equations in one, and to approximate rough functions by smooth polynomials while keeping control of averages.
- Define the Kantorovich polynomial operators and see why they converge in mean.
- Apply the Kantorovich method of reduction to ordinary differential equations.
- Solve the torsion problem for a rectangular bar approximately.
Kantorovich polynomials (1930)
Sergei Bernstein’s polynomials \(B_nf(x)=\sum_{k=0}^n f(k/n)\,p_{n,k}(x)\), with \(p_{n,k}(x)=\binom nk x^k(1-x)^{n-k}\), converge uniformly to any continuous \(f\) on \([0,1]\). They use point values and so make no sense for an integrable function. Kantorovich replaced the values by local averages:
\[ K_nf(x)=(n+1)\sum_{k=0}^{n}\Big(\int_{k/(n+1)}^{(k+1)/(n+1)} f(t)\,dt\Big)\,p_{n,k}(x). \]If \(f\) is continuous on \([0,1]\), then \(K_nf\to f\) uniformly. If \(f\in L^p[0,1]\) with \(1\le p\lt \infty\), then \(\|K_nf-f\|_p\to0\).
The operators are positive, preserve constants, and satisfy \(\|K_n f\|_1\le\|f\|_1\) because \(\int_0^1p_{n,k}=\frac{1}{n+1}\). Uniform convergence for continuous \(f\) follows from the Korovkin-type argument used for Bernstein polynomials. Mean convergence then follows by density of continuous functions and the uniform bound on \(\|K_n\|\).
The Kantorovich method
Consider a functional \(J[u]=\iint_\Omega F(x,y,u,u_x,u_y)\,dx\,dy\) on a region \(a\le x\le b\), \(\varphi_1(x)\le y\le\varphi_2(x)\). The Ritz method fixes all coordinate functions and varies only numerical coefficients. Kantorovich fixed functions of \(y\) only and let the coefficients be unknown functions of \(x\):
\[ u(x,y)\approx\sum_{k=1}^{N} f_k(x)\,g_k(x,y), \]where each \(g_k\) already satisfies the boundary condition in \(y\). Integrating out \(y\) turns \(J\) into a one-dimensional functional of \(f_1,\dots,f_N\), whose Euler–Lagrange equations are ordinary differential equations. The method keeps more of the solution exact than Ritz does, since the \(x\)-dependence is found, not guessed.
Example: torsion of a rectangular bar
The Prandtl stress function solves \(\Delta u=-2\) in \(|x|\le a,\ |y|\le b\) with \(u=0\) on the boundary; equivalently it minimizes \(J[u]=\iint(u_x^2+u_y^2-4u)\). Take \(u=(b^2-y^2)f(x)\). Integrating in \(y\),
\[ J=\int_{-a}^{a}\Big(\tfrac{16b^5}{15}f'^2+\tfrac{8b^3}{3}f^2-\tfrac{16b^3}{3}f\Big)dx, \qquad\text{so}\qquad f''-k^2f=-k^2,\quad k^2=\frac{5}{2b^2}. \]With \(f(\pm a)=0\), the solution is \(f(x)=1-\dfrac{\cosh kx}{\cosh ka}\).
Exercise 4.1Show that \(K_n1=1\) and \(K_n t=\frac{nx}{n+1}+\frac{1}{2(n+1)}\). Deduce that \(K_n t\to t\) uniformly.
Show solution
Exercise 4.2Verify the reduced functional and the ODE in the torsion example. What happens as \(a\to\infty\), and why is this the right answer?
Show solution
Ordered vector spaces
For the general reader. Most spaces of functions have a natural order: one function is below another if it is smaller at every point. Kantorovich made this order a structural axiom and showed that much of analysis, including limits, extensions of linear maps and even parts of numerical methods, can be done using the order alone. The spaces he singled out are now called Kantorovich spaces.
A vector lattice (or Riesz space) is a real vector space \(E\) with a partial order compatible with addition and positive scaling, in which every pair \(x,y\) has a supremum \(x\vee y\). Then \(x^+=x\vee0\), \(x^-=(-x)\vee0\) and \(|x|=x^++x^-\) are defined. \(E\) is Dedekind complete, a K-space in Kantorovich’s terminology, if every non-empty subset that is bounded above has a supremum. Examples are \(\mathbb R^n\) with the coordinatewise order, \(L^p(\mu)\) for \(1\le p\le\infty\) with the almost-everywhere order, and \(\ell^\infty\). The space \(C[0,1]\) is a vector lattice but not Dedekind complete. Kantorovich introduced the notion in 1935, at the same time as related work by Garrett Birkhoff in the United States, Hans Freudenthal in the Netherlands and Hidegorô Nakano in Japan, following an earlier lecture of Frigyes Riesz (1928).
Let \(X\) be a real vector space, \(E\) a Dedekind complete vector lattice, and \(p:X\to E\) sublinear: \(p(x+y)\le p(x)+p(y)\) and \(p(tx)=tp(x)\) for \(t\ge0\). If \(T_0:X_0\to E\) is linear on a subspace \(X_0\) with \(T_0\le p\) on \(X_0\), then \(T_0\) extends to a linear \(T:X\to E\) with \(T\le p\) on \(X\).
Proof sketch. As in the scalar case, extend one dimension at a time. To extend from \(X_0\) to \(X_0+\mathbb Rz\) one needs \(c\in E\) with \(T_0y-p(y-z)\le c\le p(y'+z)-T_0y'\) for all \(y,y'\in X_0\). Sublinearity shows every left-hand element is below every right-hand element, and Dedekind completeness supplies \(c=\sup_y(T_0y-p(y-z))\). Zorn’s lemma finishes the argument. Conversely, it was later shown that the extension property for all such \(X,p\) forces \(E\) to be Dedekind complete, so Kantorovich’s condition is exactly the right one.
Let \(E,F\) be vector lattices with \(F\) Dedekind complete. The order-bounded linear operators \(L_b(E,F)\) form a Dedekind complete vector lattice, and for \(x\ge0\)
\[ T^+(x)=\sup\{Ty: 0\le y\le x\},\qquad |T|(x)=\sup\{|Ty|: |y|\le x\}. \]Exercise 5.1Show that \(C[0,1]\) is not Dedekind complete by finding a bounded set of continuous functions with no least upper bound in \(C[0,1]\).
Show solution
Exercise 5.2On \(\mathbb R^2\) with the coordinatewise order, let \(T(x_1,x_2)=(x_1-x_2,\ 0)\). Compute \(T^+\), \(T^-\) and \(|T|\) using the Riesz–Kantorovich formula.
Show solution
Functional analysis and numerical mathematics
For the general reader. In 1948 Kantorovich published a long paper arguing that the abstract language of Banach spaces was the right setting for computation. An iterative method is an operator; its convergence is a statement about norms. Two results from this program are used every day in scientific computing: a guarantee that Newton’s method converges, checkable at the starting point, and a sharp bound for the speed of steepest descent.
- State and apply the Newton–Kantorovich theorem.
- Prove the Kantorovich inequality.
- Derive the convergence rate of steepest descent.
Let \(X,Y\) be Banach spaces, \(D\subseteq X\) open and convex, and \(F:D\to Y\) Fréchet differentiable with \(\|F'(x)-F'(y)\|\le K\|x-y\|\) on \(D\). Let \(x_0\in D\) with \(F'(x_0)\) invertible, \(\|F'(x_0)^{-1}\|\le\beta\) and \(\|F'(x_0)^{-1}F(x_0)\|\le\eta\). If
\[ h=\beta K\eta\le\tfrac12 \quad\text{and}\quad \overline B(x_0,t^*)\subseteq D,\qquad t^*=\frac{1-\sqrt{1-2h}}{\beta K}, \]then the Newton iterates \(x_{n+1}=x_n-F'(x_n)^{-1}F(x_n)\) are defined, stay in \(\overline B(x_0,t^*)\), and converge to a zero \(x^*\) of \(F\), with \(\|x^*-x_n\|\le 2^{1-n}(2h)^{2^n-1}\eta\). If \(h\lt \tfrac12\), the zero is unique in \(D\cap B(x_0,t^{**})\), \(t^{**}=\frac{1+\sqrt{1-2h}}{\beta K}\).
Proof idea. Compare the iterates with Newton’s method for the scalar quadratic \(\phi(t)=\tfrac{\beta K}{2}t^2-t+\eta\), whose smaller root is \(t^*\). Starting at \(t_0=0\), its Newton iterates \(t_n\) increase to \(t^*\). Induction using the Lipschitz bound on \(F'\) and the Banach lemma on perturbed inverses gives \(\|x_{n+1}-x_n\|\le t_{n+1}-t_n\). So \((x_n)\) is Cauchy, and the quadratic \(\phi\) is a majorant. This method of majorants is Kantorovich’s main idea.
Let \(A\) be a symmetric positive definite \(n\times n\) matrix with eigenvalues in \([m,M]\), \(0\lt m\). Then for every \(x\ne0\)
\[ \frac{(x^{\mathsf T}Ax)(x^{\mathsf T}A^{-1}x)}{(x^{\mathsf T}x)^2}\le\frac{(M+m)^2}{4Mm}. \]Proof. Diagonalize, so with \(w_i=x_i^2/\|x\|^2\) the left side is \(\big(\sum w_i\lambda_i\big)\big(\sum w_i\lambda_i^{-1}\big)\). On \([m,M]\) the convex function \(1/\lambda\) lies below its chord: \(\lambda^{-1}\le\frac{m+M-\lambda}{mM}\). Hence \(\sum w_i\lambda_i^{-1}\le\frac{m+M-\bar\lambda}{mM}\) with \(\bar\lambda=\sum w_i\lambda_i\), and \(\bar\lambda(m+M-\bar\lambda)\le\frac{(m+M)^2}{4}\).
Consequence. For \(f(x)=\tfrac12x^{\mathsf T}Ax-b^{\mathsf T}x\), steepest descent with exact line search satisfies
\[ f(x_{k+1})-f(x^*)\le\Big(\frac{M-m}{M+m}\Big)^2\big(f(x_k)-f(x^*)\big). \]Exercise 6.1Apply the Newton–Kantorovich theorem to \(F(x)=x^2-2\) on \(D=\mathbb R\) with \(x_0=1.5\). Compute \(\beta,\eta,K,h\) and the ball that must contain \(\sqrt2\).
Show solution
Exercise 6.2Show the Kantorovich bound is attained for \(A=\operatorname{diag}(m,M)\) and a suitable \(x\). For \(m=1\), \(M=100\), how many steepest-descent steps guarantee a reduction of the error in \(f\) by a factor \(10^{-6}\)?
Show solution
Linear programming, 1939
For the general reader. A factory has machines, materials and orders, and it wants the best plan. When every constraint and the goal are linear, the possible plans form a polyhedron, and the best one sits at a corner. Kantorovich showed in 1939 how to find it, and how to certify that a plan is optimal by attaching a number, a multiplier, to every constraint.
- Formulate Kantorovich’s production problem as a linear program.
- State the duality theorem and complementary slackness.
- Interpret resolving multipliers as dual variables.
The basic problem of 1939
In modern notation, the central problem of the booklet reads: machines \(i=1,\dots,n\) can each produce parts \(k=1,\dots,m\); machine \(i\) working a full day on part \(k\) makes \(\alpha_{ik}\) units. Parts are needed in proportions \(p_1:\dots:p_m\) to assemble finished products. Choose fractions of time \(h_{ik}\ge0\) with \(\sum_k h_{ik}=1\) for each machine to maximize the number \(z\) of complete sets:
\[ \max z\quad\text{subject to}\quad \sum_{i=1}^n\alpha_{ik}h_{ik}=p_kz\ \ (k=1,\dots,m),\qquad \sum_{k=1}^m h_{ik}=1,\qquad h_{ik}\ge0. \]Kantorovich’s criterion: a plan is optimal if and only if there are multipliers \(\lambda_k\ge0\), one per part, such that each machine works only on parts where \(\lambda_k\alpha_{ik}\) is maximal over \(k\). The \(\lambda_k\) value the parts; each machine then does what is most valuable for it. This is the complementary slackness condition of linear programming duality.
For \(A\in\mathbb R^{m\times n}\), \(b\in\mathbb R^m\), \(c\in\mathbb R^n\), consider
\[ \text{(P)}\ \ \max\{c^{\mathsf T}x: Ax\le b,\ x\ge0\},\qquad \text{(D)}\ \ \min\{b^{\mathsf T}y: A^{\mathsf T}y\ge c,\ y\ge0\}. \]If both are feasible, both have optimal solutions and the optimal values are equal. Feasible \(x,y\) are both optimal if and only if \(y_i(b-Ax)_i=0\) for all \(i\) and \(x_j(A^{\mathsf T}y-c)_j=0\) for all \(j\).
Proof sketch. Weak duality is one line: \(c^{\mathsf T}x\le(A^{\mathsf T}y)^{\mathsf T}x=y^{\mathsf T}Ax\le y^{\mathsf T}b\). Equality is the hard part. It follows from Farkas’s lemma, a separating-hyperplane statement for the closed convex cone generated by the columns of \(A\). The general duality theorem was stated by John von Neumann in 1947 and proved by David Gale, Harold Kuhn and Albert Tucker in 1951. Kantorovich had the optimality criterion for his class of problems in 1939 and developed the method further in a 1940 note in the Doklady.
Exercise 7.1Two machines make parts A and B, needed in equal numbers. Machine 1 makes 10 A or 20 B per day; machine 2 makes 30 A or 30 B per day. Find the optimal plan and verify it with multipliers.
Show solution
Exercise 7.2Write the dual of \(\max\{3x_1+2x_2: x_1+x_2\le4,\ x_1+3x_2\le6,\ x\ge0\}\), solve both, and check complementary slackness.
Show solution
The translocation of masses
For the general reader. In 1781 Gaspard Monge asked how to move a pile of earth into a hole of the same volume with least total effort. His problem resisted solution for more than a century because each grain had to go to a single destination. Kantorovich, in a four-page note of 1942, let mass be split. The problem then became linear, had a dual problem, and acquired a simple optimality test. Modern optimal transport grows from this note.
- State the Monge and Kantorovich formulations.
- State Kantorovich duality and the 1942 potential criterion.
- Solve one-dimensional and discrete examples.
Let \(\mu,\nu\) be Borel probability measures on spaces \(X,Y\) and \(c:X\times Y\to[0,\infty)\) a cost. Monge’s problem: minimize \(\int c(x,T(x))\,d\mu(x)\) over maps \(T\) with \(T_\#\mu=\nu\). Kantorovich’s problem: minimize over couplings
\[ \min_{\pi\in\Pi(\mu,\nu)}\int_{X\times Y}c(x,y)\,d\pi(x,y),\qquad \Pi(\mu,\nu)=\{\pi\ge0:\ \pi(A\times Y)=\mu(A),\ \pi(X\times B)=\nu(B)\}. \]\(\Pi(\mu,\nu)\) is convex and never empty (it contains \(\mu\otimes\nu\)), the objective is linear, and when \(X,Y\) are compact and \(c\) continuous the set is weakly compact, so a minimizer exists. Kantorovich worked on a compact metric space in exactly this setting.
Let \(X,Y\) be Polish spaces and \(c\) lower semicontinuous and bounded below. Then
\[ \min_{\pi\in\Pi(\mu,\nu)}\int c\,d\pi=\sup\Big\{\int\varphi\,d\mu+\int\psi\,d\nu:\ \varphi(x)+\psi(y)\le c(x,y)\ \text{for all }x,y\Big\}. \]In the case \(X=Y\), \(c=d\) a metric on a compact space (Kantorovich’s 1942 setting): a coupling \(\pi\) is optimal if and only if there is a potential \(U\) with \(|U(x)-U(y)|\le d(x,y)\) for all \(x,y\) and \(U(y)-U(x)=d(x,y)\) whenever \((x,y)\) lies in the support of \(\pi\).
Why the criterion works. If such \(U\) exists, then for any coupling \(\pi'\), \(\int d\,d\pi'\ge\int(U(y)-U(x))\,d\pi'=\int U\,d\nu-\int U\,d\mu\), and this lower bound is achieved by \(\pi\). So \(\pi\) is optimal. The converse, building \(U\) from an optimal plan, is the substance of the 1942 note. Kantorovich read \(U\) as a price: moving a unit of mass from \(x\) to \(y\) raises its value by exactly the transport cost along used routes, and never by more.
One dimension
On \(\mathbb R\) with \(c(x,y)=h(x-y)\), \(h\) convex, the monotone rearrangement \(T=F_\nu^{-1}\circ F_\mu\) is optimal, where \(F\) are distribution functions. For \(c=|x-y|\), the optimal cost is \(\int_{\mathbb R}|F_\mu(t)-F_\nu(t)|\,dt\).
Exercise 8.1Let \(\mu=\tfrac12(\delta_0+\delta_1)\) and \(\nu=\tfrac12(\delta_1+\delta_2)\) on \(\mathbb R\) with \(c=|x-y|\). Find all optimal couplings and a potential \(U\).
Show solution
Exercise 8.2Show that Monge’s problem can have no admissible map while Kantorovich’s always has a coupling: take \(\mu=\delta_0\) and \(\nu=\tfrac12(\delta_{-1}+\delta_1)\).
Show solution
Exercise 8.3For discrete \(\mu=\sum_i a_i\delta_{x_i}\) and \(\nu=\sum_j b_j\delta_{y_j}\), write the Kantorovich problem as a linear program and identify its dual.
Show solution
The Kantorovich–Rubinstein metric
For the general reader. How far apart are two probability distributions? Kantorovich proposed in 1942 to answer with the least work needed to move one into the other. With Gennady Rubinstein in 1958 he turned this into a norm on spaces of measures. The same distance, under the names earth mover’s distance and Wasserstein distance, now compares images, trains generative neural networks and measures curvature of metric spaces.
- Define \(W_1\) and prove it is a metric.
- State the Kantorovich–Rubinstein duality.
- Compute \(W_1\) in examples.
For a metric space \((X,d)\) and probability measures with finite first moment, define
\[ W_1(\mu,\nu)=\min_{\pi\in\Pi(\mu,\nu)}\int_{X\times X}d(x,y)\,d\pi(x,y). \]\(W_1\) is a metric on probability measures with finite first moment, and
\[ W_1(\mu,\nu)=\sup\Big\{\int f\,d\mu-\int f\,d\nu:\ |f(x)-f(y)|\le d(x,y)\ \text{for all }x,y\Big\}. \]On a compact space, \(W_1\) metrizes weak convergence of measures.
From the general duality. When \(c=d\), any admissible pair \((\varphi,\psi)\) can be improved to one of the form \((f,-f)\) with \(f\) 1-Lipschitz. Set \(f(x)=\inf_y\big(d(x,y)-\psi(y)\big)\). Then \(f\ge\varphi\), \(f\) is 1-Lipschitz as an infimum of 1-Lipschitz functions, and taking \(x=y\) gives \(f(y)\le-\psi(y)\), so \(\psi\le-f\). Hence \(\int\varphi\,d\mu+\int\psi\,d\nu\le\int f\,d\mu-\int f\,d\nu\). The two families of constraints collapse into one Lipschitz condition. The triangle inequality for \(W_1\) follows from the dual formula at once, since a supremum of seminorms is a seminorm; or from the primal formula by gluing couplings. Rubinstein and Kantorovich used the dual to extend \(W_1\) to a norm on all signed measures of total mass zero.
Exercise 9.1On \(\mathbb R\), show \(W_1(\delta_a,\delta_b)=|a-b|\) using both the primal and the dual formula.
Show solution
Exercise 9.2Compute \(W_1\) between the uniform distribution on \([0,1]\) and \(\delta_0\), and between the uniform distributions on \([0,1]\) and \([0,2]\).
Show solution
Prices from optimal plans
For the general reader. Kantorovich’s most controversial idea was economic. An optimal plan for a whole economy, he argued, automatically produces a set of valuations for every resource: the extra output one more unit would allow. He called them objectively determined valuations. Economists in the West recognized them as shadow prices; Soviet critics heard an echo of market economics. The mathematics is the sensitivity theory of linear programming.
- Interpret dual variables as marginal values of resources.
- See how valuations price scarce and abundant resources.
- Formulate the cutting-stock problem studied with Zalgaller.
Let \(v(b)=\max\{c^{\mathsf T}x:Ax\le b,\ x\ge0\}\). Then \(v\) is concave and piecewise linear on its domain. If the dual problem at \(b\) has a unique optimal solution \(y^*\), then \(v\) is differentiable at \(b\) and \(\nabla v(b)=y^*\). In general the set of dual optima is the superdifferential of \(v\) at \(b\).
Proof sketch. By duality \(v(b)=\min\{b^{\mathsf T}y: y\in P\}\), where \(P=\{y\ge0: A^{\mathsf T}y\ge c\}\) does not depend on \(b\). A minimum of linear functions of \(b\) is concave and, with finitely many vertices of \(P\), piecewise linear. For any dual-optimal \(y^*\) at \(b\) and any \(b'\), \(v(b')\le b'^{\mathsf T}y^*=v(b)+(b'-b)^{\mathsf T}y^*\): that is the supergradient inequality.
Complementary slackness carries the economic message. A resource not fully used in the optimal plan has valuation 0. An activity that is used breaks even at these valuations, and an activity that would make a loss at them is not used. In The Best Use of Economic Resources (1959) Kantorovich applied this to rent for land and natural resources, to a normative efficiency coefficient for capital investment and to the pricing of transport, each a question on which official Soviet doctrine was sensitive.
Rational cutting
From 1948 to 1950 Kantorovich and the geometer Victor Zalgaller applied these methods at the Leningrad Carriage-Building Works to cut steel sheet with least waste. Given cutting patterns \(j\), each yielding \(a_{ij}\) pieces of type \(i\), and demands \(d_i\), one solves
\[ \min\sum_j x_j\quad\text{subject to}\quad \sum_j a_{ij}x_j\ge d_i,\ \ x_j\ge0. \]The dual variables price each piece type, and a new pattern is worth adding exactly when its pieces are worth more than one sheet at those prices. This test is the idea behind column generation, rediscovered for cutting stock by Paul Gilmore and Ralph Gomory in 1961. Their 1951 book also combined linear programming with stage-by-stage optimization, independently of Richard Bellman’s dynamic programming.
Exercise 10.1In Exercise 7.2, compute \(v(b)\) for \(b=(4+\varepsilon,6)\) and \(b=(4,6+\varepsilon)\) with small \(\varepsilon>0\), and compare with the dual solution \(y^*=(3,0)\).
Show solution
Exercise 10.2Rolls of width 10 must be cut into pieces of width 3 (demand 30) and 4 (demand 20). List the efficient patterns, solve the LP relaxation, and give the dual prices.
Show solution
Recognition, 1944–1975
Kantorovich returned to Leningrad in 1944 and combined his chair at the university with the leadership of the Department of Approximate Methods at the Academy’s Mathematical Institute in Leningrad. His 1948 paper on functional analysis and applied mathematics won him the Stalin (State) Prize in 1949. In 1948 he was assigned to computations for the Soviet atomic project. With Mark Gavurin and Vera Faddeeva he worked on computer design and programming automation. Functional Analysis in Semi-ordered Spaces (1950, with Boris Vulikh and Aron Pinsker) and Functional Analysis in Normed Spaces (1959, with Gleb Akilov) became standard texts.
His economics remained under a cloud until after Stalin’s death. In 1956, when Tjalling Koopmans asked about his work, the 1939 booklet was sent abroad; it appeared in English in Management Science in 1960, and the 1942 transport note had been reprinted there in 1958. Kantorovich learned only in the mid-1950s that Koopmans, George Dantzig and Frank Hitchcock had reached similar results independently in the West. At home, The Best Use of Economic Resources finally appeared in 1959. A national conference on mathematical methods in economics in 1960 backed the new direction, though attacks from orthodox economists continued.
He was elected a Corresponding Member of the USSR Academy of Sciences in 1958 and moved in 1960 to the new Siberian Branch at Novosibirsk, where he headed a mathematical economics group with Vasily Nemchinov and founded the Department of Computational Mathematics at Novosibirsk State University. He became a full Academician in 1964. In 1965 he shared the Lenin Prize with Nemchinov and Viktor Novozhilov for the application of mathematical methods to economics. Among his doctoral students across these decades were Gleb Akilov, Victor Havin, Mikhail Solomyak, Gennady Rubinstein, Iosif Romanovskii and Valery Makarov.
DiscussionLinear programming was discovered independently in Leningrad (1939) and in the United States (1947). Why did the second discovery spread so much faster? Consider publication language, military sponsorship, computing machines and ideology.
Last years and legacy, 1971–1986
In 1971 Kantorovich moved to Moscow to lead research at the Institute of National Economic Management and to advise government planning bodies. In 1975 the Royal Swedish Academy of Sciences awarded him and Tjalling Koopmans the Sveriges Riksbank Prize in Economic Sciences in Memory of Alfred Nobel, for their contributions to the theory of optimum allocation of resources. His Nobel lecture, given on 11 December 1975, was titled Mathematics in Economics: Achievements, Difficulties, Perspectives. He received honorary doctorates from Glasgow, Grenoble, Nice, Helsinki, Paris and Cambridge, among others, and was elected to the Hungarian Academy of Sciences (1967) and the American Academy of Arts and Sciences (1969).
He died of cancer in Moscow on 7 April 1986, aged 74, and is buried at Novodevichy Cemetery. His wife, the physician Natalya Ilyina, and their two children survived him; both children became mathematical economists.
Legacy
- Optimization: linear programming, duality and shadow prices, the cutting-stock problem and the method of majorants are standard in operations research.
- Optimal transport: the Kantorovich formulation and duality underlie the work of Yann Brenier, Robert McCann, Cédric Villani and Alessio Figalli. Villani (2010) and Figalli (2018) received Fields Medals for work in this area.
- Analysis: Kantorovich spaces are a branch of functional analysis, with applications to Boolean-valued analysis and to economics with infinitely many goods.
- Numerical analysis: the Newton–Kantorovich theorem and the Kantorovich inequality are textbook results in scientific computing.
- Probability and data science: the Kantorovich–Rubinstein (Wasserstein) distance is used in limit theorems, statistics, image analysis and machine learning.
Final projectTake a transportation problem with three sources and three destinations. Solve it as a linear program, find the dual potentials, and write a page explaining how the potentials answer an economic question (what would a planner pay for one more unit at each source?) and a geometric one (Kantorovich’s 1942 criterion).
Connections
People linked to Leonid Kantorovich across this collection. A name without a link has no page here yet. See Leonid Kantorovich in the network of all 337 lives →
Teachers and advisors
- Grigorii Fichtenholz
- Vladimir Smirnov
Students
- Gennady Rubinstein 1956
- Gleb Akilov 1947
- Iosif Romanovskii 1962
- Mikhail Solomyak 1959
- Svetlozar Rachev 1979
- Valery Makarov 1965
- Victor Havin 1958
Influenced by
- Gaspard Monge
- Nikolai Luzin 1929-1930
Influenced
- Cédric Villani
Colleagues and collaborators
- John von Neumann
- Aron Pinsker
- Boris Vulikh
- George Dantzig
- Mark Gavurin
- Tjalling Koopmans
- Vasily Nemchinov
- Victor Zalgaller 1948-1951
- Viktor Novozhilov
- Vladimir Krylov 1936
Friends
- Sergei Sobolev
Also mentioned in this life
- Pavel Alexandrov
- Stefan Banach
- Augustin-Louis Cauchy
- Richard Dedekind
- Leonhard Euler
- Maurice Fréchet
- Jacques Hadamard
- Andrey Kolmogorov
- Joseph-Louis Lagrange
- Henri Lebesgue
- Isaac Newton
- Ludwig Prandtl
- Frigyes Riesz
Fields
- functional analysis 11 lives
- numerical analysis 6 lives
On CoursesHub
Sources and further reading
Biography
- L. V. Kantorovich, Biographical and Curriculum vitae, NobelPrize.org.
- J. J. O’Connor and E. F. Robertson, Leonid Vitalyevich Kantorovich, MacTutor History of Mathematics.
- Leonid Kantorovich, Wikipedia.
- INFORMS, Leonid V. Kantorovich: biographical profile.
- A. M. Vershik, L. V. Kantorovich and linear programming, arXiv:0707.0491.
- Leonid Vitalievich Kantorovich, Mathematics Genealogy Project.
Mathematics
- L. V. Kantorovich, Mathematical methods of organizing and planning production (1939), English translation in Management Science 6 (1960), 366–422.
- L. Kantorovitch, On the translocation of masses, C. R. (Doklady) Acad. Sci. URSS 37 (1942); reprinted in Management Science 5 (1958), 1–4.
- L. V. Kantorovich, Functional analysis and applied mathematics, Uspekhi Mat. Nauk 3 (1948).
- L. V. Kantorovich and G. Sh. Rubinstein, On a space of completely additive functions, Vestnik Leningrad Univ. 13 (1958).
- L. V. Kantorovich and V. I. Krylov, Approximate Methods of Higher Analysis, 1936; English translation, Interscience/Noordhoff, 1958.
- L. V. Kantorovich and G. P. Akilov, Functional Analysis, 2nd ed., Pergamon, 1982.
- L. V. Kantorovich, The Best Use of Economic Resources, Pergamon, 1965.
- L. V. Kantorovich, Mathematics in economics: achievements, difficulties, perspectives, Nobel lecture, 1975.
- C. Villani, Optimal Transport, Old and New, Springer, 2009.
- C. D. Aliprantis and O. Burkinshaw, Positive Operators, Springer, 2006.
- A. Schrijver, Theory of Linear and Integer Programming, Wiley, 1986.