显然答案等于
ϕ(m!)∗(n!/m!)
n!很容易求。
ϕ(x)=x∗p1−1p1∗p2−1p2∗...
所以
ϕ(m!)m!=p1−1p1∗p2−1p2∗...
pi为<=m的质数。
定义其为f(i).
则有
{f(i)=f(i−1) x为合数f(i)=f(i−1)∗(i−1)/m x为质数
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
char ch=getchar();int f=0;
while(ch<'0'||ch>'9')ch=getchar();
while(ch>='0'&&ch<='9'){f=(f<<1)+(f<<3)+ch-'0';ch=getchar();}
return f;
}
bool check[10000005];long long prime[3000005],fac[10000005],n,mod,tot,f[10000005];
inline void exgcd(int a,int b,int &x,int &y)
{
if (!b){x=1;y=0;return;}
exgcd(b,a%b,x,y);
int t=x;x=y;y=t-a/b*x;
}
inline int getinv(int a)
{
int x=0,y=0;
exgcd(a,mod,x,y);
return (x%mod+mod)%mod;
}
int main()
{
n=read(),mod=read();
for(int i=2;i<=10000000;i++)
{
if(!check[i])
{
prime[tot++]=i;
}
for(int j=0;j<tot;j++)
{
if(i*prime[j]>10000000)
break;
check[i*prime[j]]=1;
if(i%prime[j]==0)
break;
}
}
f[1]=1;fac[1]=1;
for(int i=2;i<=10000000;i++)
{
if(!check[i])
f[i]=1LL*(1LL*f[i-1]*(i-1))%mod*getinv(i)%mod;
else f[i]=f[i-1];
fac[i]=(1LL*fac[i-1]*i)%mod;
}
for(int i=1;i<=n;i++)
{
int x=read(),y=read();
printf("%lld\n",1LL*f[y]*fac[x]%mod);
}
}