Polarized superspecial simple abelian surfaces with real Weil numbers

Let $q$ be an odd power of a prime $p\in \mathbb{N}$, and $\mathrm{PPSP}(\sqrt{q})$ be the finite set of isomorphism classes of principally polarized superspecial abelian surfaces in the simple isogeny class over $\mathbb{F}_q$ corresponding to the real Weil $q$-numbers $\pm \sqrt{q}$. We produce explicit formulas for $\mathrm{PPSP}(\sqrt{q})$ of the following kinds: (i) the class number formula, i.e.~the cardinality of $\mathrm{PPSP}(\sqrt{q})$; (ii) the type number formula, i.e. the number of endomorphism rings up to isomorphism of the underlying abelian surfaces of $\mathrm{PPSP}(\sqrt{q})$. Similar formulas are obtained for other collections of polarized superspecial members of this isogeny class grouped together according to their polarization modules. We observe several surprising identities involving the arithmetic genus of certain Hilbert modular surface on one side and the class number or type number of $(P, P_+)$-polarized superspecial abelian surfaces in this isogeny class on the other side.