import math
N = int(input())
A = list(map(int,input().split()))
b = [[0]*60 for i in range(2)]
max_a = math.ceil(math.log(max(A),2))
for i in range(N):
    x = A[i]
    j = 0
    for j in range(max_a):
        b[x % 2][j] += 1
        x //= 2

x = 0
for i in range(max_a):
    x += b[0][i] * b[1][i] * 2**i
    x %= 10 ** 9 + 7

print (x)