I am trying to digest the proof for Vandermonde's identity, but due to my lack of experience with manipulating sigma notation, I am unable to understand how they got the RHS in the first step (expressing the product of two polynomials in sigma notation): $$\left(\sum\limits_{i=0}^ma_ix^i\right)\left(\sum\limits_{j=0}^nb_jx^j\right)=\sum\limits_{r=0}^{m+n}\left(\sum\limits_{k=0}^r a_kb_{r-k}\right)x^r$$
Can someone please give me a step-by-step run through for how they deduced this result, making it clear what Summation rules/identities have been utilised.
For example, why have they introduced the new index variable $r$? Why not just apply the sigma notation rule (listed on Wikipedia): $$\left(\sum\limits_{i=0}^n a_i\right)\left(\sum\limits_{j=0}^n b_j\right)=\sum\limits_{i=0}^n\sum\limits_{j=0}^na_ib_j$$ Does anyone also have any book recommendations for practising manipulating Capital Sigma notation (i.e. using the identities)? https://en.wikipedia.org/wiki/Summation