#include using namespace std; using ll = long long; const int INF = 1e9 + 10; const ll INFL = 4e18; /* ガウスの消去法によって、F_2 上の連立線形方程式 Ax = 0 を解く。 計算量: O(n^3) 戻り値: 解が存在するかどうか x0: 特殊解 xs: 解空間の基底 */ int gaussJordan(vector> &a, bool extended = false) { int rank = 0; int row = a.size(), col = a[0].size(); for (int c = 0; c < col; c++) { if (extended && c == col - 1) break; int pivot = -1; for (int r = rank; r < row; r++) { if (a[r][c]) { pivot = r; break; } } if (pivot == -1) continue; swap(a[pivot], a[rank]); for (int r = 0; r < row; r++) { if (r == rank) continue; if (a[r][c]) { for (int cc = 0; cc < col; cc++) { a[r][cc] = a[r][cc] ^ a[rank][cc]; } } } rank++; } return rank; } bool solveLinearEquation(vector> a, vector b, vector &x0, vector> &xs) { int row = a.size(), col = a[0].size(); for (int i = 0; i < row; i++) { a[i].push_back(b[i]); } int rank = gaussJordan(a, true); for (int i = 0; i < row; i++) { if (a[i][col]) return false; } x0 = vector(col); for (int i = 0; i < rank; i++) x0[i] = a[i][col]; xs = vector>(col - rank, vector(col)); for (int i = 0; i < col - rank; i++) { xs[i][i + rank] = 1; for (int j = 0; j < rank; j++) { xs[i][j] = a[j][i + rank]; } } return true; } int main() { int N; cin >> N; vector A(N); for (int i = 0; i < N; i++) cin >> A[i]; N = min(N, 60); vector> a(60, vector(N)); for (int i = 0; i < N; i++) { for (int j = 0; j < 60; j++) { a[j][i] = (A[i] >> j) & 1; } } vector x0; vector> xs; if (solveLinearEquation(a, vector(N), x0, xs) && xs.size() > 0) { vector ans; for (int i = 0; i < N; i++) { if (xs[0][i]) { ans.push_back(i); } } if (ans.size() == 0) { cout << -1 << endl; return 0; } cout << ans.size() << endl; for (int x : ans) { cout << x + 1 << ' '; } cout << endl; } else { cout << -1 << endl; } }