数学知识：乘法逆元

最新推荐文章于 2022-10-11 17:20:00 发布

zheng.ys

最新推荐文章于 2022-10-11 17:20:00 发布

阅读量307

点赞数 1

文章标签：模运算乘法逆元高效计算大整数 C++

本文链接：https://blog.csdn.net/weixin_66182206/article/details/124219916

版权

例如求 a1*a2*a3*...*an/b%mod 时，除法不能同时取模(既不能同时(a%mod) / (b%mod) )，这时可以通过乘法逆元来进行转换，既 $a/b\equiv a*x(mod)$ ,x为b的乘法逆元，因为a/b和a*x对mod取模的值一样，所以只需要计算(a%mod*x%mod)%mod即可

当编号是 $1^{2}+2^{2}+3^{2}+...+t^{2}$ 操作后灯会亮着;

令t=sqrt(n);

所以ans=t*(t+1)*(2*t+1)/6%(1e9+7);

显然没有很好的办法直接求答案(用__int128也能直接求或者把6拆成2*3也行，不过个人感觉乘法逆元方法要好点)

这时就可以求6的乘法逆元来求答案；

#include<iostream>
#include<cmath>
using namespace std;
typedef long long LL;

const LL mod=1e9+7;

LL inv(int t)//求乘法逆元
{
	if(t==1)
		return 1;
	else
		return (mod-mod/t)*inv(mod%t)%mod;
}

int main()
{
	LL n;
	cin>>n;
	LL t=sqrt(n);
	LL ans=t%mod*(t+1)%mod*(2*t+1)%mod*inv(6)%mod;
	cout<<ans<<"\n";
	
	return 0;
}