#include <bits/stdc++.h>
using namespace std;
int dx[4] = {1, 0, -1, 0};
int dy[4] = {0, 1, 0, -1};
const int long long MOD = 1e9 + 7;
static const int MAX = 100;
static const int INF = (1 << 23);
template <class T>
T gcd(T a, T b) {
  return b ? gcd(b, a % b) : a;
}
template <class T>
T lcm(T a, T b) {
  return a / gcd(a, b) * b;
}
int long long p(int long long a, int long long b) {
  int long long k = 1;
  for (int i = 0; i < a; i++) {
    k *= b;
    k %= MOD;
  }
  return k;
}
int main() {
  int long long n;
  cin >> n;
  vector<int long long> v(n);
  for (int i = 0; i < n; i++) {
    cin >> v[i];
  }
  int long long ans = 0;
  for (int i = 0; i < 60; i++) {
    int long long one = 0;
    int long long zero = 0;
    for (int j = 0; j < n; j++) {
      if ((v[j] >> i) & 1ll) {
        one++;
      } else {
        zero++;
      }
    }
    int long long k = 0;
    if (one != 0) {
      k = zero * one;
      k = k * p(i, 2);
      k %= MOD;
      ans += k;
      ans %= MOD;
    }
  }
  cout << ans % MOD << endl;
  return 0;
}
