I'm trying to solve the following question:
If $s_n$ is the sum of first $n$ natural numbers, then prove that
$$2(s_1s_{2n}+s_2s_{2n-1}+\dots+s_ns_{n+1})=\frac{(2n+4)!}{5!(2n-1)!}$$
This is where I've come so far:
The general term of $(1-x)^{-2n}$ happens to be $t_{r+1}=\frac{(2n+r-1)!}{(2n-1)!r!}x^r$ and therefore the 6th term, i.e, $t_{5+1}=t_6=\frac{(2n+4)!}{(2n-1)!5!}x^5$, which looks like the RHS of the question I'm trying to solve.
Second, I found that $$(1-x)^{-3}=s_1+s_2x+\dots+s_nx^{n-1}+\dots$$ It's clear that the LHS must be something like, $(1-x)^{-a}(1-x)^{-b}$. But I'm not able to find one such combination which would give me $2(s_1s_{2n}+s_2s_{2n-1}+\dots+s_ns_{n+1})$.
Help is greatly appreciated.
More generally, $$\sum_{k=3}^{m-2}\binom{k-1}{2}\binom{m-k}{2}=\binom{m}{5},$$ which you can prove combinatorially by counting $5$-subsets of $\{1,\dots,m\}$. The RHS is clear. For the LHS, condition on the third largest (median) element $k$ chosen. Then $\binom{k-1}{2}$ counts the selections of $2$ elements that are less than $k$, and $\binom{m-k}{2}$ counts the selections of $2$ elements that are greater than $k$.
For the OP's identity, rewrite $$2\sum_{k=1}^{n} s_k s_{2n+1-k} =\sum_{k=1}^{2n} s_k s_{2n+1-k} =\sum_{k=1}^{2n} \binom{k+1}{2} \binom{2n+2-k}{2} =\sum_{k=3}^{2n+2} \binom{k-1}{2} \binom{2n+4-k}{2} $$ and take $m=2n+4$.