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

Time Limit: 10 Sec   Memory Limit: 259 MB 

Submit: 2616   Solved: 880 

[ ​​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

HINT

Source

[ ​​Submit​​​][ ​​Status​​​][ ​​Discuss​​] 

​​HOME​​ Back


大视野2186: 沙拉公主的困惑(求逆元）_i++


代码：


#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;
}


接下来还有一个关于逆元的有意思的题目，描述如下


     


大视野2186: 沙拉公主的困惑(求逆元）_#include_02


      证明完毕！