luogu P6570 [NOI Online #3 提高组] 优秀子序列

题面传送门
现在来看这种东西真的很屑啊当时我是真的菜。
首先这个欧拉函数很难搞，我们考虑设\(F_i\)为和为\(i\)的方案数。
然后发现其实这个也是所有数或的和为\(i\)的方案数。
所以枚举一下子集不就可以dp了。
然后因为要保证顺序所以我们强制加入的数最高位递增。
时间复杂度\(O(3^{loga}+n)\)，注意特判0
code:

#include<bits/stdc++.h>
#define I inline
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define re register
#define ll long long
#define db double
#define N 1000000
#define W (1<<18)
#define mod 1000000007
#define eps (1e-5)
#define U unsigned int
#define it iterator
#define Gc() getchar() 
#define Me(x,y) memset(x,y,sizeof(x))
#define d(x,y) (n*(x-1)+(y))
using namespace std;
int n,A[N+5],fl[W+5],phi[W+5],pr[W+5],ph,F[N+5],ToT;ll dp[W+5],Ans;
int main(){
	freopen("1.in","r",stdin);
	re int i,j;scanf("%d",&n);for(i=1;i<=n;i++) scanf("%d",&A[i]),A[i]?(F[A[i]]++,dp[A[i]]++):(ToT++);
	for(i=2;i<=W;i++){
		!fl[i]&&(pr[++ph]=i,phi[i]=i-1);for(j=1;j<=ph&&i*pr[j]<=W;j++){
			fl[pr[j]*i]=1;if(i%pr[j]==0) {phi[i*pr[j]]=phi[i]*pr[j];break;}phi[i*pr[j]]=phi[i]*phi[pr[j]];
		}
	} 
	for(i=0;i<W;i++){
		for(j=i&(i-1);j>(i^j);j=i&(j-1)) dp[i]=(dp[i]+F[j]*dp[i^j])%mod;
	}for(i=0;i<W;i++) Ans=(Ans+dp[i]*phi[i+1])%mod;Ans++;while(ToT--) Ans=Ans*2%mod;printf("%lld\n",Ans);
}