An extension of a congruence by Tauraso (Q1952711)

Summary: For a positive integer \(n\) let \(H_n = \sum^n_{k=1} 1/n\) be the \(n\)th harmonic number. In this paper we prove that, for any prime \(p \geq 7\), \[ \sum^{p-1}_{k=1} \frac{H_k}{k^2} \equiv \sum^{p-1}_{k=1}\frac{H^2_k}{k} \equiv \frac 3{2p} \sum^{p-1}_{k=1} \frac 1{k^2} \pmod{p^2}. \] Notice that the first part of this congruence is proposed in 2008 by \textit{R. Tauraso} [J. Number Theory 130, No. 12, 2639--2649 (2010; Zbl 1208.11027)]. In our elementary proof of the second part of the above congruence we use certain classical congruences modulo a prime and the square of a prime, some congruences involving harmonic numbers, and a combinatorial identity due to \textit{V. Hernández} [``Solution IV of problem 10490 (a reciprocal summation identity),'' Am. Math. Mon. 106, 589--590 (1999)]. Our auxiliary results contain many interesting combinatorial congruences involving harmonic numbers.