4695. 【GDOI2017模拟8.14】佐助的难题

Description

由于佐助整天叫鸣人做“吊车尾的”，鸣人很不爽，于是他决定向佐助提出挑战，佐助可是个要复仇的男人啊，他才没空跟个小屁孩较劲，于是他随便给了一道题，如果鸣人没做出来，那他就要自认是“吊车尾的”，如果做出来了。。。到时再说吧
题目是这样的：
求在1到n!范围内，与m!互质的数的数量，由于答案太大，只需计算答案对R取模之后的答案即可，保证R是一个质数
鸣人作为一个整天上课睡觉的人，当然。。不会了。。于是他就向你求助了哈。

Input

第一行为两个整数T,R。表示该数据中的数据组数，R为模数，后面T行每行一对n,m

Output

共T行，对于每一对n,m，输出在1到n!范围内，与m!互质的数的数量

Sample Input

1 11
4 2

Sample Output

1

Data Constraint

对于30%的数据，n<=10，T<=30
对于60%的数据n<=3000000，T<=10000
对于100%的数据n<=10000000，T<=10000
保证R为质数，m<=n,n < R

Solution

题意相当于求出在n的阶乘中不含<=m的质因子的数的个数

由于m较小，我们可以把这些质数全筛出来

不含某些质因子的题，是容斥的经典模型，但是这里质数数量较大，普通容斥显然不行，我们要用线性容斥。

设m以内的质数为p1.p2....pk.

则答案为.

注意可以求逆元时用exgcd,或用常数优化。

Code

#include<cstdio> 
#include<algorithm>
#include<cstring>
#include<cstdlib>
#define I int
#define F(i,a,b) for(I i=a;i<=b;i++)
#define Fd(i,a,b) for(I i=a;i>=b;i--)
#define mem(a,b) memset(a,b,sizeof(a))
#define N 10000000
#define ll long long
#define P(x) freopen(x".in","r",stdin);freopen(x".out","w",stdout)
#define rt return
using namespace std;
void rd(I &x){
	x=0;char ch=getchar();
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
}
I T,R,n,m,f[N+10],inv[N+10],ans[N+10];
I p[N/2],bz[N+10];
I main(){
	P("difficult problem");
	rd(T),rd(R);
	f[0]=f[1]=inv[1]=ans[1]=1;
	F(i,2,N){
		ans[i]=ans[i-1];
		f[i]=((ll)f[i-1]*i)-((ll)f[i-1]*i)/R*R;
		inv[i]=(R-R/i)*(ll)inv[R-R/i*i]%R;
		if(!bz[i]){
			p[++p[0]]=i;
			ans[i]=((ll)ans[i]*(i-1))%R*(ll)inv[i]%R;
		}
		F(j,1,p[0]){
			if((ll)i*p[j]>N) break;
			bz[i*p[j]]=1;
			if(i-i/p[j]*p[j]==0) break;
		}
	}
	while(T--){
		rd(n),rd(m);
		printf("%d\n",((ll)ans[m]*f[n])%R);
	}
	rt 0;
}

---

---

作者：zsjzliziyang 
QQ:1634151125 
转载及修改请注明 
本文地址:https://blog.csdn.net/zsjzliziyang/article/details/99717744