bzoj2186 [Sdoi2008]沙拉公主的诱惑

题目链接:bzoj2186
题目大意:
给定N,M,R。输出1至N!中与M!互质的数的数量对R取模后的值

题解:
数论、线性求逆元
可能我傻了吧
首先要知道x与y互质,则x+y也与y互质。
对于比 m ! m! m!小的 n ! n! n!中的数来说共有 ϕ ( m ! ) \phi(m!) ϕ(m!)个与 m ! m! m!互质,那么 n ! n! n!中就有 n ! m ! ∗ ϕ ( m ! ) \frac{n!}{m!}*\phi(m!) m!n!ϕ(m!)个与 m ! m! m!互质。预处理阶乘,把 ϕ ( m ! ) m ! \frac{\phi(m!)}{m!} m!ϕ(m!)一起搞,我脑子一定没有了。把 ϕ ( m ! ) \phi(m!) ϕ(m!)拆开就是 m ! ∗ ( 1 − 1 p 1 ) ∗ ( 1 − 1 p 2 ) ∗ … ( 1 − 1 p t ) m!*(1-\frac{1}{p_1})*(1-\frac{1}{p_2})*…(1-\frac{1}{p_t}) m!(1p11)(1p21)(1pt1)那个 m ! m! m!是可以约掉的啊!
所以就可以把 ( p 1 − 1 p 1 ) ∗ ( p 2 − 1 p 2 ) ∗ … ( p t − 1 p t ) (\frac{p_1-1}{p_1})*(\frac{p_2-1}{p_2})*…(\frac{p_t-1}{p_t}) (p1p11)(p2p21)(ptpt1) 放进线性筛里一起搞出来啊。
一直想着直接 n ! ∗ m ! − 1 ∗ ϕ ( m ! ) n!*{m!}^{-1}*\phi(m!) n!m!1ϕ(m!)【黑人问号.jpg
哪来那么大的数组给我存逆元???
线性求逆元的做法戳这里。【等我写完总结才有链接233

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
#define N 10000100

int ny[N],jc[N],jm[N],mod;
bool ispri[N];int cnt,pri[N/10];
void pre(int lim)
{
	cnt=0;jm[1]=1;
	for (int i=2;i<=lim;i++)
	{
		jm[i]=jm[i-1];
		if (!ispri[i]) {pri[++cnt]=i;jm[i]=(LL)jm[i]*(LL)(i-1)%mod*(LL)ny[i]%mod;}
		for (int j=1;j<=cnt && i*pri[j]<=lim;j++)
		{
			ispri[i*pri[j]]=true;
			if (i%pri[j]==0) break;
		}
	}
}
int main()
{
	//freopen("a.in","r",stdin);
	//freopen("a.out","w",stdout);
	int T,i,j,n,m,lim,ans;
	scanf("%d%d",&T,&mod);lim=1e7;
	jc[1]=1;for (i=2;i<=lim;i++) jc[i]=(LL)i*(LL)jc[i-1]%mod;
	ny[1]=1;for (i=2;i<=lim;i++) ny[i]=(LL)(mod-(mod/i))%mod*(LL)ny[mod%i]%mod;
	pre(lim);
	while (T--)
	{
		scanf("%d%d",&n,&m);
		ans=(LL)jc[n]*(LL)jm[m]%mod;
		printf("%d\n",ans);
	}
	return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值