记个数论题(15.6.25Test )

求$\Sigma _{i=1} ^{n}$$\Sigma_{j=1}^{m}$$(gcd(i,j))^k$

$\Sigma _{i=1} ^{n}$$\Sigma_{j=1}^{m}$$(gcd(i,j))^k$
$=\Sigma_{i=1}^{n}$$\Sigma_{j=1}^{m}$$\Sigma_{p/i,p/j}p^k$$e(gcd(\frac{i}{p},\frac{j}{p}))$
$=\Sigma_{p=1}^{n}p^k$$\Sigma_{i=1}^{\frac{n}{p}}$$\Sigma_{j=1}^{\frac{m}{p}}$$e(gcd(i,j))$
$=\Sigma_{p=1}^{n}p^k$$\Sigma_{i=1}^{\frac{n}{p}}$$\Sigma_{j=1}^{\frac{m}{p}}$$\Sigma_{d/i,d/j}μ(n)$
$=\Sigma_{p=1}^{n}p^k$$\Sigma_{d=1}^{\frac{n}{p}}μ(n)$$⌊$$\frac{n}{pd}$$⌋$$⌊$$\frac{n}{pd}$$⌋$

因为对于n，除p只有$\sqrt{n}$种，同理除d也是，所以我们维护前缀和，对于相同情况的p和相同情况的d同时统计，效率就变成$O(n)$的了。