ans=(phi(m))*(n!/m!)
phi(m!)=m!*(p-1)/p p为m!的质因数
ans=n!*(p-1)/p
然后就是求逆元
用费马小定理求得 巨慢如雷电 10秒卡过 ........————————
这里说一下费马小定理求逆元
R是质数
求a关于R的逆元可以用快速幂
a^(R-2)%R即为a的逆元
因为 a*a^(R-2)=a^(R-1)=1(mod R)
好暴力啊……
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<algorithm>
#include<iostream>
#define P 10000000
#define ll long long
using namespace std;
bool prime[P+10];
int fac[P+10],ine[555555],pri[555555],ans[P+10];
int T,R,m,n,tot;
int cal(ll x,ll y,ll z)
{
ll ans=1;
while(y)
{
if(y&1) ans=ans*x%z;
x=x*x%z;
y>>=1;
}
return ans;
}
void get_pri()
{
fac[1]=1;
for(int i=2; i<=P; i++) fac[i]=(ll)fac[i-1]*i%R;
for(int i=2; i<=P; i++)
{
if(!prime[i]) pri[++tot]=i;
for(int j=1; i*pri[j]<=P; j++)
{
prime[i*pri[j]]=1;
if(i%pri[j]==0) break;
}
}
ans[1]=1;
for(int i=2; i<=P; i++)
{
ans[i]=ans[i-1];
if(!prime[i]) ans[i]=(ll)ans[i]*(i-1)%R*cal(i,R-2,R)%R;
}
}
int main()
{
scanf("%d%d",&T,&R);
get_pri();
while(T--)
{
scanf("%d%d",&n,&m);
printf("%lld\n",(ll)fac[n]*ans[m]%R);
}
return 0;
}