Abstract
The aim of this work is to optimize the existing formula based on Wilson's theorem to reduce the magnitude of the computation results.
Wilson's theorem states: if p is a prime number, then (p-1)! + 1 is divisible by p (p-1)! ≡ -1 (mod p).
The function (p-1)! increases very rapidly and reaches huge values.
When the values of p are large, the calculations become resource-intensive, so it is necessary to reduce the upper limit of the calculation results.



![Author ORCID: We display the ORCID iD icon alongside authors names on our website to acknowledge that the ORCiD has been authenticated when entered by the user. To view the users ORCiD record click the icon. [opens in a new tab]](https://www.cambridge.org/engage/assets/public/coe/logo/orcid.png)