Respuesta parcial: el cálculo al final resulta ser menos trivial de lo que esperaba, pero sí da una forma equivalente en términos de congruencia de los tres números de Fibonacci mod $p^2$.
Idea: $F_n$ es una combinación lineal de las progresiones geométricas, y por un clásico truco (la escritura como una integral y hacer un cambio de variables) tenemos:
$$\sum_{k=1}^n\frac{x^k}k = \sum_{k=1}^n\binom nk\frac{(x-1)^k-(-1)^k}k$$
Tomando $n=p$ el coeficiente binomial casi desaparece y podemos obtener
$$H(x)=\sum_{k=1}^{p-1}\frac{x^k}k \equiv\frac{(x-1)^p+1-x^p}p\pmod p$$
(para todos los $x\in \mathbb Z$, o como polinomios en $x$) $p$ impar.
Para filtrar par o impar de términos, tome $H(x)\pm H(-x)$.
De vuelta al problema: Tenemos $F_n = C \cdot (\phi^n-(\phi')^{n})$ para algunas constantes $C=1/(\phi-\phi')$. Aquí $\phi'=-\phi^{-1}$.
Algunos detalles para hacer todo legal: para asegurarnos de no dividir por $0$ y todas las congruencias sentido, a partir de ahora trabajamos en el ring $A=\mathbb Z_{(p)}[(\sqrt 5)^{-1}]\subset\mathbb Q(\sqrt5)$ donde $1,2, \ldots,p-1,\sqrt5=1/C,\phi$ es invertible. La elección de un ideal maximal $\mathfrak m$ contiene $p$ da un morfismos de $\mathbb Z$-álgebras de a $\mathbb F_p$ o $\mathbb F_{p^2}$, que envía los números de Fibonacci en $\mathbb Z$ a los números de Fibonacci mod $p$. Un entero es en $\mathfrak m^k$ fib es divisible por $p^k$, ya que el $p\neq5$ es unramified.
La aplicación de la excesiva $A$, sólo tenemos que calcular
$$\phi (H(\phi)-H(-\phi))\\ - \phi'( H(\phi')- H(-\phi'))$$
resp.
$$\phi' (H(\phi)-H(-\phi))\\ - \phi( H(\phi')- H(-\phi'))$$
modulo $p$, por lo que ahora tenemos formas cerradas. Equivalentemente, podemos deshacernos de el denominador $p$ y demostrar que $p$ veces que es $0$ mod $p^2$.
En realidad no es obvio que es $0\pmod {p^2}$. Para el primero, escrito todo en términos de $\pm$ poderes de $\phi$ (con $-\phi'-1=-\phi^2$, etc), obtenemos
$$[-p+1]-[p+1]+[2p+1]-[p+1]\\
-[p-1]+[p-1]+[-2p-1]+[p-1]$$
donde $[n]$ denota $\phi^n$. El uso de $[n]=\phi F_n+F_{n-1}$ $F_{-n}=(-1)^nF_n$ este es:
$$(2\phi-1)(F_{2p+1}-2F_{p+1}-F_{p-1})$$
por lo que es equivalente a mostrar
$$\bbox[5px,border:2px solid red]{F_{2p+1}-2F_{p+1}-F_{p-1} \equiv0\pmod{p^2}}\qquad p\equiv\pm1\pmod5$$
Asimismo, para el segundo:
$$[-p-1]-[p-1]+[2p-1]-[p-1]\\
- (p+1]+[-p+1]+[-2p+1]+[-p+1]$$
que le da el equivalente (veces $2\phi-1$)
$$\bbox[5px,border:2px solid red]{F_{2p-1}-2F_{p-1}-F_{p+1} \equiv0\pmod{p^2}}\qquad p\equiv\pm2\pmod5$$