Mathematics & Statistics

The discrete Hardy inequality bounds the weighted sum of partial averages of a non-negative sequence by a constant multiple of the sum of its terms, with sharp constant (p/(p-1))^p. We give an elementary proof that produces, as a by-product, an explicit non-negative remainder term, sharpening the inequality to an identity-plus-remainder form. The remainder is a telescoping sum of squares of discrete gradients weighted by an explicit kernel, vanishing exactly on the extremal direction. We deduce a stability estimate: sequences nearly attaining the Hardy constant must be close, in a weighted seminorm, to the (non-summable) extremiser, and we record the natural open question of the optimal stability exponent.