传送门
Description
大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在,请你帮助沙拉公主解决这个问题,由于可能张数非常大,你只需计算出对R取模后的答案即可。R是一个质数。
Input
第一行为两个整数T,R。R<=10^9+10,T<=10000,表示该组中测试数据数目,R为模后面T行,每行一对整数N,M,见题目描述 m<=n
Output
共T行,对于每一对N,M,输出1至N!中与M!素质的数的数量对R取模后的值
Sample Input
1 11
4 2
Sample Output
1
数据范围:
对于100%的数据,1 < = N , M < = 10000000
题解
看到数据范围T≤10000,就肯定不能在线处理了
所以要离线计算
题目要求的是1~N!中与m!互质的数的个数,M≤N,所以N!肯定是M!的倍数
所以我们可以将N!分成
N!M!
个M!大小的块,那么每一个块内都有
phi(M!)
个与M!互质的数,那么答案就等于
N!M!∗phi(M!)
继续化式子
N!M!∗M!∗∏phi(M!)i=1(pi−1)∏phi(M!)i=1pi(pi为所有与M!互质的数)
化简
N!∗∏phi(M!)i=1(pi−1)∏phi(M!)i=1pi
这样我们就可以暴力计算了
由于取模数r是一个质数,所以我们可以O(n)计算逆元
用ans数组记录式子中右边的分数就可以了
CODE:
#include<cstdio>
const int N=1e7+1;
int fac[N],inv[N],ans[N];
int prime[1000001];
bool b[N];
int n,m,x,y,num;
inline void calc()
{
for(int i=2;i<N;i++)
{
if(!b[i]) prime[++num]=i;
for(int j=1;j<=num&&prime[j]*i<N;j++)
{
b[prime[j]*i]=1;
if(i%prime[j]==0) break;
}
}
fac[1]=inv[1]=ans[1]=1;
for(int i=2;i<N;i++)
{
fac[i]=1ll*fac[i-1]*i%m;
inv[i]=1ll*(m-m/i)*inv[m%i]%m;
if(!b[i]) ans[i]=1ll*ans[i-1]*inv[i]%m*(i-1)%m;
else ans[i]=ans[i-1];
}
}
int main()
{
scanf("%d%d",&n,&m);
calc();
for(int i=1;i<=n;i++)
{
scanf("%d%d",&x,&y);
printf("%lld\n",1ll*fac[x]*ans[y]%m);
}
return 0;
}
总结
要学会灵活地使用题目中给的各种有用信息,化简式子,合理储存答案。