[SDOI2008]沙拉公主的困惑----公式推导+快速幂+逆元+欧拉函数

题目链接

所以对于任何一个正整数x和y，想求小于等于x里和y互质的个数，就用x乘y的欧拉函数再除y就行了吧
因为y真正有效的是它的所有质因子，容斥再二项式推一下就是欧拉函数里分数连乘的那一坨

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
#define int long long
const int N = 1e7+5;
int mod;
int qmi(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1)res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
//筛质数
int primes[N],idx=0,id[N];
bool st[N];
void get_prime(int n){
    for(int i=2;i<=n;i++){
        if(!st[i])primes[++idx]=i;
        id[i]=idx;
        for(int j=1;primes[j]*i<=n;j++){
            st[primes[j]*i]=true;
            if(i%primes[j]==0)break;
        }
    }
}
int phi[N];// pi-1/pi相乘
int jc[N]; //阶乘
void get_phi(int n)
{
    phi[0]=1;
    for(int i=1;i<=idx;i++)
        phi[i]=phi[i-1]*(primes[i]-1)%mod*qmi(primes[i]%mod,mod-2)%mod;
    
    jc[0]=1;
    for(int i=1;i<=n;i++)jc[i]=jc[i-1]*i%mod;
}
signed main()
{
	int T,n,m;
    cin>>T>>mod;
    get_prime(N-1);
    get_phi(N-1);
    while(T--)
    {
        cin>>n>>m;
        cout<<phi[id[m]]*jc[n]%mod<<endl;
    }
}