sequences and series - Generalized Euler sum $sum_{n=1}^infty frac{H_n}{n^q}$

I found the following formula

$$\sum_{n=1}^\infty \frac{H_n}{n^q}= \left(1+\frac{q}{2} \right)\zeta(q+1)-\frac{1}{2}\sum_{k=1}^{q-2}\zeta(k+1)\zeta(q-k)$$

and it is cited that Euler proved the formula above , but how ?

Do there exist other proofs ?

Can we have a general formula for the alternating form

$$\sum_{n=1}^\infty (-1)^{n+1}\frac{H_n}{n^q}$$

Answer

$$
\begin{align}
&\sum_{j=0}^k\zeta(k+2-j)\zeta(j+2)\\

&=\sum_{m=1}^\infty\sum_{n=1}^\infty\sum_{j=0}^k\frac1{m^{k+2-j}n^{j+2}}\tag{1}\\
&=(k+1)\zeta(k+4)
+\sum_{\substack{m,n=1\\m\ne n}}^\infty\frac1{m^2n^2}
\frac{\frac1{m^{k+1}}-\frac1{n^{k+1}}}{\frac1m-\frac1n}\tag{2}\\
&=(k+1)\zeta(k+4)
+\sum_{\substack{m,n=1\\m\ne n}}^\infty\frac1{nm^{k+2}(n-m)}-\frac1{mn^{k+2}(n-m)}\tag{3}\\
&=(k+1)\zeta(k+4)
+2\sum_{m=1}^\infty\sum_{n=m+1}^\infty\frac1{nm^{k+2}(n-m)}-\frac1{mn^{k+2}(n-m)}\tag{4}\\
&=(k+1)\zeta(k+4)
+2\sum_{m=1}^\infty\sum_{n=1}^\infty\frac1{(n+m)m^{k+2}n}-\frac1{m(n+m)^{k+2}n}\tag{5}\\

&=(k+1)\zeta(k+4)\\
&+2\sum_{m=1}^\infty\sum_{n=1}^\infty\frac1{m^{k+3}n}-\frac1{(m+n)m^{k+3}}\\
&-2\sum_{m=1}^\infty\sum_{n=1}^\infty\frac1{m(n+m)^{k+3}}+\frac1{n(n+m)^{k+3}}\tag{6}\\
&=(k+1)\zeta(k+4)
+2\sum_{m=1}^\infty\frac{H_m}{m^{k+3}}
-4\sum_{n=1}^\infty\sum_{m=1}^\infty\frac1{n(n+m)^{k+3}}\tag{7}\\
&=(k+1)\zeta(k+4)
+2\sum_{m=1}^\infty\frac{H_m}{m^{k+3}}
-4\sum_{n=1}^\infty\sum_{m=n+1}^\infty\frac1{nm^{k+3}}\tag{8}\\
&=(k+1)\zeta(k+4)

+2\sum_{m=1}^\infty\frac{H_m}{m^{k+3}}
-4\sum_{n=1}^\infty\sum_{m=n}^\infty\frac1{nm^{k+3}}+4\zeta(k+4)\tag{9}\\
&=(k+5)\zeta(k+4)
+2\sum_{m=1}^\infty\frac{H_m}{m^{k+3}}
-4\sum_{m=1}^\infty\sum_{n=1}^m\frac1{nm^{k+3}}\tag{10}\\
&=(k+5)\zeta(k+4)
+2\sum_{m=1}^\infty\frac{H_m}{m^{k+3}}
-4\sum_{m=1}^\infty\frac{H_m}{m^{k+3}}\tag{11}\\
&=(k+5)\zeta(k+4)
-2\sum_{m=1}^\infty\frac{H_m}{m^{k+3}}\tag{12}

\end{align}
$$
Letting $q=k+3$ and reindexing $j\mapsto j-1$ yields
$$
\sum_{j=1}^{q-2}\zeta(q-j)\zeta(j+1)
=(q+2)\zeta(q+1)-2\sum_{m=1}^\infty\frac{H_m}{m^q}\tag{13}
$$
and finally
$$
\sum_{m=1}^\infty\frac{H_m}{m^q}

=\frac{q+2}{2}\zeta(q+1)-\frac12\sum_{j=1}^{q-2}\zeta(q-j)\zeta(j+1)\tag{14}
$$

---

Explanation

$\hphantom{0}(1)$ expand $\zeta$
$\hphantom{0}(2)$ pull out the terms for $m=n$ and use the formula for finite geometric sums on the rest
$\hphantom{0}(3)$ simplify terms
$\hphantom{0}(4)$ utilize the symmetry of $\frac1{nm^{k+2}(n-m)}+\frac1{mn^{k+2}(m-n)}$
$\hphantom{0}(5)$ $n\mapsto n+m$ and change the order of summation
$\hphantom{0}(6)$ $\frac1{mn}=\frac1{m(m+n)}+\frac1{n(m+n)}$
$\hphantom{0}(7)$ $H_m=\sum_{n=1}^\infty\frac1n-\frac1{n+m}$ and use the symmetry of $\frac1{m(n+m)^{k+3}}+\frac1{n(n+m)^{k+3}}$
$\hphantom{0}(8)$ $m\mapsto m-n$
$\hphantom{0}(9)$ subtract and add the terms for $m=n$
$(10)$ combine $\zeta(k+4)$ and change the order of summation
$(11)$ $H_m=\sum_{n=1}^m\frac1n$
$(12)$ combine sums