大富翁国因为通货膨胀，以及假钞泛滥，政府决定推出一项新的政策：现有钞票编号范围为1到N的阶乘，但是，政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在，请你帮助沙拉公主解决这个问题，由于可能张数非常大，你只需计算出对R取模后的答案即可。R是一个质数。
标签：
大富翁国因为通货膨胀，以及假钞泛滥，政府决定推出一项新的政策：现有钞票编号范围为1到N的阶乘，但是，政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在，请你帮助沙拉公主解决这个问题，由于可能张数非常大，你只需计算出对R取模后的答案即可。R是一个质数。
第一行为两个整数T，R。R<=10^9+10，T<=10000，表示该组中测试数据数目，R为模后面T行，每行一对整数N，M，见题目描述 m<=n
共T行，对于每一对N，M，输出1至N！中与M！素质的数的数量对R取模后的值
题意：求
分析：因为
对于两个正整数
本结论是很好证明的，因为
结论成立。那么对于本题，答案就是
其中
求逆元就有技巧了，用刚刚介绍的递推法预处理，否则会TLE的。
代码：
#include <iostream>
#include <string.h>
#include <stdio.h>
#include <bitset>
using namespace std;
typedef long long LL;
const int N = 10000005;
bitset<N> prime;
void isprime()
{
prime.set();
for(int i=2; i<N; i++)
{
if(prime[i])
{
for(int j=i+i; j<N; j+=i)
prime[j] = false;
}
}
}
LL ans1[N],ans2[N];
LL inv[N];
int main()
{
isprime();
int MOD,m,n,T;
scanf("%d%d",&T,&MOD);
ans1[0] = 1;
for(int i=1; i<N; i++)
ans1[i] = ans1[i-1] * i % MOD;
inv[1] = 1;
for(int i=2;i<N;i++)
{
if(i >= MOD) break;
inv[i] = (MOD - MOD / i) * inv[MOD % i] % MOD;
}
ans2[1] = 1;
for(int i=2; i<N; i++)
{
if(prime[i])
{
ans2[i] = ans2[i-1] * (i - 1) % MOD;
ans2[i] = ans2[i] * inv[i % MOD] % MOD;
}
else
{
ans2[i] = ans2[i-1];
}
}
while(T--)
{
scanf("%d%d",&n,&m);
LL ans = ans1[n] * ans2[m] % MOD;
printf("%lld\n",ans);
}
return 0;
}
接下来还有一个关于逆元的有意思的题目，描述如下
证明：由
其中
所以只需要证明
所以进一步得到
证明完毕！
标签：
原文地址：http://blog.csdn.net/liangzhaoyang1/article/details/51104532
