#include <bits/stdc++.h>
using namespace std;
int dx[] = {1, -1, 0, 0, 1, -1, 1, -1};
int dy[] = {0, 0, 1, -1, 1, -1, -1, 1};
const long double PI = acos(-1);
int r() {
  int x;
  cin >> x;
  return x;
}
void solve() {
  int n;
  cin >> n;
  long long a[n];
  for (int64_t i = (0); i < (n); i++) cin >> a[i];
  long long ans = 0;
  for (int64_t i = (0); i < (60); i++) {
    int num = 0;
    for (int64_t j = (0); j < (n); j++)
      if (a[j] >> i & 1) num++;
    num = num * (n - num) % (1000000007);
    for (int64_t j = (0); j < (i); j++) num = num * 2 % (1000000007);
    ans += num;
    ans %= (1000000007);
  }
  cout << ans << endl;
}
int main() {
  cin.tie(0);
  ios::sync_with_stdio(false);
  solve();
  return 0;
}
