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

传送门

题意：求1-n!中与m!互质的数的个数，其中m小于n。

思路：由于m小于n，所以n!一定是m!的倍数，小于m!与m!互质的数的个数为phi(m!)，所以1-n!中与m!互质的数的个数为n!/m!*phi(m!)=n!/m!*m!*(1-1/p1)*(1-1/p2)*...(1-1/pk) =n!(1-1/p1)*(1-1/p2)*...(1-1/pk)其中p1-pk是m!的素因子，其实也就是小于m的所有素数，最后结果再对一个素数取余。

这个题时间卡的确实很紧，在求逆元的时候用的是o(p)复杂度来求：

void init()
{
inv[0]=inv[1]=1;
for(int i=2;i<maxn;i++)
{
if(i>=mod) break;
inv[i]=((mod-mod/i)*inv[mod%i])%mod;
}
}

而且在求(1-1/p1)*(1-1/p2)*...(1-1/pk)时先预处理，否则也是超时

完整代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
const int maxn=1e7+5;
bool prime[maxn];
int cnt;
int mod;
ll inv[maxn];
ll jc[maxn];
ll res[maxn];
void isprime()
{
prime[0]=false;
prime[1]=false;
prime[2]=true;
for(int i=3;i<maxn;i+=2)
{
prime[i]=true;
prime[i+1]=false;
}
for(int i=3;i<=sqrt((double)maxn);i+=2)
if(prime[i])
for(int j=i+i;j<maxn;j+=i)
prime[j]=false;
}
void init()
{
inv[0]=inv[1]=1;
for(int i=2;i<maxn;i++)
{
if(i>=mod) break;
inv[i]=((mod-mod/i)*inv[mod%i])%mod;
}
}
void factorial()
{
jc[0]=jc[1]=1;
for(int i=2;i<maxn;i++)
{
jc[i]=jc[i-1]*i%mod;
}
}
int main()
{
int t;
scanf("%d%d",&t,&mod);
isprime();
init();
factorial();

res[1]=1;
for(int i=2;i<maxn;i++)
{
if(prime[i])
{
res[i]=res[i-1]*(i-1)%mod;
res[i]=res[i]*inv[i%mod]%mod;
}
else
res[i]=res[i-1];
}

while(t--)
{
int n,m;
scanf("%d%d",&n,&m);
ll ans=jc
*res[m]%mod;
printf("%lld\n",ans);
}
return 0;
}