Coprime Subsequences(数论gcd)

题目链接：https://ac.nowcoder.com/acm/problem/112055

题目大意

给你n个正整数，求gcd=1的子序列的个数（mod 1e9+7）

思路

gcd = i的序列是由一些i的倍数的数构成的，假设a[i]表示给的数中i的倍数的个数，那么这些数一共有(2^a[i] - 1)个组合，但是这些组合并不一定全都满足gcd = i，还要筛去gcd = i*2, i*3....这么多，剩下的才是满足条件的数量。而gcd = i*2 的组合一共有(2^a[i*2] - 1)这么多，所以只需要求解a数组即可，这里可以对每个数直接暴力O()存桶。

ac代码

#include<bits/stdc++.h>
using namespace std;
#define io cin.tie(0);ios::sync_with_stdio(false);
#define ok(x, y) x >= 1 && x <= n && y >= 1 && y <= m
#define debug(x) cout<<#x<<"="<<x<<endl
#define lowbit(x) x&(-x)
#define pii pair<int,int>
#define mk make_pair
#define ll long long
#define ull unsigned long long
#define lb long double
#define rs p<<1|1
#define ls p<<1
#define eps 1e-6
#define pi acos(-1)
const int maxn = 1e5 + 5;
const int mod = 1e9+7;
const int inf = 0x3f3f3f3f;
const int bas = 2333;
inline ll read(){
    ll p=0,f=1;char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){p=(p<<1)+(p<<3)+(c^48),c=getchar();}
    return f*p;
}
void print(__int128 x){
    if(x<0) {putchar('-'); x=-x;}
    if (x>9) print(x/10);
    putchar('0'+x%10);
}

ll a[maxn], b[maxn];
ll _pow(ll a, ll b){
    ll ans = 1;
    while(b){
        if(b & 1) ans = ans * a % mod;
        a=  a * a % mod;
        b >>= 1;
    }
    return ans;
}
void solve(){
    int n; cin >> n;
    for(int i = 1; i <= n; i ++){
        int x; cin >> x;
        for(int j = 1; j * j <= x; j ++){
            if(x % j == 0){ //存两个因子
                a[j] ++; 
                if(x / j == j) continue;
                a[x/j]++;
            }
        }
    }
    for(int i = 1; i < maxn; i ++){ //a[i]个数一共有(2^a[i] - 1)个组合
        b[i] = (_pow(2, a[i]) - 1 + mod) % mod;
    }
    for(int i = maxn - 1; i >= 1; i --){ //这里是从大到小，因为求解b[i]就是要减去b[i*2],减去b[i*3]......，类似于递归的过程
        for(int j = i + i; j < maxn; j += i){
            b[i] = (b[i] - b[j] + mod) % mod;
        }
    }
    cout << b[1] << endl;
}

int main(){
    io;
    solve();
    return 0;
}