【codevs2301】【BZOJ2186】沙拉公主的困惑，数论练习之逆元与φ

传送门1
传送门2
写在前面：数论！数论！数论！
思路：
1.分析出题目本意就是求phi(m!)*n!/m!%r（这一步大概挺难想的，想出来这个后面就简单多了）
2.n!%r可以O(10^7)预处理
3.由phi(x)=x×(p1-1)/p1×(p2-1)/p2×……×(pi-1)/pi，其中p1,p2,p3……pi为x的质因数且各不相同我们可以得到phi(m!)/m!=(p1-1)/p1×(p2-1)/p2×……×(pi-1)/pi，而且由于是阶乘，那么m!的质因数一定小于等于m,即小于10^7，所以同样可以预处理1-10^7中的质数并求出形同上式的答案，对于除法运算采取求逆元操作,即求出每个质因数关于r的逆元即可
代码：

#include<cstdio>#include<cstring>#include<iostream>#include<cmath> #define LL long long#define MAXN 10000000using namespace std;int t,r,n,m;bool pd[10000001];int phi[10000001],fac[10000001],prime[10000001];//这里的phi[i]并不是指i的φ，而是i!的φinline void exgcd(int a,int b,int &x,int &y){    if (b==0) {x=1;y=0;return;}    exgcd(b,a%b,x,y);    int x1=y,y1=x-a/b*y;    x=x1;y=y1;}inline int getinv(int p){    int x,y;    exgcd(p,r,x,y);    return (x%r+r)%r;}main(){    scanf("%d%d",&t,&r);    fac[0]=1;    phi[1]=1;    for (int i=1;i<=MAXN;i++)    fac[i]=(LL)fac[i-1]*i%r;//加强制类型转换防止乘法爆int，下同    for (int i=2;i<=MAXN;i++)    if (!pd[i])    {        phi[i]=(LL)phi[i-1]*(i-1)%r*getinv(i)%r;        for (int j=2;i*j<=MAXN;j++)        pd[i*j]=1;//最朴素的筛法，可以用欧拉筛的说    }    else phi[i]=phi[i-1];    while (t--)    {        scanf("%d%d",&n,&m);        printf("%lld\n",(LL)fac[n]*phi[m]%r);    }}