HDU6265狄利克雷卷积

下面就做一道欧拉函数|积性函数|狄利克雷卷积的题吧，相关原理在上一篇博文：点我跳转

题意：

​ 翻译成汉语，大概意思就是首先输入一个数T代表有T组样例，每组样例第一行m代表有m对<p,q>（其中p为素数，每对中的p各不相同），计算得到全部pi^qi的乘积n，最后求一下∑d|n φ(d) × n / d的值取模998244353

思路：

​ 由于p和q数据范围是1e8，我们没办法求出全部p^q的乘积n（爆long long），只能想办法化简公式和通过中间过程来得到答案，又因为给出的∑d|n φ(d) × n / d包含求和，我们遍历辣么大的数的求和铁定超时，所以还应该想办法把公式化简成O(1)的，也就是把求和给去掉

​ 我们接着按照上面的需求想问题：

​ 令函数f(n) = n，令h(n) = ∑d|n φ(d) × n / d，是不是发现了点什么！！！对，h(n）是φ(n)与f(n)的狄利克雷卷积。由于φ(n)是欧拉函数，所以它是积性函数，f(n) = n是单位函数，所以它是完全积性函数，那么h(n)一定也是积性函数，它满足a,b互质，则h(a * b) = h(a) * h(b)的性质

​ 又因为每对中的p都是不同的素数！！！因此每对pi^qi之间都是互质的，所以,我们就可以解决爆long long 问题：

h(n) = h(p1^q1 * p2^q2 * ... * pm^qm) = ∏ h(pi ^ qi)  //i为从1～m

​ 这样我们从求h(n)变成了求m个h(pi^qi)的乘积了

​ 此时pi^qi也会爆long long，我们先别管继续往下走

​ 这时候∑d|n φ(d) × n / d就变成了求∑d|p^q φ(d) × p^q / d .我们需要通过化简∑d|p^q φ(d) × p^q / d为O(1)来解决超时问题，怎么化简呢？手写了一份化简过程：

我们就可以把h(p^q) = ∑d|p^q φ(d) × p^q / d 化简为h(p^q) = p^q + q * p^(q - 1) * (p - 1)

​ 这样我们就可以通过h(n) = ∏ h(pi ^ qi)来求h(n)了（通过不断的取余998244353来防治爆longlong）

​ Over～

代码：

#include <stdio.h>
#include <string.h>
#include <inttypes.h>
#define ll int64_t
#define d int32_t
#define f double
#define r return
#define mod 998244353
#define mem(a) memset(a, 0, sizeof(a))
#define scanfd(a) scanf("%d", &a)
#define scanfl(a) scanf("%lld", &a)
#define printfl(a) printf("%lld\n", a)
#define For(i, star, endd) for (d i = star; i <= endd; i++)
#define Forr(i, endd, star) for (d i = endd; i >= star; i--)

//快速幂
ll quick(ll a, ll b) {
    ll ans = 1;
    while (b) {
        if (b & 1) {
            ans = ans * a % mod;
        }
        b >>= 1;
        a = a * a % mod;
    }
    r ans;
}


d main () {
    d T, m;
    ll p, q;
    scanfd(T);
    while (T--) {
        ll ans = 1;
        scanfd(m);
        For (i, 1, m) {
            ll res = 0;
            scanfl(p);
            scanfl(q);
            res = quick(p, q - 1) * q % mod * (p - 1) % mod;
            res = (res + quick(p, q)) % mod;
            ans = ans * res % mod;
        }
        printfl(ans);
    }
    r 0;
}

转载请注明出处！！！

如果有写的不对或者不全面的地方 可通过主页的联系方式进行指正，谢谢