lst=(expand(x^1100*(1-x^851)^4*(1-x^1501))).coefficients()
f(x)=0
coins=3000
for l in lst:
    if l[1] <= coins:
        f += l[0]*binomial(x-l[1]+4,4)
print f(coins)