费马求 3 的逆元

运行下面这段程序：

本节模型：inv(a, p) 用费马小定理求逆元：p 是质数时 a^(p-2) mod p。

#include <array>
#include <iostream>
#include <string>
#include <vector>
using namespace std;

long long inv(long long a, long long p) {
    // 费马小定理求逆元：p 是质数时 a^(p-2) 就是 a 的逆元
    long long r = 1;
    a %= p;
    long long e = p - 2;
    while (e > 0) {
        if (e & 1) r = r * a % p;
        a = a * a % p;
        e >>= 1;
    }
    return r;
}

int main() {
    cout << inv(3, 7) << endl;
}

（本题用 g++ -std=c++17 -O0 编译。）