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 个整数 Ai 表示CA的每个数。 1≤T≤50, 1≤N≤1000, 1≤Ai≤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;
}