18_ccpc_杭州_B Master of Phi

当代大学生的第一道数论题

2018年CCPC杭州站B题,Master of Phi 题目地址

---

题目大意：

给你 $F(x)={\sum }_{d|n}^n\phi (d)(\frac{n}{d})$ module 998244353, n 以唯一素数分解的形式给出，求 $F(x)$ 的值。

需要用到的知识：

1. $\phi (x)$是积性函数，积性函数的狄利克雷卷积同样也是积性函数。积性函数都满足若$gcd(n,m)=1$, 则 $f(nm)=f(n)f(m)$
2. $\phi (x)$性质：$n=p^k$ 时 $\phi (p^k)=(p-1)p^{k-1}$
3. $\phi (x)$性质：$n$为素数是 $\phi (n)=n-1$,特殊的$\phi (1)=1$
4. 整数的素数唯一分解定理，$n={\prod }_{i=1}^k{p}_i^{e_i}$ ，代入积性函数$g(n)={\sum }_{d|n}f(d)$时可以分解为$g(n)={\prod }_{i=1}^k{\sum }_{j=0}^{e_i}f(p_i^j)$

---

解题过程：

$$f(n)=\sum _{d|n}\phi (d)(\frac{n}{d})$$
由上 4 对 $n$进行素数分解，得
$$f(n)=\prod _{i=1}^k\sum _{j=1}^{e_i}\phi (p_i^j)(\frac{p_i^{e_i}}{p_i^j})$$
由 2 对上式化简，得
$$f(n)=\prod _{i=1}^k\sum _{j=1}^{e_i}(p_i-1)p_i^{j-1}(\frac{p_i^{e_i}}{p_i^j})$$
化简，得
$$f(n)=\prod _{i=1}^k(p_i-1)p_i^{e_i-1}\sum _{j=1}^{e_i}$$
上面的式子中没有$j=0$的情况，所以每个连乘中每项还要加一项 $p_i^{e_i}$,然后再乘，得
$$f(n)=\prod _{i=1}^k(p_i-1)p_i^{e_i-1}{e}_i+p_i^{e_i}$$

---

代码：

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn = 1e5+500;
const ll mod = 998244353;
ll ans,p,e;
int n;
ll po(ll a,ll b){
    ll an = 1;
    while(b){
        if(b&1) an = (an*a)%mod;
        a = (a*a)%mod;
        b>>=1;
    }
    return an;
}
int main(){
    int _;scanf("%d",&_);
    while(_--){
        scanf("%d",&n);
        ans = 1;
        while(n--){
            scanf("%lld%lld",&p,&e);
            ll t = po(p,e-1);
            ll temp = (((((p-1)*t)%mod)*e)%mod+po(p,e))%mod;
            ans = (ans * temp)%mod;
        }
        printf("%lld\n",ans);
    }
    return 0;
}