#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <limits>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

using namespace std;

using Int = long long;

template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> ostream &operator<<(ostream &os, const vector<T> &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; }
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }
#define COLOR(s) ("\x1b[" s "m")


struct Op {
  int t, x, y, z;
  friend ostream &operator<<(ostream &os, const Op &op) {
    return os << "(" << op.t << "; " << op.x << "," << op.y << "," << op.z << ")";
  }
};

int N, M, Q, C;

vector<Op> ops;
void oper(int x) {
  assert(0 <= x); assert(x < M);
  ops.push_back(Op{1, x, -1, -1});
}
void oper(int x, int y) {
  assert(0 <= x); assert(x < M);
  assert(0 <= y); assert(y < M);
  assert(x != y);
  ops.push_back(Op{2, x, y, -1});
}
void oper(int x, int y, int z) {
  assert(0 <= x); assert(x < M);
  assert(0 <= y); assert(y < M);
  assert(0 <= z); assert(z < M);
  assert(x != y); assert(x != z); assert(y != z);
  ops.push_back(Op{3, x, y, z});
}

void judge(int t) {
  for (int p = 0; p < 1 << (t ? M : (N + 1)); ++p) {
    int q = p, r = p;
    if (!(~q & ((1 << N) - 1))) q ^= 1 << N;
    for (const auto &op : ops) {
      switch (op.t) {
        case 1: r ^= 1 << op.x; break;
        case 2: if (r >> op.x & 1) r ^= 1 << op.y; break;
        case 3: if ((r >> op.x & 1) && (r >> op.y & 1)) r ^= 1 << op.z; break;
        default: assert(false);
      }
    }
    if (q != r) {
      cerr << "FAIL" << endl;
      r = p;
      for (const auto &op : ops) {
        for (int i = 0; i < M; ++i) cerr << (r >> i & 1);
        cerr << " " << op << endl;
        switch (op.t) {
          case 1: r ^= 1 << op.x; break;
          case 2: r ^= (r >> op.x & 1) << op.y; break;
          case 3: r ^= (r >> op.x & r >> op.y & 1) << op.z; break;
          default: assert(false);
        }
      }
      for (int i = 0; i < M; ++i) cerr << (r >> i & 1);
      cerr << endl;
    }
    assert(q == r);
  }
}

void solve(const vector<int> &us, int tar) {
  const int usLen = us.size();
  if (usLen == 0) {
    oper(tar);
  } else if (usLen == 1) {
    oper(us[0], tar);
  } else if (usLen == 2) {
    oper(us[0], us[1], tar);
  } else if (usLen == 4) {
    // ?
    __int128 used = 0;
    for (const int u : us) used |= (__int128)1 << u;
    used |= (__int128)1 << tar;
    int lot = N + 1;
    for (; used >> lot & 1; --lot) {}
    oper(us[2], us[3], tar);
    oper(us[0], us[1], us[3]);
    oper(us[3], lot, tar);
    oper(us[0], us[1], us[3]);
    oper(us[2], us[3], lot);
    oper(us[0], us[1], us[3]);
    oper(us[3], lot, tar);
    oper(us[0], us[1], us[3]);
    oper(us[2], us[3], lot);
  } else {
    __int128 used = 0;
    for (const int u : us) used |= (__int128)1 << u;
    used |= (__int128)1 << tar;
    int lot = N + 1;
    for (; used >> lot & 1; --lot) {}
    const int half = (usLen + 1) / 2;
    vector<int> vs0, vs1;
    for (int i = 0; i < half; ++i) vs0.push_back(us[i]);
    for (int i = half; i < usLen; ++i) vs1.push_back(us[i]);
    vs1.push_back(lot);
    solve(vs0, lot);
    solve(vs1, tar);
    solve(vs0, lot);
    solve(vs1, tar);
  }
}

int main() {
  for (; ~scanf("%d%d%d%d", &N, &M, &Q, &C); ) {
    ops.clear();
    
    if (C == 1 || C == 2) {
      if (N == 0) {
        oper(0);
      } else if (N == 1) {
        oper(0, 1);
      } else if (N == 2) {
        oper(0, 1, 2);
      } else {
        oper(0, 1, N + 1);
        for (int i = 2; i < N - 1; ++i) oper(N + i - 1, i, N + i);
        oper(N + N - 2, N - 1, N);
        for (int i = N - 1; --i >= 2; ) oper(N + i - 1, i, N + i);
        oper(0, 1, N + 1);
      }
    } else {
      vector<int> us(N);
      for (int u = 0; u < N; ++u) us[u] = u;
      solve(us, N);
    }
    
    printf("%d\n", (int)ops.size());
    for (const auto &op : ops) {
      switch (op.t) {
        case 1: printf("1 %d\n", op.x); break;
        case 2: printf("2 %d %d\n", op.x, op.y); break;
        case 3: printf("3 %d %d %d\n", op.x, op.y, op.z); break;
        default: assert(false);
      }
    }
    
#ifdef LOCAL
    judge(1);
#endif
  }
  return 0;
}