[cf803F] Coprime Subsequences(组合数学,容斥)

题目链接:http://codeforces.com/contest/803/problem/F

题意:n个数,统计所有满足gcd为1的子序列的个数。

统计各个数字出现的次数,然后容斥。

从10000扫到1,枚举所有倍数,把x的倍数出现的次数和统计出来,然后再统计x的倍数出现的次数和-1(减掉空集),然后作差。

倒着扫是因为,假如先前遇到了同倍数的,要减掉,如果正着扫就会容斥不到。

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 
 4 typedef long long LL;
 5 const int maxn = 100100;
 6 const LL mod = (LL)1e9+7;
 7 int n, a;
 8 LL vis[maxn], tot[maxn];
 9 
10 LL mul(LL x, LL n) {
11     LL tot = 1;
12     for(; n; x=(x*x)%mod, n>>=1) if(n & 1) tot = (tot * x) % mod;
13     return tot;
14 }
15 
16 int main() {
17     // freopen("in", "r", stdin);
18     while(~scanf("%d",&n)) {
19         memset(vis, 0, sizeof(vis));
20         memset(tot, 0, sizeof(tot));
21         for(int i = 1; i <= n; i++) {
22             scanf("%d", &a);
23             vis[a]++;
24         }
25         for(int i = 100000; i >= 1; i--) {
26             LL X = 0, Y = 0;
27             for(int j = 1; i * j <= 100000; j++) {
28                 X = (X + vis[i*j]) % mod;
29                 Y = (Y + tot[i*j]) % mod;
30             }
31             tot[i] = (mul(2, X) - 1 - Y + mod) % mod;
32         }
33         printf("%I64d\n", (tot[1]+mod)%mod);
34     }
35     return 0;
36 }

 

posted @ 2017-05-11 22:05  Kirai  阅读(201)  评论(0编辑  收藏  举报