确定性米勒-拉宾实现

我正在尝试使用确定性的米勒-拉宾算法实现素数检查函数，但结果并不总是正确的：当检查前 1,000,000 个数字时，它只找到 78,495 而不是 78,498。

这是使用 [2， 7， 61] 作为基础获得的，根据维基百科，对于高达 4,759,123,141 的值应该始终正确。
有趣的是，缺少的 3 个素数正是构成基数的素数(2、7 和 61(。

为什么会这样？我使用的代码如下：

T modular_power(T base, T exponent, T modulo) {
    base %= modulo;
    T result = 1;
    while (exponent > 0) {
        if (exponent % 2 == 1)
            result = (result * base) % modulo;
        base = (base * base) % modulo;
        exponent /= 2;
    }
    return result;
}
bool miller_rabin(const T& n, const vector<T>& witnesses) {
    unsigned int s = 0;
    T d = n - 1;
    while (d % 2 == 0) {
        s++;
        d /= 2;
    }
    for (const auto& a : witnesses) {
        if (modular_power<T>(a, d, n) == 1)
            continue;
        bool composite = true;
        for (unsigned int r = 0; r < s; r++) {
            if (modular_power<T>(a, (T) pow(2, r) * d, n) == n - 1) {
                composite = false;
                break;
            }
        }
        if (composite)
            return false;
    }
    return true;
}
bool is_prime(const T& n) {
    if (n < 4759123141)
        return miller_rabin(n, {2, 7, 61});
    return false; // will use different base
}