1. 程式人生 > >[51Nod 1238] 最小公倍數之和 V3【數論】【杜教篩】【未完成】

[51Nod 1238] 最小公倍數之和 V3【數論】【杜教篩】【未完成】

Description

i=1nj=1nlcm(i,j)\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{n} lcm(i,j) n1010n\leq 10^{10}

Solution

lcmlcm拆開 原式化為 i=1nj=1nijgcd(i,j)\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{n} {ij\over \gcd(i,j)}d=gcd(i,j),i=i/d,j=j/dd=\gcd(i,j),i=i/d,j=j/d

方便起見以下的除法預設為下取整 =d=1ndi=1ndj=1ndij[gcd(i,j)=1]=\sum\limits_{d=1}^{n}d\sum\limits_{i=1}^{n\over d}\sum\limits_{j=1}^{n\over d} ij[\gcd(i,j)=1]

考慮對稱性,有 =d=1nd(2(i=1ndij=1ij[gcd(i,j)=1])1)=\sum\limits_{d=1}^{n}d*\left(2\left(\sum\limits_{i=1}^{n\over d}i\sum\limits_{j=1}^{i} j[\gcd(i,j)=1]\right)-1\right)

(因為1被計算了2次,因此最後要-1)

我們知道,小於等於n且與n互質的數的和是φ(n)n2φ(n)*n\over 2(考慮一個小於等於n且與n互質的數i,顯然n-i也與n互質,這樣總共有φ(n)2φ(n)\over 2組) 1除外,小於等於1且與1互質的數的和是1,因此還要加上一個121\over 2 因此又有 =d=1nd(2(i=1ndφ(i)i22)1+122)=\sum\limits_{d=1}^{n}d*\left(2\left(\sum\limits_{i=1}^{n\over d}{φ(i)*i^2\over 2}\right)-1+{1\over 2}*2\right)

=d=1ndi=1ndφ(i)i2=\sum\limits_{d=1}^{n}d*\sum\limits_{i=1}^{n\over d}φ(i)*i^2