2
1 31
 2
1 0
 10
450661 128600 993228 725823 293549 33490 843121 903634 556169 448234
 1
1
 5
1 3 5 7 9
 6
524 529 5249 524 529 529
 3
1 2 3
 2
0 1
 55
0 1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 9 2 3 4 5 6 7 8 9 3 4 5 6 7 8 9 4 5 6 7 8 9 5 6 7 8 9 6 7 8 9 7 8 9 8 9 9
 16
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
 3
128 1024 2048
 10
3 3 3 5 5 3 3 3 3 5
 1
0
 2
1 2
 2
0 0
 2
1 1

0
 2
 632
 0
 0
 24
 2
 2
 621247139
 64594
 4
 0
 1
 1
 3
 0

#include <bits/stdc++.h>
using namespace std;
template <typename T>
ostream& operator<<(ostream& out, vector<T> const& a) {
  out << "[";
  for (int i = 0; i < int((a).size()); ++i) {
    if (i != 0) {
      out << ", ";
    }
    out << a[i];
  }
  out << "]";
  return out;
}
template <typename T1, typename T2>
ostream& operator<<(ostream& out, pair<T1, T2> const& p) {
  out << "(" << p.first << ", " << p.second << ")";
  return out;
}
const int N = 1e6;
struct Data {
  int n;
  int a[N];
  bool read() {
    if (scanf("%d", &n) != 1) {
      return false;
    }
    for (int i = 0; i < n; ++i) {
      scanf("%d", &a[i]);
    }
    return true;
  }
  long long ans;
  void write() { printf("%d\n", (int)ans); }
  virtual void solve() {}
  virtual void clear() { *this = Data(); }
  void gen() {
    n = 1e6;
    for (int i = 0; i < n; ++i) {
      a[i] = i;
    }
  }
};
const int mod = 1e9 + 7;
int sum(int x, int y) {
  x += y;
  if (x >= mod) {
    x -= mod;
  }
  return x;
}
const int L = 20;
const int C = 1 << L;
int cnt[C + 1], cnt2[C + 1];
int pw2[N + 1];
struct Solution : Data {
  void solve() {
    int mx = 0;
    for (int i = 0; i < n; ++i) {
      ++cnt[a[i]];
    }
    for (int i = 0; i < L; ++i) {
      for (int j = 0; j < C; ++j) {
        if (j & (1 << i)) {
          cnt[j ^ (1 << i)] += cnt[j];
        }
      }
    }
    pw2[0] = 1;
    for (int i = 1; i <= n; ++i) {
      pw2[i] = sum(pw2[i - 1], pw2[i - 1]);
    }
    ans = 0;
    for (int i = 1; i < C; ++i) {
      long long e = -1;
      for (int j = 0; j < L; ++j) {
        if (i & (1 << j)) {
          e *= -1;
        }
      }
      ans += e * (pw2[cnt[i]] - 1);
      ans %= mod;
      ans += mod;
      ans %= mod;
    };
    ans = (pw2[n] - 1 - ans) % mod;
    ans += mod;
    ans %= mod;
  }
  Solution() {}
  Solution(Data d) : Data(d) {}
  void clear() {
    *this = Solution();
    memset(cnt, 0, sizeof(cnt));
    memset(cnt2, 0, sizeof(cnt2));
  }
};
Solution sol;
void test() {
  const int C = 1e6;
  long long cnt = 0;
  for (int i = 0; i <= C; ++i) {
    for (int j = i; j > 0; j = (j - 1) & i) {
      ++cnt;
    }
  };
}
int main() {
  sol.read();
  sol.solve();
  sol.write();
  return 0;
}
