[step:Construct a Borel set $B \supset A$ with $\mathcal{H}^s(B) = \mathcal{H}^s(A)$ (Borel regularity)]
Let $A \subset \mathbb{R}^n$. If $\mathcal{H}^s(A) = +\infty$, take $B = \mathbb{R}^n$, which is Borel and satisfies $\mathcal{H}^s(B) \geq \mathcal{H}^s(A) = +\infty$.
Assume $\mathcal{H}^s(A) < +\infty$. For each $k \geq 1$, set $\delta_k = 1/k$. For each pair $(k, m)$ with $k, m \geq 1$, choose a countable $\delta_k$-cover $\{C_{k,m,j}\}_{j=1}^\infty$ of $A$ such that
\begin{align*}
\sum_{j=1}^\infty \alpha(s)\left(\frac{\operatorname{diam}(C_{k,m,j})}{2}\right)^s \leq \mathcal{H}^s_{\delta_k}(A) + \frac{1}{m}.
\end{align*}
Define the open set
\begin{align*}
U_{k,m} = \bigcup_{j=1}^\infty \left\{x \in \mathbb{R}^n : \operatorname{dist}(x, C_{k,m,j}) < \frac{1}{m}\right\}.
\end{align*}
Each set in the union is open (as an open $1/m$-neighborhood of $C_{k,m,j}$), so $U_{k,m}$ is open, and $A \subset \bigcup_j C_{k,m,j} \subset U_{k,m}$.
Define the $G_\delta$ set
\begin{align*}
B = \bigcap_{k=1}^\infty \bigcap_{m=1}^\infty U_{k,m}.
\end{align*}
Then $B$ is a countable intersection of open sets, hence a Borel set ($G_\delta$), and $A \subset B$.
We claim $\mathcal{H}^s(B) = \mathcal{H}^s(A)$. By monotonicity, $\mathcal{H}^s(B) \geq \mathcal{H}^s(A)$. For the reverse, fix $k \geq 1$. For every $m \geq 1$, the set $B \subset U_{k,m}$ is covered by the open $1/m$-neighborhoods of the sets $C_{k,m,j}$. Each neighborhood satisfies
\begin{align*}
\operatorname{diam}\!\left(\{x : \operatorname{dist}(x, C_{k,m,j}) < 1/m\}\right) \leq \operatorname{diam}(C_{k,m,j}) + 2/m.
\end{align*}
For any $\delta > 0$, choosing $m$ large enough so that $\delta_k + 2/m \leq \delta + \delta_k$ and $1/m < \delta/2$, these neighborhoods form a cover of $B$ by sets of diameter at most $\delta_k + 2/m$. As $m \to \infty$, the diameters of the covering sets approach those of the original $C_{k,m,j}$, and the covering sums approach $\mathcal{H}^s_{\delta_k}(A)$. Specifically, for each fixed $k$:
\begin{align*}
\mathcal{H}^s_{\delta_k + 2/m}(B) \leq \sum_{j=1}^\infty \alpha(s)\left(\frac{\operatorname{diam}(C_{k,m,j}) + 2/m}{2}\right)^s.
\end{align*}
As $m \to \infty$, the right-hand side approaches $\mathcal{H}^s_{\delta_k}(A)$, because $\operatorname{diam}(C_{k,m,j}) \leq \delta_k$ and the extra $2/m$ vanishes. Therefore $\mathcal{H}^s_{\delta_k}(B) \leq \mathcal{H}^s_{\delta_k}(A)$ (since $\delta_k + 2/m \to \delta_k$). Taking $k \to \infty$ (equivalently $\delta_k \to 0$) gives $\mathcal{H}^s(B) \leq \mathcal{H}^s(A)$.
Combining, $\mathcal{H}^s(B) = \mathcal{H}^s(A)$.
[/step]