Function!

$f_a(x)=a^x(a>0,a=1)$，我們要求$\sum _{a=2}^n\left(a\sum _{b=a}^n?f_a^{?1}(b)??f_b^{?1}(a)?\right)$。
$$\begin{gathered} \sum _{a=2}^n\left(a\sum _{b=a}^n?f_a^{?1}(b)??f_b^{?1}(a)?\right) \\ \sum _{a=2}^n\left(a\sum _{b=a}^n?{\log }_{a}b??{\log }_{b}a?\right) \\ b\ge a，則有?{\log }_{b}a?=1 \\ \sum _{a=2}^n\left(a\sum _{b=a}^n?{\log }_{a}b?\right) \end{gathered}$$

且容易發現，當$a\times a>n$時，有$?{\log }_{a}b?$恒為$1$，所以可以單獨用公式計算。

#include <bits/stdc++.h>using namespace std;const int mod = 998244353, inv2 = mod + 1 >> 1, inv6 = (mod + 1) / 6;typedef long long ll;ll calc1(ll l, ll r) {return (l + r) % mod * ((r - l + 1) % mod) % mod * inv2 % mod; }ll calc2(ll n) {n %= mod;return n * (n + 1) % mod * (2 * n + 1) % mod * inv6 % mod; }int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);ll a, n, ans = 0;cin >> n;for (a = 2; a * a <= n; a++) {for (ll l = a, r, k = 1; l <= n; l = r + 1, k++) {r = min(l * a - 1, n);int tot = (r - l + 1) % mod;ans = (ans + a * tot % mod * k % mod) % mod;}}// for (; a <= n; a++) {//原本我們是這樣算的，當時這里可以變成，自然冪次求和的形式，所以可以快速算出來。// ans = (ans + a * (n - a + 1) % mod) % mod;// }ans = (ans + (n + 1) % mod * calc1(a, n) % mod) % mod;ans = ((ans - calc2(n) + calc2(a - 1)) % mod + mod) % mod;cout << ans << "\n";return 0; }

總結

以上是生活随笔為你收集整理的Function!（计蒜客 - 42386）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。