luogu P3811 【模板】乘法逆元

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\n",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\n",(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\n");
    for(int i = 2; i <= n; i++) {
        inv[i] = (ll)(p - p / i) * inv[p % i] % p;
        printf("%lld\n",inv[i]);
    }
    return 0;
}

原文地址:https://www.cnblogs.com/sevenyuanluo/p/10307363.html

时间: 2024-08-30 15:06:10

luogu P3811 【模板】乘法逆元的相关文章

Luogu P3811 [模板]乘法逆元 题解报告

题目传送门 [题目大意] 给定$n$,求$1-n$在膜$p$意义下的乘法逆元. [思路分析] 好的原本我只会求单个数的逆元,然后被告知了这道题之后发现自己不会做(我果然还是太弱了),于是就学了一下递推求逆元. 设$p=k*i+r$,则可得$k*i+r\equiv0(mod\ p)$,然后乘上$i^{-1},r^{-1}$即可得到$k*r^{-1}+i^{-1}\equiv0(mod\ p)$ 由于$k=\lfloor \frac{p}{i}\rfloor,r=p\ mod\ i$,所以$i^{-

[模板]乘法逆元

本博客所有代码基于题目 luogu_P3811 逆元: 一般用于求 (a/b) mod p  定义: 若 a*x ≡ 1 (mod p) ,且 a 与 p 互质,那么我们就能定义: x 为 a 的逆元,记为 a^-1 ,所以我们也可以称 x 为 a 的倒数(mod p意义下). 所以对于 (a/b) mod p ,我们就可以求出 b 在 mod p 意义下的逆元,然后乘上 a ,再 mod p ,就是这个乘法逆元的值了. 求法: First:费马小定理 定理内容:如果 a , p 互质,那么 a

luogu P3811线性求逆元

首先扩O:T了一个点(因为上界松),83分. #include <cstdio> using namespace std; int n, p; void exgcd(int a, int p, int &b, int &x){ if (p==0){ b=1, x=0; return; } exgcd(p, a%p, b, x); int tmp=b; b=x; x=tmp-a/p*x; return; } int main(){ scanf("%d%d",

P3811 【模板】乘法逆元

P3811 [模板]乘法逆元 题目背景 这是一道模板题 题目描述 给定n,p求1~n中所有整数在模p意义下的乘法逆元. 输入输出格式 输入格式: 一行n,p 输出格式: n行,第i行表示i在模p意义下的逆元. 输入输出样例 输入样例#1: 10 13 输出样例#1: 1 7 9 10 8 11 2 5 3 4 说明 1≤n≤3×10?6??,n<p<20000528 输入保证 p 为质数. 我们有三种办法求逆元 由欧拉定理可知 当gcd(a,n)==1 时 我们有 Aφ(n-1)≡ 1(mod

洛谷 P3811 【模板】乘法逆元 如题

P3811 [模板]乘法逆元 时空限制1s / 256MB 题目背景 这是一道模板题 题目描述 给定n,p求1~n中所有整数在模p意义下的乘法逆元. 输入输出格式 输入格式: 一行n,p 输出格式: n行,第i行表示i在模p意义下的逆元. 输入输出样例 输入样例#1: 10 13 输出样例#1: 1 7 9 10 8 11 2 5 3 4 说明 1 \leq n \leq 3 \times 10 ^ 6, n < p < 200005281≤n≤3×106,n<p<20000528

题解 P3811 【模板】乘法逆元

题意求\(i\)在模\(p\)意义下的逆元\(\frac{1}{i}\)即\(inv(i)\).题目数据范围很明显规定了要求一个线性求逆元的算法. 令\(p=ai+b\),则有: \[ai+b\equiv 0(\mod p)\] 移项得: \[ai\equiv -b(\mod p)\] 系数化简得: \[i\equiv -\frac{b}{a}(\mod p)\] 取倒数得: \[\frac{1}{i} \equiv -\frac{a}{b}(\mod p)\] 即 \[inv(i) \equi

[洛谷P3811]【模板】乘法逆元

题目大意:给你n和质数p,求1~n在模p意义下的乘法逆元(n<p). 解题思路:由于$n<p<20000528$,所以扩展欧几里得是会超时的.这儿就要用到线性推逆元大法辣!→不懂戳这里← 注意乘法可能会超过int,所以计算时先转化为long long即可. C++ Code: #include<cstdio> int n,p,inv[20000529]; int main(){ scanf("%d%d",&n,&p); puts("

【模板】求1~n的整数的乘法逆元

洛谷3811 先用n!p-2求出n!的乘法逆元 因为有(i-1)!-1=i!-1*i (mod p),于是我们可以O(n)求出i!-1 再用i!-1*(i-1)!=i-1 (mod p)即是答案 1 #include<iostream> 2 #include<cstring> 3 #include<cstdlib> 4 #include<cstdio> 5 #include<cmath> 6 #include<algorithm> 7

T106021 【模板】乘法逆元(快速幂)

题目地址 注意点: 使用exgcd求乘法逆元需额外进行(相对)较多操作,建议使用快速幂求乘法逆元. #include<cstdio> #include<iostream> #define ll long long using namespace std; int n,p; int poww(int a,int b){ ll ans=1,tmp=a; while(b){ if(b&1){ ans*=tmp; ans%=p; } tmp=tmp*tmp; tmp%=p; b&g

乘法逆元(模板)

乘法逆元定义: b存在乘法逆元的充要条件是b与模数m互质 原因:b * x ≡ 1 (mod m) 如果b和m不互质,则 b * x肯定是m的倍数,b * x%m=0  所以b%m==0 ,b不存在乘法逆元 1.当n为质数时,可以用快速幂求逆元: a / b(整除) ≡ a * x (mod m) 两边同乘b可得 a ≡ a * b * x (mod m) -> 1 ≡ b * x (mod m) 同 b * x ≡ 1 (mod m) 由费马小定理可知,当m为质数时:  b ^ (m - 1)