Czy istnieje sposób na zbudowanie np.(853467 * 21660421200929) % 100000000000007 bez bibliotek BigInteger (zauważ, że każda liczba pasuje do 64-bitowej liczby całkowitej, ale wynik mnożenia nie)?

To rozwiązanie wydaje się nieefektywne:

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;
}