ll calc_binom(ll n, ll k, ll p) { if (n >= p || k > n) return 0; ll x = k % 2 != 0 ? p-1 : 1; for (ll c = p-1; c > n; c--) { x *= c - k; x %= p; x *= multInv(c, p); x %= p; } return x; }