I have been playing with Pollard's rho algorithm for elliptic curves over finite fields. I have noticed after some experimenting, that the algorithm almost always 'fails' for supersingular elliptic curves over a finite field with characteristic equal to a Mersenne prime. I cannot seem to find any logical explanation for it, so I wonder whether anyone has an idea why this happens.
Let me elaborate on the meaning of 'fail'. So I use the standard Pollard rho method which you find in chapter 11 in 'The Arithmetic of Elliptic Curves' by Joseph Silverman. In the implementation I used a partition function, which basically looks at the sum of the coefficients of each point on the curve modulo 3 (so that the points on the curve are just about equally divided into three sets). I don't think actually this influences the result.
Using the notation in the book, I consider the algorithm to fail if $n \mid \delta_i - \beta_i$. In this case as Silverman says, we get no information on $m$ (where we want to find $m$ such that $m \cdot x = y$, for $x,y \in E(\mathbf{F}_p)$).
Apparently, this almost always happens when $p$ is a Mersenne prime and $E(\mathbf{F}_p)$ is supersingular (so the amount of rational points is exactly $p+1$). If anyone has any suggestions why this happens, I'm all ears!