bzoj2186: [Sdoi2008]沙拉公主的困惑

链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2186

题意：中文题。。。

分析：题目要求1~n!内与m!互质的数的个数。。我们回顾一下gcd(a,b)我们发现gcd(a,b)=gcd(a,a+b)这是什么呢？设x<=m!在这个题意中就有gcd(x,m!)=gcd(x+m!,m!)=gcd(x+2*m!,m!)....那么这个题目就变成了求phi(m!)*(n!/m!)，且phi(m!)=m!*(Pi-1)/Pi，Pi为m!的质因子即小于m的所有质数，那么就有ans=n!*(Pi-1)/Pi，题目保证R为质数，所有我们只要求Pi的逆元即可，但是这题出题人估计就是这种做法但是缺忽略了Pi>=R这种情况，这个时候R是不存在关于模R的逆元的，这个时候要用其他方法处理了，根据之前的分析可以知道是n!/m!这里的处理出了问题，我们其实只要对(m+1)*(m+2)*...*n<R且m>=R这种情况进行特殊处理就好了，而这种复杂度也很低，暴力处理掉就行了。后半段特殊处理是瞎bb的，没有实现，如果有什么问题，请多指教。

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=10000010;
const int MAX=151;
const int MOD1=100000007;
const int MOD2=100000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=5000011;
const ll INF=10000000010;
typedef unsigned long long ull;
ll h[N],f[N];
int R,a[N],bo[N];
void deal() {
    int i,j,k=0,n=10000000;
    memset(bo,0,sizeof(bo));
    bo[1]=h[1]=f[1]=1;
    for (i=2;i<=n;i++) {
        h[i]=h[i-1]*i%R;
        if (!bo[i]) a[++k]=i;
        for (j=1;j<=k;j++) {
            if (a[j]*i>n) break ;
            bo[a[j]*i]=1;
            if (i%a[j]==0) break ;
        }
    }
    a[1]=1;
    for (i=2;i<=n;i++) a[i]=(ll)(R-R/i)*a[R%i]%R;
    for (i=2;i<=n;i++)
    if (bo[i]) f[i]=f[i-1];
    else f[i]=(f[i-1]*(i-1)%R)*a[i]%R;
}
int main()
{
    int n,m,t;
    scanf("%d%d", &t, &R);
    deal();
    while (t--) {
        scanf("%d%d", &n, &m);
        printf("%lld\n", (h[n]*f[m]%R+R)%R);
    }
    return 0;
}