HDU 5656：CA Loves GCD

CA Loves GCD

 

 Accepts: 135

 

 Submissions: 586

 Time Limit: 6000/3000 MS (Java/Others)

 

 Memory Limit: 262144/262144 K (Java/Others)

问题描述

CA喜欢是一个热爱党和人民的优秀同♂志，所以他也非常喜欢GCD（请在输入法中输入GCD得到CA喜欢GCD的原因）。
现在他有N个不同的数，每次他会从中选出若干个（至少一个数），求出所有数的GCD然后放回去。
为了使自己不会无聊，CA会把每种不同的选法都选一遍，CA想知道他得到的所有GCD的和是多少。
我们认为两种选法不同，当且仅当有一个数在其中一种选法中被选中了，而在另外一种选法中没有被选中。

输入描述

第一行 $T$，表示有 $T$ 组数据。
接下来 $T$ 组数据，每组数据第一行一个整数 $N$，表示CA的数的个数，接下来一行 $N$ 个整数 A_iA​i​​ 表示CA的每个数。
1 \le T \le 50,~1 \le N \le 1000,~1 \le A_i \le 10001≤T≤50, 1≤N≤1000, 1≤A​i​​≤1000

输出描述

对于每组数据输出一行一个整数表示CA所有的选法的GCD的和对 $100000007$ 取模的结果。

输入样例

2
2
2 4
3
1 2 3

输出样例

8
10

下午刚刚做浙江大学校赛的题，被这个容斥折磨的不行。结果晚上看BC发现出了一道原题。。。

dp[i]表示gcd为i的集合数，然后就查有多少个含有该因子的元素个数，减去之前算过的gcd个数dp[j](j=i*k，因为想一想，这些个数组成在一起的集合，不是比自己大的那些因子，就剩下i了)。

代码：

#pragma warning(disable:4996)
#include <iostream>
#include <functional>
#include <algorithm>
#include <cstring>
#include <vector>
#include <string>
#include <cstdio>
#include <cmath>
#include <queue>
#include <stack>
#include <deque>
#include <set>
#include <map>
using namespace std;
typedef long long ll;

#define INF 0x333f3f3f
#define repp(i, n, m) for (int i = n; i <= m; i++)
#define rep(i, n, m) for (int i = n; i < m; i++)
#define sa(n) scanf("%d", &(n))

const ll mod = 100000007;
const int maxn = 1005;
const double PI = acos(-1.0);

ll n;
ll dp[maxn], cnt[maxn], val[maxn];

ll po(ll a, ll k)
{
	ll ans = 1;
	while (k)
	{
		if (k & 1)
			ans = ans*a%mod;
		a = a*a%mod;
		k = k >> 1;
	}
	return ans;
}

void solve()
{
	ll i, j;
	ll maxx = 0;
	scanf("%lld", &n);
	
	memset(dp, 0, sizeof(dp));
	memset(cnt, 0, sizeof(cnt));

	for (i = 1; i <= n; i++)
	{
		scanf("%lld", &val[i]);
		cnt[val[i]]++;
		maxx = max(maxx, val[i]);
	}
	ll ans = 0;
	for (i = maxx; i >= 1; i--)
	{
		ll c = 0;
		for (j = i; j <= maxx; j += i)
		{
			c += cnt[j];
			dp[i] = (dp[i] - dp[j] + mod) % mod;
		}
		dp[i] = ((dp[i] + po(2, c) - 1) % mod + mod) % mod;
		ans = (ans + dp[i] * i) % mod;
	}
	printf("%lld\n", ans);
}

int main()
{
#ifndef ONLINE_JUDGE  
	freopen("i.txt", "r", stdin);
	freopen("o.txt", "w", stdout);
#endif

	int t;
	scanf("%d", &t);
	while (t--)
	{
		solve();
	}
	return 0;
}