1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
#include <bits/stdc++.h>
using namespace std;
#define fwd(i, a, n) for (int i = (a); i < (n); i++)
#define rep(i, n) fwd(i, 0, n)
#define all(X) X.begin(), X.end()
#define sz(X) int(size(X))
#define pb push_back
#define eb emplace_back
#define st first
#define nd second
using pii = pair<int, int>; using vi = vector<int>;
using ll = long long; using ld = long double;
#ifdef LOC
auto SS = signal(6, [](int) { *(int *)0 = 0; });
#define DTP(x, y) auto operator << (auto &o, auto a) -> decltype(y, o) { o << "("; x; return o << ")"; }
DTP(o << a.st << ", " << a.nd, a.nd);
DTP(for (auto i : a) o << i << ", ", all(a));
void dump(auto... x) { (( cerr << x << ", " ), ...) << '\n'; }
#define deb(x...) cerr << setw(4) << __LINE__ << ":[" #x "]: ", dump(x)
#else
#define deb(...) 0
#endif

#define S 1000007

ll mod;

ll pot(ll x, ll y) {
	if(y == 0)
		return 1;
	if(y == 1)
		return x;
	if(y % 2 == 0) {
		ll p = pot(x,y/2);
		return (p*p)%mod;
	} 
	return (pot(x,y-1) * x) % mod;
}

ll sil[S];

int32_t main() {
	cin.tie(0)->sync_with_stdio(0);
	cout << fixed << setprecision(10);
	int a,b,c;
	cin >> a >> b >> c >> mod;
	if(c == a+b-1) {
		int n = a*b;
		sil[0] = 1;
		for(int i = 1; i <= n; i++) {
			sil[i] = sil[i-1] * (ll)i;
			sil[i] %= mod;
		}
		ll p = 1;
		for(int i = 1; i <= a; i++) {
			for(int j = 1; j <= b; j++) {
				p *= (ll)(a+b-i-j+1);
				p %= mod;
			}
		}
		p = pot(p,mod-2);
		p *= sil[n];
		p %= mod;
		p *= p;
		p %= mod;
		cout << n << " " << p << "\n";
	} else {
		cout << "dunno\n";
	}
}