[NOI Online #3 提高组]优秀子序列欧拉函数，SOS DP

This way

题意：

题解：

首先要了解的一个概念就是：枚举$[1,2^x]$的所有状态的子状态的时间复杂度是$3^x$左右。那么我们这道题肯定不能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;
}