#include <bits/stdc++.h>
using namespace std;
const int MAXN = 11234, MAXM = 32;
inline void swp(int i, int j, vector<pair<int, int> > &seq) {
  seq.push_back(make_pair(i, j));
  seq.push_back(make_pair(j, i));
  seq.push_back(make_pair(i, j));
}
struct gauss {
  bitset<32> e[MAXN];
  int pos[MAXM];
  vector<pair<int, int> > seq;
  void print(int n, int m = 32) const {
    for (int i = 0; i < int(n); i++)
      for (int j = 0; j < int(m); j++)
        printf("%d%c", int(e[i][j]), " \n"[j + 1 == m]);
    puts("-----");
  }
  int reduce(int n, int m = 32) {
    int r = 0;
    for (int c = 0, i, j, k; c < m && r < n; ++c) {
      for (k = r; k < n && !e[k][c]; ++k)
        ;
      if (k >= n)
        continue;
      else if (k != r) {
        swap(e[k], e[r]);
        swp(k, r, seq);
      }
      pos[r] = c;
      for (i = 0; i < n; ++i)
        if (i != r && e[i][c]) {
          e[i] ^= e[r];
          seq.push_back(make_pair(i, r));
        }
      r++;
    }
    return r;
  }
} x, y;
int xx[11234], yy[11234];
int main() {
  int n;
  scanf("%d", &n);
  for (int i = 0; i < int(n); i++) scanf("%d", xx + i);
  for (int i = 0; i < int(n); i++) scanf("%d", yy + i);
  for (int i = 0; i < int(n); i++) {
    x.e[i] = xx[i];
    y.e[i] = yy[i];
  }
  int rx = x.reduce(n);
  int ry = y.reduce(n);
  if (rx < ry) {
    puts("-1");
    return 0;
  }
  int(&pos)[MAXM] = x.pos;
  for (int i = 0; i < int(ry); i++)
    if (y.e[i] != x.e[i]) {
      for (int j = i; j < rx; ++j) {
        if (x.e[i][pos[j]] != y.e[i][pos[j]]) {
          x.e[i] ^= x.e[j];
          x.seq.push_back(make_pair(i, j));
        }
      }
      if (x.e[i] != y.e[i]) {
        puts("-1");
        return 0;
      }
    }
  for (int i = ry; i < rx; ++i) {
    x.e[i] = 0;
    x.seq.push_back(make_pair(i, i));
  }
  int M = x.seq.size() + y.seq.size();
  printf("%d\n", M);
  for (int i = 0; i < int(x.seq.size()); i++) {
    pair<int, int> pr = x.seq[i];
    printf("%d %d\n", pr.first + 1, pr.second + 1);
    xx[pr.first] ^= xx[pr.second];
  }
  for (int j = y.seq.size() - 1; j >= 0; --j) {
    pair<int, int> pr = y.seq[j];
    printf("%d %d\n", pr.first + 1, pr.second + 1);
    xx[pr.first] ^= xx[pr.second];
  }
  for (int i = 0; i < int(n); i++) assert(xx[i] == yy[i]);
}
