luogu P3811 模板乘法逆元

Posted 2021-02-05 sevenyuanluo

solution 1：

费马小定理

若p为素数，a为正整数，且a、p互质，则有a^（p - 1）≡ 1（mod p）

那么a的逆元就是a^(p - 2)

用一个快速幂即可

//t两个点

#include<cstdio>
using namespace std;
#define int long long

int p;

int quickpow(int n,int k) {
    int ans = 1;
    while(k) {
        if(k & 1)
            ans = ans * n % p;
        n = n * n % p;
        k >>= 1;
    }
    return ans % p;
}

 main() {
    int n;
    scanf("%lld%lld",&n,&p);
    for(int i = 1; i <= n; i++)
        printf("%lld
",quickpow(i,p - 2));
    return 0;
}

solution 2：

exgcd

a * x = 1（mod p）

则 a * x + p * y = 1

//t一个点

#include<cstdio>
using namespace std;

int x,y;

void exgcd(int a,int b) {
    if(b == 0) {
        x = 1;
        y = 0;
        return ;
    }
    exgcd(b,a % b);
    int tmp = x;
    x = y;
    y = tmp - a / b * y;
}

int main() {
    int n,p;
    scanf("%d%d",&n,&p);
    for(int i = 1; i <= n; i++) {
        exgcd(i,p);
        printf("%d
",(x % p + p)% p);
    }
    return 0;
}

solution 3：

线性递推

根据费马小定理的公式进行推导得出

#include<cstdio>
using namespace std;
#define maxn 3000010
#define ll long long

ll inv[maxn];

int main() {
    int n,p;
    scanf("%d%d",&n,&p);
    inv[1] = 1;
    printf("1
");
    for(int i = 2; i <= n; i++) {
        inv[i] = (ll)(p - p / i) * inv[p % i] % p;
        printf("%lld
",inv[i]);
    }
    return 0;
}