除法取模意义下的运算

除法取模意义下的运算

逆元

如果整数$b$, $m$互质，并且对于任意的整数$a$，如果满足$a$能被$b$整除，记为$b|a$，则存在一个整数$x$，使得$a÷b\equiv a\times x(modm)$。则称$x$为$b$模$m$乘法逆元，记为$b^{-1}（modm）$

如何求乘法逆元：

当$b$存在乘法逆元的条件是$b$与$m$互质。当模数$m$为质数的时候，$b^{m-2}$为$b$的乘法逆元

另一种定义：

$a\times x\equiv 1(modp)$，且$gcd(a,p)=1$，则称$a$关于模$p$的乘法逆元为$x$

推导过程

记$b$的逆元为$b^{-1}$，故有：
$$\frac{a}{b}\times bmodp=a\times b^{-1}\times bmodp$$

左边消去分母有：
$$amodp=a\times b^{-1}\times bmodp$$
左右两边消去$a$，有：
$$b^{-1}bmodp=1$$
费马小定理：
$$b^{p-1}modp=1$$
代入得：
$$b^{-1}b\equiv b^{p-1}modp$$
即：
$$b^{-1}=b^{p-2}modp$$
所以求$b$模$p$的逆元可以用快速幂求$b^{p-2}modp$

快速幂

通过二分的方式在$logn$的复杂度内求快速幂，由于可能溢出，所以每一次计算的时候要$modp$

int quick(int a, int b, int p)
{
    int res = 1;
    while(b)
    {
        if(b & 1) res = (res * a) % p;
        b >>= 1;
        a = (a * a) % p;
    }
    return res;
}

如果是Python的话，可以用求$x$在模$m$下的逆元。

pow(x, -1, mod)

例题

力扣2514 统计同位异构字符串数目

排列组合知识：

• 长度为$len(s)$，不包含重复字符的字符串，将所有字符全排列构成的字符串数目是$len(s)!$

• 长度为$len(s)$，包含重复字符的字符串。重复字符的个数分别为$a,b,...,n$。将所有字符全排列构成的字符串数目是$\frac{len(s)!}{a!\times b!\times ...\times n!}$

分别计算分子和分母，在对分母做除法的时候用到上面的除法下的取模运算。

代码：

k, mod = 10 ** 5, 10 ** 9 + 7
fac = [1]
for i in range(1, k+1):
    fac.append(fac[-1] * i % mod)

class Solution:
    def countAnagrams(self, s: str) -> int:
        words = s.split()
        ans = 1
        for word in words:
            cnt = Counter(word)
            ans *= fac[len(word)]
            ans %= mod
            for x in cnt.values():
                ans *= pow(fac[x], -1, mod)
                ans %= mod
        return ans

class Solution {
    const int MOD = 1e9 + 7;

    long pow(long x, int n) {
        long res = 1L;
        for (; n; n /= 2) {
            if (n % 2) res = res * x % MOD;
            x = x * x % MOD;
        }
        return res;
    }

public:
    int countAnagrams(string &s) {
        long ans = 1L, mul = 1L;
        int cnt[26]{};
        for (int i = 0, j = 0; i < s.length(); ++i) {
            if (s[i] == ' ') {
                memset(cnt, 0, sizeof(cnt));
                j = 0;
            } else {
                mul = mul * ++cnt[s[i] - 'a'] % MOD;
                ans = ans * ++j % MOD;
            }
        }
        return ans * pow(mul, MOD - 2) % MOD;
    }
};

作者：灵茶山艾府
链接：https://leetcode.cn/problems/count-anagrams/solutions/2031837/zu-he-ji-shu-by-endlesscheng-leem/
来源：力扣（LeetCode）