import sys

input = lambda: sys.stdin.readline().rstrip()
mod = 10**9 + 7


def solve():
    N = int(input())
    A = list(map(int, input().split()))

    a_max = max(A)
    a_len = len(A)
    keta = 0

    while True:
        a_max = a_max // 2
        keta += 1
        if a_max == 0:
            break

    # print(keta)
    t = [[0 for _ in range(keta)] for _ in range(2)]

    for k in range(keta):
        for a in A:
            a = a & (2**k)
            a = a >> k
            t[a][k] += 1

    # print(t)

    ans = 0
    for k in range(keta):
        for a in A:
            tmp = a & (2**k)
            tmp = tmp >> k
            tmp = not tmp
            ans += t[tmp][k] * (2**k)
            ans %= mod

    ans = (ans // 2) % mod
    print(ans)


if __name__ == '__main__':
    solve()
