题面
解法
其实很久之前就想写这道题了……
- 首先我们先考虑一下 d ( i j ) d(ij) d(ij)怎么处理。
- 有一个结论: d ( i j ) = ∑ x ∣ i ∑ y ∣ j [ g c d ( x , y ) = = 1 ] d(ij)=\sum_{x|i}\sum_{y|j}[gcd(x,y)==1] d(ij)=∑x∣i∑y∣j[gcd(x,y)==1],具体证明可以看popoqqq大爷的证明:链接。
- 然后我们就可以把式子写成这样: ∑ i = 1 n ∑ j = 1 m ∑ x ∣ i ∑ y ∣ j [ g c d ( x , y ) = = 1 ] \sum_{i=1}^n\sum_{j=1}^m\sum_{x|i}\sum_{y|j}[gcd(x,y)==1] ∑i=1n∑j=1m∑x∣i∑y∣j[gcd(x,y)==1]
- 感觉后面的整除并不是那么方便,不妨把 x , y x,y x,y提前到前面枚举,就变成: ∑ x = 1 n ∑ y = 1 m ⌊ n x ⌋ ⌊ m y ⌋ [ g c d ( x , y ) = = 1 ] \sum_{x=1}^n\sum_{y=1}^m\lfloor\frac{n}{x}\rfloor\lfloor\frac{m}{y}\rfloor[gcd(x,y)==1] ∑x=1n∑