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;
long long N, sol;
long long A[1100100], tmp[1100100], F[1100100];
long long pw[1100100], freq[1100100];
void calc(long long b, long long pF) {
  if (b == 0) {
    F[pF] = freq[pF];
    F[pF + 1] = freq[pF + 1];
    F[pF] += F[pF + 1];
    return;
  }
  long long n = (1 << (b));
  calc(b - 1, pF);
  calc(b - 1, pF + n);
  for (long long I = 0; I < n; I++) F[pF + I] += F[pF + I + n];
}
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  pw[0] = 1;
  for (long long I = 1; I < 1100100; I++) pw[I] = (pw[I - 1] * 2) % 1000000007;
  cin >> N;
  for (long long I = 0; I < N; I++) {
    cin >> A[I];
    freq[A[I]]++;
  }
  calc(2, 0);
  for (long long I = 0; I < (1 << 20); I++)
    if (__builtin_parity(I))
      sol -= pw[F[I]];
    else
      sol += pw[F[I]];
  cout << sol << endl;
}
