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