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

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

HINT

Source

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

欧拉函数+逆元+思路~

易知答案=phi[m!]*(n!/m!)，

因为phi[m!]=m!*(k1-1)/k1*(k2-1)/k2*...，其中ki是m!的质因子，所以带入原式，得n!*(k1-1)/k1*(k2-1)/k2*...

ki相当于是1~m的所有质因数，所以我们预处理出1~1e7的所有质因数以及它们的逆元（用扩展欧几里得），然后用num[i]记录下m=i时的n！要乘的东西。

又遇到了昨天的那种情况，于是问了游神，被告知是因为a数组开小了，改大了就真的A了，真是神奇啊~

#include<cstdio>
#include<iostream>
using namespace std;
#define ll long long

int t,mod,n,m,jie[10000001],a[800001],num[10000001],ni[10000001];
bool b[10000001];

int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0' && ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}

void exgcd(int u,int v,int &x,int &y)
{
	if(!v)
	{
		x=1;y=0;return;
	}
	exgcd(v,u%v,x,y);
	int t=x;x=y;y=t-u/v*y;
}

int cal(int u)
{
	int x,y;
	exgcd(u,mod,x,y);
	return (x%mod+mod)%mod;
}

void init()
{
	jie[1]=1;
	for(int i=2;i<=10000000;i++)
	{
		jie[i]=(ll)jie[i-1]*i%mod;
		if(!b[i]) a[++a[0]]=i,ni[i]=cal(i);
		for(int j=1;i*a[j]<=10000000 && j<=a[0];j++)
		{
			b[i*a[j]]=1;
			if(!(i%a[j])) break;
		}
	}
	num[1]=1;
	for(int i=2;i<=10000000;i++)
	{
		num[i]=num[i-1];
		if(!b[i]) num[i]=(ll)num[i]*(i-1)%mod*ni[i]%mod;
	}
}

int main()
{
	t=read();mod=read();init();
	while(t--)
	{
		n=read();m=read();
		printf("%d\n",(ll)jie[n]*num[m]%mod);
	}
	return 0;
}