Deje $p$ ser un extraño prime, y denotan $$f(x)=\sum_{k=0}^{p-1}\binom{2k}{k}^2x^k.$$ Demostrar que para cada $x\in \mathbf Z$,$$(-1)^\frac{p-1}2f(x)\equiv f\left(\frac{1}{16}-x\right)\pmod{p^2}.$$ Este es un concurso que pregunta, no sé cómo demostrarlo. Gracias.
Además: me parece que es equivalente a probar que: $$\sum _{k=r}^{p-1} \frac{(-1)^r \binom{2 k}{k}^2 \binom{k}{r}}{16^{k-r}}\equiv(-1)^{\frac{p-1}{2}} \binom{2 r}{r}^2 \pmod {p^2}\tag1$$ para $r=0,1,2,\cdots p-1.$
Y $(1)$ es equivalente a $$\sum _{k=r}^{p-1} \binom{2 k}{k}^2 \binom{k}{r}16^{-k} \equiv(-1)^{\frac{p-1}{2}} (-16)^{-r}\binom{2 r}{r}^2 \pmod {p^2}.$$