#include <bits/stdc++.h>
using namespace std;
 
using ll = long long;
using ld  = long double;
using pint = pair<int, int>;
using pll = pair<ll, ll>;
using pld = pair<ld, ld>;
const int INF=1e9+7;
const ll LINF=1LL<<60;
const ll MOD=1e9+7;
const ld PI=acos(-1);
const ld EPS = 1e-9; //微調整用（EPSより小さいと0と判定など）

//MODINT
//変数名.nでint型の数値を受け取る
struct mint {
  int n;
  mint(int n_ = 0) : n(n_) {}
};
 
mint operator+(mint a, mint b) { a.n += b.n; if (a.n >= MOD) a.n -= MOD; return a; }
mint operator-(mint a, mint b) { a.n -= b.n; if (a.n < 0) a.n += MOD; return a; }
mint operator*(mint a, mint b) { return (long long)a.n * b.n % MOD; }
mint &operator+=(mint &a, mint b) { return a = a + b; }
mint &operator-=(mint &a, mint b) { return a = a - b; }
mint &operator*=(mint &a, mint b) { return a = a * b; }

int ii() { int x; scanf("%d", &x); return x; }
long long il() { long long x; scanf("%lld", &x); return x; }
string is() { string x; cin >> x; return x; }
char ic() { char x; cin >> x; return x; }
void oi(int x) { printf("%d ", x); }
void ol(long long x) { printf("%lld ", x); }
void od(double x) { printf("%.15f ", x); }
void os(const string &s) { printf("%s ", s.c_str()); }
void oc(const char &c) { printf("%c ", c); }
auto op = [&](auto p) -> void{ cout << p; };
auto ov = [&](auto vec) -> void{ cout << vec; };
#define o_map(v){cerr << #v << endl; for(const auto& xxx: v){cout << xxx.first << " " << xxx.second << "\n";}} //動作未確認
void br() { putchar('\n'); }
 
#define gcd __gcd //llは受け取ってくれない
int lcm(int a, int b){return a / gcd(a, b) * b;}
#define ALL(a) a.begin(),a.end() //sort(ALL(vec));
#define REP(i,m,n) for(int i=(int)(m) ; i < (int) (n) ; ++i )
#define rep(i,n) REP(i,0,n)
#define MP(a,b) make_pair(a,b)
#define SORT_UNIQUE(c) (sort(c.begin(),c.end()), c.resize(distance(c.begin(),unique(c.begin(),c.end()))))
#define GET_POS(c,x) (lower_bound(c.begin(),c.end(),x)-c.begin())
#define PB push_back
#define SZ(x) ((int)(x).size) //size()がunsignedなのでエラー避けに

//4近傍（上下左右）
vector<int> dx_4 = {1, 0, -1, 0};
vector<int> dy_4 = {0, 1, 0, -1};

// coutによるpairの出力（空白区切り）
template<typename T1, typename T2> ostream& operator<<(ostream& s, const pair<T1, T2>& p) {return s << "(" << p.first << " " << p.second << ")";}
// 空白区切りだけ --> return s << "(" << p.first << " " << p.second << ")";
// 見やすくしたいとき --> return s << "(" << p.first << " " << p.second << ")";
// coutによるvectorの出力（空白区切り）
template<typename T> ostream& operator<<(ostream& s, const vector<T>& v) {
	int len = v.size();
	for (int i = 0; i < len; ++i) {
		s << v[i]; if (i < len - 1) s << " "; //"\t"に変えるとTabで見やすく区切る
	}
	return s;
}
// coutによる多次元vectorの出力（空白区切り）
template<typename T> ostream& operator<<(ostream& s, const vector< vector<T> >& vv) {
	int len = vv.size();
	for (int i = 0; i < len; ++i) {
		s << vv[i] << endl;
	}
	return s;
}

//最大値、最小値を更新する。aよりbのが大きい（小さい）か等しければaを更新してtrueを返す。そうでなければ何もせずfalseを返す chmax(nowmax,x);
template<typename T>
bool chmax(T& a, T b){return (a = max(a, b)) == b;}
template<typename T>
bool chmin(T& a, T b){return (a = min(a, b)) == b;}

// ------------ template - //


// - global -------------- //
 

// -------------- global - //


// - library ------------- //

//約数と素因数分解
//divisor(n): 約数が入ったvectorを返す。
//返される約数はソートされていない
//factorization(n): 素因数分解。mapの形で、{素因数, 指数}の形式で返す。
vector<int> divisor(int n){
  vector<int> ans;
  REP(i, 1, sqrt(n) + 1){
    if(!(n % i)){
      ans.PB(i);
      if(i * i < n){
        ans.PB(n / i);
      }
    }
  }
  return ans;
}

// ------------- library - //


int main(){

  int N = ii();

  vector<int> As;
  As.PB(0); //1-indexedにするため、ダミーを入れる
  rep(i, N){
    As.PB(ii());
  }

  vector<int> stats(N+1, 0); //今、ボールが入っている数
  int ans_num = 0;
  vector<int> ans_box;

  rep(i, N){
    // j : N ~ 1 についてループ
    int j = N - i;

    if (As.at(j) == 1){
      if (stats.at(j) % 2 == 0){
        //jの約数全てのstatsに+1
        for ( auto each : divisor(j)){
          stats.at(each)++;
        }
        //ボールを入れる
        ans_num++;
        ans_box.PB(j);
      }
      //既に奇数なら何もしない
    }
    else{
      //目標が偶数なら
      if (stats.at(j) % 2 == 1){
        //jの約数全てのstatsに+1
        for ( auto each : divisor(j)){
          stats.at(each)++;
        }
        //ボールを入れる
        ans_num++;
        ans_box.PB(j);
      }
      //既に偶数なら何もしない
    }
  }

  oi(ans_num); br();
  reverse(ALL(ans_box));
  ov(ans_box);

  


// -- main() end -- //
}


Measured memory_bytes: 3444