前置知识:
算法一
对于 30 % 30 \% 30% 的数据, n ≤ 3 × 1 0 3 n \leq 3 \times 10^3 n≤3×103.
直接模拟就好了。
时间复杂度: O ( n 2 log n ) O(n^2 \log n) O(n2logn). 实际得分: 30 p t s 30pts 30pts.
for(i from 1 to n)
for(j from 1 to n)
s+=gcd(i,j)
print(s)
算法二
对于 60 % 60 \% 60% 的数据, 7000 ≤ n ≤ 7100 7000 \leq n \leq 7100 7000≤n≤7100.
注意到 O ( n 2 log n ) O(n^2 \log n) O(n2logn) 无法通过了。但是,我们可以计算出 n = 7000 n = 7000 n=7000 的答案为 275797760 275797760 275797760. 然后,对于 7001 7001 7001 ~ n n n 区间的答案,实际上就是 两倍的该区间与 1 ≤ i ≤ 7000 1 \leq i \leq 7000 1≤i≤7000 的两两 gcd \gcd gcd 之和,再加上自己和自己区间答案。
???
就比方说, a a a 和 b b b 表示区间,定义 a × b a \times b a×b 为 a a a 区间与 b b b 区间两两 gcd \gcd gcd 和。那么, x = a + b x = a + b x=a+b,则 x 2 = ( a + b ) 2 = a 2 + 2 × a × b + b 2 x^2 = (a+b)^2 = a^2 + 2 \times a \times b + b^2 x2=(a+b)2=a2+2×a×b+b2.
而 a 2 a^2 a2 的答案被提前算出,然后 b 2 b^2 b2 区间长度 ≤ 100 \leq 100 ≤100 结束, a × b a \times b a×b 计算也只需要 7000 × 100 = 7 × 1 0 5 7000 \times 100 = 7 \times 10^5 </