#include <cstdio>

constexpr int maxn = 100 + 5;

int n, m, p;

inline void add(int& x, int y) {
    (x += y) >= p && (x -= p);
}

int quick_power(int x, int k) {
    int ret = 1;
    for (;;) {
        if (k & 1) {
            ret = (long long) ret * x % p;
        }
        if (!(k >>= 1)) {
            break;
        }
        x = (long long) x * x % p;
    }
    return ret;
}

int inverse(int x) {
    return quick_power(x, p - 2);
}

namespace Case1 {
    int f[maxn];

    void solve() {
        f[0] = 1;
        for (int i = 1; i <= m; i ++) {
            long long inv = inverse(i);
            for (int j = n; j; j --) {
                add(f[j], f[j - 1] * inv % p);
            }
        }
        printf("%d\n", f[n]);
    }
}

namespace Case2 {
    int f[maxn * maxn][maxn];

    void solve() {
        int ans = 0;

        for (int a1 = 1; a1 <= m; a1 ++) {
            for (int i = 0; i <= m * (m + 1) / 2; i ++) {
                for (int j = 0; j < n; j ++) {
                    f[i][j] = 0;
                }
            }
            f[0][0] = 1;

            for (int i = 1; i <= m; i ++) {
                if (a1 == i) {
                    continue;
                }

                long long inv = inverse(i);
                for (int j = i * (i + 1) / 2; j >= 0; j --) {
                    for (int k = n; k; k --) {
                        add(f[j][k], (f[j][k - 1] - (i > j ? 0 :
                            f[j - i][k - 1]) + p) * inv % p);
                    }
                }
            }

            for (int i = 0; i <= m * (m + 1) / 2; i ++) {
                add(ans, (long long) a1 * inverse(i + a1) % p * f[i][n - 1] % p);
            }
        }

        printf("%d\n", ans);
    }
}

int main() {
    int k;
    scanf("%d%d%d%d", &n, &m, &k, &p);

    if (k == 1) {
        Case1::solve();
    } else {
        Case2::solve();
    }

    return 0;
}