
of size \(n_{s}=\left|\Upsilon \right|\left|\mathcal {N}_{s}\right|+\left|\mathcal {E}_{s}\right|\) .
If the pmf \(p_{\eta }\)  of the number \(N_{s,\ell }\)  of users is known,
the value of \(\gamma _{s,\ell }\)  such that \(p_{s,\ell }\left(\gamma _{s,\ell }\right)=\underline{p}_{s}^{\text{sp}}\) 
can be obtained by bisection search methods, see, e.g., [1]}.
The multidimensional integral in (REF )
can be evaluated, e.g., using a quasi-Monte Carlo integration
algorithm presented in [2]}.
