题目大意:求出1~N!中与M!互质的数的个数,保证N>M。
前言:最讨厌数学了。。。
思路:因为保证了n>m,所以n!一定是m!的倍数。如果找到了一个x使得gcd(x,m!)==1,那么gcd(x + m!,m!) == 1一定成立,gcd(x + k * m!,m!) == 1(k >= 1)也一定成立。x的个数就是φ(m!),那么总的个数就是φ(m!) * n! / m!。之后各种推公式,就可以用线性筛来解决这个问题。(我太弱,讲不明白。。。
CODE:
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define MAX 10000100
using namespace std;
int T,P,n,m,primes;
long long prime[500100],ans[MAX],_link[MAX],rev[MAX];
bool notp[MAX];
void Pretreatment()
{
for(int i = 2; i < MAX; ++i) {
if(!notp[i]) prime[++primes] = i;
for(int j = 1;j <= primes && prime[j] * i < MAX; ++j) {
notp[prime[j] * i] = true;
if(i % prime[j] == 0)
break;
}
}
_link[0] = 1;
for(int i = 1; i < MAX; ++i)
_link[i] = _link[i - 1] * i % P;
rev[1] = 1;
for(int i = 2; i < MAX && i < P; ++i)
rev[i] = (P - P / i) * rev[P % i] % P;
ans[1] = 1;
for(int i = 2; i < MAX; ++i) {
if(!notp[i])
ans[i] = ans[i - 1] * (i - 1) % P * rev[i % P] % P;
else ans[i] = ans[i - 1];
}
}
int main()
{
cin >> T >> P;
Pretreatment();
for(int i = 1; i <= T; ++i) {
scanf("%d%d",&n,&m);
printf("%d\n",_link[n] * ans[m] % P);
}
return 0;
}