#pragma GCC optimize ("Ofast")
#include<bits/stdc++.h>
using namespace std;
#define MD (1000000007U)
struct Rand{
unsigned x;
unsigned y;
unsigned z;
unsigned w;
Rand(void){
x=123456789;
y=362436069;
z=521288629;
w=(unsigned)time(NULL);
}
Rand(unsigned seed){
x=123456789;
y=362436069;
z=521288629;
w=seed;
}
inline unsigned get(void){
unsigned t;
t = (x^(x<<11));
x=y;
y=z;
z=w;
w = (w^(w>>19))^(t^(t>>8));
return w;
}
inline double getUni(void){
return get()/4294967296.0;
}
inline int get(int a){
return (int)(a*getUni());
}
inline int get(int a, int b){
return a+(int)((b-a+1)*getUni());
}
inline long long get(long long a){
return(long long)(a*getUni());
}
inline long long get(long long a, long long b){
return a+(long long)((b-a+1)*getUni());
}
inline double get(double a, double b){
return a+(b-a)*getUni();
}
inline int getExp(int a){
return(int)(exp(getUni()*log(a+1.0))-1.0);
}
inline int getExp(int a, int b){
return a+(int)(exp(getUni()*log((b-a+1)+1.0))-1.0);
}
}
;
unsigned long long HashMap_ullP_L[4];
template<class KEY, class VAL> struct HashMap{
char*used;
KEY*key;
VAL*val;
int mem;
int n;
int mask;
HashMap(){
mem = 0;
}
~HashMap(){
free();
}
void expand(int nn){
if(mem >= nn){
return;
}
if(mem){
free();
}
mem = nn;
used = new char[nn];
key = new KEY[nn];
val = new VAL[nn];
}
void free(){
if(mem){
mem = 0;
delete[] used;
delete[] key;
delete[] val;
}
}
void init(int nn){
int i;
n = 1;
nn = nn + (nn + 1) / 2;
while(n < nn){
n *= 2;
}
mask = n - 1;
expand(n);
for(i=(0);i<(n);i++){
used[i] = 0;
}
}
inline int getHash(const int a){
unsigned long long d = a;
d = (((d * HashMap_ullP_L[0]) >> 32) * HashMap_ullP_L[1]) & mask;
return d;
}
inline int getHash(const unsigned a){
unsigned long long d = a;
d = (((d * HashMap_ullP_L[0]) >> 32) * HashMap_ullP_L[1]) & mask;
return d;
}
inline int getHash(const long long a){
unsigned long long d = a;
d = (((((d * HashMap_ullP_L[0]) >> 32) * HashMap_ullP_L[1]) >> 32) * HashMap_ullP_L[2]) & mask;
return d;
}
inline int getHash(const unsigned long long a){
unsigned long long d = a;
d = (((((d * HashMap_ullP_L[0]) >> 32) * HashMap_ullP_L[1]) >> 32) * HashMap_ullP_L[2]) & mask;
return d;
}
inline int getHash(const pair<int,int> a){
unsigned long long d = (((unsigned long long)a.first) << 32) + ((unsigned long long)a.second);
d = (((((d * HashMap_ullP_L[0]) >> 32) * HashMap_ullP_L[1]) >> 32) * HashMap_ullP_L[2]) & mask;
return d;
}
inline VAL& operator[](const KEY a){
int k = getHash(a);
for(;;){
if(used[k]==1 && key[k]==a){
break;
}
if(used[k]==0){
used[k] = 1;
key[k] = a;
break;
}
k = (k+1) & mask;
}
return val[k];
}
inline bool exist(const KEY a){
int k = getHash(a);
for(;;){
if(used[k]==1 && key[k]==a){
return true;
}
if(used[k]==0){
break;
}
k = (k+1) & mask;
}
return false;
}
template<class S> inline bool exist(const KEY a, S &res){
int k = getHash(a);
for(;;){
if(used[k]==1 && key[k]==a){
res = val[k];
return true;
}
if(used[k]==0){
break;
}
k = (k+1) & mask;
}
return false;
}
}
;
#define main dummy_main
int main(){
{
int i;
int j;
int k;
Rand rnd;
for(i=(0);i<(20);i++){
rnd.get(2);
}
for(i=(0);i<(4);i++){
for(j=(0);j<(32);j++){
k = rnd.get(1,62);
HashMap_ullP_L[i] |= (1ULL << k);
}
HashMap_ullP_L[i] |= (1ULL << 0);
HashMap_ullP_L[i] |= (1ULL << 63);
}
}
return 0;
}
#undef main
HashMap<int,int> hs;
class Solution{
public:
int countPairs(vector<int>& A){
int i;
dummy_main();
int N = A.size();
int nx = 1;
long long res = 0;
sort(A.begin(), A.end());
hs.init(N);
for(i=(0);i<(N);i++){
hs[A[i]] = 0;
}
for(i=(0);i<(N);i++){
while(A[i] > nx){
nx *= 2;
}
if(A[i] == nx){
res += hs[A[i]];
}
if(hs.exist(nx-A[i])){
res += hs[nx-A[i]];
}
hs[A[i]]++;
}
return res % MD;
}
}
;
// cLay version 20210103-1
// --- original code ---
// #define main dummy_main
// {}
// #undef main
//
// HashMap<int,int> hs;
//
// class Solution {
// public:
// int countPairs(vector<int>& A) {
// dummy_main();
// int N = A.size(), nx = 1;
// ll res = 0;
// sort(A.begin(), A.end());
// hs.init(N);
// rep(i,N) hs[A[i]] = 0;
// rep(i,N){
// while(A[i] > nx) nx *= 2;
// if(A[i] == nx) res += hs[A[i]];
// if(hs.exist(nx-A[i])) res += hs[nx-A[i]];
// hs[A[i]]++;
// }
// return res % MD;
// }
// };