use scipy nchoosek function