Compute the value of the following expression
$$\left (1+\frac{1}{2}+\cdots+\frac{1}{n} \right )^2+\left ( \frac{1}{2}+\cdots + \frac{1}{n}\right )^2+\cdots+\left (\frac{1}{n-1}+\frac{1}{n} \right )^2+\left (\frac{1}{n} \right )^2$$
The answer is $\boxed{2n-\left (1+\frac{1}{2}+\cdots+\frac{1}{n} \right )}$.
I've been trying to do it but I've been failed. Any ideas ?
Using Iverson brackets and omitting the limits of sums when their index runs from $1$ to $n$, this is $$ S_n=\sum_k\left(\sum_i\frac1i[k\leqslant i]\right)^2=\sum_{k}\sum_{i,j}\frac1{ij}[k\leqslant i,k\leqslant j]=\sum_{i,j}\frac1{ij}\sum_{k}[k\leqslant\min(i,j)] $$ hence $$ S_n=\sum_{i,j}\frac1{ij}\min(i,j)=\sum_{i,j}\frac1{\max(i,j)}=\sum_{m}\frac1m\sum_{i,j}[\max(i,j)=m] $$ that is, as claimed by the OP, $$ S_n=\sum_m\frac1m(2m-1)=2n-\sum_m\frac1m=2n-H_n $$