[洛谷1390]公约数的和

题目大意：

求$displaystyle{sum_{1leq i<jleq n}}gcd(i,j)$的值。

思路：

由于数据水，可以直接用动态规划做。
用$f_k$表示在n以内以$k$为$gcd$的整数对个数，那么可以得到状态转移方程：
$f_i=lfloorfrac{n}{i} floor-displaystyle{sum_{j=2}^{lfloorfrac{n}{i} floor}}f_{ij}$
因为要减去$gcd(d,d)=d$的和$gcd(i,j)=gcd(j,i)$重复的，答案为：
$frac{displaystyle{sum_{i=1}^n}f_i-frac{n imes n+1}{2}}{2}$

#include<cstdio>
const long long N=2000001;
long long f[N]={0};
inline long long sqr(const long long x) {
    return x*x;
}
int main() {
    long long n;
    scanf("%lld",&n);
    long long ans=0;
    for(long long i=n;i;i--) {
        f[i]=sqr(n/i);
        for(long long j=2;j<=n/i;j++) {
            f[i]-=f[i*j];
        }
        ans+=f[i]*i;
    }
    printf("%lld
",(ans-n*(n+1)/2)/2);
    return 0;
}