G. Problem G. 你说的对,但是不如莫反

    传统题 2000ms 256MiB

Problem G. 你说的对,但是不如莫反

该比赛已结束,您无法在比赛模式下递交该题目。您可以点击“在题库中打开”以普通模式查看和递交本题。

Problem G. 你说的对,但是不如莫反

时间限制 : 2000 ms

空间限制 : 256 MB

题目背景

5dde3f5ef77d3b9430458c41bd0ce43

image-20231107130300270

以下文字仅为娱乐,与题目本身无关。

你说的对,但是感觉不如莫反。莫比乌斯反演,是一种数论技巧。设 τ\tau 是狄利克雷卷积单位元 τ(x)=[x=1]\tau(x)=[x=1]uu 是乘法单位元函数 u(x)=1u(x)=1,则莫比乌斯函数 μ\mu 定义为 uu 在狄利克雷卷积下的逆元,即 μu=τ\mu * u=\tau。莫比乌斯反演归根到底就是对于数论函数f,gf,g,有 f=gug=fμf=g*u\Longleftrightarrow g=f*\mu(这里 * 是狄利克雷卷积运算)。你的数学很差,我现在每天用莫反都能求 1e5 次数据规模 1e12 的积性函数前缀和,每个月差不多 3e6 次反演变换,也就是现实生活中 3e18 次加法运算,换算过来最少也要算 1000 年。虽然我只有 14 岁,但是已经超越了中国大多数人(包括你)的水平,这便是莫反给我的骄傲的资本。

题目描述

pzr 刚刚学会了莫比乌斯反演,使用莫比乌斯反演的技巧,可以很快地对于一个数组求出

i=1nj=1ngcd(ai,aj)\sum_{i=1}^n\sum_{j=1}^n \gcd(a_i,a_j) $$ \sum_{i=1}^n\sum_{j=1}^n \operatorname{lcm}(a_i,a_j) $$

等信息。

现在,他想对于长度为 nnaa 数组求出:

$$ \sum_{i=1}^n\sum_{j=1}^n [\gcd(a_i,a_j)\ge 300]\times a_i\times a_j \times \gcd(a_i,a_j) $$

其中,gcd(ai,aj)\gcd(a_i,a_j) 表示 aia_iaja_j 的最大公约数。公式 [s][s] 在命题 ss 为真时取 11,其余时刻取 00

换句话说,我们每次选出数组中任意两个下标 (i,j)(i, j)ii 可以等于 jj),且满足 gcd(ai,aj)300\gcd(a_i,a_j)\ge 300 。 对于任意这样的 (i,j)(i, j),都可以计算出 ai×aj×gcd(ai,aj)a_i\times a_j\times \gcd(a_i,a_j) 的值,请求出它们的和。

由于答案过大,请对 998244353 取模后输出

输入格式

第一行一个整数 nn,表示数组的长度。

接下来一行 nn 个整数 a1,a2,,ana_1,a_2,\cdots ,a_n,表示这个数组。

输出格式

仅包含一个整数,表示公式对 998244353 取模的结果

样例输入1

8
5 4 3 2 1 500 2000 9997

样例输出1

996717796

样例 1 解释

最大公约数大于 300300 的取值方案有以下几种:

(i,j)=(6,6)(ai,aj)=(500,500):gcd=500(i,j)=(6,6)\to (a_i,a_j)=(500, 500):\gcd=500,对答案的贡献为 5003500^3

(i,j)=(6,7)(ai,aj)=(500,2000):gcd=500(i,j)=(6,7)\to (a_i,a_j)=(500,2000):\gcd=500,对答案的贡献为 500×2000×500500\times 2000\times 500

(i,j)=(7,6)(ai,aj)=(2000,500):gcd=500(i,j)=(7,6)\to (a_i,a_j)=(2000,500):\gcd=500,对答案的贡献为 500×2000×500500\times 2000\times 500

(i,j)=(7,7)(ai,aj)=(2000,2000):gcd=2000(i,j)=(7,7)\to (a_i,a_j)=(2000,2000):\gcd=2000,对答案的贡献为 200032000^3

(i,j)=(8,8)(ai,aj)=(9997,9997):gcd=9997(i,j)=(8,8)\to (a_i,a_j)=(9997,9997):\gcd=9997,对答案的贡献为 999739997^3

答案的总和为 1008225269973,对 998244353998244353 取模的结果是 996717796。

样例输入2

10
5000 2000 1000 2000 9999 1000 19998 29997 39996 49995

样例输出2

750967481

数据范围与约定

1n,ai1051\le n,a_i\le 10^5

南京师范大学第九届互联网创新创业科技节计算机程序设计大赛

未参加
状态
已结束
规则
ACM/ICPC
题目
13
开始于
2024-3-20 17:40
结束于
2024-3-20 20:10
持续时间
2.5 小时
主持人
参赛人数
133