乘法逆元模板题题解 #cpp

最新推荐文章于 2024-08-01 16:50:47 发布

原创最新推荐文章于 2024-08-01 16:50:47 发布 · 313 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#c++ #线性代数

c++暴零生涯同时被 3 个专栏收录

57 篇文章

订阅专栏

数论

14 篇文章

订阅专栏

题解

12 篇文章

订阅专栏

本文介绍了如何解决洛谷P3811题目，包括80分暴力枚举的解决方案，以及通过添加inline关键字优化后的满分策略。探讨了求解一连串数mod p逆元的公式，并提供了相关代码。

题目洛谷 P3811

传送门：题目
80分做法
6个点 TLE了一个点
思路：暴力枚举，求每一个数的乘法逆元，用我博客里的代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll n,p;
ll exgcd(ll a,ll b,ll& x,ll& y)
{
	if(b==0)
	{
		x=1,y=0;
		return a;
	}
	ll g=exgcd(b,a%b,x,y);
	ll t;
	t=x;
	x=y;
	y=t-a/b*y;
	return g;
}
ll niyuan(ll a,ll b)//求a对b取模的逆元
{
	ll x,y;
	ll aa=exgcd(a,b,x,y);
	return (x + b) % b;
}
int main()
{
	scanf("%ld%ld",&n,&p);
	for(ll i=1;i<=n;i++)
	{
		ll answer=niyuan(i,p);
		printf("%ld\n",answer);
	}
	return 0;
}

满分做法:
后来我加了个 inline
~~满分了~~
还多 TLE 了一个点

要知道：
可以用公式求一连串数 mod p 的逆元
公式在代码里，自己看看吧
或者我再发篇博客
代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll n,p,inv[3000001];
int main()
{
	scanf("%ld%ld",&n,&p);
    inv[1] = 1;
    printf("1\n");
    for(int i=2;i<=n;++i)
    {
	    inv[i]=(p - p / i)*inv[p%i]%p;
	    printf("%ld\n",inv[i]);
    }
	return 0;
}