import numpy as np
n=int(input())
a=np.array(list(map(lambda x:list(map(int,list(('0'*62+bin(int(x))[2:])[-62:]))),input().split())),np.int64)
mod=pow(10,9)+7
ans=0
for i in range(61):# 2^iの位
  bit1=np.sum(a[:,-i-1])
  b=pow(2,i,mod)
  b*=bit1
  b%=mod
  b*=n-bit1
  b%=mod
  ans+=b
  ans%=mod
print(ans)