N=int(input())
A=list(map(int,input().split()))
l=len(format(max(A),'b'))
Z=[0]*l
for a in A:
  a=format(a,'b').zfill(l)
  for i,ai in enumerate(a):
    if ai=='0':
      Z[i]+=1
x=1
ans=0
for z in Z[::-1]:
  ans+=x*(z*(N-z))
  x*=2
print(ans%(10**9+7))