Codeforces 585E 莫比乌斯反演 + 容斥

       给定 $n(2\le n\le 5e5)$ 个数，$2\le a_i\le 1e7$。

       一次操作中，从数组 $a$ 中选定一个 $x$，再从剩下的数中选出一个 $\gcd >1$ 的非空集合 $Y$，使得所有选中的数的 $\gcd =1$，即 $\gcd (x,Y)=1$。

       问有多少种这样的方案，答案 $mod1e9+7$。

       我们尝试用公式表达出我们的所求：
$$\sum _{Y\subseteq U}\sum _{x⊈Y}[\gcd (Y)\ne 1][\gcd (x,Y)=1]$$

       当 $\gcd (Y)=1$ 时，$\gcd (x,Y)=1$ 必定满足，所以我们考虑用 $\gcd (x,Y)=1$ 的所有情况减去 $\gcd (Y)=1$ 的所有情况。

       即上式变为 :
$$\sum _{Y\subseteq U}\sum _{x⊈Y}[\gcd (x,Y)=1]-\sum _{Y\subseteq U}\sum _{x⊈Y}[\gcd (Y)=1]$$

       先考虑前半部分，当我们从一个 $\gcd =1$ 的集合 $Z(|Z|>1)$ 取出任何一个元素 $x$，都满足 $\gcd (x,Z-x)=1$，也就是说 $Z-x$ 就是我们要找的集合 $Y$。

       所以所有的集合产生的贡献就可以表述成:
$$\sum _{Y\subseteq U}\sum _{x⊈Y}[\gcd (x,Y)=1]=\sum _{|Z|>1,Z\subseteq U}[\gcd (Z)=1]|Z|$$

       对于一个集合 $Z$，贡献总共 $|Z|$ 种构造 $x$ 和 $Y$ 的方法。

       我们令 $m$ 为 $a_i$ 中为 $d$ 的倍数的数的数量。

$$\begin{array}{rl}\sum _{|Z|>1,Z\subseteq U}[\gcd (Z)=1]|Z|& =\sum _{|Z|>1,Z\subseteq U}\sum _{d|\gcd (Z)}\mu (d)|Z|\\ & =\sum _{d=1}^{max(a_i)}\mu (d)\sum _{|Z|=2}^m{C}_m^{|Z|}\ast |Z|\end{array}$$

       ${\sum }_{|Z|=2}^m{C}_m^{|Z|}\ast |Z|$ 可以表述成 ${\sum }_{i=2}^n{C}_n^i\ast i$，我们令 $f(n)={\sum }_{i=1}^n{C}_n^i\ast i$。

       则上述式子写成:
$$\sum _{d=1}^{max(a_i)}\mu (d)(f(m)-m)$$

       根据组合数学公式 $(1+x{)}^n={\sum }_{i=0}^n{C}_n^i\ast x^i$，两边求导可得 ${\sum }_{i=0}^n{C}_n^i\ast i=n\ast 2^{n-1}$ (而 $i=0$ 时无贡献可省去，即为上式)。

       则有 $f(x)=2f(x-1)+2^{x-1}$

       再考虑后半部分，道理类似，任意集合 $Y$，产生的贡献为 $n-|Y|$，$n-|Y|$ 即为 $x$ 的数量:

$$\begin{array}{rl}\sum _{Y\subseteq U}\sum _{x⊈Y}[\gcd (Y)=1]& =\sum _{Y\subseteq U}[\gcd (Y)=1](n-|Y|)\\ & =\sum _{Y\subseteq U}\sum _{d|\gcd (Y)}\mu (d)\ast (n-|Y|)\\ & =\sum _{d=1}^{max(a_i)}\mu (d)((n\ast \sum _{|Y|=1}^m{C}_m^i)-(\sum _{|Y|=1}^m{C}_m^i\ast |Y|))\\ & =\sum _{d=1}^{max(a_i)}\mu (d)(n\ast (2^m-1)-f(m))\end{array}$$

       合并起来可得:

$$\sum _{d=1}^{max(a_i)}\mu (d)(f(m)-m)-\sum _{d=1}^{max(a_i)}\mu (d)(n\ast (2^m-1)-f(m))=\sum _{d=1}^{max(a_i)}\mu (d)(2f(m)-m-n\ast (2^m-1))$$

       其中 $f(1)=1,f(x)=2f(x-1)+2^{x-1}$。

       这里有一个小 $trick$，我们不需要枚举 $a_i$ 的因子来计算贡献。而是通过统计每个 $a_i$ 的贡献，然后利用枚举 $d$ 的同时，统计它的倍数的贡献（调和级数复杂度 $nlogn$，如果枚举 $a_i$ 的因子则是 $n\sqrt{a_i}$ 复杂度）

#include<bits/stdc++.h>
using namespace std;
#define IOS ios::sync_with_stdio(false); cin.tie(0); cout.tie(0)
#define endl '\n'
typedef long long LL;
const int maxn = 1e7 + 1;
const LL mod = 1e9 + 7;

int bit[maxn], f[maxn], d[maxn], mu[maxn], prime[maxn], cnt;
bool judge[maxn];

inline void init(){
    mu[1] = 1;
    for(int i = 2; i < maxn; i++){
        if(!judge[i]){
            prime[++cnt] = i;
            mu[i] = -1;
        }
        for(int j = 1, v; j <= cnt && (v = i * prime[j]) < maxn; j++){
            judge[v] = true;
            if(!(i % prime[j])) break;
            mu[v] = -mu[i];
        }
    }
    f[1] = 1;
    bit[0] = 1;
    bit[1] = 2;
    for(int i = 2; i < maxn; i++){
        f[i] = (2LL * f[i - 1] + bit[i - 1]) % mod;
        bit[i] = bit[i - 1] * 2LL % mod;
    }
}

int main(){
    IOS;
    init();
    int n, a;
    cin >> n;
    for(int i = 1; i <= n; i++){
        cin >> a;
        d[a]++;
    }
    int ans = mu[1] * (2LL * f[n] % mod - n - n * (bit[n] - 1LL) % mod) % mod;
    for(int i = 2; i < maxn; i++){
        if(mu[i]){
            int tol = 0;
            for(int j = i; j < maxn; j += i) tol += d[j];
            ans = (ans + mu[i] * (2LL * f[tol] % mod - tol - n * (bit[tol] - 1LL) % mod)) % mod;
        }
    }
    cout << (ans + mod) % mod << endl;
}