 Solution

Work first with essential bounds. On a Euclidean ball $B_R(x_0)⋐B_1$, put $M={\mathrm{esssup}}_{B_R}u$, $m={\mathrm{essinf}}_{B_R}u$ and $\omega =M-m$. The functions $M-u$ and $u-m$ are nonnegative weak solutions, hence weak supersolutions. In $B_{R/2}$, at least one of the sets $\{u\le (M+m)/2\}$ and $\{u\ge (M+m)/2\}$ has at least half the measure.

In the first case, the Weak Harnack inequality applied to $M-u$ givesIn the second case, applying it to $u-m$ gives the same lower bound for the improvement of the minimum. Thus the oscillation decay estimate isThis measure argument works even when $q<1$, avoiding a false use of the triangle inequality in that range.

Iteration yields ${\mathrm{osc}}_{B_r(x_0)}u\le C(r/R{)}^{\mu }{\mathrm{osc}}_{B_R(x_0)}u$, with, for example, $\mu =\min \{1/2,-\log \tau /\log 2\}\in (0,1)$. At Lebesgue points the shrinking oscillations define a unique continuous representative, and the estimate extends to every point by continuity. For $x,y\in B_{1/4}$ at distance less than $1/4$, apply the estimate to a Euclidean ball centred at $x$ with initial radius $1/2$ and final radius comparable to $|x-y|$. This Euclidean ball stays in $B_1$. Larger distances use the trivial bound $2\Vert u{\Vert }_{\infty }$. We obtainThe constants depend only on $n,\lambda ,\Lambda$. This derives the relevant conclusion of the De Giorgi-Nash-Moser theorem directly from part (a).