这种题肯定是考虑按位计算的跑不了

因为选三个
要么3个0、2个0和1个1、1个0和2个1、3个1.
容易发现0 0 1和1 1 1才有贡献

那么遍历二进制的64位，根据组合数选取方案数乘上该二进制位的数即可。。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
ll a[1<<18],cnt[65];
ll f[65];
int main(){

    /*cout<<(0^0^0)<<endl;
    cout<<(0^0^1)<<endl;
    cout<<(0^1^1)<<endl;
    cout<<(1^1^1)<<endl;*/
    f[0]=1;
    for(int i=1;i<=63;i++) f[i]=(f[i-1]*2)%mod;
    int n;cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        for(int j=0;j<64;j++){
            cnt[j]+=a[i]>>j&1;
        }
    }
    ll ans=0;
    for(ll i=0;i<64;i++){
        ll one=cnt[i],zero=n-cnt[i];
        if(zero>=2 && one>=1){ //0 0 1
            ll k=zero*(zero-1)/2%mod*one%mod;
            (ans+=k*f[i]%mod)%=mod;
        }
        /// 2e5 4e10  8e15
        if(one>=3){//1 1 1
            ll k=one*(one-1)*(one-2)/6%mod;
            (ans+=k*f[i]%mod)%=mod;
        }
    }
    cout<<ans;

    return 0;
}