#include<bits/stdc++.h>
using namespace std;
#define N 3020
typedef long long ll;
int n,m,mod,fac[N*N],ifac[N*N],dp[N][N];
inline void add(int &x,int y){(x+=y)>=mod?x-=mod:233;}
int qpow(int a,int b){
	int ans=1;
	while(b){
		if(b&1)ans=1LL*ans*a%mod;
		a=1LL*a*a%mod;
		b>>=1;
	}
	return ans;
}
void init(int n){
	fac[0]=1;
	for(int i=1;i<=n;++i){
		fac[i]=1LL*fac[i-1]*i%mod;
	}
	ifac[n]=qpow(fac[n],mod-2);
	for(int i=n-1;i>=0;--i){
		ifac[i]=1LL*ifac[i+1]*(i+1)%mod;
	}
}
inline int C(int n,int m){
	return 1LL*fac[n]*ifac[m]%mod*ifac[n-m]%mod;
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	cin>>n>>m>>mod;
	init(n*m);
	dp[1][1]=n*m;
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j){
			add(dp[i+1][j],1LL*dp[i][j]*C(n*m-i*j-1,j-1)%mod*fac[j]%mod*(n-i)%mod);
			add(dp[i][j+1],1LL*dp[i][j]*C(n*m-i*j-1,i-1)%mod*fac[i]%mod*(m-j)%mod);
		}
	}
	cout<<dp[n][m]<<'\n';
	return 0;
}