CodeForces839D[莫比乌斯反演] Codeforces Round #428 (Div. 2)

/*CodeForces839D[莫比乌斯反演]*/
#include <bits/stdc++.h>
typedef long long LL;
const LL MOD = 1000000007LL;
using namespace std;
int n, maxa = , mina = 0x3f3f3f3f;
LL m[], F[], f[];
int sum[], a[];
int  prime[], vis[], mu[];
void mo_init(int N) {
    memset(vis, , sizeof(vis));
    mu[] = ;
    int cnt = ;
    for (int i = ; i < N; i++)
    {
        if (!vis[i])
        {
            prime[cnt++] = i;
            mu[i] = -;
        }
        for (int j = ; j < cnt && i * prime[j] < N; j++)
        {
            vis[i * prime[j]] = ;
            if (i % prime[j]) mu[i * prime[j]] = -mu[i];
            else
            {
                mu[i * prime[j]] = ;
                break;
            }
        }
    }
}
void solve() {
    for (int i = ; i < n; i++) {
        for (int j = ; j * j <= a[i]; j++) {
            if (a[i] % j == ) {
                if (j * j == a[i]) {
                    sum[j]++;
                }
                else {
                    sum[j]++, sum[a[i] / j]++;
                }
            }
        }
    }
    LL ans = ;
    for (int i = ; i <= maxa; i++) {
        F[i] = sum[i] * m[sum[i] - ] % MOD;
    }
    for (int i = ; i <= maxa; i++) {
        for (int j = i; j <= maxa; j += i) {
            f[i] = (f[i] + mu[j / i] * F[j] % MOD) % MOD;
        }
        ans = (ans + f[i] * i % MOD) % MOD;
    }
    printf("%lld\n", ans);
}
void init() {
    mo_init();
    m[] = 1LL;
    for (int i = ; i <= n; i++) {
        m[i] = m[i - ] * 2LL % MOD;
    }
}
int main() {
    scanf("%d", &n);
    for (int i = ; i < n; i++) {
        scanf("%d", &a[i]);
        maxa = max(maxa, a[i]);
        mina = min(mina, a[i]);
    } init(), solve();
    return ;
}