Let 𝑝 be an odd prime and $ℎ(𝑥)=𝑥^2+𝑎𝑥+1$ be an irreducible polynomial over the field $\mathbb{Z}_p$. The polynomial function

$$\Psi:\mathbb{Z}_p^2⟶\mathbb{Z}_p,\quad (𝑥,𝑦)\mapsto 𝑥^2+𝑦^2−𝑥+𝑦−𝑎𝑥𝑦$$

is surjective, as proved here: https://mathoverflow.net/questions/356935/image-of-a-polynomial-function-x2y2-xy-axy-over-mathbbf-p.

I would like to compute a set of representatives of the classes of the kernel of Ψ (i.e. the relation $𝑘𝑒𝑟(\Psi)={((𝑥,𝑦),(𝑡,𝑤))\in \mathbb{Z}_p^4, \, \Psi(𝑥,𝑦)=\Psi(𝑡,𝑤)})$. So basically I would like to provide a set of $𝑝−1$ elements of $\mathbb{ℤ}_p^2$ that have all different values when you applies $\Psi$ (I am not interested to the case $\Psi(𝑥,𝑦)=0$).

Is there a way to do that in general, regardless who is $p$?