BZOJ 2186: [Sdoi2008]沙拉公主的困惑

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

Time Limit:10 Sec  Memory Limit: 259 MB
Submit: 3945  Solved: 1370
[Submit][Status][Discuss]

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

ans=phi(m!)*n!/m!

多组数据，所以预处理阶乘，素数逆元

#include<cmath>
#include<ctime>
#include<cstdio> 
#include<climits>
#include<cstring> 
#include<cstdlib> 
#include<iostream> 
#include<algorithm> 
#include<iomanip> 
#include<vector> 
#include<string>
#include<queue>  
#include<map> 
#include<set>
using namespace std; 
typedef long long ll;
inline int read() 
{ 
	int x=0,f=1;char ch=getchar(); 
	while(ch<'0'|ch>'9'){if(ch=='-')f=-1;ch=getchar();} 
	while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();} 
	return f*x; 
}
const int N=10000010; 
int cnt,n,m,mod,k,prime[N>>3],ans[N],x,y,inverse[N],factorial[N];
bool book[N];
inline void exgcd(ll a,ll b)
{
	if(b==0){x=1;y=0;return ;}
	exgcd(b,a%b);
	ll t=x;x=y;y=t-a/b*y;
}
inline void initial()
{
	factorial[1]=1;for(int i=2;i<=N;i++)factorial[i]=(ll)factorial[i-1]*i%mod;
	inverse[1]=1;for(int i=2;i<=N;i++)
	{
		if(!book[i])prime[++cnt]=i,exgcd(i,mod),inverse[i]=(x%mod+mod)%mod;
		for(int j=1;j<=cnt&&prime[j]*i<=N;j++)
		{
			book[i*prime[j]]=1;
			if(i%prime[j]==0)break;
		}
	}ans[1]=1;
	for(int i=2;i<=N;i++){ans[i]=ans[i-1];if(!book[i])ans[i]=(ll)ans[i]*inverse[i]%mod*(i-1)%mod;}
}
int main()
{
	k=read();mod=read();
	initial();
	while(k--)
	{
		n=read();m=read();
		printf("%d\n",(ll)ans[m]*factorial[n]%mod);
	}
}