【51Nod1585】买邮票送邮票

商店里面有n张邮票，现在去买一张，然后老板会送若干张（至少一张）邮票。如果老板送的邮票的面值的最大公约数不是1，并且老板送的邮票和我们购买的邮票的面值最大公约数是1，那么就是一组好的邮票组合。问有多少种好的邮票组合。
样例解释：
· 买第1张，送第2张；
· 买第3张，送第2张；
· 买第2张，送第1张；
· 买第2张，送第3张；
· 买第2张，送第1和第3张；

Input
单组测试数据。
第一行有一个整数n (2≤n≤5*10^5)，表示商店中邮票的张数。
第二行有n个整数 a1,a2,…,an (2≤ai≤10^7)，表示第i张邮票的面值。
Output
由于结果比较大，对 10^9+7取余后输出即可。
Input示例
样例输入1
3
2 3 2
Output示例
样例输出1
5

题解
枚举gcd，则其倍数都是可以选的，发现有重复，用莫比乌斯容斥即可。

代码

#include<bits/stdc++.h>
typedef long long ll;
const int N=10000000;
const int mod=1000000007;
const double eps=0.00000001;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,tot,p[5000000],mu[N+5],ans,fac[500005],a[N+5];
bool flag[N+5];
void getmu()
{
    mu[1]=1;
    for (int i=2;i<=N;i++)
    {
        if (!flag[i]) p[++tot]=i,mu[i]=-1;
        for (int j=1;j<=tot&&p[j]*i<=N;j++)
        {
            flag[p[j]*i]=1;
            if (i%p[j]==0) break;
            mu[i*p[j]]=-mu[i];
        }
    }
}
int main()
{
    n=read();getmu();fac[0]=1;
    for (int i=1;i<=n;i++)fac[i]=(fac[i-1]<<1)%mod;
    for (int i=1;i<=n;i++) a[read()]++;
    for (int i=2;i<=N;i++)if (mu[i])
    {
        int sum=0;
        for (int j=1;j*i<=N;j++) sum+=a[j*i];
        (ans+=((ll)fac[sum]-1)*(n-sum)%mod*(mu[i]==1?-1:1))%=mod;
    }
    cout<<(ans+mod)%mod;
    return 0;
}