1 | initial version |
There's always brute force (there is only a small finite number of possibilities for $y$):
sage: var('z')
sage: rhs = sum(z*binomial((1001-z),950),z,1,51)
sage: next(y for y in range(1, 1002) if binomial(1002, y) == rhs)
50
sage: [y for y in range(1, 1002) if binomial(1002, y) == rhs]
[50, 952]