bzoj2186【SDOI2008】沙拉公主的困惑

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

Time Limit: 10 Sec Memory Limit: 259 MB

Submit: 2363 Solved: 779

[

id=2186" style="color:blue; text-decoration:none">Submit][

id=2186" style="color:blue; text-decoration:none">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

HINT

Source

数论

欧拉函数+线性筛法+
乘法逆元

数论题的做法简直不能再6，感觉自己智商严重不够用…

首先答案为phi(m!)*n!/m!%p。由于全部小于m!且与m!互质的数加上m!的整数倍都与m!互质，而其它数都不与m!互质。(正确性显然)

那么这个式子怎么求呢？？？

我们能够分成两部分来求，phi(m!)/mi和n!。

n!%p是非常easy预处理的。这里的主要问题是怎样求phi(m!)/m!。

令f(m)=phi(m!)/m!，依据phi(x)=x*(p1-1)/p1*(p2-1)/p2*…

可得f(m)=(p1-1)/p1*(p2-1)/p2*…当中pi为不大于m的质数

所以对于f(i)，假设i是质数f(i)=f(i-1)*(i-1)/m。否则f(i)=f(i-1)。

依据以上关系式能够预处理f(1)-f(10^7)。

每次询问仅仅须要输出f(m)*n!%p就可以。

#include<iostream>

#include<cstdio>

#include<cstdlib>

#include<cstring>

#include<cmath>

#include<algorithm>

#define F(i,j,n) for(int i=j;i<=n;i++)

#define D(i,j,n) for(int i=j;i>=n;i--)

#define ll long long

#define maxn 10000005

using namespace std;

int n,m,p,t;

ll fac[maxn],ans[maxn];

bool f[maxn];

inline 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*10+ch-'0';ch=getchar();}

	return x*f;

}

inline void exgcd(int a,int b,int &x,int &y)

{

	if (!b){x=1;y=0;return;}

	exgcd(b,a%b,x,y);

	int t=x;x=y;y=t-a/b*x;

}

inline int getinv(int a)

{

	int x=0,y=0;

	exgcd(a,p,x,y);

	return (x%p+p)%p;

}

int main()

{

	t=read();p=read();

	int x=10000000;

	fac[1]=1;

	F(i,2,x) fac[i]=fac[i-1]*i%p;

	ans[1]=1;

	F(i,2,x)

	{

		if (!f[i])

		{

			ans[i]=ans[i-1]*(i-1)%p*getinv(i)%p;

			F(j,2,x/i) f[i*j]=true;

		}

		else ans[i]=ans[i-1];

	}

	while (t--)

	{

		n=read();m=read();

		printf("%lld\n",ans[m]*fac[n]%p);

	}

}

巴特西