【学术篇】SDOI2008 沙拉公主的困惑

最新推荐文章于 2021-08-17 08:51:53 发布

Enzymii

最新推荐文章于 2021-08-17 08:51:53 发布

阅读量411

点赞数 1

分类专栏：数论什么都去死吧文章标签：数论逆元欧拉函数

本文链接：https://blog.csdn.net/Enzymii/article/details/78823789

版权

数论什么都去死吧专栏收录该内容

8 篇文章 0 订阅

订阅专栏

传送门!

题目在这里…

题目大意?

~~难道不是说的很清楚了么OvO~~
求n!中与m!互质的数的个数..

题目分析.

显然的数论… 所以就是化式子呗..
一个很显然的性质就是如果 $gcd(a,b)=1$ ,那么 $gcd(a+kb,b)=1$ …
而题目中说了 $m\leqslant n$ , ∴ $m!|n!$
于是我们只需要计算 $m!$ 中与 $m!$ 互质的数的个数,然后乘以 $\frac{n!}{m!}$ 即可..
我们发现上面加粗的这一坨就是 $\varphi(m!)$ 嘛…
所以 $ans=\varphi(m!)*\frac{n!}{m!}$
又有 $\varphi(x)=x*\prod_{i}^{n}(1-\frac{1}{p_i})$ 其中 $p_i$ 表示x的质因数…
而 $m!=1*2*...*m$ , 所以 $m!$ 的质因数很显然就是不大于 $m$ 的质数…
然后带入上式约掉 $m!$ 就有了 $ans=n!*\prod_{i}^{n}\frac{p_i-1}{p_i}$ (其中 $p_i\leqslant m$ 且 $p_i$ 为质数)…
由于多组询问, ~~而且内存开了256MB不是~~ 所以我们要预处理… 不然会T…
由于上式, 我们要预处理的东西有:
- 筛素数(简单欧拉筛)
- 阶乘(顺着乘一遍取模就行了)
- 逆元(要递推求出所有数的哦) (所以最好用 $O(n)$ 的, 不会的话直接看代码就行了百度一下一堆详细讲解OvO)
- $mul_i=\prod_{i}^{n}\frac{p_i-1}{p_i}$ 这一坨东西…(不大于 $m$ 的质数 $p_i$ 们的 $(1-\frac{1}{p_i})$ 的乘积…)
然后处理这一坨的时候也很容易…递推即可.. 显然, 我们有
1. 当 $i$ 是质数时, $mul_i=mul_{i-1}*\frac{i-1}{i}$
2. 否则 $mul_i=mul_{i-1}$ 即可…
这样就做完了.

实现代码:

#include <cstdio>
typedef long long LL;
const int X=1e7+3;
inline int gn(int a=0,char c=0){
    for(;c<48||c>57;c=getchar());
    for(;c>47&&c<58;c=getchar())
        a=a*10+c-48; return a;
}
int inv[X],fac[X],eu[X],mul[X],pri[X/10],tot;
bool notp[X]; int T,R,M,N;
void prime(){
    notp[1]=1;
    for(int i=2;i<X;++i){
        if(!notp[i])pri[++tot]=i;
        for(int j=1;j<=tot&&i*pri[j]<=1e7;++j){
            notp[i*pri[j]]=1; if(i%pri[j]==0) break;
        }
    }
}
void calcinv(){
    inv[1]=1;
    for(int i=2;i<X;++i){
        inv[i]=(LL)(R-R/i)*inv[R%i]%R;
        if(inv[i]<0) inv[i]+=R;
    }
}
void calcfac(){
    fac[1]=1;
    for(int i=2;i<X;++i)
        fac[i]=(LL)fac[i-1]*i%R;
}
void calcmul(){
    mul[1]=1;
    for(int i=2;i<X;++i)
        if(!notp[i]) mul[i]=(LL)mul[i-1]*(i-1)%R*inv[i]%R;
        else mul[i]=mul[i-1];
}
int main(){
    T=gn(),R=gn();
    prime(); calcinv(); calcfac(); calcmul();
    while(T--){
        N=gn(),M=gn();
        printf("%d\n",(int)((LL)fac[N]*mul[M]%R));
    }
}