Existe uma maneira de construir, e.(853467 * 21660421200929) % 100000000000007 sem BigInteger bibliotecas (note que cada número se encaixa em um inteiro de 64 bits, mas o resultado da multiplicação não)?

Esta solução parece ineficiente:

int64_t mulmod(int64_t a, int64_t b, int64_t m) {
    if (b < a)
        std::swap(a, b);
    int64_t res = 0;
    for (int64_t i = 0; i < a; i++) {
        res += b;
        res %= m;
    }
    return res;
}