CodeForces 585D Present for Vitalik the Philatelist（容斥原理+gcd）

题意：给出一个数字序列，求有多少组满足原来子集不互质，先加一个元素后新集合互质的集合数量。

思路：首先解决子问题，如何求出一个集合中互质的子集个数，首先求出I的倍数的元素的个数，记为cnt[I]，

然后运用容斥原理来做，当I为某些不重复的素数的组合时，当质因数的个数为奇时，减去2^cnt[I]-1，反之则加。

如果I的质因数分解的结果中有重复的质因数的情况，那么不加也不减，因为I已经在之前被考虑过了。

然后求出来互质的子集的数量后，固定每个元素，求出剩下的不互质的子集中与这个元素互质的集合数量就是答案，

因为这道题卡空间，即使预处理出每个数的所有质因数也会被卡掉，所以可以用筛法先求出每个数的最小的质因数，

求每个数的所有不同质因数组合成的因数时，先求出每个数的所有质因数

然后用状压求出每个数的所有不同质因数组合成的因数。

#include <bits/stdc++.h>
#define eps 1e-6
#define LL long long
#define pii pair<int, int>
#define pb push_back
#define mp make_pair
//#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;

const int MAXN = 1000000;
const int MAXM = 10001000;
const int MOD = 1e9+7;
//const int INF = 0x3f3f3f3f;
int n;
int a[MAXN], cnt[MAXM], mul[1000], c[MAXM];
int sieve[MAXM];
LL po[MAXN];
void init(int m) {
	for (int i = 1; i <= m; i++) sieve[i] = i;
	for (int i = 2; i*i <= m; i++) if (sieve[i] == i) {
		for (int j = i*i; j <= m; j+=i) sieve[j] = i; 
	}
	po[0] = 1;
	for (int i = 1; i <= n+10; i++) { 
		po[i] = po[i-1]<<1;
		if (po[i] >= MOD) po[i] -= MOD;
	}
}
vector<pii> get_prime(int n) {
	vector<int> prime;
	while (n > 1) {
		int p = sieve[n];
		while (n%p == 0) {
			n /= p;
		}
		prime.push_back(p);
	}
	vector<pii> ans;
	int sz = prime.size();
	int maxS = 1 << sz;
	mul[0] = 1;
	for (int j = 0; j < sz; j++) 
		mul[1<<j] = prime[j];
	for (int j = 1; j < maxS; j++) {
		mul[j] = mul[j&-j] * mul[j^(j&-j)];
		ans.push_back(mp(mul[j], __builtin_popcount(j)));
	}
	return ans;
}


int main()
{
	//freopen("input.txt", "r", stdin);
	scanf("%d", &n);
	int maxA = 0;
	for (int i = 1; i <= n; i++) {
		scanf("%d", &a[i]);
		maxA = max(a[i], maxA);
	} 
	init(maxA+10);
	for (int i = 1; i <= n; i++) {
		vector<pii> d = get_prime(a[i]);
		for (int j = 0; j < d.size(); j++) {
			cnt[d[j].first]++;
			c[d[j].first] = d[j].second;
		}
	}
	LL cop = po[n] - 1;
	for (int i = 2; i <= maxA; i++) {
		if (c[i]&1)
			cop = (cop-po[cnt[i]]+1) % MOD;
		else if (c[i] != 0)
			cop = (cop+po[cnt[i]]-1) % MOD;
		//cout << i << " " << d.size() << " " << cop << endl;
	}
	LL noCop = po[n] - 1 - cop;
	//cout << noCop << endl;
	LL ans = 0;
	for (int i = 1; i <= n; i++) {
		vector<pii> d = get_prime(a[i]);
		LL tmp = noCop;
		for (int j = 0; j < d.size(); j++) {
			if (d[j].second&1)
				tmp = (tmp-po[cnt[d[j].first]]+1) % MOD;
			else 
				tmp = (tmp+po[cnt[d[j].first]]-1) % MOD;
		}
		ans = (ans + tmp) % MOD;
		//cout << ans << endl;
	}
	//cout << sieve[9] << endl;
	ans = (ans+MOD) % MOD;
	printf("%I64d", ans);
	return 0;
}