本文共 1762 字，大约阅读时间需要 5 分钟。
为了解决这个问题，我们需要计算一个数学表达式，它涉及到求和和模运算。这个表达式可以分解为两个部分：首先是双重求和，然后是减去相等的情况。我们将使用C++编写代码来高效地计算结果。
我们需要计算以下表达式：∑_{n}^i ∑_{m}^j (n mod i) * (m mod j)，其中i ≠ j。
为了高效计算，我们首先将问题分解为两个部分：
solve(n)和solve(m)来处理。solve2(n, m)来处理两个变量的联合求和，同时考虑它们的大小关系。#include#include #include #include using namespace std;const LL mod = 19940417;const LL inv = 3323403;LL sum(LL l, LL r) { return ((LL)(r - l + 1) * (LL)(l + r) / 2) % mod;}LL sum_2(LL n) { return ((LL)n * (LL)(n + 1) % mod * (LL)(2 * n + 1) % mod * inv % mod;}LL solve(LL n) { LL ans = 0; LL pos; for (LL i = 1; i <= n; i = pos + 1) { pos = (n / (n / i)); LL term = (n * (pos - i + 1) - (n / i) * sum(i, pos)) % mod; ans = (ans + term) % mod; } return ans;}LL solve2(LL n, LL m) { if (n > m) swap(n, m); LL ans = 0; LL pos; for (LL i = 1; i <= n; i = pos + 1) { pos = min(n / (n / i), m / (m / i)); LL t1 = (n * m % mod) * (pos - i + 1) % mod; LL t2 = ((n / i) % mod) * sum(i, pos) % mod; LL t3 = ((m / i) % mod) * ((sum_2(pos) - sum_2(i - 1)) % mod) % mod; LL t3 = (t3 + mod) % mod; LL term = (t1 - t2 + t3) % mod; ans = (ans + term) % mod; } return ans;}LL n, m;int main() { scanf("%lld %lld", &n, &m); LL res = (solve(n) * solve(m) % mod - solve2(n, m)) % mod; res = (res + mod) % mod; printf("%lld", res);}
通过这个方法，我们可以高效地计算所需的数学表达式，并确保结果在合理范围内。
转载地址：http://rdza.baihongyu.com/