def binomcoeff(n,k):
nom = 1
denum = 1
for i in range(n-k):
nom *= (n-i)
denum *= (i+1)
v = nom//denum
return v
def E(n,k):
return binomcoeff(2*n+k,n)*binomcoeff(n+k,n)//2
print(E(1,2), E(2,1), E(2,2), E(3,1))