Codeforces 803F Coprime Subsequences DP+GCD

点击打开链接

题意:n个数a[i],问有多少个子序列的gcd为1? (子序列下标不同就算不同) n,a[i]<=1e5 

令dp[i]:gcd为i的子序列个数 
i的倍数的子序列有2^c[i]-1个,这些子序列的gcd可能为(i,2i,3i...) 
则gcd为i的子序列个数:dp[i]=2^(c[i])-dp[2i]-dp[3i]...

#include <bits/stdc++.h>using namespace std;typedef long long ll;const ll mod=1e9+7;const int N=2e5+20;const int M=1e5;ll cnt[N],c[N];//c[i] i的倍数的个数 ll dp[N];//dp[i]:gcd为i的子序列个数 int n;ll pw2[N]; int main(){ while(cin>>n){memset(dp,0,sizeof(dp));memset(cnt,0,sizeof(cnt));int x;pw2[0]=1;for(int i=1;i<=n;i++){scanf("%I64d",&x);cnt[x]++;pw2[i]=(2ll*pw2[i-1])%mod;}for(int i=M;i>=1;i--){ll tot=0;for(int k=i;k<=M;k+=i)tot+=cnt[k];//tot:i的倍数个数 dp[i]=pw2[tot]-1;for(int k=i+i;k<=M;k+=i)dp[i]=(dp[i]-dp[k]+mod)%mod;}cout<<dp[1]<<endl;}return 0;}