Revision history [back]

"Brute-force" solution :

sage: var("y, z")
(y, z)

The right-hand size involves a constant ; compute it once...

sage: R = sum(z*binomial((1001-z),950),z,1,51) ; R
10480853106371895870377052872259391385423261386799458513700836377945732866106336546780

Brute force check of all possible (integer) solution candidates :

sage: [u for u in range(1003) if binomial(1002,u)==R]
[50, 952]

If this is homework, your professor may or may not appreciate this answer : there might be subtler points in his/her course...

HTH nevertheless...