Computation of Hecke eigenvalues (mod $p$) via quaternions

In a 1987 letter, Serre proves that the systems of Hecke eigenvalues arising from mod $p$ modular forms (of fixed level $\mathrm{\Gamma }(N)$ coprime to $p$, and any weight $k$) are the same as those arising from functions $\mathrm{\Omega }(N)\to {\bar{\mathbb{F}}}_p$, where $\mathrm{\Omega }(N)$ is some double quotient of $D^{\times }({\mathbb{A}}_f)$ and $D$ is the unique quaternion algebra over $\mathbb{Q}$ ramified at $\{p,\mathrm{\infty }\}$. We present an algorithm which then computes these Hecke eigenvalues on the quaternion side in a combinatorial manner.