题目链接: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!∗(1−p11)∗(1−p21)∗…(1−pt1)那个
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})
(p1p1−1)∗(p2p2−1)∗…(ptpt−1) 放进线性筛里一起搞出来啊。
一直想着直接
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;
}