【C++】浮点数运算误差与逆元

问题的出现

使用double与int型数据进行相同的运算过程：

int main() {
	int mod = 1e9 + 7;
	double a = pow(13094024580916, 9);
	double b = 2;
	cout << int(fmod(a / b, mod)) << endl;
	long long r = 1;
	for (int i = 0; i < 9; i++)
	{
		r *= 13094024580916 % mod;
		r %= mod;
	}
	cout << r / 2 << endl;
	return 0;
}

会得到不同的运算结果：

算法题中常常遇到大数运算，取模作为结果的情况，这种情况采用浮点数类型计算就很容易得到错误的结果，而中间计算结果又会超出整型数据的存储范围，要怎么解决呢？
参考leetcode No1830的题解，使用乘法逆元，将问题进行转换：
两数相除+取模 → 两数相乘+取模 → 两数取模+相乘

这样就可以避免出现较大的中间结果了。

乘法逆元

定义

计算

当 mm 为质数时，一种简单的方法是使用「费马小定理」，即

那么有

至于如何快速求解 幂运算+取模：

// 快速幂，用来计算 x^y mod m
    int quickmul(int x, int y) {
        int ret = 1, mul = x;
        while (y) {
            if (y & 1) {
                ret = (LL)ret * mul % mod;
            }
            mul = (LL)mul * mul % mod;
            y >>= 1;
        }
        return ret;
    }

作者：LeetCode-Solution
链接：https://leetcode-cn.com/problems/minimum-number-of-operations-to-make-string-sorted/solution/shi-zi-fu-chuan-you-xu-de-zui-shao-cao-z-qgra/
来源：力扣（LeetCode）
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

当 m 不为质数时，我们可以使用「扩展欧几里得算法」求出乘法逆元。

应用

使用逆元，可以将除法的取模变成取模后相乘，避免大结果的出现：