#include <bits/stdc++.h>
using namespace std;
int m[2][2], mm[2][2];
char c[10000], a[10000];
int main() {
  int n;
  scanf("%d", &n);
  scanf("%s%s", c + 1, a + 1);
  for (int i = 1; i <= n; i++) m[c[i] - '0'][a[i] - '0']++;
  for (int i = 0; i <= m[0][1]; i++)
    for (int j = 0; j <= m[1][0]; j++)
      if (i + j <= n / 2) {
        int dif = abs(m[0][1] - i - j);
        if ((dif & 1) != (m[1][1] & 1)) continue;
        int n11 = (m[0][1] - i - j + m[1][1]) / 2;
        if (n11 < 0 || n11 > m[1][1]) continue;
        if (i + j + n11 <= n / 2 && i + j + n11 + m[0][0] >= n / 2) {
          mm[0][1] = i;
          mm[1][0] = j;
          mm[1][1] = n11;
          mm[0][0] = n / 2 - i - j - n11;
          for (int k = 1; k <= n; k++)
            if (mm[c[k] - '0'][a[k] - '0']) {
              mm[c[k] - '0'][a[k] - '0']--;
              printf("%d ", k);
            }
          return 0;
        }
      }
  puts("-1");
}
