By applying the identity
X+YXY−Y2=Y−X+Y2Y2
with X=ai and Y=bi to the terms in the LHS and
X=i=1∑nai and
Y=i=1∑nbi to the RHS,
LHS=i=1∑nai+biaibi−bi2=i=1∑n(bi−ai+bi2bi2)=i=1∑nbi−2i=1∑nai+bibi2,
RHS=i=1∑nai+i=1∑nbii=1∑nai⋅i=1∑nbi−(i=1∑nbi)2=i=1∑nbi−2i=1∑n(ai+bi)(i=1∑nbi)2.
Therefore, the statement is equivalent with
i=1∑nai+bibi2≥i=1∑n(ai+bi)(i=1∑nbi)2,
which is the same as the well-known variant of the Cauchy-Schwarz
inequality,
i=1∑nYiXi2≥Y1+⋯+Yn(X1+⋯+Xn)2(Y1,…,Yn>0)
with Xi=bi and Yi=ai+bi.