[guided]We use the partition by first continued-fraction digit. For each $k\in\mathbb N$, define
\begin{align*}
I_k:=\left(\frac{1}{k+1},\frac{1}{k}\right].
\end{align*}
If $x\in I_k$, then $k\leq 1/x<k+1$, so $\lfloor 1/x\rfloor=k$. Hence on this interval the Gauss map has the explicit branch formula
\begin{align*}
G|_{I_k}(x)=\frac{1}{x}-k.
\end{align*}
This branch maps $I_k$ onto $[0,1)$ up to endpoints, and the corresponding inverse branch is
\begin{align*}
\psi_k:[0,1)\to I_k,\quad y\mapsto \frac{1}{k+y}.
\end{align*}
The relevant entropy theorem is the countable full-branch version of Rokhlin's entropy formula. It requires an invariant nonsingular probability measure, a countable generating branch partition of finite entropy, and integrability of the logarithmic derivative. We verify these requirements explicitly.
First, define the Gauss density $\rho:(0,1]\to(0,\infty)$ by
\begin{align*}
\rho(x):=\frac{1}{\log 2}\frac{1}{1+x}.
\end{align*}
To prove invariance, take a non-negative Borel function $\varphi:(0,1]\to[0,\infty]$. Changing variables on the branch $I_k$ with $x=\psi_k(y)=1/(k+y)$ and $|\psi_k'(y)|=(k+y)^{-2}$ gives
\begin{align*}
\int_{(0,1]}\varphi(G(x))\rho(x)\,d\mathcal L^1(x)=\sum_{k=1}^{\infty}\int_0^1\varphi(y)\rho(\psi_k(y))|\psi_k'(y)|\,d\mathcal L^1(y).
\end{align*}
The branch densities add up to the original density:
\begin{align*}
\sum_{k=1}^{\infty}\rho(\psi_k(y))|\psi_k'(y)|=\frac{1}{\log 2}\sum_{k=1}^{\infty}\frac{1}{(k+y)(k+y+1)}=\frac{1}{\log 2}\frac{1}{1+y}=\rho(y).
\end{align*}
The middle equality is a telescoping sum, since
\begin{align*}
\frac{1}{(k+y)(k+y+1)}=\frac{1}{k+y}-\frac{1}{k+y+1}.
\end{align*}
Hence
\begin{align*}
\int_{(0,1]}\varphi\circ G\,d\mu_G=\int_{(0,1]}\varphi\,d\mu_G,
\end{align*}
so $\mu_G$ is $G$-invariant.
Second, the partition has finite entropy. We compute
\begin{align*}
\mu_G(I_k)=\frac{1}{\log 2}\int_{1/(k+1)}^{1/k}\frac{1}{1+x}\,d\mathcal L^1(x)=\frac{1}{\log 2}\log\left(\frac{(k+1)^2}{k(k+2)}\right).
\end{align*}
Since $\log(1+t)\leq t$ for $t>-1$, this implies $\mu_G(I_k)\leq 1/(k^2\log 2)$. Therefore $-\sum_{k=1}^{\infty}\mu_G(I_k)\log\mu_G(I_k)$ converges by comparison with $\sum_{k=1}^{\infty}k^{-2}\log k$. Define the partition entropy $H_{\mu_G}(\mathcal P)$ by
\begin{align*}
H_{\mu_G}(\mathcal P):=-\sum_{k=1}^{\infty}\mu_G(I_k)\log\mu_G(I_k).
\end{align*}
The convergence just proved gives $H_{\mu_G}(\mathcal P)<\infty$.
Third, the partition is generating modulo a null set. If two points outside the countable endpoint orbit have the same itinerary with respect to the intervals $I_k$, then they have the same continued-fraction digits and hence the same continued-fraction expansion. The endpoint orbit is countable, so it has $\mu_G$-measure zero because $\mu_G$ is absolutely continuous with respect to $\mathcal L^1$.
Finally, we verify integrability of the logarithmic derivative. On every branch,
\begin{align*}
G'(x)=-\frac{1}{x^2},
\end{align*}
so $\log |G'(x)|=2\log(1/x)$. Using the density of $\mu_G$ and the bound $1/(1+x)\leq 1$ for $x\in(0,1]$, we get
\begin{align*}
\int_{(0,1]}\log |G'(x)|\,d\mu_G(x)=\frac{2}{\log 2}\int_0^1\log\left(\frac{1}{x}\right)\frac{1}{1+x}\,d\mathcal L^1(x)\leq\frac{2}{\log 2}\int_0^1\log\left(\frac{1}{x}\right)\,d\mathcal L^1(x)=\frac{2}{\log 2}.
\end{align*}
Thus $\log |G'|\in L^1(\mu_G)$. All hypotheses of the countable full-branch Rokhlin entropy formula are now verified, and the formula gives
\begin{align*}
h_{\mu_G}(G)=\int_{(0,1]} \log |G'(x)|\,d\mu_G(x).
\end{align*}
Substituting $|G'(x)|=x^{-2}$ and the Gauss density gives
\begin{align*}
h_{\mu_G}(G)
=
\frac{1}{\log 2}
\int_0^1
\log\left(\frac{1}{x^2}\right)\frac{1}{1+x}\,d\mathcal L^1(x).
\end{align*}[/guided]