Example

Input

2 2 2 1
0 0 0


Output

24
#include <bits/stdc++.h>
using namespace std;
#define SZ(x) (int)(x.size())
#define REP(i, n) for(int i=0;i<(n);++i)
#define FOR(i, a, b) for(int i=(a);i<(b);++i)
#define RREP(i, n) for(int i=(int)(n)-1;i>=0;--i)
#define RFOR(i, a, b) for(int i=(int)(b)-1;i>=(a);--i)
#define ALL(a) a.begin(),a.end()
#define DUMP(x) cerr<<#x<<" = "<<(x)<<endl
#define DEBUG(x) cerr<<#x<<" = "<<(x)<<" (L"<<__LINE__<<")"<< endl;
using ll = long long;
using vi = vector<int>;
using vvi = vector<vi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using P = pair<int, int>;
const double eps = 1e-8;
const ll MOD = 1000000007;
const int INF = INT_MAX / 2;
const ll LINF = LLONG_MAX / 2;
template <typename T1, typename T2>
bool chmax(T1 &a, const T2 &b) {
    if (a < b) { a = b; return true; }
    return false;
}
template <typename T1, typename T2>
bool chmin(T1 &a, const T2 &b) {
    if (a > b) { a = b; return true; }
    return false;
}
template<typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p) {
    os << p.first << ":" << p.second;
    return os;
}
template<typename T1, typename T2>
ostream &operator<<(ostream &os, const map<T1, T2> &mp) {
    os << "{";
    int a = 0;
    for (auto &tp : mp) {
        if (a) os << ", "; a = 1;
        os << tp;
    }
    return os << "}";
}
template<typename T>
ostream &operator<<(ostream &os, const vector<T> &v) {
    os << "[";
    REP(i, SZ(v)) {
        if (i) os << ", ";
        os << v[i];
    }
    return os << "]";
}


int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(10);

    ll A, B, C, N; cin >> A >> B >> C >> N;
    ll area = 2 * (A * B + B * C + C * A);

    map<tuple<ll,ll,ll>,bool> mp;

    REP(i, N) {
        ll x, y, z; cin >> x >> y >> z;
        area += 6;
        if (x == 0) area -= 2;
        if (x == A-1) area -= 2;
        if (y == 0) area -= 2;
        if (y == B-1) area -= 2;
        if (z == 0) area -= 2;
        if (z == C-1) area -= 2;
        if (mp.count(make_tuple(x-1, y, z))) area -= 2;
        if (mp.count(make_tuple(x+1, y, z))) area -= 2;
        if (mp.count(make_tuple(x, y-1, z))) area -= 2;
        if (mp.count(make_tuple(x, y+1, z))) area -= 2;
        if (mp.count(make_tuple(x, y, z-1))) area -= 2;
        if (mp.count(make_tuple(x, y, z+1))) area -= 2;
        mp[make_tuple(x, y, z)] = 1;
    }

    cout << area << endl;

    return 0;
}

