题意:
题解:
首先要了解的一个概念就是:枚举
[
1
,
2
x
]
[1,2^x]
[1,2x]的所有状态的子状态的时间复杂度是
3
x
3^x
3x左右。那么我们这道题肯定不能for一遍n然后枚举所有当前点的补集。然后 发现a的数据范围是2e5,那么就可以先记录所有状态有多少种,然后再枚举所有状态的同时枚举补集进行转移。最后0要特别的运算一下因为它每一个取不取都不会对结果造成影响。
#include<bits/stdc++.h>
using namespace std;
const int maxn=(1<<18)+5;
bitset<maxn>notprime;
int phi[maxn],prime[maxn],cnt=0;
void pre(){
phi[1]=1;
for(int i=2;i<=maxn;i++){
if(!notprime[i]){
prime[++cnt]=i;
phi[i]=i-1;
}
for(int j=1;j<=cnt&&prime[j]*i<=maxn;j++){
notprime[i*prime[j]]=1;
if(i%prime[j]==0){
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
else phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
}
const int N=(1<<18)+5;
int a;
#define ll long long
ll dp[N],num[N];
const ll mod=1e9+7;
int main(){
pre();
int n;
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d",&a),num[a]++;
dp[0]=1;
for(int i=1;i<(1<<18);i++){
for(int j=i;;j=(j-1)&i){
int res=i^j;
if(res>j)break;
dp[i]=(dp[i]+dp[res]*num[j])%mod;
}
}
ll ans=1;
for(int i=1;i<N-1;i++)
ans=(ans+dp[i]*phi[i+1])%mod;
for(int i=1;i<=num[0];i++)ans=ans*2%mod;
printf("%lld\n",ans);
return 0;
}