！！逆元详解！！————Luogu P3811 【模板】乘法逆元【逆元】（未完成）

最新推荐文章于 2025-01-09 14:56:16 发布

原创最新推荐文章于 2025-01-09 14:56:16 发布 · 447 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#逆元 #模板 #数论 #Luogu #题解

题解（较高质量）同时被 2 个专栏收录

60 篇文章

订阅专栏

模板题

15 篇文章

订阅专栏

这篇博客详细介绍了模运算中的逆元概念，包括逆元的定义，并提供了三种求解逆元的方法：扩展欧几里得算法、费马小定理和线性求逆元。博主特别指出在Luogu P3811题目中，线性求逆元是唯一能在时限内通过的解法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目

给定 $n, p$ 求 $1\sim n$ 中所有整数在模 $p$ 意义下的乘法逆元。

思路

定义

首先看逆元的定义：
若 $a*x\equiv1 \pmod {b}$ ，且 $a$ 与 $b$ 互质，
那么我们就能定义: $x$ 为 $a$ 的逆元，记为 $a^{-1}$ ，
所以我们也可以称 $x$ 为 $a$ 在 $\bmod b$ 意义下的倒数。

逆元有很多应用，比如在模意义下的乘法中。
对于 $\displaystyle\frac{a}{b} \pmod {p}$ ，我们就可以求出 $b$ 在 $\bmod {p}$ 下的逆元，然后乘上 $a$ ，再 $\bmod {p}$ ，就是这个分数的值了。

那怎么求逆元呢？
求逆元呢，有三种做法,我会一一解释。

求法

1.exgcd求逆元

……

exgcd求逆元代码

//已卡常，已开O2
#include<iostream>
#include<cstdio>
using namespace std;
int n,p,x,y;
inline void exgcd(int a,int b,int &x,int &y)
{
	if(b==0)
	 {
	 	x=1;
	 	y=0;
	 	return;
	 }
	exgcd(b,a%b,x,y);
    int t=x;
	x=y;
	y=t-(a/b)*y;
}
inline void print(int x)
{
    if(x<0)
     {
    	putchar('-');
    	x=-x;
     }
    if(x>9)
      print(x/10);
    putchar(x%10+'0');
}
int main()
{
	scanf("%d%d",&n,&p);
	for(register int i=1; i<=n; ++i)
	 {
	 	exgcd(i,p,x,y);
	 	x=(x%p+p)%p;
	 	print(x);
	 	puts("");
	 }
	return 0;
}

2.费马小定理求逆元

……

3.线性求逆元

线性求逆元呢，其实不需要用到什么算法或定理，直接推式子就可以了。

本题只有线性求逆元能AC，其余两种卡常也只能跑到550ms。（时限500ms）。

线性求逆元代码

#include<iostream>
#include<cstdio>
using namespace std;
long long n,p,inv[30000010];
int main()
{
   scanf("%lld%lld",&n,&p);
   printf("1\n");
   inv[1]=1;
   for(long long i=2; i<=n; i++)
    {
       inv[i]=p-(p/i)*inv[p%i]%p;
	   printf("%lld\n",inv[i]);
    } 
	return 0;
}