【乘法逆元】好题选讲

最新推荐文章于 2022-07-15 19:10:04 发布

ygmjsjdboy

最新推荐文章于 2022-07-15 19:10:04 发布

阅读量294

点赞数

分类专栏： ---------数论---------

本文链接：https://blog.csdn.net/ygmjsjdboy/article/details/102756897

版权

---------数论--------- 专栏收录该内容

8 篇文章 0 订阅

订阅专栏

一、【模板】有理数取余

题目：

$求 a / b m o d 19260817$

题解：

~~板子题没什么好说的~~
$\frac{a}{b}=a*inv(b)$
输入太大，输的时候边读边取模就可以辣

#include<bits/stdc++.h>
#define ll long long 
using namespace std;
ll a,b,mod=19260817;
inline ll read()
{
	ll ans=0;char s=getchar();
	while(s>'9'||s<'0')s=getchar();
	while(s>='0'&&s<='9')
	{
		ans=((ans<<1)%mod+(ans<<3)%mod+(s&15))%mod;
		s=getchar();
	}return ans;
}
ll quick_m(ll i,ll n)
{
	i%=mod;
	ll ans=1,ds=i;
	while(n)
	{
		if(n&1)
		{
			ans=ans*ds%mod;
		}
		ds=ds*ds%mod;
		n>>=1;
	}return ans%mod;
}
signed main(){
	a=read();b=read();
	if(b==0){printf("Angry!");return 0;}
	printf("%lld",a*quick_m(b,mod-2)%mod);
	return 0;
}

二、 Sumdiv

题目：

求 $a^b的约数和mod \ 9901$

题解：

若a的唯一质因数分解为 $a=p_1^{k_1}p_2^{k_2}p_3^{k_3}...p_n^{k_n}$
则 $a^b$ 的唯一质因数分解为 $a^b=p_1^{k_1b}p_2^{k_2b}p_3^{k_3b}...p_n^{k_nb}$
约数和数学表达
$s=(1+p_1+p_1^2...+p_1^{k_1b})(1+p_2+p_2^2...+p_2^{k_2b})...(1+p_n+p_n^2...+p_n^{k_nb})$
等比数列求解
$s=(\frac{p_1^{k_1b+1}-1}{p_1-1})(\frac{p_2^{k_2b+1}-1}{p_2-1})...(\frac{p_n^{k_nb+1}-1}{p_n-1})$
我们逆元+快速幂？
不， $p i - 1 可能为 9901 的倍数，不与 9901 互质，此时没有逆元$
那怎么办？
这里有个非常nice 的idea,上图
$若求\frac a b (mod \ p ),且刚好有a为b的倍数$
在这里插入图片描述
注意：有多组数据。。。还有快速幂时可能会爆long long，用龟速乘。
~~卡在RE的应该不止我一个吧~~（现已解决，a>1卡了几乎两页）

#include<iostream>
#include<string.h>
#include<stdio.h>
#define ll long long
using namespace std;
ll qk_c(ll a,ll b,ll mod)
{
	ll ans=0;a%=mod;
	while(b)
	{
		if(b&1)
		{
			ans=(ans+a)%mod;
		}
		a=(a+a)%mod;
		b>>=1;
	}
	return ans;
}
ll qk_m(ll a,ll k,ll mod)
{
	ll ans=1;a%=mod;
	while(k)
	{
		if(k&1)
		{
			ans=qk_c(ans,a,mod);
		}
		a=qk_c(a,a,mod);
		k>>=1;
	}
	return ans;
}
const int N=10005;
int z[N],p;
bool vis[N];
void init()
{
	for(int i=2;i<N;i++)
	{
		if(!vis[i])z[++p]=i;
		for(int j=1;j<=p&&i*z[j]<N;j++)
		{
			vis[i*z[j]]=1;
			if(i%z[j]==0)break;
		}
	}
}
const int mod=9901;
ll a,b,ans;
int k;
int main()
{ 
	init();
	while(cin>>a>>b)
	{
		ans=1;
		for(int i=1;i<=p&&z[i]*z[i]<=a;i++)
		{
			if(a%z[i]==0)
			{
				k=0;while(a%z[i]==0)a/=z[i],k++;
				ll M=mod*(z[i]-1);
				ans*=(qk_m(z[i],k*b+1,M)+M-1)/(z[i]-1);
				ans%=mod;
			}
		}
		if(a>1)
		{
			ll M=mod*(a-1);
			ans*=(qk_m(a,b+1,M)+M-1)/(a-1);
			ans%=mod;
		}
		cout<<ans<<endl;
	}
    return 0;
}

三、沙拉公主的困惑

题目：

求 $1 - > N!$ 中与 $M!$ 互质的个数

题解：

推演一下，最后得 $\frac {N!}{M!} \phi{(M!)}$
=> $N!(\frac {1-p_1}{p_1})(\frac {1-p_2}{p_2})...(\frac {1-p_k}{p_k})$
一堆预处理，质数+线性求逆+阶乘。。。

#include<bits/stdc++.h>
#define ll long long 
using namespace std;
inline ll read()
{
	ll ans=0;char s=getchar();
	while(!isdigit(s))s=getchar();
	while(isdigit(s))
	{
		ans=(ans*10)+(s&15);
		s=getchar();
	}return ans;
}
const int N=1e7+10;
int T;
ll n,m,mod;
int z[N],p;
bool vis[N];
void get_prime()
{
	for(int i=2;i<N;i++)
	{
		if(!vis[i])z[++p]=i;
		for(int j=1;j<=p&&i*z[j]<N;j++)
		{
			vis[i*z[j]]=1;
			if(i%z[j]==0)break;
		}
	}
}
ll inv[N],mul[N],pt[N];
int main()
{
	scanf("%d%lld",&T,&mod);
	get_prime();
	inv[1]=1;inv[0]=1;
	for(int i=2;i<N;i++)
		inv[i]=inv[mod%i]*(mod-mod/i)%mod;
	mul[1]=1;
	for(int i=2;i<N;i++)
		mul[i]=mul[i-1]*i%mod;
	pt[1]=1;
	for(int i=2;i<N;i++)
		if(!vis[i])
			 pt[i]=pt[i-1]*(i-1)%mod*inv[i%mod]%mod;
		else pt[i]=pt[i-1];
	while(T--)
	{
		n=read();m=read();
		printf("%lld\n",mul[n]*pt[m]%mod);
	}
}//假装没开O2就过了