#include <iostream>
#include <vector>
#include <algorithm>
#include <string>
#include <ctype.h>
#include <deque>
#include <queue>
#include <cstring>
#include <set>
#include <list>
#include <map>
#include <random>
#include <unordered_map>
#include <stdio.h>

using namespace std;

typedef long long ll;
typedef std::vector<int> vi;
typedef std::vector<bool> vb;
typedef std::vector<string> vs;
typedef std::vector<double> vd;
typedef std::vector<long long> vll;
typedef std::vector<std::vector<int> > vvi;
typedef vector<vvi> vvvi;
typedef vector<vll> vvll;
typedef std::vector<std::pair<int, int> > vpi;
typedef vector<vpi> vvpi;
typedef std::pair<int, int> pi;
typedef std::pair<ll, ll> pll;
typedef std::vector<pll> vpll;

const long long mod = 1000000007;

#define all(c) (c).begin(),(c).end()
#define sz(c) (int)(c).size()
#define forn(i, a, b) for(int i = a; i < b; i++)

#define pb push_back
#define mp make_pair
vvi ans;
void go(ll n, ll b, ll f, ll c, ll add) {
    if(n==1) return;
    if(c >= n-2) {
        int b1 = b;
        if(b>0) {
            b1--;
        }
        go(n-1, b1, f, c-n+2, add);
        ans[add].pb(add+n-1);
        forn(i,1,n-1) ans[add+n-1].pb(add+i);
        if(b1<b) ans[add+n-1].pb(add);
    }
    else {
        ans[add].pb(add+1);
        go(n-1, max(0ll, b-n+1), max(0ll, f-n+2), c, add+1);
        forn(i,0,min(n-1,b)) {
            ans[add+i+1].pb(add);
        }
        forn(i,0,min(n-2,f)) {
            ans[add].pb(add+i+2);
        }
    }
}

int main()
{

    ll t,b,f,c;
    cin>>t>>b>>f>>c;
    ll n = t+1;
    ans.resize(n);
//    if(b+f+c > (n*(n-1))/2 || f+c > ((n-1)*(n-2))/2) {
//        cout<<-1;
//        return 0;
//    }
    ll in = max(t+f, b);
    if(in + c > (n*(n-1))/2) {
        cout<<-1;
        return 0;
    }
    go(n,b,f,c,0);
//    vi g;
//    ll left = ((n-1)*(n-2))/2 - c;
//    int v = n-1;
//    while(v>0) {
//        g.pb(1);
//        v--;
//        while(left >= g.back() && v >0) {
//            v--;
//            left -= g.back();
//            g.back()++;
//        }
//    }
////    for(auto x : g) cout<<x<<' ';
//    vvi ans(n);
//    int l = g.size();
//    int num = 0;
//    forn(i,0,l) {
//        forn(j,0,g[i]) {
//            if(j==0) ans[0].pb(num+j+1);
//            else ans[num+j].pb(num+j+1);
//        }
//        forn(j,0,g[i]) {
//            forn(k,j+1,g[i]) {
//                if(f==0) break;
//                f--;
//                if(j==0) ans[0].pb(num+k+1);
//                else ans[num+j].pb(num+k+1);
//            }
//            if(f==0) break;
//        }
//        forn(j,0,g[i]) {
//            forn(k,0,j) {
//                if(b==0) break;
//                b--;
//                if(k==0) ans[num+j].pb(0);
//                else ans[num+j].pb(num+k);
//            }
//            if(b==0) break;
//        }
//        forn(j,1,num+1) {
//            forn(k,1,g[i]+1) {
//                if(c==0) break;
//                c--;
//                ans[num+k].pb(j);
//            }
//            if(c==0) break;
//        }
//        num += g[i];
//    }
    printf("%lld\n", n);
    forn(i,0,n) {
        printf("%d", (int)ans[i].size());
        for(auto x : ans[i]) printf(" %d", x+1);
        printf("\n");
    }
}